Search

Supernova explosion and the associated neutron star natal kicks are important events on a pathway of a binary to become a gravitational wave source, an X-ray binary or a millisecond radio pulsar. Weak natal kicks often lead to binary survival, while strong kicks frequently disrupt the binary. In this article, we aim to further constrain neutron star natal kicks in binaries. We explore binary population synthesis models by varying prescription for natal kick, remnant mass and mass accretion efficiency. We introduce a robust statistical technique to analyse combined observations of different nature. Using this technique, we further test different models using parallax and proper motion measurements for young isolated radio pulsars and similar measurements for Galactic Be X-ray binaries. Our best model for natal kicks is consistent with both measurements and contains a fraction of $w=0.2\pm 0.1$ weak natal kicks with $\sigma_1 = 45^{+25}_{-15}$ km/s, the remaining natal kicks are drawn from the high-velocity component, same as in previous works: $\sigma_2 = 336$ km/s. We found that currently used models for natal kicks of neutron stars produced by electron capture supernova (combination of maxwellian $\sigma=265$ km/s and $\sigma = 30$ km/s for electron capture) are inconsistent or marginally consistent with parallaxes and proper motions measured for isolated radio pulsars. We suggest a new model for natal kicks of ecSN, which satisfy both observations of isolated radio pulsars and Be X-ray binaries., Comment: 21 page, 36 figures, accepted to MNRAS on September 21

We search for merger products among the 25 most massive white dwarfs in the Montreal White Dwarf Database 100 pc sample through follow-up spectroscopy and high-cadence photometry. We find an unusually high fraction, 40%, of magnetic white dwarfs among this population. In addition, we identify four outliers in transverse velocity and detect rapid rotation in five objects. Our results show that $56^{+9}_{-10}$\% of the $M\approx1.3~M_{\odot}$ ultramassive white dwarfs form through mergers. This fraction is significantly higher than expected from the default binary population synthesis calculations using the $\alpha$-prescription (with $\alpha \lambda = 2$), and provides further support for efficient orbital shrinkage, such as with low values of the common envelope efficiency., Comment: MNRAS, in press

In this paper, we consider triple systems composed of main-sequence (MS) stars, and their internal evolution due to stellar and binary evolution. Our focus is on triples that produce white dwarfs (WDs), where Roche lobe overflow of an evolving tertiary triggers accretion onto the inner binary via a circumbinary disk (CBD) driving it toward a mass ratio of unity. We present a combination of analytic- and population synthesis-based calculations performed using the \texttt{SeBa} code to constrain the expected frequency of such systems, given a realistic initial population of MS triples, and provide the predicted distributions of orbital periods. We identify the parameter space for triples that can accommodate a CBD, to inform future numerical simulations of suitable initial conditions. We find that $\lesssim$ 10\% of all MS triples should be able to accommodate a CBD around the inner binary, and compute lower limits for the production rates. This scenario broadly predicts mergers of near equal-mass binaries, producing blue stragglers (BSs), Type Ia supernovae, gamma ray bursts and gravitational wave-induced mergers, along with the presence of an outer WD tertiary companion. We compare our predicted distributions to a sample of field BS binaries, and argue that our proposed mechanism explains the observed range of orbital periods. Finally, the mechanism considered here could produce hypervelocity MS stars, WDs and even millisecond pulsars with masses close to the Chandrasekhar mass limit, and be used to constrain the maximum remnant masses at the time of any supernova explosion., Comment: 15 pages, 7 figures; accepted for publication in MNRAS

Recent observational evidence has demonstrated that white dwarf (WD) mergers are a highly efficient mechanism for mass accretion onto WDs in the galaxy. In this paper, we show that WD mergers naturally produce highly-magnetized, uniformly-rotating WDs, including a substantial population within a narrow mass range close to the Chandrasekhar mass ($M_{\rm Ch}$). These near-$M_{\rm Ch}$ WD mergers subsequently undergo rapid spin up and compression on a $\sim 10^2$ yr timescale, either leading to central ignition and a normal SN Ia via the DDT mechanism, or alternatively to a failed detonation and SN Iax through pure deflagration. The resulting SNe Ia and SNe Iax will have spectra, light curves, polarimetry, and nucleosynthetic yields similar to those predicted to arise through the canonical near-$M_{\rm Ch}$ single degenerate (SD) channel, but with a $t^{-1}$ delay time distribution characteristic of the double-degenerate (DD) channel. Furthermore, in contrast to the SD channel, WD merger near-$M_{\rm Ch}$ SNe Ia and SNe Iax will not produce observable companion signatures. We discuss a range of implications of these findings, from SNe Ia explosion mechanisms, to galactic nucleosynthesis of iron peak elements including manganese., Comment: 20 pages, 4 figures, nucleosynthetic yields included as ancillary data files. Accepted to The Astrophysical Journal. Comments welcome

The trajectory of the center of light of an unresolved binary is different from that of its center of mass. Binary-induced stellar centroid wobbling can therefore be detected as an excess in the goodness-of-fit of the single-star astrometric model. We use reduced $\chi^2$ of the astrometric fit in the {\it Gaia} Early Data Release 3 to detect the likely unresolved double white dwarfs (DWDs). Using parallax-based distances we convert the excess of reduced $\chi^2$ into the amplitude of the centroid wobble $\delta a$, which is proportional to the binary separation $a$. The measured $\delta a$ distribution drops towards larger wobble amplitudes and shows a break around $\delta a \approx 0.2$ where it steepens. The integral of the distribution yields DWD fraction of $6.5 \pm 3.7$ per cent in the range $0.01 < a (\text{au}) < 2$. Using synthetic models of the Galactic DWDs we demonstrate that the break in the $\delta a$ distribution corresponds to one side of a deep gap in the DWD separation distribution at around $a\approx 1$ au. Model DWDs with separations less than several au shrink dramatically due to (al least one) common envelope phase, reshaping the original separation distribution, clearing a gap and creating a pile-up of systems with $a\approx 0.01$ au and $\delta a < 0.01$. Our models reproduce the overall shape of the observed $\delta a$ distribution and its normalisation, however the predicted drop in the numbers of DWDs beyond the break is steeper than in the data., Comment: submitted to MNRAS main journal, comments are welcome

Barium stars have been studied extensively over the past few decades, yet our current understanding of how these intriguing objects formed leaves much to be desired. Many trends observed in systems containing barium stars cannot be satisfactorily explained by classical binary evolution models, naturally raising the question of whether triples and other higher-order multiples can give rise to such exotic objects. In this paper, we study the possibility that a Roche Lobe overflow from a tertiary in a hierarchical triple system can potentially lead to surface barium enrichment within the inner binary, while at the same time causing the inner binary to merge, thereby producing a barium star. This possibility has the potential to form a large proportion of Barium stars, as Roche Lobe overflow from a tertiary is typically much more stable for close orbits than that from a binary companion. Various formation channels and mechanisms by which this can be achieved are considered, and constraints on relative formation rates are placed on each scenario. We conclude that a significant, if not dominant, proportion of barium stars are formed from hierarchical triple systems, and that further studies are required in this area before a complete understanding of Barium star populations can be achieved., Comment: Submitted to MNRAS

Stellar-mass black holes are the final remnants of stars born with more than 15 solar masses. Billions are expected to reside in the Local Group, yet only a few are known, mostly detected through X-rays emitted as they accrete material from a companion star. Here, we report on VFTS 243: a massive X-ray-faint binary in the Large Magellanic Cloud. With an orbital period of 10.4 d, it comprises an O-type star of 25 solar masses and an unseen companion of at least nine solar masses. Our spectral analysis excludes a non-degenerate companion at a 5σ confidence level. The minimum companion mass implies that it is a black hole. No other X-ray-quiet black hole is unambiguously known outside our Galaxy. The (near-)circular orbit and kinematics of VFTS 243 imply that the collapse of the progenitor into a black hole was associated with little or no ejected material or black-hole kick. Identifying such unique binaries substantially impacts the predicted rates of gravitational-wave detections and properties of core-collapse supernovae across the cosmos., This research has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (T.S. and H.S., grant agreement 772225: MULTIPLES). T.S. acknowledges support from the European Union’s Horizon 2020 programme under Marie Skłodowska-Curie grant agreement 101024605. This work is based on observations collected at the ESO under programme IDs 182.D-0222, 090.D-0323 and 092.D-0136. L.M. thanks the ESA and the Belgian Federal Science Policy Office (BELSPO) for their support in the framework of the PRODEX programme. P. Marchant acknowledges support from FWO junior postdoctoral fellowship 12ZY520N. We thank J. Hillier for making the CMFGEN code available. C.H. acknowledges support from the KU Leuven Research Council (grant C16/17/007: MAESTRO). P.A.C. acknowledges support from UK Science and Technology Facilities Council research grant ST/V000853/1. V.H.-B. acknowledges the support of the Natural Sciences and Engineering Research Council of Canada (NSERC) through grant RGPIN-2020-05990. M. Gieles acknowledges support from the Ministry of Science and Innovation through a Europa Excelencia grant (EUR2020-112157). The PoWR code, as well as the associated post-processing and visualization tool WRplot, has been developed under the guidance of W.-R. Hamann with substantial contributions from L. Koesterke, G. Gräfener, A. Sander, T.S. and other co-workers and students. This work has received funding from the ERC under the European Union’s Horizon 2020 research and innovation programme (grant agreement 945806). It is also supported by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany’s Excellence Strategy EXC 2181/1-390900948 (the Heidelberg STRUCTURES Excellence Cluster). S.T. acknowledges support from the Netherlands Research Council NWO (grants VENI 639.041.645, VIDI 203.061). A.H. and D.J.L. acknowledge support by the Spanish MCI through grant PGC-2018-0913741-B-C22 and the Severo Ochoa Programme through CEX2019-000920-S. J.M.A. acknowledges support from the Spanish Government Ministerio de Ciencia, Innovación y Universidades through grant PGC2018-095-049-B-C22.

Dwarf carbon (dC) stars, main sequence stars showing carbon molecular bands, are enriched by mass transfer from a previous asymptotic-giant-branch (AGB) companion, which has since evolved to a white dwarf. While previous studies have found radial-velocity variations for large samples of dCs, there are still relatively few dC orbital periods in the literature and no dC eclipsing binaries have yet been found. Here, we analyze photometric light curves from DR5 of the Zwicky Transient Facility for a sample of 944 dC stars. From these light curves, we identify 34 periodically variable dC stars. Remarkably, of the periodic dCs, 82\% have periods less than two days. We also provide spectroscopic follow-up for four of these periodic systems, measuring radial velocity variations in three of them. Short-period dCs are almost certainly post-common-envelope binary systems, since the periodicity is most likely related to the orbital period, with tidally locked rotation and photometric modulation on the dC either from spots or from ellipsoidal variations. We discuss evolutionary scenarios that these binaries may have taken to accrete sufficient C-rich material while avoiding truncation of the thermally pulsing AGB phase needed to provide such material in the first place. We compare these dCs to common-envelope models to show that dC stars probably cannot accrete enough C-rich material during the common-envelope phase, suggesting another mechanism like wind-Roche lobe overflow is necessary. The periodic dCs in this paper represent a prime sample for spectroscopic follow-up and for comparison to future models of wind-Roche lobe overflow mass transfer., Comment: 23 pages, 5 figures, 5 tables, accepted to ApJ

As main-sequence stars with C$>$O, dwarf carbon (dC) stars are never born alone but inherit carbon-enriched material from a former asymptotic giant branch (AGB) companion. In contrast to M dwarfs in post-mass transfer binaries, C$_2$ and/or CN molecular bands allow dCs to be identified with modest-resolution optical spectroscopy, even after the AGB remnant has cooled beyond detectability. Accretion of substantial material from the AGB stars should spin up the dCs, potentially causing a rejuvenation of activity detectable in X-rays. Indeed, a few dozen dCs have recently been found to have photometric variability with periods under a day. However, most of those are likely post-common-envelope binaries (PCEBs), spin-orbit locked by tidal forces, rather than solely spun-up by accretion. Here, we study the X-ray properties of a sample of the five nearest known dCs with $Chandra$. Two are detected in X-rays, the only two for which we also detected short-period photometric variability. We suggest that the coronal activity detected so far in dCs is attributable to rapid rotation due to tidal locking in short binary orbits after a common-envelope phase, late in the thermally pulsing (TP) phase of the former C-AGB primary (TP-AGB)., 18 pages, 9 figures, 4 tables, Accepted to ApJ, title changed based on referee recommendation

Many white dwarf stars show signs of having accreted smaller bodies, implying that they may host planetary systems. A small number of these systems contain gaseous debris discs, visible through emission lines. We report a stable 123.4min periodic variation in the strength and shape of the CaII emission line profiles originating from the debris disc around the white dwarf SDSSJ122859.93+104032.9. We interpret this short-period signal as the signature of a solid body held together by its internal strength., Comment: 11 Figures, 3 Tables

Double white dwarf binaries with merger timescales smaller than the Hubble time and with a total mass near the Chandrasekhar limit (i.e. classical Chandrasekhar population) or with high-mass primaries (i.e. sub-Chandrasekhar population) are potential supernova type Ia (SNIa) progenitors. However, we have not yet unambiguously confirmed the existence of these objects observationally, a fact that has been often used to criticise the relevance of double white dwarfs for producing SNIa. We analyse whether this lack of detections is due to observational effects. To that end we simulate the double white dwarf binary population in the Galaxy and obtain synthetic spectra for the SNIa progenitors. We demonstrate that their identification, based on the detection of Halpha double-lined profiles arising from the two white dwarfs in the synthetic spectra, is extremely challenging due to their intrinsic faintness. This translates into an observational probability of finding double white dwarf SNIa progenitors in the Galaxy of (2.1+-1.0)x10^{-5} and (0.8+-0.4)x10^{-5} for the classical Chandrasekhar and the sub-Chandrasekhar progenitor populations, respectively. Eclipsing double white dwarf SNIa progenitors are found to suffer from the same observational effect. The next generation of large-aperture telescopes are expected to help in increasing the probability for detection by ~1 order of magnitude. However, it is only with forthcoming observations such as those provided by LISA that we expect to unambiguously confirm or disprove the existence of double white dwarf SNIa progenitors and to test their importance for producing SNIa., Comment: Re-submitted to MNRAS after applying all referee's suggestions

The population of Milky Way satellite galaxies is of great interest for cosmology, fundamental physics, and astrophysics. They represent the faint end of the galaxy luminosity function, are the most dark-matter dominated objects in the local Universe, and contain the oldest and most metal-poor stellar populations. Recent surveys have revealed around 60 satellites, but this could represent less than half of the total. Characterization of these systems remains a challenge due to their low luminosity. We consider the gravitational wave observatory LISA as a potential tool for studying these satellites through observations of their short-period double white dwarf populations. LISA will observe the entire sky without selection effects due to dust extinction, complementing optical surveys, and could potentially discover massive satellites hidden behind the disk of the galaxy., 7 pages, 2 figures

Normal type Ia supernovae (SNe) are thought to arise from the thermonuclear explosion of massive ($>0.8$ M$_\odot$) carbon-oxygen white dwarfs (WDs), although the exact mechanism is debated. In some models helium accretion onto a carbon-oxygen (CO) WD from a companion was suggested to dynamically trigger a detonation of the accreted helium shell. The helium detonation then produces a shock that after converging on itself close to the core of the CO-WD, triggers a secondary carbon detonation and gives rise to an energetic explosion. However, most studies of such scenarios have been done in one or two dimensions, and/or did not consider self-consistent models for the accretion and the He-donor. Here we make use of detailed 3D simulation to study the interaction of a He-rich hybrid $0.69\,\mathrm{M_\odot}$ HeCO WD with a more massive $0.8\,\mathrm{M_\odot}$ CO~WD. We find that accretion from the hybrid WD onto the CO~WD gives rise to a helium detonation. However, the helium detonation does not trigger a carbon detonation in the CO~WD. Instead, the helium detonation burns through the accretion stream to also burn the helium shell of the donor hybrid HeCO-WD. The detonation of its massive helium shell then compresses its CO core, and triggers its detonation and full destruction. The explosion gives rise to a faint, likely highly reddened transient, potentially observable by the Vera Rubin survey, and the high-velocity ($\sim 1000\,\mathrm{km s^{-1}}$) ejection of the heated surviving CO~WD companion. Pending on uncertainties in stellar evolution we estimate the rate of such transient to be up to $\sim10\%$ of the rate of type Ia SNe., 14 pages, 10 figures, accepted by MNRAS, comments welcome

We present a spectroscopic survey of 230 white dwarf candidates within 40 pc of the Sun from the William Herschel Telescope and Gran Telescopio Canarias. All candidates were selected from Gaia Data Release 2 (DR2) and in almost all cases, had no prior spectroscopic classifications. We find a total of 191 confirmed white dwarfs and 39 main-sequence star contaminants. The majority of stellar remnants in the sample are relatively cool (99 per cent) has now been reached for Gaia DR2 sources within 40-pc sample, in the Northern hemisphere (δ > 0) and located on the white dwarf cooling track in the Hertzsprung-Russell diagram. A statistical study of the full northern sample is presented in a companion paper., The research leading to these results has received funding from the European Research Council under the European Union’s Horizon 2020 research and innovation programme no. 677706 (WD3D). This article is based on observations made in the Observatorios de Canarias del IAC with the WHT operated on the island of La Palma by the Isaac Newton Group of Telescopes in the Observatorio del Roque de los Muchachos. Based on observations made with the Gran Telescopio Canarias (GTC), installed at the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofísica de Canarias, in the island of La Palma. This work presents results from the European Space Agency (ESA) space mission Gaia. Gaia data are being processed by the Gaia Data Processing and Analysis Consortium (DPAC). Funding for the DPAC is provided by national institutions, in particular the institutions participating in the Gaia MultiLateral Agreement (MLA). BTG was supported by the UK STFC grants ST/P000495 and ST/T000406/1. ARM acknowledges support from the MINECO under the Ramón y Cajal programme (RYC-2016-20254), the AYA2017-86274-P grant, and the AGAUR grant SGR-661/2017. MRS thanks for support from Fondecyt (grant 1181404). DdM acknowledges financial support from 〈0:funding-source 3:href="http://dx.doi.org/10.13039/501100 003981"〉ASI〈/0:funding-source〉-INAF I/037/12/0 and ASI-INAF no. 2017-14-H.0 and from INAF “INAF main streams”, Presidential Decree 43/2018. RR has received funding from the postdoctoral fellowship programme Beatriu de Pinós, funded by the Secretary of Universities and Research (Government of Catalonia) and by the Horizon 2020 programme of research and innovation of the European Union under the Maria Skłodowska-Curie grant agreement no. 801370.

We present an overview of the sample of northern hemisphere white dwarfs within 40 pc of the Sun detected from $Gaia$ Data Release 2 (DR2). We find that 521 sources are spectroscopically confirmed degenerate stars, 111 of which were first identified as white dwarf candidates from $Gaia$ DR2 and followed-up recently with the William Herschel Telescope and Gran Telescopio Canarias. Three additional white dwarf candidates remain spectroscopically unobserved and six unresolved binaries are known to include a white dwarf but were not in our initial selection of white dwarfs in the $Gaia$ DR2 Hertzsprung-Russell diagram (HRD). Atmospheric parameters are calculated from $Gaia$ and Pan-STARRS photometry for all objects in the sample, confirming most of the trends previously observed in the much smaller 20 pc sample. Local white dwarfs are overwhelmingly consistent with Galactic disc kinematics, with only four halo candidates. We find that DAZ white dwarfs are significantly less massive than the overall DA population ($\overline{M}_\mathrm{DAZ} = 0.59\,\mathrm{M}_\odot$, $\overline{M}_\mathrm{DA} = 0.66\,\mathrm{M}_\odot$). It may suggest that planet formation is less efficient at higher mass stars, producing more massive white dwarfs. We detect a sequence of crystallised white dwarfs in the mass range from $0.6\ \lesssim M/\mathrm{M}_\odot \lesssim\ 1.0\,$ and find that the vast majority of objects on the sequence have standard kinematic properties that correspond to the average of the sample, suggesting that their nature can be explained by crystallisation alone. We also detect 56 wide binaries including a white dwarf and 26 double degenerates., 20 pages, 14 figures. Accepted for publication in MNRAS. Online tables available at http://deneb.astro.warwick.ac.uk/phrgwr/40pcTables/index.html

Context. Milky Way dwarf satellites are unique objects that encode the early structure formation and therefore represent a window into the high redshift Universe. So far, their study has been conducted using electromagnetic waves only. The future Laser Interferometer Space Antenna (LISA) has the potential to reveal Milky Way satellites through gravitational waves emitted by double white dwarf (DWD) binaries. Aims. We investigate gravitational wave signals that will be detectable by LISA as a possible tool for the identification and characterisation of the Milky Way satellites. Methods. We used the binary population synthesis technique to model the population of DWDs in dwarf satellites and we assessed the impact on the number of LISA detections when making changes to the total stellar mass, distance, star formation history, and metallicity of satellites. We calibrated predictions for the known Milky Way satellites on their observed properties. Results. We find that DWDs emitting at frequencies ≳3 mHz can be detected in Milky Way satellites at large galactocentric distances. The number of these high frequency DWDs per satellite primarily depends on its mass, distance, age, and star formation history, and only mildly depends on the other assumptions regarding their evolution such as metallicity. We find that dwarf galaxies with M⋆ > 106 M⊙ can host detectable LISA sources; the number of detections scales linearly with the satellite’s mass. We forecast that out of the known satellites, Sagittarius, Fornax, Sculptor, and the Magellanic Clouds can be detected with LISA. Conclusions. As an all-sky survey that does not suffer from contamination and dust extinction, LISA will provide observations of the Milky Way and dwarf satellites galaxies, which will be valuable for Galactic archaeology and near-field cosmology.

HR 6819 was recently claimed to be a hierarchical triple system of a Be star in a wide orbit around an inner binary system of a black hole (BH) and a B III type star. We argue that this system is unlikely to be a hierarchical triple due to three reasons: (i) Given that this system is discovered in a magnitude limited Bright Star Catalog, the expected number of such systems in the Milky Way amounts to about $10^4$ while the estimate for the MW budget for such systems is between $10^2-10^3$ systems under generous assumptions. Such a large gap cannot be reconciled as it would otherwise likely overflow the MW budget for BHs; (ii) The dynamical stability of this system sets lower bounds on the orbital separation of the outer Be star, while it not being resolved by Gaia places an upper limit on its projected sky separation. We show that these two constraints would imply a narrow range for the outer orbit without resorting to geometrical fine-tuning; (iii) The triple system should have survived the stellar evolution prior to the formation of the BH in the inner binary. We perform numerical simulations starting with conservative initial conditions of this system and show that a small parameter space for BH progenitor star's mass loss, BH natal kicks, and initial orbital separation can reproduce HR 6819. Therefore, we propose this system is a chance superposition of a Be star with a binary., Matches the published version at ApJ Letters

Many stars do not live alone, but instead have one or more stellar companions. Observations show that these binaries, triples and higher-order multiples are common. Whereas the evolution of single stars and binaries have been studied extensively, the same is not true for the evolution of stellar triples. To fill this gap in our general understanding of stellar lives, we aim to systematically explore the long-term evolution of triples and to map out the most common evolutionary pathways that triples go through. We quantitatively study how triples evolve, which processes are most relevant, and how this differs from binary evolution. We simulate the evolution of several large populations of triples with a population synthesis approach. We make use of the triple evolution code TRES to simulate the evolution of each triple in a consistent way; including three-body dynamics (based on the secular approach), stellar evolution and their mutual influences. We simulate the evolution of the system up until mass transfer starts, the system becomes dynamically unstable, or a Hubble time has passed. We find that stellar interactions are common in triples. Compared to a binary population, we find that the fraction of systems that can undergo mass transfer is about 2 to 3 times larger in triples. Moreover, whereas in binaries the orbits typically reach circularisation before Roche-lobe overflow, this is not true anymore in triples. In our simulations, about 40% of systems retain an eccentric orbit. Additionally, we discuss various channels of triple evolution in detail such as those where the secondary or the tertiary is the first star to initiate a mass transfer event., Comment: updated version: accepted for publication in A&A 18 pages, 16 figures, 2 tables

Stellar triples with massive stellar components are common, and can lead to sequential binary black-hole mergers. Here, we outline the evolution towards these sequential mergers, and explore these events in the context of gravitational-wave astronomy and the pair-instability mass gap. We find that binary black-hole mergers in the pair-instability mass gap can be of triple origin and therefore are not exclusively formed in dense dynamical environments. We discuss the sequential merger scenario in the context of the most massive gravitational-wave sources detected to date: GW170729 and GW190521. We propose that the progenitor of GW170729 is a low-metallicity field triple. We support the premise that GW190521 could not have been formed in the field. We conclude that triple stellar evolution is fundamental in the understanding of gravitational-wave sources, and likely, other energetic transientsas well., Comment: 13 main pages, 4 figures, 1 table plus Appendices. Accepted for publication in ApJ Letters

Using Gaia DR2 data, we present an up-to-date sample of white dwarfs within 20 pc of the Sun. In total we identified 139 systems in Gaia DR2, nine of which are new detections, with the closest of these located at a distance of 13.05 pc. We estimated atmospheric parameters for all stellar remnants based on the Gaia parallaxes and photometry. The high precision and completeness of the Gaia astrometry allowed us to search for wide binary companions. We re-identified all known binaries where both components have accurate DR2 astrometry, and established the binarity of one of the nine newly identified white dwarfs. No new companions were found to previously known 20 pc white dwarfs. Finally, we estimated the local white dwarf space-density to be $(4.49\pm0.38)\times10^{-3}$ pc$^{-3}$, having given careful consideration to the distance-dependent Gaia completeness, which misses known objects at short distances, but is close to complete for white dwarfs near 20 pc., Accepted for publication in MNRAS

Hydrogen-rich supernovae, known as Type II (SNe II), are the most common class of explosions observed following the collapse of the core of massive stars. We use analytical estimates and population synthesis simulations to assess the fraction of SNe II progenitors that are expected to have exchanged mass with a companion prior to explosion. We estimate that 1/3 to 1/2 of SN II progenitors have a history of mass exchange with a binary companion before exploding. The dominant binary channels leading to SN II progenitors involve the merger of binary stars. Mergers are expected to produce a diversity of SN II progenitor characteristics, depending on the evolutionary timing and properties of the merger. Alternatively, SN II progenitors from interacting binaries may have accreted mass from their companion, and subsequently been ejected from the binary system after their companion exploded. We show that the overall fraction of SN II progenitors that are predicted to have experienced binary interaction is robust against the main physical uncertainties in our models. However, the relative importance of different binary evolutionary channels is affected by changing physical assumptions. We further discuss ways in which binarity might contribute to the observed diversity of SNe II by considering potential observational signatures arising from each binary channel. For supernovae which have a substantial H-rich envelope at explosion (i.e., excluding Type IIb SNe), a surviving non-compact companion would typically indicate that the supernova progenitor star was in a wide, non-interacting binary. We argue that a significant fraction of even Type II-P SNe are expected to have gained mass from a companion prior to explosion., The paper has been accepted for publication in Astronomy and Astrophysics

Double white dwarf (double-WD) binaries may merge within a Hubble time and produce high-mass WDs. Compared to other high-mass WDs, the double-WD merger products have higher velocity dispersion because they are older. With the power of Gaia data, we show strong evidence for double-WD merger products among high-mass WDs by analyzing the transverse-velocity distribution of more than a thousand high-mass WDs (0.8--1.3 $M_\odot$). We estimate that the fraction of double-WD merger products in our sample is about 20 %. We also obtain a precise double-WD merger rate and its mass dependence. Our merger rate estimates are close to binary population synthesis results and support the idea that double-WD mergers may contribute to a significant fraction of type Ia supernovae., Comment: 9 pages, 5 figures; published in ApJ

Context: White dwarfs (WDs) are important and abundant tools to study the structure and evolution of the Galactic environment. However, the multiplicity of WD progenitors is generally neglected. Specifically, a merger in a binary system can lead to a single WD, which could result in wrongly inferred quantities if only single stellar evolution (SSE) is considered. These mergers are linked to transients such as luminous red novae and Type Ia supernovae. Aims: We investigate the impact of binary evolution (BE) upon observable single WDs, and compare their properties to WDs formed through SSE. We assess the evolutionary channels and the age and mass distributions of the resulting single Wds. Methods: We employed \texttt{SeBa} to model the evolution of single star and binary populations. We synthesised the observable single WD population within $100$ pc, including cooling and observational selection effects. Additionally, we constructed models with different evolution and primordial population properties to study the effects on the properties of the resulting single WDs. Results: White dwarfs from binary mergers make up about $10-30\%$ of all observable single WDs and $30-50\%$ of massive WDs. On average, individual WDs take $3.1-5$ times longer to form through BE than SE, and so appear ${\sim} 1$ Gyr younger than they are if BE is ignored. In all models, the effect of mergers on the age distribution is clearly noticeable. The median age typically increases by $85-430$ Myr and $200-390$ Myr for massive WDs. Although abundant, we do not find evidence that WDs from mergers significantly alter the shape of the WD mass distribution. Conclusions: Assuming SSE for inferring properties of single WDs gives rise to intrinsic errors as single WDs can also be formed following a binary merger. Strategies for mitigating the effect of mergers on the WD age distributions are discussed., 16 pages, 14 figures. Published in A&A

Double detonations in double white dwarf (WD) binaries undergoing unstable mass transfer have emerged in recent years as one of the most promising Type Ia supernova (SN Ia) progenitor scenarios. One potential outcome of this "dynamically driven double-degenerate double-detonation" (D^6) scenario is that the companion WD survives the explosion and is flung away with a velocity equal to its > 1000 km/s pre-SN orbital velocity. We perform a search for these hypervelocity runaway WDs using Gaia's second data release. In this paper, we discuss seven candidates followed up with ground-based instruments. Three sources are likely to be some of the fastest known stars in the Milky Way, with total Galactocentric velocities between 1000 and 3000 km/s, and are consistent with having previously been companion WDs in pre-SN Ia systems. However, although the radial velocity of one of the stars is > 1000 km/s, the radial velocities of the other two stars are puzzlingly consistent with 0. The combined five-parameter astrometric solutions from Gaia and radial velocities from follow-up spectra yield tentative 6D confirmation of the D^6 scenario. The past position of one of these stars places it within a faint, old SN remnant, further strengthening the interpretation of these candidates as hypervelocity runaways from binary systems that underwent SNe Ia., Accepted for publication in ApJ. Minor corrections for clarity. D6 spectra are available as ancillary data files

We present a catalogue of white dwarf candidates selected from the second data release of $Gaia$ (DR2). We used a sample of spectroscopically confirmed white dwarfs from the Sloan Digital Sky Survey (SDSS) to map the entire space spanned by these objects in the $Gaia$ Hertzsprung-Russell diagram. We then defined a set of cuts in absolute magnitude, colour, and a number of $Gaia$ quality flags to remove the majority of contaminating objects. Finally, we adopt a method analogous to the one presented in our earlier SDSS photometric catalogues to calculate a probability of being a white dwarf ($P_{\mathrm{WD}}$) for all $Gaia$ sources which passed the initial selection. The final catalogue is composed of $486\,641$ stars with calculated $P_{\mathrm{WD}}$ from which it is possible to select a sample of $\simeq 260\,000$ high-confidence white dwarf candidates in the magnitude range $8 7000$ K, at high Galactic latitudes ($|b|>20^{\circ}$). However, the completeness drops at low Galactic latitudes, and the magnitude limit of the catalogue varies significantly across the sky as a function of $Gaia$'s scanning law. We also provide the list of objects within our sample with available SDSS spectroscopy. We use this spectroscopic sample to characterise the observed structure of the white dwarf distribution in the H-R diagram., Comment: 23 pages, 21 figures, 4 tables. Accepted for publication in MNRAS. Version 3 accepted. Gaia DR2 catalogue of white dwarf candidates available . SDSS spectra sample available

White dwarfs (WDs) are the stellar core remnants of low mass stars. They are typically divided into three main composition groups: Oxygen Neon (ONe), Carbon Oxygen (CO) and Helium (He) WDs. The evolution of binary systems can significantly change the evolution of the binary stellar components. In particular, striping the envelope of an evolved star can give rise to a core remnant, which can later evolve into a WD with significantly different composition. Here we focus on the formation and evolution of hybrid HeCO WDs. We follow the formation and stellar evolution of such WDs for a range of initial conditions and provide their detailed structure, mass-radius relation and luminosity temperature evolution. We find that both low mass WDs (< 0.45M , typically thought to be He WDs) and intermediate-mass WDs (0.45 < MWD < 0.7, typically thought to be CO WDs) could in fact be hybrid HeCO WDs, with 5-25 (75 -95)% of their mass in He (CO). We use population synthesis calculations to infer the birth rate and properties of such WDs. We find that hybrid HeCO WD comprise the majority of young (< 2Gyr) WDs in binaries, but are more rare among older WDs in binaries. The high frequency and large He content of such WDs could have an important role in WD WD mergers, and may give rise to sub Chandrasekhar thermonuclear supernova explosions., 9 pages, 12 figures,, 1 table, Submitted to MNRAS

White dwarfs (WDs) are powerful tools to study the evolutionary history of stars and binaries in the Galaxy. But do we understand their multiplicity from a theoretical point of view? This can be tested by a comparison with the sample of WDs within 20 pc, which is minimally affected by selection biases. From the literature, we compile the available information of the local WD sample with a particular emphasis on their multiplicity, and compare this to synthetic models of WD formation in single stars and binaries. As part of our population synthesis approach, we also study the effect of different assumptions concerning the star formation history, binary evolution, and the initial distributions of binary parameters. We find that the observed space densities of single and binary WDs are well reproduced by the models. The space densities of the most common WD systems (single WDs and unresolved WD-MS binaries) are consistent within a factor two with the observed value. We find a discrepancy only for the space density of resolved double WDs. We exclude that observational selection effects, fast stellar winds, or dynamical interactions with other objects in the Milky Way explain this discrepancy. We find that either the initial mass ratio distribution in the Solar neighbourhood is biased towards low mass-ratios, or more than ten resolved DWDs have been missed observationally in the 20 pc sample. Furthermore, we show that the low binary fraction of WD systems (~25%) compared to Solartype MS-MS binaries (~50%) is consistent with theory, and is mainly caused by mergers in binary systems, and to a lesser degree by WDs hiding in the glare of their companion stars. Lastly, Gaia will dramatically increase the size of the volume-limited WD sample, detecting the coolest and oldest WDs out to 50 pc. We provide a detailed estimate of the number of single and binary WDs in the Gaia sample., 8 figures,published in A&A, typos & affiliations updated

White dwarfs (WDs) can increase their mass by accretion from companion stars, provided the mass-accretion rate is high enough to avoid nova eruptions. The accretion regimes that allow growth of the WDs are usually calculated assuming constant mass-transfer rates. However, it is possible that these systems are influenced by effects that cause the rate to fluctuate on various timescales. We investigate how long-term mass-transfer variability affects accreting WDs systems. We show that, if such variability is present, it expands the parameter space of binaries where the WD can effectively increase its mass. Furthermore, we find that the supernova type Ia (SNIa) rate is enhanced by a factor 2-2.5 to a rate that is comparable with the lower limit of the observed rates. The changes in the delay-time distribution allow for more SNIae in stellar populations with ages of a few Gyr. Thus, mass-transfer variability gives rise to a new formation channel of SNIa events that can significantly contribute to the SNIa rate. Mass-transfer variability is also likely to affect other binary populations through enhanced WD growth. For example, it may explain why WDs in cataclysmic variables are observed to be more massive than single WDs, on average., 10 pages, 6 figures, accepted for publication in MNRAS

Binary mass transfer is at the forefront of some of the most exciting puzzles of modern astrophysics, including Type Ia supernovae, gamma-ray bursts, and the formation of most observed exotic stellar populations. Typically, the evolution is assumed to proceed in isolation, even in dense stellar environments such as star clusters. In this paper, we test the validity of this assumption via the analysis of a large grid of binary evolution models simulated with the SeBa code. For every binary, we calculate analytically the mean time until another single or binary star comes within the mean separation of the mass-transferring binary, and compare this time-scale to the mean time for stable mass transfer to occur. We then derive the probability for each respective binary to experience a direct dynamical interruption. The resulting probability distribution can be integrated to give an estimate for the fraction of binaries undergoing mass transfer that are expected to be disrupted as a function of the host cluster properties. We find that for lower-mass clusters ($\lesssim 10^4$ M$_{\odot}$), on the order of a few to a few tens of percent of binaries undergoing mass-transfer are expected to be interrupted by an interloping single, or more often binary, star, over the course of the cluster lifetime, whereas in more massive globular clusters we expect $\ll$ 1% to be interrupted. Furthermore, using numerical scattering experiments performed with the FEWBODY code, we show that the probability of interruption increases if perturbative fly-bys are considered as well, by a factor $\sim 2$., 6 pages, 4 figures, 1 table; accepted for publication in ApJ

The theoretical and observed populations of pre-cataclysmic variables (pre-CVs) are dominated by systems with low-mass white dwarfs (WDs), while the WD masses in CVs are typically high. In addition, the space density of CVs is found to be significantly lower than theoretical models. We investigate the influence of nova outbursts on the formation and (initial) evolution of CVs. In particular, we calculate the stability of the mass transfer in case all the material accreted on the WD is lost in classical novae, and part of the energy to eject the material comes from a common-envelope like interaction with the companion. In addition, we study the effect of an asymmetry in the mass ejection, that may lead to small eccentricities in the orbit. We find that a common-envelope like ejection significantly decreases the stability of the mass transfer, in particular for low-mass WD. Similarly, the influence of asymmetric mass loss can be important for short-period systems and even more so for low-mass WD, but likely disappears long before the next nova outburst due to orbital circularization. In both cases the mass-transfer rates increase, which may lead to observable (and perhaps already observed) consequences for systems that do survive to become CVs. However, a more detailed investigation of the interaction between nova ejecta and the companion and the evolution of slightly eccentric CVs is needed before definite conclusions can be drawn., Accepted for publication in ApJ

Aims. We model the present-day population of 'classical' low-mass X-ray binaries (LMXBs) with neutron star accretors, which have hydrogen-rich donor stars. Their population is compared with that of hydrogen-deficient LMXBs, known as ultracompact X-ray binaries (UCXBs). We model the observable LMXB population and compare it to observations. Methods. We combine the binary population synthesis code SeBa with detailed LMXB evolutionary tracks to model the size and properties of the present-day LMXB population in the Galactic Bulge. Whether sources are persistent or transient, and what their instantaneous X-ray luminosities are, is predicted using the thermal-viscous disk instability model. Results. We find a population of ~2.1 x 10^3 LMXBs with neutron star accretors. Of these about 15 - 40 are expected to be persistent (depending on model assumptions), with luminosities higher than 10^35 erg s^-1. About 7 - 20 transient sources are expected to be in outburst at any given time. Within a factor of two these numbers are consistent with the observed population of bright LMXBs in the Bulge. This gives credence to our prediction of the existence of a population of ~1.6 x 10^3 LMXBs with low donor masses that have gone through the period minimum, and have present-day mass transfer rates below 10^-11 Msun yr^-1. Conclusions. Even though the observed population of hydrogen-rich LMXBs in the Bulge is larger than the observed population of (hydrogen-deficient) UCXBs, the latter have a higher formation rate. While UCXBs may dominate the total LMXB population at the present, the majority would be very faint, or may have become detached and produced millisecond radio pulsars. In that case UCXBs would contribute significantly more to the formation of millisecond radio pulsars than hydrogen-rich LMXBs. [abridged], Comment: 8 pages, 10 figures. Accepted for publication in Astronomy and Astrophysics. v2: minor language corrections

Binary population synthesis (BPS) modelling is a very effective tool to study the evolution and properties of close binary systems. The uncertainty in the parameters of the model and their effect on a population can be tested in a statistical way, which then leads to a deeper understanding of the underlying physical processes involved. To understand the predictive power of BPS codes, we study the similarities and differences in the predicted populations of four different BPS codes for low- and intermediate-mass binaries. We investigate whether the differences are caused by different assumptions made in the BPS codes or by numerical effects. To simplify the complex problem of comparing BPS codes, we equalise the inherent assumptions as much as possible. We find that the simulated populations are similar between the codes. Regarding the population of binaries with one WD, there is very good agreement between the physical characteristics, the evolutionary channels that lead to the birth of these systems, and their birthrates. Regarding the double WD population, there is a good agreement on which evolutionary channels exist to create double WDs and a rough agreement on the characteristics of the double WD population. Regarding which progenitor systems lead to a single and double WD system and which systems do not, the four codes agree well. Most importantly, we find that for these two populations, the differences in the predictions from the four codes are not due to numerical differences, but because of different inherent assumptions. We identify critical assumptions for BPS studies that need to be studied in more detail., Comment: 13 pages, +21 pages appendix, 35 figures, accepted for publishing in A&A, Minor change to match published version, most important the added link to the website http://www.astro.ru.nl/~silviato/popcorn for more detailed figures and information

Context: The nature of the progenitors of type Ia supernova (SN) progenitors remains unclear. While it is usually agreed that single-degenerate progenitor systems would be luminous supersoft X-ray sources (SSSs), it was recently suggested that double-degenerate progenitors might also go through a SSS phase. Aims: We aim to examine the possibility of double-degenerate progenitor systems being SSSs, and place stringent upper limits on the maximally possible durations of any SSS phases and expected number of these systems in a galactic population. Method: We employ the binary population synthesis code SeBa to examine the mass-transfer characteristics of a possible SSS phase of double-degenerate type Ia SN progenitor systems for 1) the standard SeBa assumptions, and 2) an optimistic best-case scenario. The latter case establishes firm upper limits on the possible population of supersoft source double-degenerate type Ia SN progenitor systems. Results: Our results indicate that unlike what is expected for single-degenerate progenitor systems, the vast majority of the material accreted by either pure wind mass transfer or a combination of wind and RLOF mass transfer is helium rather than hydrogen. Even with extremely optimistic assumptions concerning the mass-transfer and retention efficiencies, the average mass accreted by systems that eventually become double-degenerate type Ia SNe is small. Consequently, the lengths of time that these systems may be SSSs are short, even under optimal conditions, and the expected number of such systems in a galactic population is negligible. Conclusions: The population of double-degenerate type Ia SN progenitors that are SSSs is at least an order of magnitude smaller than the population of single-degenerate progenitors expected to be SSSs, and the supersoft X-ray behaviour of double-degenerate systems typically ceases long before the SN explosions., Comment: 10 pages, 3 figures, A&A accepted

[abridged] Aims. We model the number and properties of ultracompact X-ray binaries (UCXBs) in the Galactic Bulge. The objective is to learn about the formation of UCXBs and their evolution, such as the onset of mass transfer and late-time behavior. Methods. The binary population synthesis code SeBa and stellar evolutionary tracks are used to model the UCXB population. The luminosity behavior of UCXBs is predicted using long-term X-ray observations of the known UCXBs and the thermal-viscous disk instability model. Results. In our model, the majority of UCXBs initially have a helium burning star donor. In the absence of a mechanism that destroys old UCXBs, we predict (0.2 - 1.9) x 10^5 UCXBs in the Galactic Bulge, mostly at orbital periods longer than 60 min (a large number of long-period systems also follows from the observed short-period UCXB population). About 5 - 50 UCXBs should be brighter than 10^35 erg/s, mostly persistent sources with orbital periods shorter than 30 min and with degenerate helium and carbon-oxygen donors. This is one order of magnitude more than the observed number of (probably) three. Conclusions. This overprediction of short-period UCXBs implies that fewer systems are formed, or that a super-Eddington mass transfer rate is more difficult to survive than assumed. The very small number of observed long-period UCXBs with respect to short-period UCXBs, the surprisingly high luminosity of the observed UCXBs with orbital periods around 50 min, and the properties of the PSR J1719-1438 system all point to much faster UCXB evolution than expected from angular momentum loss via gravitational wave radiation alone. Old UCXBs, if they still exist, probably have orbital periods longer than 2 h and have become very faint due to either reduced accretion or quiescence, or have become detached. UCXBs are promising candidate progenitors of isolated millisecond radio pulsars., 16 pages, 14 figures, 1 table. Accepted for publication in Astronomy and Astrophysics

The study of Type Ia supernovae (SNIa) has lead to greatly improved insights into many fields in astrophysics, however a theoretical explanation of the origin of these events is still lacking. We investigate the potential contribution to the SNIa rate from the population of merging double carbon-oxygen white dwarfs. We aim to develope a model that fits the observed SNIa progenitors as well as the observed close double white dwarf population. We differentiate between two scenarios for the common envelope (CE) evolution; the alpha-formalism based on the energy equation and the gamma-formalism that is based on the angular momentum equation. In one model we apply the alpha-formalism always. In the second model the gamma-formalism is applied, unless the binary contains a compact object or the CE is triggered by a tidal instability for which the alpha-formalism is used. The binary population synthesis code SeBa was used to evolve binary systems from the zero-age main sequence to the formation of double white dwarfs and subsequent mergers. SeBa has been thoroughly updated since the last publication of the content of the code. The limited sample of observed double white dwarfs is better represented by the simulated population using the gamma-formalism than the alpha-formalism. For both CE formalisms, we find that although the morphology of the simulated delay time distribution matches that of the observations within the errors, the normalisation and time-integrated rate per stellar mass are a factor 7-12 lower than observed. Furthermore, the characteristics of the simulated populations of merging double carbon-oxygen white dwarfs are discussed and put in the context of alternative SNIa models for merging double white dwarfs., 16 pages (including 4 pages appendix), 15 figures

An update is presented to the software package SeBa for simulating single star and binary evolution in which new stellar evolution tracks have been implemented. SeBa is applied to study the population of close double white dwarf and the delay time distribution of double white dwarf mergers that may lead to Supernovae Type Ia., 4 pages, 4 figures, appeared in proceedings for "The 17th European White Dwarf Workshop", shorter version appeared in proceedings for "Binary Star Evolution: Mass Loss, Accretion and Mergers"

Field stars are frequently formed in pairs, and many of these binaries are part of triples or even higher-order systems. Even though, the principles of single stellar evolution and binary evolution, have been accepted for a long time, the long-term evolution of stellar triples is poorly understood. The presence of a third star in an orbit around a binary system can significantly alter the evolution of those stars and the binary system. The rich dynamical behavior in three-body systems can give rise to Lidov-Kozai cycles, in which the eccentricity of the inner orbit and the inclination between the inner and outer orbit vary periodically. In turn, this can lead to an enhancement of tidal effects (tidal friction), gravitational-wave emission and stellar interactions such as mass transfer and collisions. The lack of a self-consistent treatment of triple evolution, including both three-body dynamics as well as stellar evolution, hinders the systematic study and general understanding of the long-term evolution of triple systems. In this paper, we aim to address some of these hiatus, by discussing the dominant physical processes of hierarchical triple evolution, and presenting heuristic recipes for these processes. To improve our understanding on hierarchical stellar triples, these descriptions are implemented in a public source code TrES which combines three-body dynamics (based on the secular approach) with stellar evolution and their mutual influences. Note that modeling through a phase of stable mass transfer in an eccentric orbit is currently not implemented in TrES , but can be implemented with the appropriate methodology at a later stage., Comment: 40 pages, 14 figures, accepted for publication