-
A spectroscopic and kinematic survey of fast hot subdwarfs
Authors:
S. Geier,
U. Heber,
A. Irrgang,
M. Dorsch,
A. Bastian,
P. Neunteufel,
T. Kupfer,
S. Bloemen,
S. Kreuzer,
L. Möller,
M. Schindewolf,
D. Schneider,
E. Ziegerer,
I. Pelisoli,
V. Schaffenroth,
B. N. Barlow,
R. Raddi,
S. J. Geier,
N. Reindl,
T. Rauch,
P. Nemeth,
B. T. Gänsicke
Abstract:
Hot subdwarfs (sdO/B) are the stripped helium cores of red giants formed by binary interactions. Close hot subdwarf binaries with massive white dwarf companions have been proposed as possible progenitors of thermonuclear supernovae type Ia (SN Ia). If the supernova is triggered by stable mass transfer from the helium star, the companion should survive the explosion and should be accelerated to hig…
▽ More
Hot subdwarfs (sdO/B) are the stripped helium cores of red giants formed by binary interactions. Close hot subdwarf binaries with massive white dwarf companions have been proposed as possible progenitors of thermonuclear supernovae type Ia (SN Ia). If the supernova is triggered by stable mass transfer from the helium star, the companion should survive the explosion and should be accelerated to high velocities. The hypervelocity star US 708 is regarded as the prototype for such an ejected companion. To find more of those objects we conducted an extensive spectroscopic survey. Candidates for such fast stars have been selected from the spectroscopic database of the Sloan Digital Sky Survey (SDSS) and several ground-based proper motion surveys. Follow-up spectroscopy has been obtained with several 4m- to 10m-class telescopes. Combining the results from quantitative spectroscopic analyses with space-based astrometry from \textit{Gaia} Early Data Release 3 (EDR3) we determined the atmospheric and kinematic parameters of 53 fast hot subdwarf stars. None of these stars is unbound to the Galaxy, although some have Galactic restframe velocities close to the Galactic escape velocity. 21 stars are apparently single objects, which crossed the Galactic disc within their lifetimes in the sdO/B stage and could be regarded as potential candidates for the SN Ia ejection scenario. However, the properties of the full sample are more consistent with a pure old Galactic halo population. We therefore conclude that the fast sdO/B stars we found are likely to be extreme halo stars.
△ Less
Submitted 5 July, 2024;
originally announced July 2024.
-
SCATTER: A New Common Envelope Formalism
Authors:
Rosanne Di Stefano,
Matthias U. Kruckow,
Yan Gao,
Patrick G. Neunteufel,
Chiaki Kobayashi
Abstract:
One of the most mysterious astrophysical states is the common envelope (CE) phase of binary evolution, in which two stars are enshrouded by the envelope shed by one of them. Interactions between the stars and the envelope shrinks the orbit. The CE can lead to mergers or to a subsequent phase of interactions. Mergers may involve any combination of two compact objects and/or stars. Some involving wh…
▽ More
One of the most mysterious astrophysical states is the common envelope (CE) phase of binary evolution, in which two stars are enshrouded by the envelope shed by one of them. Interactions between the stars and the envelope shrinks the orbit. The CE can lead to mergers or to a subsequent phase of interactions. Mergers may involve any combination of two compact objects and/or stars. Some involving white dwarfs, may produce Type Ia supernovae, while merging neutron stars may yield gamma-ray bursts, and merging compact objects of all kinds produce gravitational radiation. Since CEs can arise from a variety of different initial conditions, and due to the complexity of the processes involved, it is difficult to predict their end states. When many systems are being considered, as in population synthesis calculations, conservation principles are generally employed. Here we use angular momentum in a new way to derive a simple expression for the final orbital separation. This method provides advantages for the study of binaries and is particularly well-suited to higher order multiples, now considered to be important in the genesis of potential mergers. Here we focus on CEs in binaries, and the follow-up paper extends our formalism to multiple star systems within which a CE occurs.
△ Less
Submitted 23 December, 2022; v1 submitted 13 December, 2022;
originally announced December 2022.
-
Triple evolution: an important channel in the formation of type Ia supernovae
Authors:
Abinaya Swaruba Rajamuthukumar,
Adrian Hamers,
Patrick Neunteufel,
Ruediger Pakmor,
Selma E. de mink
Abstract:
Type Ia supernovae (SNe Ia) are thought to be the result of thermonuclear explosions in white dwarfs (WDs). Commonly considered formation pathways include two merging WDs (the double degenerate channel), and a single WD accreting material from a H or He donor (the single degenerate channel). Since the predicted SNe Ia rates from WD in binaries are thought to be insufficient to explain the observed…
▽ More
Type Ia supernovae (SNe Ia) are thought to be the result of thermonuclear explosions in white dwarfs (WDs). Commonly considered formation pathways include two merging WDs (the double degenerate channel), and a single WD accreting material from a H or He donor (the single degenerate channel). Since the predicted SNe Ia rates from WD in binaries are thought to be insufficient to explain the observed SNe Ia rate, it is important to study similar interactions in higher-order multiple star systems such as triple systems. We use the evolutionary population synthesis code Multiple Stellar Evolution (MSE) to study stellar evolution, binary interactions and gravitational dynamics of the triple-star systems. Also, unlike previous studies, prescriptions are included to simultaneously take into account the single and double degenerate channels, and we consider triples across the entire parameter space (including those with tight inner binaries). We explore the impact of typically ignored or uncertain physics such as fly-bys and CE prescription parameters on our results. The majority of systems undergo circular mergers to explode as SNe Ia, while eccentric collisions contribute to $0.4-4$ per cent of SNe Ia events. The time-integrated SNe Ia rate from the triple channel is found to be $(3.60 \pm 0.04) \times {10^{-4}}\,\mathrm{M_{\odot}}^{-1}$ which is, surprisingly, similar to that of the isolated binary channel where the SNe Ia rate is $(3.2 \pm 0.1) \times {10^{-4}}\,\mathrm{M_{\odot}}^{-1}$. This implies that triples, when considering their entire parameter space, yield an important contribution to the overall SNe Ia rate.
△ Less
Submitted 24 March, 2023; v1 submitted 8 November, 2022;
originally announced November 2022.
-
The Equilibrium Tide: An Updated Prescription for Population Synthesis Codes
Authors:
Holly P. Preece,
Adrian S. Hamers,
Patrick G. Neunteufel,
Adam L. Schafer,
Christopher A. Tout
Abstract:
We present an updated prescription for the equilibrium tides suitable for population synthesis codes. A grid of 1D evolutionary models was created and the viscous time-scale was calculated for each detailed model. A metallicity dependent power-law relation was fitted to both the convective cores and convective envelopes of the models. The prescription was implemented into the population synthesis…
▽ More
We present an updated prescription for the equilibrium tides suitable for population synthesis codes. A grid of 1D evolutionary models was created and the viscous time-scale was calculated for each detailed model. A metallicity dependent power-law relation was fitted to both the convective cores and convective envelopes of the models. The prescription was implemented into the population synthesis code BSE and predicts an 16.5% reduction in the overall number of merges, with those involving main-sequence stars most affected. The new prescription also reduces the overall supernova rate by 3.6% with individual channels being differently affected. The single degenerate Ia supernova occurrence is reduced by 12.8%. The merging of two Carbon Oxygen white dwarfs to cause a Ia supernova occurs 16% less frequently. The number of sub-synchronously rotating stars in close binaries is substantially increased with our prescription, as is the number of non-circularized systems at the start of common-envelope evolution.
△ Less
Submitted 23 June, 2022; v1 submitted 13 June, 2022;
originally announced June 2022.
-
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.
-
Properties and applications of a predicted population of runaway He-sdO/B stars ejected from single degenerate He-donor SNe
Authors:
P. Neunteufel,
H. Preece,
M. Kruckow,
S. Geier,
A. S. Hamers,
S. Justham,
Ph. Podsiadlowski
Abstract:
This study builds on previous works, producing the most extensive prediction of the properties of such a hypothetical population to date, taking into account both Chandrasekhar and non-Chandrasekhar mass events. These results are then used to define criteria for membership of this population and characterise putative subpopulations
This study contains 6x10^6 individual ejection trajectories out…
▽ More
This study builds on previous works, producing the most extensive prediction of the properties of such a hypothetical population to date, taking into account both Chandrasekhar and non-Chandrasekhar mass events. These results are then used to define criteria for membership of this population and characterise putative subpopulations
This study contains 6x10^6 individual ejection trajectories out of the Galactic plane calculated with the stellar kinematics framework SHyRT, which are analysed with regard to their bulk observational properties. These are then put into context with the only previously identified population member US\,708 and applied to a number of other possible candidate objects.
We find that two additional previously observed objects possess properties to warrant a designation as candidate objects. Characterisation of these object with respect to the predicted population finds all of them to be extreme in at least one astrometric observable. We find that current observations support a Galactic SN rate on the order of ~3x10^-7/yr to ~2x10^-6/yr, three orders of magnitude below the inferred Galactic SN Ia rate and two orders of magnitude below the formation rate of predicted He-donor progenitors. The number of currently observed population members suggests that the He-donor scenariois not a dominant contributor to the number of observed SNe Ia. However, we find that, even at the low event rate suggested, the majority of possibly detectable population members is still undetected. The extreme nature of current population members suggests that a still larger number of objects has simply evaded detection up to this point, hinting at a higher contribution than is currently supported by observation. - abridged -
△ Less
Submitted 28 March, 2022; v1 submitted 14 December, 2021;
originally announced December 2021.
-
A statistical view on stable and unstable Roche lobe overflow of a tertiary star onto the inner binary in triple systems
Authors:
Adrian S. Hamers,
Hila Glanz,
Patrick Neunteufel
Abstract:
In compact stellar triple systems, an evolved tertiary star can overflow its Roche lobe around the inner binary. Subsequently, the tertiary star can transfer mass to the inner binary in a stable manner, or Roche lobe overflow (RLOF) can be unstable and lead to common-envelope (CE) evolution. In the latter case, the inner binary enters the extended envelope of the tertiary star and spirals in towar…
▽ More
In compact stellar triple systems, an evolved tertiary star can overflow its Roche lobe around the inner binary. Subsequently, the tertiary star can transfer mass to the inner binary in a stable manner, or Roche lobe overflow (RLOF) can be unstable and lead to common-envelope (CE) evolution. In the latter case, the inner binary enters the extended envelope of the tertiary star and spirals in towards the donor's core, potentially leading to mergers or ejections. Although studied in detail for individual systems, a comprehensive statistical view on the various outcomes of triple RLOF is lacking. Here, we carry out 10^5 population synthesis simulations of tight triples, self-consistently taking into account stellar evolution, binary interactions, and gravitational dynamics. Also included are prescriptions for the long-term evolution of stable triple mass transfer, and triple CE evolution. Although simple and ignoring hydrodynamic effects, these prescriptions allow for a qualitative statistical study. We find that triple RLOF occurs in ~0.06% of triples, with ~64% leading to stable mass transfer, and ~36% to triple CE evolution. Triple CE is most often (~76%) followed by one or multiple mergers in short succession, most likely an inner binary merger of two main-sequence stars. Other outcomes of triple CE are a binary+single system (~23%, most of which not involving exchange interactions), and a stable triple (~1%). We also estimate the rate of Type Ia supernovae involving white dwarf mergers following triple RLOF, but find only a negligible contribution.
△ Less
Submitted 5 January, 2022; v1 submitted 30 September, 2021;
originally announced October 2021.
-
Return of the TEDI: revisiting the Triple Evolution Dynamical Instability channel in triple stars
Authors:
Adrian S. Hamers,
Hagai B. Perets,
Todd A. Thompson,
Patrick Neunteufel
Abstract:
Triple-star systems exhibit a phenomenon known as the Triple Evolution Dynamical Instability (TEDI), in which mass loss in evolving triples triggers short-term dynamical instabilities, potentially leading to collisions of stars, exchanges, and ejections. Previous work has shown that the TEDI is an important pathway to head-on stellar collisions in the Galaxy, significantly exceeding the rate of co…
▽ More
Triple-star systems exhibit a phenomenon known as the Triple Evolution Dynamical Instability (TEDI), in which mass loss in evolving triples triggers short-term dynamical instabilities, potentially leading to collisions of stars, exchanges, and ejections. Previous work has shown that the TEDI is an important pathway to head-on stellar collisions in the Galaxy, significantly exceeding the rate of collisions due to random encounters in globular clusters. Here, we revisit the TEDI evolutionary pathway using state-of-the-art population synthesis methods that self-consistently take into account stellar evolution and binary interactions, as well as gravitational dynamics and perturbations from passing stars in the field. We find Galactic TEDI-induced collision rates on the order of 1e-4/yr, consistent with previous studies which were based on more simplified methods. The majority of TEDI-induced collisions involve main sequence stars, potentially producing blue straggler stars. Collisions are also possible involving more evolved stars, potentially producing eccentric post-common-envelope systems, and white dwarfs collisions leading to Type Ia supernovae (although the latter with a negligible contribution to the Galactic rate). In our simulations, the TEDI is not only triggered by adiabatic wind mass loss, but also by Roche lobe overflow in the inner binary: when the donor star becomes less massive than the accretor, the inner binary orbit widens, triggering triple dynamical instability. We find that collision rates are increased by ~17% when fly-bys in the field are taken into account. In addition to collisions, we find that the TEDI produces ~1e-4/yr of unbound stars, although none with escape speeds in excess of 1e3 km/s.
△ Less
Submitted 2 December, 2021; v1 submitted 28 July, 2021;
originally announced July 2021.
-
A hot subdwarf-white dwarf super-Chandrasekhar candidate supernova Ia progenitor
Authors:
Ingrid Pelisoli,
P. Neunteufel,
S. Geier,
T. Kupfer,
U. Heber,
A. Irrgang,
D. Schneider,
A. Bastian,
J. van Roestel,
V. Schaffenroth,
B. N. Barlow
Abstract:
Supernova Ia are bright explosive events that can be used to estimate cosmological distances, allowing us to study the expansion of the Universe. They are understood to result from a thermonuclear detonation in a white dwarf that formed from the exhausted core of a star more massive than the Sun. However, the possible progenitor channels leading to an explosion are a long-standing debate, limiting…
▽ More
Supernova Ia are bright explosive events that can be used to estimate cosmological distances, allowing us to study the expansion of the Universe. They are understood to result from a thermonuclear detonation in a white dwarf that formed from the exhausted core of a star more massive than the Sun. However, the possible progenitor channels leading to an explosion are a long-standing debate, limiting the precision and accuracy of supernova Ia as distance indicators. Here we present HD265435, a binary system with an orbital period of less than a hundred minutes, consisting of a white dwarf and a hot subdwarf -- a stripped core-helium burning star. The total mass of the system is 1.65+/-0.25 solar-masses, exceeding the Chandrasekhar limit (the maximum mass of a stable white dwarf). The system will merge due to gravitational wave emission in 70 million years, likely triggering a supernova Ia event. We use this detection to place constraints on the contribution of hot subdwarf-white dwarf binaries to supernova Ia progenitors.
△ Less
Submitted 19 July, 2021;
originally announced July 2021.
-
A Catalogue of Potential Post Common Envelope Binaries
Authors:
Matthias U. Kruckow,
Patrick G. Neunteufel,
Rosanne Di Stefano,
Yan Gao,
Chiaki Kobayashi
Abstract:
We present a catalogue containing 839 candidate post common envelope systems. Common envelope evolution is very important in stellar astrophysics, particularly in the context of very compact and short-period binaries, including cataclysmic variables, as progenitors of e.g. supernovae type Ia or mergers of black holes and/or neutron stars. At the same time it is a barely understood process in binar…
▽ More
We present a catalogue containing 839 candidate post common envelope systems. Common envelope evolution is very important in stellar astrophysics, particularly in the context of very compact and short-period binaries, including cataclysmic variables, as progenitors of e.g. supernovae type Ia or mergers of black holes and/or neutron stars. At the same time it is a barely understood process in binary evolution. Due to limitations, since partially remedied, on direct simulation, early investigations were mainly focused on providing analytic prescriptions of the outcome of common envelope evolution. In recent years, detailed hydrodynamical calculations have produced deeper insight into the previously elusive process of envelope ejection. However, a direct link between observations and theory of this relatively short-lived phase in binary evolution has not been forthcoming. Therefore, the main insight to be gained from observations has to be derived from the current state of systems likely to have gone through a common envelope. Here we present an extensive catalogue of such observations as found in the literature. The aim of this paper is to provide a reliable set of data, obtained from observations, to be used in the theoretical modelling of common envelope evolution. In this catalogue, the former common envelope donor star is commonly observed as a white dwarf star or as a hot sub-dwarf star. This catalogue includes period and mass estimates, wherever obtainable. Some binaries are border line cases to allow an investigation of the transition between a common envelope formation and other mass-transfer processes.
△ Less
Submitted 12 July, 2021;
originally announced July 2021.
-
First and second-generation black hole and neutron star mergers in 2+2 quadruples: population statistics
Authors:
Adrian S. Hamers,
Giacomo Fragione,
Patrick Neunteufel,
Bence Kocsis
Abstract:
Recent detections of gravitational waves from mergers of neutron stars (NSs) and black holes (BHs) in the low and high-end mass gap regimes pose a puzzle to standard stellar and binary evolution theory. Mass-gap mergers may originate from successive mergers in hierarchical systems such as quadruples. Here, we consider repeated mergers of NSs and BHs in stellar 2+2 quadruple systems, in which secul…
▽ More
Recent detections of gravitational waves from mergers of neutron stars (NSs) and black holes (BHs) in the low and high-end mass gap regimes pose a puzzle to standard stellar and binary evolution theory. Mass-gap mergers may originate from successive mergers in hierarchical systems such as quadruples. Here, we consider repeated mergers of NSs and BHs in stellar 2+2 quadruple systems, in which secular evolution can accelerate the merger of one of the inner binaries. Subsequently, the merger remnant may interact with the companion binary, yielding a second-generation merger. We model the initial stellar and binary evolution of the inner binaries as isolated systems. In the case of successful compact object formation, we subsequently follow the secular dynamical evolution of the quadruple system. When a merger occurs, we take into account merger recoil, and model subsequent evolution using direct N-body integration. With different assumptions on the initial properties, we find that the majority of first-generation mergers are not much affected by secular evolution, with their observational properties mostly consistent with isolated binaries. A small subset shows imprints of secular evolution through residual eccentricity in the LIGO band, and retrograde spin-orbit orientations. Second-generation mergers are ~10^7 times less common than first-generation mergers, and can be strongly affected by scattering (i.e., three-body interactions) induced by the first-generation merger. In particular, scattering can account for mergers within the low-end mass gap, although not the high-end mass gap. Also, in a few cases, scattering could explain highly eccentric LIGO sources and negative effective spin parameters.
△ Less
Submitted 6 August, 2021; v1 submitted 5 March, 2021;
originally announced March 2021.
-
Predicted spatial and velocity distributions of ejected companion stars of helium accretion-induced thermonuclear supernovae
Authors:
P. Neunteufel,
M. Kruckow,
S. Geier,
A. S. Hamers
Abstract:
Thermonuclear supernovae (SNe), a subset of which are the highly important SNe Type\,Ia, remain one of the more poorly understood phenomena known to modern astrophysics. In recent years, the single degenerate helium (He) donor channel, where a white dwarf star (WD) accretes He-rich matter from a hydrogen-depleted companion, has emerged as a promising candidate progenitor scenario for these events.…
▽ More
Thermonuclear supernovae (SNe), a subset of which are the highly important SNe Type\,Ia, remain one of the more poorly understood phenomena known to modern astrophysics. In recent years, the single degenerate helium (He) donor channel, where a white dwarf star (WD) accretes He-rich matter from a hydrogen-depleted companion, has emerged as a promising candidate progenitor scenario for these events. An unresolved question in this scenario is the fate of the companion star, which would be evident as a runaway hot subdwarf (He sdO/B) in the aftermath of the SN event. Previous studies have shown that the kinematic properties of an ejected companion provide an opportunity to closer examination of the properties of an SN progenitor system. However, with the number of observed objects not matching predictions by theory, the viability of this mechanism is called into question. In this study, we first synthesize a population of companion stars ejected by the aforementioned mechanism, taking into account predicted ejection velocities, inferred population density in the Galactic (Gal.) mass distribution and subsequent kinematics in the Gal. potential. We then discuss the astrometric properties of this population. We present $10^{6}$ individual ejection trajectories, numerically computed with a newly developed, lightweight simulation framework. A peak in the density distribution for close objects is expected in the direction of the Gal. center. If the entire considered mass range is realized, the radial velocity distribution should show a peak at 500\kms. If only close US\,708 analogues are considered, there should be a peak at ($\sim750-850$)\kms. We show that the puzzling lack of confirmed surviving companion stars of thermonuclear SNe, though possibly an observation-related selection effect, may indicate a selection against high mass donors in the SD He donor channel. (-abridged-)
△ Less
Submitted 20 January, 2021; v1 submitted 30 November, 2020;
originally announced November 2020.
-
Multiple Stellar Evolution: a population synthesis algorithm to model the stellar, binary, and dynamical evolution of multiple-star systems
Authors:
Adrian S. Hamers,
Antti Rantala,
Patrick Neunteufel,
Holly Preece,
Pavan Vynatheya
Abstract:
In recent years, observations have shown that multiple-star systems such as hierarchical triple and quadruple-star systems are common, especially among massive stars. They are potential sources of interesting astrophysical phenomena such as compact object mergers, leading to supernovae, and gravitational wave events. However, many uncertainties remain in their often complex evolution. Here, we pre…
▽ More
In recent years, observations have shown that multiple-star systems such as hierarchical triple and quadruple-star systems are common, especially among massive stars. They are potential sources of interesting astrophysical phenomena such as compact object mergers, leading to supernovae, and gravitational wave events. However, many uncertainties remain in their often complex evolution. Here, we present the population synthesis code Multiple Stellar Evolution (MSE), designed to rapidly model the stellar, binary, and dynamical evolution of multiple-star systems. MSE includes a number of new features not present in previous population synthesis codes: (1) an arbitrary number of stars, as long as the initial system is hierarchical, (2) dynamic switching between secular and direct N-body integration for efficient computation of the gravitational dynamics, (3) treatment of mass transfer in eccentric orbits, which occurs commonly in multiple-star systems, (4) a simple treatment of tidal, common-envelope, and mass transfer evolution in which the accretor is a binary instead of a single star, (5) taking into account planets within the stellar system, and (6) including gravitational perturbations from passing field stars. MSE, written primarily in the C++ language, will be made publicly available and has few prerequisites; a convenient Python interface is provided. We give a detailed description of MSE and illustrate how to use the code in practice. We demonstrate its operation in a number of examples.
△ Less
Submitted 18 February, 2021; v1 submitted 9 November, 2020;
originally announced November 2020.
-
Velocity limits in the thermonuclear supernova ejection scenario for hypervelocity stars and the origin of US 708
Authors:
Patrick Neunteufel
Abstract:
Hypervelocity stars (HVS) are a class of stars moving at high enough velocities to be gravitationally unbound from the Galaxy. Ejection from a close binary system in which one of the components undergoes a thermonuclear supernova (SN) has emerged as a promising candidate production mechanism for the least massive specimens of this class. This study presents a thorough theoretical analysis of candi…
▽ More
Hypervelocity stars (HVS) are a class of stars moving at high enough velocities to be gravitationally unbound from the Galaxy. Ejection from a close binary system in which one of the components undergoes a thermonuclear supernova (SN) has emerged as a promising candidate production mechanism for the least massive specimens of this class. This study presents a thorough theoretical analysis of candidate progenitor systems of thermonuclear SNe in the single degenerate helium donor scenario in the relevant parameter space leading to the ejection of HVS. The primary goal is investigation of the, previously unclear, characteristics of the velocity spectra of the ejected component. Presented are the results of 390 binary model sequences computed with the MESA framework, investigating the evolution of supernova progenitors composed of a helium-rich hot subdwarf and a accreting white dwarf. Results are then correlated with an idealized kinematic analysis of the observed object US 708. It is seen that the ejection velocity spectra reach a maximum in the range $0.19~M_\odot < M_{HVS} < 0.25~M_\odot$. Depending on the local Galactic potential, all donors below $0.4~\text{M}_\odot$ are expected to become HVS. This channel is able to account for runaway velocities up to $\sim1150~\text{km s}^{-1}$ with a Chandrasekhar mass accretor, exceeding $1200~\text{km s}^{-1}$ if super-Chandrasekhar mass detonations are taken into account. It is found that the previously assumed mass of $0.3~M_\odot$ for US 708, combined with more recently obtained proper motions, favor a sub-Chandrasekhar mass explosion with a terminal WD mass between $1.1~M_\odot$ and $1.2~M_\odot$. The presence of clear ejection velocity maxima provides constraints on the terminal state of a supernova progenitor. It is possible to discern certain types of explosion mechanisms from the inferred ejection velocities alone.
△ Less
Submitted 24 June, 2020; v1 submitted 19 June, 2020;
originally announced June 2020.
-
TW Hya: an old protoplanetary disc revived by its planet
Authors:
Sergei Nayakshin,
Takashi Tsukagoshi,
Cassandra Hall,
Allona Vazan,
Ravit Helled,
Jack Humphries,
Farzana Meru,
Patrick Neunteufel,
Olja Panic
Abstract:
Dark rings with bright rims are the indirect signposts of planets embedded in protoplanetary discs. In a recent first, an azimuthally elongated AU-scale blob, possibly a planet, was resolved with ALMA in TW Hya. The blob is at the edge of a cliff-like rollover in the dust disc rather than inside a dark ring. Here we build time-dependent models of TW Hya disc. We find that the classical paradigm ca…
▽ More
Dark rings with bright rims are the indirect signposts of planets embedded in protoplanetary discs. In a recent first, an azimuthally elongated AU-scale blob, possibly a planet, was resolved with ALMA in TW Hya. The blob is at the edge of a cliff-like rollover in the dust disc rather than inside a dark ring. Here we build time-dependent models of TW Hya disc. We find that the classical paradigm cannot account for the morphology of the disc and the blob. We propose that ALMA-discovered blob hides a Neptune mass planet losing gas and dust. We show that radial drift of mm-sized dust particles naturally explains why the blob is located on the edge of the dust disc. Dust particles leaving the planet perform a characteristic U-turn relative to it, producing an azimuthally elongated blob-like emission feature. This scenario also explains why a 10 Myr old disc is so bright in dust continuum. Two scenarios for the dust-losing planet are presented. In the first, a dusty pre-runaway gas envelope of about 40 Earth mass Core Accretion planet is disrupted, e.g., as a result of a catastrophic encounter. In the second, a massive dusty pre-collapse gas giant planet formed by Gravitational Instability is disrupted by the energy released in its massive core. Future modelling may discriminate between these scenarios and allow us to study planet formation in an entirely new way -- by analysing the flows of dust and gas recently belonging to planets, informing us about the structure of pre-disruption planetary envelopes.
△ Less
Submitted 14 May, 2020; v1 submitted 21 April, 2020;
originally announced April 2020.
-
The Evolution of Helium Star Plus Carbon-Oxygen White Dwarf Binary Systems and Implications for Diverse Stellar Transients and Hypervelocity Stars
Authors:
P. Neunteufel,
S. -C. Yoon,
N. Langer
Abstract:
Helium accretion induced explosions in CO white dwarfs (WDs) are considered promising candidates for a number of observed types of stellar transients, including supernovae (SNe) of Type Ia and Type Iax. However, a clear favorite outcome has not yet emerged. We explore the conditions of helium ignition in the white dwarf and the final fates of helium star-WD binaries as function of their initial or…
▽ More
Helium accretion induced explosions in CO white dwarfs (WDs) are considered promising candidates for a number of observed types of stellar transients, including supernovae (SNe) of Type Ia and Type Iax. However, a clear favorite outcome has not yet emerged. We explore the conditions of helium ignition in the white dwarf and the final fates of helium star-WD binaries as function of their initial orbital periods and component masses. We compute 274 model binary systems with the Binary Evolution Code (BEC), where both components are fully resolved. Stellar and orbital evolution is computed simultaneously, including mass and angular momentum transfer, tides, and gravitational wave emission, as well as differential rotation and internal hydrodynamic and magnetic angular momentum transport. We find that helium detonations are expected only in systems with the shortest initial orbital periods, and for initially massive white dwarfs (MWD > 1.0 MSun ) and lower mass donors (Mdonor < 0.8 MSun), with accumulated helium layers mostly exceeding 0.1 MSun. Upon detonation, these systems would release the donor as a hypervelocity pre-WD runaway star, for which we predict the expected range of kinematic and stellar properties. Systems with more massive donors or initial periods exceeding 1.5 h will likely undergo helium deflagrations after accumulating 0.1 - 0.001 MSun of helium. Helium ignition in the white dwarf is avoided in systems with helium donor stars below - 0.6 MSun, and lead to three distinctly different groups of double white dwarf systems. The size of the parameter space open to helium detonation corresponds to only about 3 % of the galactic SN Ia rate, and to 10 % of the SN Iax rate, while the predicted large amounts of helium (>0.1 MSun) in progenitors cannot easily be reconciled with observations of archetypical SN Ia. ...
△ Less
Submitted 28 April, 2019;
originally announced April 2019.
-
Helium ignition in rotating magnetized CO white dwarfs leading to fast and faint rather than classical type Ia supernovae
Authors:
P. Neunteufel,
S. -C. Yoon,
N. Langer
Abstract:
Based mostly on stellar models which do not include rotation, CO white dwarfs which accrete helium at rates of about $\sim 10^{-8}~\mathrm{M}_\odot/\mathrm{yr}$ have been put forward as candidate progenitors for a number of transient astrophysical phenomena, including supernovae of Type Ia, and the peculiar and fainter Type Iax supernovae. Here we study the impact of accretion-induced spin-up incl…
▽ More
Based mostly on stellar models which do not include rotation, CO white dwarfs which accrete helium at rates of about $\sim 10^{-8}~\mathrm{M}_\odot/\mathrm{yr}$ have been put forward as candidate progenitors for a number of transient astrophysical phenomena, including supernovae of Type Ia, and the peculiar and fainter Type Iax supernovae. Here we study the impact of accretion-induced spin-up including the subsequent magnetic field generation, angular momentum transport, and viscous heating on the white dwarf evolution up to the point of helium ignition. We resolve the structure of the helium accreting white dwarf models with a one dimensional Langrangian hydrodynamic code, modified to include rotational and magnetic effects. We find magnetic angular momentum transport, which leads to quasi solid-body rotation, profoundly impacts the evolution of the white dwarf models. Our rotating lower mass ($0.54$ and $0.82~\mathrm{M}_\odot$) models accrete up to 50\% more mass up to ignition compared to the non-rotating case and the opposite for our more massive models. Furthermore, we find that rotation leads to up to 10-times smaller helium ignition densities, except for the lowest adopted initial white dwarf mass. Ignition densities of order $10^6\,\mathrm{g/cm}^3$ are only found for the lowest accretion rates and for large amounts of accreted helium $\gtrsim 0.4\,\mathrm{M}_\odot$. However, correspondingly massive donor stars would transfer mass at much higher rates. We expect explosive He-shell burning to mostly occur as deflagrations and at $\dot{M}>2\cdot10^{-8}~\mathrm{M}_\odot/\mathrm{yr}$, regardless of white dwarf mass. Our results imply that helium accretion onto CO white dwarfs at the considered rates is unlikely to lead to explosions of the CO core or to classical Type,Ia supernovae, but may instead produce events like faint and fast hydrogen-free supernovae.
△ Less
Submitted 28 February, 2017; v1 submitted 3 January, 2017;
originally announced January 2017.
-
Models for the evolution of close binaries with He-Star and WD components towards Type Ia supernova explosions
Authors:
P. Neunteufel,
S. -C. Yoon,
N. Langer
Abstract:
Type Ia supernovae (SNe Ia) have been an important tool for astronomy for quite some time; however, the nature of their progenitors remains somewhat mysterious. Recent theoretical studies indicated the possibility of producing thermonuclear detonations of carbon-oxygen white dwarfs (CO WDs) at masses less than the Chandrasekhar mass through accretion of helium-rich matter, which would, depending o…
▽ More
Type Ia supernovae (SNe Ia) have been an important tool for astronomy for quite some time; however, the nature of their progenitors remains somewhat mysterious. Recent theoretical studies indicated the possibility of producing thermonuclear detonations of carbon-oxygen white dwarfs (CO WDs) at masses less than the Chandrasekhar mass through accretion of helium-rich matter, which would, depending on mass accretion rate, mass, and initial temperature of the WD, spectrally resemble either a normal SN Ia or a peculiar one. This study aims to further resolve the state of binary systems comprised of a sub-Chandrasekhar-mass CO WD and a helium star at the point where an accretion-induced detonation occurs and constrains the part of the initial parameter space where this kind of phenomenon is possible. Preexisting data obtained through simulations of single, constantly accreting CO WDs is used as an indicator for the behavior of new binary models in which the WD is treated as a point mass and which include the non-degenerate partner as a fully resolved stellar model. We parameterize the ignition of the accumulated helium layer, changes in the WD temperature, and changes in the CO core mass depending on the mass transfer rate. The initial conditions allowing for detonation do not form a single contiguous area in the parameter space, whose shape is profoundly influenced by the behavior of the donor star. Mass loss due to Nova outbursts acts in favor of detonation. According to our criteria, about 10% of the detonations in this study can be expected to show spectra consistent with ordinary SNe Ia; the rest exhibit peculiar features.
△ Less
Submitted 2 March, 2016;
originally announced March 2016.
-
Feasibility of radar detection of extensive air showers
Authors:
J. Stasielak,
R. Engel,
S. Baur,
P. Neunteufel,
R. Šmída,
F. Werner,
H. Wilczyński
Abstract:
Reflection of radio waves off the short-lived plasma produced by the high-energy shower particles in the air is simulated, considering various radar setups and shower geometries. We show that the plasma produced by air showers has to be treated always as underdense. Therefore, we use the Thomson cross-section for scattering of radio waves corrected for molecular quenching and we sum coherently con…
▽ More
Reflection of radio waves off the short-lived plasma produced by the high-energy shower particles in the air is simulated, considering various radar setups and shower geometries. We show that the plasma produced by air showers has to be treated always as underdense. Therefore, we use the Thomson cross-section for scattering of radio waves corrected for molecular quenching and we sum coherently contributions of the reflected radio wave over the volume of the plasma disk to obtain the time evolution of the signal arriving at the receiver antenna. The received power and the spectral power density of the radar echo are analyzed. Based on the obtained results, we discuss possible modes of radar detection of extensive air showers. We conclude that the scattered signal is too weak for the radar method to provide an efficient and inexpensive method of air shower detection.
△ Less
Submitted 4 September, 2015; v1 submitted 26 November, 2014;
originally announced November 2014.
-
First Experimental Characterization of Microwave Emission from Cosmic Ray Air Showers
Authors:
R. Smida,
F. Werner,
R. Engel,
J. C. Arteaga-Velazquez,
K. Bekk,
M. Bertaina,
J. Bluemer,
H. Bozdog,
I. M. Brancus,
A. Chiavassa,
F. Cossavella,
F. Di Pierro,
P. Doll,
B. Fuchs,
D. Fuhrmann,
C. Grupen,
A. Haungs,
D. Heck,
J. R. Hoerandel,
D. Huber,
T. Huege,
K. -H. Kampert,
D. Kang,
H. Klages,
M. Kleifges
, et al. (33 additional authors not shown)
Abstract:
We report the first direct measurement of the overall characteristics of microwave radio emission from extensive air showers. Using a trigger provided by the KASCADE-Grande air shower array, the signals of the microwave antennas of the CROME (Cosmic-Ray Observation via Microwave Emission) experiment have been read out and searched for signatures of radio emission by high-energy air showers in the…
▽ More
We report the first direct measurement of the overall characteristics of microwave radio emission from extensive air showers. Using a trigger provided by the KASCADE-Grande air shower array, the signals of the microwave antennas of the CROME (Cosmic-Ray Observation via Microwave Emission) experiment have been read out and searched for signatures of radio emission by high-energy air showers in the GHz frequency range. Microwave signals have been detected for more than 30 showers with energies above 3*10^16 eV. The observations presented in this Letter are consistent with a mainly forward-directed and polarised emission process in the GHz frequency range. The measurements show that microwave radiation offers a new means of studying air showers at energies above 10^17 eV.
△ Less
Submitted 30 October, 2014;
originally announced October 2014.
-
Enhancement of the radar signal of air showers due to time compression
Authors:
J. Stasielak,
S. Baur,
R. Engel,
P. Neunteufel,
J. Pekala,
R. Šmída,
F. Werner,
H. Wilczyński
Abstract:
We investigate the feasibility of detecting extensive air showers by the radar technique at viewing angles smaller than $\sim 25^\circ$ to the shower axis. Considering a bistatic radar setup and shower geometries in which the receiver points into the arriving shower, we simulate reflection of radio waves off the short-lived plasma produced by the high-energy shower particles in the air. Using the…
▽ More
We investigate the feasibility of detecting extensive air showers by the radar technique at viewing angles smaller than $\sim 25^\circ$ to the shower axis. Considering a bistatic radar setup and shower geometries in which the receiver points into the arriving shower, we simulate reflection of radio waves off the short-lived plasma produced by the high-energy shower particles in the air. Using the Thomson cross-section for scattering of radio waves and summing coherently contributions of the reflected radio wave over the volume of the plasma disk, we obtain the time evolution of the signal. We analyze the spectral power density of the radar echo and the received power. Based on the obtained results, we discuss possible modes of the radar detection of extensive air showers.
△ Less
Submitted 2 October, 2013;
originally announced October 2013.
-
Observation of Polarised Microwave Emission from Cosmic Ray Air Showers
Authors:
R. Smida,
F. Werner,
R. Engel,
J. C. Arteaga-Velazquez,
K. Bekk,
M. Bertaina,
J. Bluemer,
H. Bozdog,
I. M. Brancus,
A. Chiavassa,
F. Cossavella,
F. Di Pierro,
P. Doll,
B. Fuchs,
D. Fuhrmann,
C. Grupen,
A. Haungs,
D. Heck,
J. R. Hoerandel,
D. Huber,
T. Huege,
K. -H. Kampert,
D. Kang,
H. Klages,
M. Kleifges
, et al. (33 additional authors not shown)
Abstract:
We report on the first direct measurement of the basic features of microwave radio emission from extensive air showers. Using a trigger provided by the KASCADE-Grande air shower array, the signals of the microwave antennas of the CROME (Cosmic-Ray Observation via Microwave Emission) experiment have been read out and searched for signatures of radio emission by high-energy air showers. Microwave si…
▽ More
We report on the first direct measurement of the basic features of microwave radio emission from extensive air showers. Using a trigger provided by the KASCADE-Grande air shower array, the signals of the microwave antennas of the CROME (Cosmic-Ray Observation via Microwave Emission) experiment have been read out and searched for signatures of radio emission by high-energy air showers. Microwave signals have been detected for more than 30 showers with energies above $3\times10^{16}$\,eV. The observations presented in this Letter are consistent with a mainly forward-beamed, coherent and polarised emission process in the GHz frequency range. An isotropic, unpolarised radiation is disfavoured as the dominant emission model. The measurements show that microwave radiation offers a new means of studying air showers at very high energy.
△ Less
Submitted 28 June, 2013;
originally announced June 2013.
-
Radar reflection off extensive air showers
Authors:
J. Stasielak,
S. Baur,
M. Bertaina,
J. Blümer,
A. Chiavassa,
R. Engel,
A. Haungs,
T. Huege,
K. -H. Kampert,
H. Klages,
M. Kleifges,
O. Krömer,
M. Ludwig,
S. Mathys,
P. Neunteufel,
J. Pekala,
J. Rautenberg,
M. Riegel,
M. Roth,
F. Salamida,
H. Schieler,
R. Šmída,
M. Unger,
M. Weber,
F. Werner
, et al. (2 additional authors not shown)
Abstract:
We investigate the possibility of detecting extensive air showers by the radar technique. Considering a bistatic radar system and different shower geometries, we simulate reflection of radio waves off the static plasma produced by the shower in the air. Using the Thomson cross-section for radio wave reflection, we obtain the time evolution of the signal received by the antennas. The frequency upsh…
▽ More
We investigate the possibility of detecting extensive air showers by the radar technique. Considering a bistatic radar system and different shower geometries, we simulate reflection of radio waves off the static plasma produced by the shower in the air. Using the Thomson cross-section for radio wave reflection, we obtain the time evolution of the signal received by the antennas. The frequency upshift of the radar echo and the power received are studied to verify the feasibility of the radar detection technique.
△ Less
Submitted 4 October, 2012;
originally announced October 2012.
