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.

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.

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.

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.

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 $\textit{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 ($\langle T_{\rm eff} \rangle$ = 6200 K), showing either hydrogen Balmer lines or a featureless spectrum, corresponding to 89 DA and 76 DC white dwarfs, respectively. We also recover two DBA white dwarfs and 9--10 magnetic remnants. We find two carbon-bearing DQ stars and 14 new metal-rich white dwarfs. This includes the possible detection of the first ultra-cool white dwarf with metal lines. We describe three DZ stars for which we find at least four different metal species, including one which is strongly Fe- and Ni-rich, indicative of the accretion of a planetesimal with core-Earth composition. We find one extremely massive (1.31 $\pm$ 0.01 M$_{\odot}$) DA white dwarf showing weak Balmer lines, possibly indicating stellar magnetism. Another white dwarf shows strong Balmer line emission but no infrared excess, suggesting a low-mass sub-stellar companion. High spectroscopic completeness ($>$99%) has now been reached for $\textit{Gaia}$ DR2 sources within 40 pc sample, in the northern hemisphere ($\delta >$ 0 deg) 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.

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.

We present the first volume-limited sample of cataclysmic variables (CVs), selected using the accurate parallaxes provided by the second data release (DR2) of the ESA $\mathit{Gaia}$ space mission. The sample is composed of 42 CVs within $150\,$pc, including two new systems discovered using the $\mathit{Gaia}$ data, and is $(77 \pm 10)\,$per cent complete. We use this sample to study the intrinsic properties of the Galactic CV population. In particular, the CV space density we derive, $\rho=(4.8^{+0.6}_{-0.8}) \times10^{-6}\,\mathrm{pc}^{-3}$, is lower than predicted by most binary population synthesis studies. We also find a low fraction of period bounce CVs, seven per cent, and an average white dwarf mass of $\langle M_\mathrm{WD} \rangle = (0.83 \pm 0.17)\,\mathrm{M}_\odot$. Both findings confirm previous results, ruling out the presence of observational biases affecting these measurements, as has been suggested in the past. The observed fraction of period bounce CVs falls well below theoretical predictions, by at least a factor of five, and remains one of the open problems in the current understanding of CV evolution. Conversely, the average white dwarf mass supports the presence of additional mechanisms of angular momentum loss that have been accounted for in the latest evolutionary models. The fraction of magnetic CVs in the $150\,$pc sample is remarkably high at $36\,$per cent. This is in striking contrast with the absence of magnetic white dwarfs in the detached population of CV progenitors, and underlines that the evolution of magnetic systems has to be included in the next generation of population models.

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.

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.

We model the present day, observable, normal radio pulsar population of the Small Magellanic Cloud (SMC). The pulsars are generated with SeBa, a binary population synthesis code that evolves binaries and the constituent stellar objects up to remnant formation and beyond. We define radio pulsars by selecting neutron stars that satisfy a selection of criteria defined by Galactic pulsars, and apply the detection thresholds of previous and future SMC pulsar surveys.The number of synthesised and recovered pulsars are exceptionally sensitive to the assumed star formation history and applied radio luminosity model, but is not affected extensively by the assumed common envelope model, metallicity, and neutron star kick velocity distribution. We estimate that the SMC formed (1.6$\pm$0.3)$\times 10^4$ normal pulsars during the last 100 Myrs. We study which pulsars could have been observed by the Parkes multibeam survey of the SMC, by applying the survey's specific selection effects, and recover 4.0$\pm$0.8 synthetic pulsars.This is in agreement with their five observed pulsars. We also apply a proposed MeerKAT configuration for the upcoming SMC survey, and predict that the MeerKAT survey will detect 17.2$\pm$2.5 pulsars.

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 was conducted using electromagnetic waves only. The future Laser Interpreter Space Antenna (LISA) has the potential to reveal Milky Way satellites in gravitational waves emitted by double white dwarf (DWD) binaries. We investigate gravitational wave (GW) signals detectable by LISA as a possible tool for the identification and characterisation of the Milky Way satellites. We use the binary population synthesis technique to model the population of DWDs in dwarf satellites and we assess 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 calibrate predictions for the known Milky Way satellites on their observed properties. We find that DWDs emitting at frequencies $\gtrsim 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_\star>10^6\,$M$_{\odot}$ can host detectable LISA sources with a number of detections that 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. 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 valuable for Galactic archaeology and near-field cosmology.

We perform an analysis of the single white dwarf and the double degenerate binary populations in the solar neighbourhood following a population synthesis approach to investigate the effects of unresolved double degenerates in the white dwarf luminosity function. We consider all unresolved synthetic binaries to be associated with fictitious effective temperatures and surface gravities that are obtained in the same way as if these objects were observed as single point sources. We evaluate the effects of unresolved double white dwarfs assuming that the synthetic samples are observed both by the magnitude-limited SDSS and the volume-limited Gaia surveys, the latter limited to a distance of no more than 100pc. We find that, for our standard model, the impact of unresolved double degenerates in the white dwarf luminosity function derived from the Gaia sample is nearly negligible. Unresolved double degenerates are hence expected to have no effect on the age of the Galactic disc, nor on the star formation history from this population. However, for the SDSS sample, the effect of unresolved double degenerates is significant at the brighter bins (Mbol<6.5 mag), with the fraction of such systems reaching ~40% of the total white dwarf population at Mbol=6 mag. This indicates unresolved double degenerates may influence the constraints on the star formation history derived from the SDSS white dwarf sample.

Tertiary tides (TTs), or the continuous tidal distortion of the tertiary in a hierarchical triple system, can extract energy from the inner binary, inducing within it a proclivity to merge. Despite previous work on the subject, which established that it is significant for certain close triple systems, it is still not a well-understood process. A portion of our ignorance in this regard stems from our inability to integrate a simulation of this phenomenon into conventional stellar evolution codes, since full calculations of these tidal interactions are computationally expensive on stellar evolution timescales. Thus, to attain a better understanding of how these TTs act on longer timescales, an empirical expression of its effects as a function of parameters of the triple system involved is required. In our work, we evaluate the rate at which TTs extract energy from the inner binary within a series of constructed hierarchical triple systems under varying parameters, and study the rate at which the inner binary orbital separation shrinks as a function of those parameters. We find that this rate varies little with the absolute values of the masses of the three component objects, but is very sensitive to the mass ratio of the inner binary $q$, the tertiary radius $R_{\rm 3}$, the inner binary orbital separation $a_{\rm 1}$, the outer orbital separation $a_{\rm 2}$, and the viscoelastic relaxation time of the tertiary $\tau$. More specifically, we find that the percentage by which $a_{\rm 1}$ shrinks per unit time can be reasonably approximated by (1/$a_{\rm 1}$)(d$a_{\rm 1}$/d$t$)=$\left(2.22{\times}10^{-8}{\rm yrs}^{-1}\right)4q\left(1+q\right)^{-2}(R_{\rm 3}/100{\rm R}_{\odot})^{5.2}(a_{\rm 1}/0.2{\rm AU})^{4.8}(a_{\rm 2}/2{\rm AU})^{-10.2}$ $({\tau}/0.534{\rm yrs})^{-1.0}$. We also provide tests of how precise this fitting function is.

Type Ia supernovae (SNe) are thought to originate from the thermonuclear explosions of carbon-oxygen (CO) white dwarfs (WDs). The proposed progenitors of standard type Ia SNe have been studied for decades and can be, generally, divided into explosions of CO WDs accreting material from stellar non-degenerate companions (single-degenerate; SD models), and those arising from the explosive interaction of two CO WDs (double-degenerate; DD models). However, current models for the progenitors of such SNe fail to reproduce the diverse properties of the observed explosions, nor do they explain the inferred rates and the characteristics of the observed populations of type Ia SNe and their expected progenitors. Here we show that the little-studied mergers of CO-WDs with hybrid Helium-CO (He-CO) WDs can provide for a significant fraction of the normal type Ia SNe. Here we use detailed thermonuclear-hydrodynamical and radiative-transfer models to show that a wide range of mergers of CO WDs with hybrid He-CO WDs can give rise to normal type Ia SNe. We find that such He-enriched mergers give rise to explosions for which the synthetic light-curves and spectra resemble those of observed type Ia SNe, and in particular, they can produce a wide range of peak-luminosities, MB(MR)~ 18.4 to 19.2 (~ 18.5 to 19:45), consistent with those observed for normal type Ia SNe. Moreover, our population synthesis models show that, together with the contribution from mergers of massive double CO-WDs (producing the more luminous SNe), they can potentially reproduce the full range of type Ia SNe, their rate and delay-time distribution.

Mergers of white dwarfs (WDs) may lead to a variety of transient astrophysical events, SNIa being one possible outcome. Lyutikov & Toonen (2017, 2019) argued that mergers of WDs result, under various parameter regimes, in unusual central engine-powered supernova and a type of short Gamma Ray Bursts that show extended emission tails. Observations by Gvaramadze et al. (2019) of the central star and the nebula J005311 match to the details the model of Lyutikov & Toonen (2017, 2019) for the immediate product of a merger of a heavy ONeMg WD with CO WD (age, luminosity, stellar size, hydrogen deficiency and chemical composition).

In the last decade, the Kepler and CoRoT space-photometry missions have demonstrated the potential of asteroseismology as a novel, versatile and powerful tool to perform exquisite tests of stellar physics, and to enable precise and accurate characterisations of stellar properties, with impact on both exoplanetary and Galactic astrophysics. Based on our improved understanding of the strengths and limitations of such a tool, we argue for a new small/medium space mission dedicated to gathering high-precision, high-cadence, long photometric series in dense stellar fields. Such a mission will lead to breakthroughs in stellar astrophysics, especially in the metal poor regime, will elucidate the evolution and formation of open and globular clusters, and aid our understanding of the assembly history and chemodynamics of the Milky Way's bulge and few nearby dwarf galaxies.

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.

We suggest that fast-rising blue optical transients (FBOTs) and the brightest event of the class AT2018cow result from an electron-capture collapse to a \NS\ following a merger of a massive ONeMg white dwarf (WD) with another WD. Two distinct evolutionary channels lead to the disruption of the less massive WD during the merger and the formation of a shell burning non-degenerate star incorporating the ONeMg core. During the shell burning stage a large fraction of the envelope is lost to the wind, while mass and angular momentum are added to the core. As a result, the electron-capture collapse occurs with a small envelope mass, after $\sim 10^2-10^4$ years. During the formation of a neutron star as little as $\sim 10^{-2} M_\odot $ of the material is ejected at the bounce-off with mildly relativistic velocities and total energy $\sim$ few $ 10^{50}$ ergs. This ejecta becomes optically thin on a time scale of days - this is the FBOT. During the collapse, the neutron star is spun up and magnetic field is amplified. The ensuing fast magnetically-dominated relativistic wind from the newly formed neutron star shocks against the ejecta, and later against the wind. The radiation-dominated forward shock produces the long-lasting optical afterglow, while the termination shock of the relativistic wind produces the high energy emission in a manner similar to Pulsar Wind Nebulae. If the secondary WD was of the DA type, the wind will likely have $\sim 10^{-4} M_\odot$ of hydrogen; this explains the appearance of hydrogen late in the afterglow spectrum. The model explains many of the puzzling properties of FBOTs/AT2018cow: host galaxies, a fast and light anisotropic ejecta producing a bright optical peak, afterglow high energy emission of similar luminosity to the optical, and late infra-red features.

(Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more.

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.

Interacting binaries containing white dwarfs can lead to a variety of outcomes that range from powerful thermonuclear explosions, which are important in the chemical evolution of galaxies and as cosmological distance estimators, to strong sources of low frequency gravitational wave radiation, which makes them ideal calibrators for the gravitational low-frequency wave detector LISA mission. However, current theoretical evolution models still fail to explain the observed properties of the known populations of white dwarfs in both interacting and detached binaries. Major limitations are that the existing population models have generally been developed to explain the properties of sub-samples of these systems, occupying small volumes of the vast parameter space, and that the observed samples are severely biased. The overarching goal for the next decade is to assemble a large and homogeneous sample of white dwarf binaries that spans the entire range of evolutionary states, to obtain precise measurements of their physical properties, and to further develop the theory to satisfactorily reproduce the properties of the entire population. While ongoing and future all-sky high- and low-resolution optical spectroscopic surveys allow us to enlarge the sample of these systems, high-resolution ultraviolet spectroscopy is absolutely essential for the characterization of the white dwarfs in these binaries. The Hubble Space Telescope is currently the only facility that provides ultraviolet spectroscopy, and with its foreseeable demise, planning the next ultraviolet mission is of utmost urgency.

LISA will open the mHz band of gravitational waves (GWs) to the astronomy community. The strong gravity which powers the variety of GW sources in this band is also crucial in a number of important astrophysical processes at the current frontiers of astronomy. These range from the beginning of structure formation in the early universe, through the origin and cosmic evolution of massive black holes in concert with their galactic environments, to the evolution of stellar remnant binaries in the Milky Way and in nearby galaxies. These processes and their associated populations also drive current and future observations across the electromagnetic (EM) spectrum. We review opportunities for science breakthroughs, involving either direct coincident EM+GW observations, or indirect multimessenger studies. We argue that for the US community to fully capitalize on the opportunities from the LISA mission, the US efforts should be accompanied by a coordinated and sustained program of multi-disciplinary science investment, following the GW data through to its impact on broad areas of astrophysics. Support for LISA-related multimessenger observers and theorists should be sized appropriately for a flagship observatory and may be coordinated through a dedicated mHz GW research center.

The Maunakea Spectroscopic Explorer (MSE) is a planned 11.25-m aperture facility with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. A rebirth of the 3.6m Canada-France-Hawaii Telescope on Maunakea, MSE will use 4332 fibers operating at three different resolving powers (R ~ 2500, 6000, 40000) across a wavelength range of 0.36-1.8mum, with dynamical fiber positioning that allows fibers to match the exposure times of individual objects. MSE will enable spectroscopic surveys with unprecedented scale and sensitivity by collecting millions of spectra per year down to limiting magnitudes of g ~ 20-24 mag, with a nominal velocity precision of ~100 m/s in high-resolution mode. This white paper describes science cases for stellar astrophysics and exoplanet science using MSE, including the discovery and atmospheric characterization of exoplanets and substellar objects, stellar physics with star clusters, asteroseismology of solar-like oscillators and opacity-driven pulsators, studies of stellar rotation, activity, and multiplicity, as well as the chemical characterization of AGB and extremely metal-poor stars.

Anticipating the kinematic constraints from the Gaia mission, we perform an extensive numerical study of the evolution of massive binary systems to predict the peculiar velocities that stars obtain when their companion collapses and disrupts the system. Our aim is to (1) identify which predictions are robust against model uncertainties and assess their implications, (2) investigate which physical processes leave a clear imprint and may therefore be constrained observationally and (3) provide a suite of publicly available model predictions. We find that $22_{-8}^{+26}$% of all massive binary systems merge prior to the first core collapse in the system. Of the remainder, $86_{-9}^{+11}$% become unbound because of the core-collapse. Remarkably, this rarely produce runaway stars (i.e., stars with velocities above 30 km/s). These are outnumbered by more than an order of magnitude by slower unbound companions, or "walkaway stars". This is a robust outcome of our simulations and is due to the reversal of the mass ratio prior to the explosion and widening of the orbit, as we show analytically and numerically. We estimate a $10^{+5}_{-8}$% of massive stars to be walkaways and only $0.5^{+1.0}_{-0.4}$% to be runaways, nearly all of which have accreted mass from their companion. Our findings are consistent with earlier studies, however the low runaway fraction we find is in tension with observed fractions 10%. If Gaia confirms these high fractions of massive runaway stars resulting from binaries, it would imply that we are currently missing physics in the binary models. Finally, we show that high end of the mass distributions of runaway stars is very sensitive to the assumed black hole natal kicks and propose this as a potentially stringent test for the explosion mechanism. We discuss companions remaining bound which can evolve into X-ray and gravitational wave sources.

Neutron-star (NS) - white-dwarf (WD) mergers may give rise to observable explosive transients, but have been little explored. We use 2D coupled hydrodynamical-thermonuclear FLASH-code simulations to study the evolution of WD debris-disks formed following WD-disruptions by NSs. We use a 19-elements nuclear-network and a detailed equation-of-state to follow the evolution, complemented by a post-process analysis using a larger 125-isotopes nuclear-network. We consider a wide range of initial conditions and study the dependence of the results on the NS/WD masses ($1.4-2{\rm M_{\odot}}$;$\,{\rm 0.375-0.7\,M_{\odot}}$, respectively), WD-composition (CO/He/hybrid-He-CO) and the accretion-disk structure. We find that viscous inflow in the disk gives rise to continuous wind-outflow of mostly C/O material mixed with nuclear-burning products arising from a weak detonation occurring in the inner-region of the disk. We find that such transients are energetically weak ($10^{48}-10^{49}$ergs) compared with thermonuclear-supernovae (SNe), and are dominated by the (gravitational) accretion-energy. Although thermonuclear-detonations occur robustly in all of our simulations (besides the He-WD) they produce only little energy $(1-10\%$ of the kinetic energy) and $^{56}{\rm Ni}$ ejecta (few$\times10^{-4}-10^{-3}{\rm M_{\odot}})$, with overall low ejecta masses of $\sim0.01-0.1{\rm M_{\odot}}$. Such explosions may produce rapidly-evolving transients, much shorter and fainter than regular type-Ia SNe. The composition and demographics of such SNe appear to be inconsistent with those of Ca-rich type Ib SNe. Though they might be related to the various classes of rapidly evolving SNe observed in recent years, they are likely to be fainter than the typical ones, and may therefore give rise a different class of potentially observable transients.

The galactic halo likely grew over time in part by assembling smaller galaxies, the so-called building blocks. We investigate if the properties of these building blocks are reflected in the halo white dwarf (WD) population in the Solar neighborhood. Furthermore, we compute the halo WD luminosity functions (WDLFs) for four major building blocks of five cosmologically motivated stellar haloes. We couple the SeBa binary population synthesis model to the Munich-Groningen semi-analytic galaxy formation model, applied to the high-resolution Aquarius dark matter simulations. Although the semi-analytic model assumes an instantaneous recycling approximation, we model the evolution of zero-age main sequence stars to WDs, taking age and metallicity variations of the population into account. Although the majority of halo stars is old and metal-poor and therefore the WDs in the different building blocks have similar properties (including present-day luminosity), we find in our models that the WDs originating from building blocks that have young and/or metal-rich stars can be distinguished from WDs that were born in other building blocks. In practice however, it will be hard to prove that these WDs really originate from different building blocks, as the variations in the halo WD population due to binary WD mergers result in similar effects. The five joined stellar halo WD populations that we modelled result in WDLFs that are very similar to each other. We find that simple models with a Kroupa or Salpeter initial mass function (IMF) fit the observed luminosity function slightly better, since the Chabrier IMF is more top-heavy, although this result is dependent on our choice of the stellar halo mass density in the Solar neighborhood.

We report the discovery of J1953-1019, the first resolved triple white dwarf system. The triplet consists of an inner white dwarf binary and a wider companion. Using Gaia DR2 photometry and astrometry combined with our follow-up spectroscopy, we derive effective temperatures, surface gravities, masses and cooling ages of the three components. All three white dwarfs have pure-hydrogen (DA) atmospheres, masses of 0.60-0.63 Msun and cooling ages of 40-290 Myr. We adopt eight initial-to-final mass relations to estimate the main sequence progenitor masses (which we find to be similar for the three components, 1.6-2.6 Msun) and lifetimes. The differences between the derived cooling times and main sequence lifetimes agree for most of the adopted initial-to-final mass relations, hence the three white dwarfs in J1953-1019 are consistent with coeval evolution. Furthermore, we calculate the projected orbital separations of the inner white dwarf binary (303.25 +- 0.01 au) and of the centre of mass of the inner binary and the outer companion (6398.97 +- 0.09 au). From these values, and taking into account a wide range of possible configurations for the triplet to be currently dynamically stable, we analyse the future evolution of the system. We find that a collision between the two inner white dwarfs due to Lidov-Kozai oscillations is unlikely, though if it occurs it could result in a sub-Chandrasekhar Type Ia supernova explosion.

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<G<21$. By comparing this catalogue with a sample of SDSS white dwarf candidates we estimate an upper limit in completeness of $85$ per cent for white dwarfs with $G \leq 20$ mag and $T_{\mathrm{eff}} > 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.

The mergers of neutron stars (NSs) and white dwarfs (WDs) could give rise to explosive transients, potentially observable with current and future transient surveys. However, the expected properties and distribution of such events is not well understood. Here we characterize the rates of such events, their delay time distribution, their progenitors and the distribution of their properties. We use binary populations synthesis models and consider a wide range of initial conditions and physical processes. In particular we consider different common-envelope evolution models and different NS natal kick distributions. We provide detailed predictions arising from each of the models considered. We find that the majority of NS-WD mergers are born in systems in which mass-transfer played an important role, and the WD formed before the NS. For the majority of the mergers the WDs have a carbon-oxygen composition (60-80%) and most of the rest are with oxygen-neon WDs. The rates of NS-WD mergers are in the range of 3-15% of the type Ia supernovae (SNe) rate. Their delay time distribution is very similar to that of type Ia SNe, but slightly biased towards earlier times. They typically explode in young 0.1-1Gyr environments, but have a tail distribution extending to long, Gyrs-timescales. Models including significant kicks give rise to relatively wide offset distribution extending to hundreds of kpcs. The demographic and physical properties of NS-WD mergers suggest they are likely to be peculiar type Ic-like SNe, mostly exploding in late type galaxies. Their overall properties could be related to a class of rapidly evolving SNe recently observed, while they are less likely to be related to the class of Ca-rich SNe.

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.

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.

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.

A significant fraction of stars are members of gravitationally bound hierarchies containing three or more components. Almost all low mass stars in binaries with periods shorter three days are part of a hierarchical system. We therefore anticipate that a large fraction of compact galactic binaries detected by the Laser Interferometer Space Antenna (LISA) will be members of hierarchical triple or quadruple system. The acceleration imparted by the hierarchical companions can be detected in the gravitational wave signal for outer periods as large as 100 years. For systems with periods that are shorter than, or comparable to, the mission lifetime, it will be possible to measure the period and eccentricity of the outer orbit. LISA observations of hierarchical stellar systems will provide insight into stellar evolution, including the role that Kozai-Lidov oscillations play in driving systems towards merger.

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.

The formal division of the distance indicators into primary and secondary leads to difficulties in description of methods which can actually be used in two ways: with, and without the support of the other methods for scaling. Thus instead of concentrating on the scaling requirement we concentrate on all methods of distance determination to extragalactic sources which are designated, at least formally, to use for individual sources. Among those, the Supernovae Ia is clearly the leader due to its enormous success in determination of the expansion rate of the Universe. However, new methods are rapidly developing, and there is also a progress in more traditional methods. We give a general overview of the methods but we mostly concentrate on the most recent developments in each field, and future expectations.

Type Ia supernovae (SNe Ia) exhibit a wide diversity of peak luminosities and light curve shapes: the faintest SNe Ia are 10 times less luminous and evolve more rapidly than the brightest SNe Ia. Their differing characteristics also extend to their stellar age distributions, with fainter SNe Ia preferentially occurring in old stellar populations and vice versa. In this Letter, we quantify this SN Ia luminosity - stellar age connection using data from the Lick Observatory Supernova Search (LOSS). Our binary population synthesis calculations agree qualitatively with the observed trend in the >1 Gyr-old populations probed by LOSS if the majority of SNe Ia arise from prompt detonations of sub-Chandrasekhar mass white dwarfs (WDs) in double WD systems. Under appropriate assumptions, we show that double WD systems with less massive primaries, which yield fainter SNe Ia, interact and explode at older ages than those with more massive primaries. We find that prompt detonations in double WD systems are capable of reproducing the observed evolution of the SN Ia luminosity function, a constraint that any SN Ia progenitor scenario must confront.

We discuss a scenario of short Gamma-Ray Bursts (GRBs) following a merger of a massive ONeMg white dwarf (WD) with a CO WD, and an ensuing accretion induced collapse (AIC). An initial system with the primary mass $M_1 \sim 6-10 \, M_\odot$ and the secondary mass $M_2 \sim 3-6\, M_\odot$ forms, via two distinct evolutionary channels, a double degenerate CO-ONeMg WD system. For sufficiently large mass ratio $q \equiv M_2/M_1 > q_{crit} \sim 0.25$ the ensuing gravitational wave-driven mass transfer is unstable, whereby the less massive CO WD is disrupted and transfers its mass to the primary ONeMg WD on a few orbital time scales. The merger product ignites shell CO burning, adding mass to the degenerate core; at the same time mass and angular momentum is lost due to powerful winds. For an ONeMg WD sufficiently close to the Chandrasekhar mass an electron-capture accretion induced collapse (AIC) follows $\sim 10^4$ years afterwards. We associate the prompt short GRB emission with a direct collapse of an ONeMg WD to a neutron star, without formation of an accretion disk. After the collapse the accretion of the unburnt part of the shell onto the newly formed NS powers the extended emission (EE). During the collapse the neutron star is spun to millisecond periods and produce long lasting relativistic winds that shock against the material lost during the shell-burning stage, and produce afterglow emission from the wind termination shock.

Type Ia supernovae (Ia-SNe) are thought to arise from the thermonuclear explosions of white dwarfs (WDs). The progenitors of such explosions are still highly debated; in particular the conditions leading to detonations in WDs are not well understood in most of the suggested progenitor models. Nevertheless, direct head-on collisions of two WDs were shown to give rise to detonations and produce Ia-SNe - like explosions, and were suggested as possible progenitors. The rates of such collisions in dense globular clusters are far below the observed rates of type Ia SNe, but it was suggested that quasi-secular evolution of hierarchical triples could produce a high rate of such collisions. Here we used detailed triple stellar evolution populations synthesis models coupled with dynamical secular evolution to calculate the rates of WD-WD collisions in triples and their properties. We explored a range of models with different realistic initial conditions and derived the expected SNe total mass, mass-ratio and delay time distributions for each of the models. We find that the SNe rate from WD-WD collisions is of the order of 0.1% of the observed Ia-SNe rate across all our models, and the delay-time distribution is almost uniform in time, and is inconsistent with observations. We conclude that SNe from WD-WD collisions in isolated triples can at most provide for a small fraction of Ia-SNe, and can not serve as the main progenitors of such explosions.

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.

We consider the formation of binary black hole mergers through the evolution of field massive triple stars. In this scenario, favorable conditions for the inspiral of a black hole binary are initiated by its gravitational interaction with a distant companion, rather than by a common-envelope phase invoked in standard binary evolution models. We use a code that follows self-consistently the evolution of massive triple stars, combining the secular triple dynamics (Lidov-Kozai cycles) with stellar evolution. After a black hole triple is formed, its dynamical evolution is computed using either the orbit-averaged equations of motion, or a high-precision direct integrator for triples with weaker hierarchies for which the secular perturbation theory breaks down. Most black hole mergers in our models are produced in the latter non-secular dynamical regime. We derive the properties of the merging binaries and compute a black hole merger rate in the range (0.3- 1.3) Gpc^{-3}yr^{-1}, or up to ~2.5Gpc^{-3}yr^{-1} if the black hole orbital planes have initially random orientation. Finally, we show that black hole mergers from the triple channel have significantly higher eccentricities than those formed through the evolution of massive binaries or in dense star clusters. Measured eccentricities could therefore be used to uniquely identify binary mergers formed through the evolution of triple stars. While our results suggest up to ~10 detections per year with Advanced-LIGO, the high eccentricities could render the merging binaries harder to detect with planned space based interferometers such as LISA.

Even though Type Ia supernovae (SNIa) play an important role in many fields in astronomy, the nature of the progenitors of SNIa remain a mystery. One of the classical evolutionary pathways towards a SNIa explosion is the single degenerate (SD) channel, in which a carbon-oxygen white dwarf accretes matter from its non-degenerate companion until it reaches the Chandrasekhar mass. Constraints on the contribution from the SD channel to the overall SNIa rate come from a variety of methods, e.g. from abundances, from signatures of the companion star in the light curve or near the SNIa remnant, and from synthetic SNIa rates. In this proceedings, I show that when incorporating our newest understandings of binary evolution, the SNIa rate from the single degenerate channel is enhanced. I also discuss the applicability of these constraints on the evolution of SNIa progenitors.

Double white dwarf (DWD) binaries are expected to be very common in the Milky Way, but their intrinsic faintness challenges the detection of these systems. Currently, only a few tens of detached DWDs are know. Such systems offer the best chance of extracting the physical properties that would allow us to address a wealth of outstanding questions ranging from the nature of white dwarfs, over stellar and binary evolution to mapping the Galaxy. In this paper we explore the prospects for detections of ultra-compact (with binary separations of a few solar radii or less) detached DWDs in: 1) optical radiation with Gaia and the LSST and 2) gravitational wave radiation with LISA. We show that Gaia, LSST and LISA have the potential to detect respectively around a few hundreds, a thousand, and 25 thousand DWD systems. Moreover, Gaia and LSST data will extend by respectively a factor of two and seven the guaranteed sample of binaries detected in electromagnetic and gravitational wave radiation, opening the era of multi-messenger astronomy for these sources.

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.

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$.

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.

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]

Aims: We study single and binary white dwarfs in the inner halo of the Milky Way in order to learn more about the conditions under which the population of halo stars was born, such as the initial mass function (IMF), the star formation history, or the binary fraction. Methods: We simulate the evolution of low-metallicity halo stars at distances up to ~ 3 kpc using the binary population synthesis code SeBa. We use two different white dwarf cooling models to predict the present-day luminosities of halo white dwarfs. We determine the white dwarf luminosity functions (WDLFs) for eight different halo models and compare these with the observed halo WDLF of white dwarfs in the SuperCOSMOS Sky Survey. Furthermore, we predict the properties of binary white dwarfs in the halo and determine the number of halo white dwarfs that is expected to be observed with the Gaia satellite. Results: By comparing the WDLFs, we find that a standard IMF matches the observations more accurately than a top-heavy one, but the difference with a bottom-heavy IMF is small. A burst of star formation 13 Gyr ago fits slightly better than a star formation burst 10 Gyr ago and also slightly better than continuous star formation $10-13$ Gyrs ago. Gaia will be the first instrument to constrain the bright end of the field halo WDLF, where contributions from binary WDs are considerable. Many of these will have He cores, of which a handful have atypical surface gravities ($\log g < 6$) and reach luminosities $\log(L/L_\odot) > 0$ in our standard model for WD cooling. These so called pre-WDs, if observed, can help us to constrain white dwarf cooling models and might teach us something about the fraction of halo stars that reside in binaries.

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.

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.

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.

Context. An important ingredient in binary evolution is the common-envelope (CE) phase. Although this phase is believed to be responsible for the formation of many close binaries, the process is not well understood. Aims. We investigate the characteristics of the population of post-common-envelope binaries (PCEB). As the evolution of these binaries and their stellar components are relatively simple, this population can be directly used to constraint CE evolution. Methods. We use the binary population synthesis code SeBa to simulate the current-day population of PCEBs in the Galaxy. We incorporate the selection effects in our model that are inherent to the general PCEB population and that are specific to the SDSS survey, which enables a direct comparison for the first time between the synthetic and observed population of visible PCEBs. Results. We find that selection effects do not play a significant role on the period distribution of visible PCEBs. To explain the observed dearth of long-period systems, the {\alpha}-CE efficiency of the main evolutionary channel must be low. In the main channel, the CE is initiated by a red giant as it fills its Roche lobe in a dynamically unstable way. Other evolutionary paths cannot be constrained more. Additionally our model reproduces well the observed space density, the fraction of visible PCEBs amongst white dwarf (WD)- main sequence (MS) binaries, and the WD mass versus MS mass distribution, but overestimates the fraction of PCEBs with helium WD companions.

[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.

There is general agreement that supernovae Ia correspond to the thermonuclear runaway of a white dwarf that is part of a compact binary, but the details of the progenitor systems are still unknown and much debated. One of the proposed progenitor theories is the single-degenerate channel in which a white dwarf accretes from a companion, grows in mass, reaches a critical mass limit, and is then consumed after thermonuclear runaway sets in. However, there are major disagreements about the theoretical delay time distribution and the corresponding time-integrated supernova Ia rate from this channel. We investigate whether the differences are due to the uncertainty in the common envelope phase and the fraction of transferred mass that is retained by the white dwarf. This so-called retention efficiency may have a strong influence on the final amount and timing of supernovae Ia. Using the population synthesis code SeBa, we simulated large numbers of binaries for various assumptions on common envelopes and retention efficiencies. We compare the resulting supernova Ia rates and delay time distributions with each other and with those from the literature, including observational data. For the three assumed retention efficiencies, the integrated rate varies by a factor 3-4 to even more than a factor 100, so in extreme cases, the retention efficiency strongly suppresses the single-degenerate channel. Our different assumptions for the common envelope phase change the integrated rate by a factor 2-3. Although our results do recover the trend in the theoretical predictions from different binary population synthesis codes, they do not fully explain the large disagreement among them.

Despite the significance of Type Ia supernovae (SNeIa) in many fields in astrophysics, SNeIa lack a theoretical explanation. The standard scenarios involve thermonuclear explosions of carbon/oxygen white dwarfs approaching the Chandrasekhar mass; either by accretion from a companion or by a merger of two white dwarfs. We investigate the contribution from both channels to the SNIa rate with the binary population synthesis (BPS) code SeBa in order to constrain binary processes such as the mass retention efficiency of WD accretion and common envelope evolution. We determine the theoretical rates and delay time distribution of SNIa progenitors and in particular study how assumptions affect the predicted rates.

Despite the significance of type Ia supernovae (SNeIa) in many fields in astrophysics, SNeIa lack a theoretical explanation. We investigate the potential contribution to the SNeIa rate from the most common progenitor channels using the binary population synthesis (BPS) code SeBa. Using SeBa, we aim constrain binary processes such as the common envelope phase and the efficiency of mass retention of white dwarf accretion. We find that the simulated rates are not sufficient to explain the observed rates. Further, we find that the mass retention efficiency of white dwarf accretion significantly influences the rates, but does not explain all the differences between simulated rates from different BPS codes.

We compare the results of three binary population synthesis codes to understand the differences in their results. As a first result we find that when equalizing the assumptions the results are similar. The main differences arise from deviating physical input.

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.

We briefly discuss the method of population synthesis to calculate theoretical delay time distributions of type Ia supernova progenitors. We also compare the results of the different research groups and conclude that although one of the main differences in the results for single degenerate progenitors is the retention efficiency with which accreted hydrogen is added to the white dwarf core, this cannot explain all the differences.

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.

Stars on eccentric orbits around a massive black hole (MBH) emit bursts of gravitational waves (GWs) at periapse. Such events may be directly resolvable in the Galactic centre. However, if the star does not spiral in, the emitted GWs are not resolvable for extra-galactic MBHs, but constitute a source of background noise. We estimate the power spectrum of this extreme mass ratio burst background (EMBB) and compare it to the anticipated instrumental noise of the Laser Interferometer Space Antenna (LISA). To this end, we model the regions close to a MBH, accounting for mass-segregation, and for processes that limit the presence of stars close to the MBH, such as GW inspiral and hydrodynamical collisions between stars. We find that the EMBB is dominated by GW bursts from stellar mass black holes, and the magnitude of the noise spectrum (f S_GW)^{1/2} is at least a factor ~10 smaller than the instrumental noise. As an additional result of our analysis, we show that LISA is unlikely to detect relativistic bursts in the Galactic centre.

The Fabry-Perot interferometer FaNTOmM was used at the 3.6m Canada France Hawaii Telescope and the 1.6m Mont Megantic Telescope to obtain data cubes in H-alpha of 9 nearby spiral galaxies from which maps in integrated intensity, velocity, and velocity dispersion were derived. We then applied the Tremaine-Weinberg method, in which the pattern speed can be deduced from its velocity field, by finding the integrated value of the mean velocity along a slit parallel to the major axis weighted by the intensity and divided by the weighted mean distance of the velocity points from the tangent point measured along the slit. The measured variables can be used either to make separate calculations of the pattern speed and derive a mean, or in a plot of one against the other for all the points on all slits, from which a best fit value can be derived. Linear fits were found for all the galaxies in the sample. For two galaxies a clearly separate inner pattern speed with a higher value, was also identified and measured.

We study the kinematics of the barred spiral galaxy NGC 6946 by investigating the velocity field from H-alpha Fabry-Perot observations, determined the pattern speed of the bar by using the Tremaine-Weinberg method, and find a main pattern speed of 21.7 (+4.0,-0.8) km/s/kpc. Our data clearly suggest the presence of an additional pattern with a pattern speed more than twice that of the large pattern in this galaxy. We use the epicycle approximation to deduce the location of the resonance radii and subsequently determine the pattern speed between the radii, and find that inside the inner Inner Lindblad Resonance radius, a bar-like system has evolved.

We study the two-dimensional kinematics of the H-alpha-emitting gas in the nearby barred Scd galaxy, NGC 6946, in order to determine the pattern speed of the primary m=2 perturbation mode. The pattern speed is a crucial parameter for constraining the internal dynamics, estimating the impact velocities of the gravitational perturbation at the resonance radii, and to set up an evolutionary scenario for NGC 6946. Our data allows us to derive the best fitting kinematic position angle and the geometry of the underlying gaseous disk, which we use to derive the pattern speed using the Tremaine-Weinberg method. We find a main pattern speed Omega_p=22 km/s/kpc, but our data clearly reveal the presence of an additional pattern speed Omega_p=47 km/s/kpc in a zone within 1.25 kpc of the nucleus. Using the epicyclic approximation, we deduce the location of the resonance radii and confirm that inside the outer Inner Lindblad Resonance radius of the main oval, a primary bar has formed rotating at more than twice the outer pattern speed. We further confirm that a nuclear bar has formed inside the Inner Lindblad Resonance radius of the primary bar, coinciding with the inner Inner Lindblad Resonance radius of the large-scale m=2 mode oval.