Present time-domain astronomy efforts will unveil a variety of rare transients. We focus here on pulsational pair-instability evolution, which can result in signatures observable with electromagnetic and gravitational waves. We simulate grids of bare helium stars to characterize the resulting black hole (BH) masses and ejecta composition, velocity, and thermal state. The stars do not react "elastically" to the thermonuclear explosion: there is not a one-to-one correspondence between pair-instability driven ignition and mass ejections, causing ambiguity in what is an observable pulse. In agreement with previous studies, we find that for carbon-oxygen core masses 28Msun< M_CO<30.5Msun the explosions are not strong enough to affect the surface. With increasing mass, they first cause large radial expansion (30.5Msun<M_CO<31.4Msun), and finally, also mass ejection episodes (M_CO>31.4Msun). The lowest mass to be fully disrupted in a pair-instability supernova is M_CO=57Msun. Models with M_CO>121Msun reach the photodisintegration regime, resulting in BHs with M_BH>125Msun. If the pulsating models produce BHs via (weak) explosions, the previously-ejected material might be hit by the blast wave. We characterize the H-free circumstellar material from the pulsational pair-instability of helium cores assuming simply that the ejecta maintain a constant velocity after ejection. Our models produce He-rich ejecta with mass 10^{-3}Msun<M_CSM<40Msun. These ejecta are typically launched at a few thousand \kms and reach distances of ~10^{12}-10^{15} cm before core-collapse. The delays between mass ejection events and the final collapse span a wide and mass-dependent range (from sub-hour to 10^4 years), and the shells ejected can also collide with each other. The range of properties we find suggests a possible connection with (some) type Ibn supernovae.

Until recently, black holes (BHs) could be discovered only through accretion from other stars in X-ray binaries, or in merging double compact objects. Improvements in astrometric and spectroscopic measurements have made it possible to detect BHs also in non-interacting BH binaries (nBHB) through a precise analysis of the companion's motion. In this study, using an updated version of the \startrack\ binary-star population modelling code and a detailed model of the Milky Way (MW) galaxy we calculate the expected number of detections for \gaia\ and \lamost\ surveys. We develop a formalism to convolve the binary population synthesis output with a realistic stellar density distribution, star-formation history (SFH), and chemical evolution for the MW, which produces a probability distribution function of the predicted compact-binary population over the MW. This avoids the additional statistical uncertainty which is introduced by methods which Monte Carlo sample from binary population synthesis output to produce one potential specific realisation of the MW compact-binary distribution, and our method is also comparatively fast to such Monte Carlo realisations. Specifically, we predict $\sim41$--$340$ nBHBs to be observed by \gaia, although the numbers may drop to $\sim10$--$70$ if the recent ($\lesssim100\;$ Myr) star formation is low ($\sim1\;M_\odot$/yr ). For \lamost\ we predict $\lesssim14$ detectable nBHBs, which is lower partially because its field-of-view covers just $\sim6\%$ of the Galaxy.

Gravitational wave detections are starting to allow us to probe the physical processes in the evolution of very massive stars through the imprints they leave on their final remnants. Stellar evolution theory predicts the existence of a gap in the black hole mass distribution at high mass due to the effects of pair-instability. Previously, we showed that the location of the gap is robust against model uncertainties, but it does depend sensitively on the uncertain $^{12}\rm{C}\left(\alpha,\gamma\right)^{16}\!\rm{O}$ rate. This rate is of great astrophysical significance and governs the production of oxygen at the expense of carbon. We use the open source MESA stellar evolution code to evolve massive helium stars to probe the location of the mass gap. We find that the maximum black hole mass below the gap varies between $40\rm{M}_\odot$ to $90\rm{M}_\odot$, depending on the strength of the uncertain $^{12}\rm{C}\left(\alpha,\gamma\right)^{16}\!\rm{O}$ reaction rate. With the first ten gravitational-wave detections of black holes, we constrain the astrophysical S-factor for $^{12}\rm{C}\left(\alpha,\gamma\right)^{16}\!\rm{O}$, at $300\rm{keV}$, to $S_{300}>175\rm{\,keV\, barns}$ at 68% confidence. With $\mathcal{O}(50)$ detected binary black hole mergers, we expect to constrain the S-factor to within $\pm10$-$30\rm{\,keV\, barns}$. We also highlight a role for independent constraints from electromagnetic transient surveys. The unambiguous detection of pulsational pair instability supernovae would imply that $S_{300}>79\rm{\,keV\, barns}$. Degeneracies with other model uncertainties need to be investigated further, but probing nuclear stellar astrophysics poses a promising science case for the future gravitational wave detectors.

The theory for single stellar evolution predicts a gap in the mass distribution of black holes (BHs) between approximately 45-130M$_{\odot}$, the so-called "pair-instability mass gap". We examine whether BHs can pollute the gap after accreting from a stellar companion. To this end, we simulate the evolution of isolated binaries using a population synthesis code, where we allow for super-Eddington accretion. Under our most extreme assumptions, we find that at most about 2% of all merging binary BH systems contains a BH with a mass in the pair-instability mass gap, and we find that less than 0.5% of the merging systems has a total mass larger than 90M$_{\odot}$. We find no merging binary BH systems with a total mass exceeding 100M$_{\odot}$. We compare our results to predictions from several dynamical pathways to pair-instability mass gap events and discuss the distinguishable features. We conclude that the classical isolated binary formation scenario will not significantly contribute to the pollution of the pair-instability mass gap. The robustness of the predicted mass gap for the isolated binary channel is promising for the prospective of placing constraints on (i) the relative contribution of different formation channels, (ii) the physics of the progenitors including nuclear reaction rates, and (iii), tentatively, the Hubble parameter.

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.

Massive binaries that merge as compact objects are the progenitors of gravitational-wave sources. Most of these binaries experience one or more phases of mass transfer, during which one of the stars loses part or all of its outer envelope and becomes a stripped-envelope star. The evolution of the size of these stripped stars is crucial in determining whether they experience further interactions and their final fate. We present new calculations of stripped-envelope stars based on binary evolution models computed with MESA. We use these to investigate their radius evolution as a function of mass and metallicity. We further discuss their pre-supernova observable characteristics and potential consequences of their evolution on the properties of supernovae from stripped stars. At high metallicity we find that practically all of the hydrogen-rich envelope is removed, in agreement with earlier findings. Only progenitors with initial masses below 10\Msun expand to large radii (up to 100\Rsun), while more massive progenitors stay compact. At low metallicity, a substantial amount of hydrogen remains and the progenitors can, in principle, expand to giant sizes (> 400\Rsun), for all masses we consider. This implies that they can fill their Roche lobe anew. We show that the prescriptions commonly used in population synthesis models underestimate the stellar radii by up to two orders of magnitude. We expect that this has consequences for the predictions for gravitational-wave sources from double neutron star mergers, in particular for their metallicity dependence.

Gravitational-wave detections are now probing the black hole (BH) mass distribution, including the predicted pair-instability mass gap. These data require robust quantitative predictions, which are challenging to obtain. The most massive BH progenitors experience episodic mass ejections on timescales shorter than the convective turn-over timescale. This invalidates the steady-state assumption on which the classic mixing-length theory relies. We compare the final BH masses computed with two different versions of the stellar evolutionary code \texttt{MESA}: (i) using the default implementation of \cite{paxton:18} and (ii) solving an additional equation accounting for the timescale for convective deceleration. In the second grid, where stronger convection develops during the pulses and carries part of the energy, we find weaker pulses. This leads to lower amounts of mass being ejected and thus higher final BH masses of up to $\sim$\,$5\,M_\odot$. The differences are much smaller for the progenitors which determine the maximum mass of BHs below the gap. This prediction is robust at $M_{\rm BH, max}\simeq 48\,M_\odot$, at least within the idealized context of this study. This is an encouraging indication that current models are robust enough for comparison with the present-day gravitational-wave detections. However, the large differences between individual models emphasize the importance of improving the treatment of convection in stellar models, especially in the light of the data anticipated from the third generation of gravitational wave detectors.

The majority of massive stars, the progenitors of core-collapse supernovae (SNe), are found in close binary systems. Zapartas et al. (2019) modeled the fraction of hydrogen-rich, Type II SN progenitors which have their evolution affected by mass exchange with their companion, finding this to be between 1/3 and 1/2 for most assumptions. Here we study in more depth the impact of this binary history of Type II SN progenitors on their final pre-SN core mass distribution, using population synthesis simulations. We find that binary star progenitors of Type II SNe typically end their life with a larger core mass than they would have had if they had lived in isolation, because they gained mass or merged with a companion before explosion. The combination of the diverse binary evolutionary paths typically lead to a marginally shallower final core mass distribution. Discussing our results in the context of the red supergiant problem, i.e., the reported lack of detected high luminosity progenitors, we conclude that binary evolution does not seem to significantly affect the issue. This conclusion is quite robust against our variations in the assumptions of binary physics. We also predict that inferring the initial masses of Type II SN progenitors from "age-dating" its surrounding environment systematically yields lower masses compared to methods that probe the pre-SN core mass or luminosity. A robust discrepancy between the inferred initial masses of a SN progenitor from those different techniques could indicate an evolutionary history of binary mass accretion or merging.

Close Double Neutron Stars (DNSs) have been observed as Galactic radio pulsars, while their mergers have been detected as gamma-ray bursts and gravitational-wave sources. They are believed to have experienced at least one common-envelope episode (CEE) during their evolution prior to DNS formation. In the last decades there have been numerous efforts to understand the details of the common-envelope phase, but its computational modelling remains challenging. We present and discuss the properties of the donor and the binary at the onset of the Roche-lobe overflow (RLOF) leading to these CEEs as predicted by rapid binary population synthesis models. These properties can be used as initial conditions for detailed simulations of the common-envelope phase. There are three distinctive populations, classified by the evolutionary stage of the donor at the moment of the onset of the RLOF: giant donors with fully-convective envelopes, cool donors with partially-convective envelopes, and hot donors with radiative envelopes. We also estimate that, for standard assumptions, tides would not circularise a large fraction of these systems by the onset of RLOF. This makes the study and understanding of eccentric mass transferring systems relevant for DNS populations.

A class of optical transients known as Luminous Red Novae (LRNe) have recently been associated with mass ejections from binary stars undergoing common-envelope evolution. We use the population synthesis code COMPAS to explore the impact of a range of assumptions about the physics of common-envelope evolution on the properties of LRNe. In particular, we investigate the influence of various models for the energetics of LRNe on the expected event rate and light curve characteristics, and compare with the existing sample. We find that the Galactic rate of LRNe is $\sim 0.2$ yr$^{-1}$, in agreement with the observed rate. In our models, the luminosity function of Galactic LRNe covers multiple decades in luminosity and is dominated by signals from stellar mergers, consistent with observational constraints from iPTF and the Galactic sample of LRNe. We discuss how observations of the brightest LRNe may provide indirect evidence for the existence of massive ($> 40$ M$_\odot$) red supergiants. Such LRNe could be markers along the evolutionary pathway leading to the formation of double compact objects. We make predictions for the population of LRNe observable in future transient surveys with the Large Synoptic Survey Telescope and the Zwicky Transient Facility. In all plausible circumstances, we predict a selection-limited observable population dominated by bright, long-duration events caused by common envelope ejections. We show that the Large Synoptic Survey Telescope will observe $20$--$750$ LRNe per year, quickly constraining the luminosity function of LRNe and probing the physics of common-envelope events.

A 70Msun BH was discovered in Milky Way disk in a long period and almost circular detached binary system (LB-1) with a high metallicity 8Msun B star companion. Current consensus on the formation of BHs from high metallicity stars limits the black hole mass to be below 20Msun. Using simple evolutionary model, we show that the formation of a 70Msun BH in high metallicity environment is possible if stellar wind mass loss rates are reduced by factor of 5. As observations indicate, a fraction of massive stars have surface magnetic fields which may quench the wind mass-loss, independently of stellar mass and metallicity. We also computed detailed stellar evolution models and we confirm such a scenario. A non-rotating 85Msun model at Z=0.014 with decreased winds ends up as a 71Msun star prior core-collapse with a 32Msun helium core and a 28Msun CO core. Such star avoids pair-instability pulsation supernova mass loss and may form a 70Msun BH in the direct collapse. Stars that can form such BHs expand to significant size with radius of R>600Rsun, exceeding the size of LB-1 orbit. Therefore, we can explain the formation of BHs upto 70Msun at high metallicity and this result is independent from LB-1. However, if LB-1 hosts a massive BH we are unable to explain how such a binary star system could have formed without invoking some exotic scenarios.

All stellar mass black holes have hitherto been identified by X-rays emitted by gas that is accreting onto the black hole from a companion star. These systems are all binaries with black holes below 30 M$_{\odot}$$^{1-4}$. Theory predicts, however, that X-ray emitting systems form a minority of the total population of star-black hole binaries$^{5,6}$. When the black hole is not accreting gas, it can be found through radial velocity measurements of the motion of the companion star. Here we report radial velocity measurements of a Galactic star, LB-1, which is a B-type star, taken over two years. We find that the motion of the B-star and an accompanying H$\alpha$ emission line require the presence of a dark companion with a mass of $68^{+11}_{-13}$ M$_{\odot}$, which can only be a black hole. The long orbital period of 78.9 days shows that this is a wide binary system. The gravitational wave experiments have detected similarly massive black holes$^{7,8}$, but forming such massive ones in a high-metallicity environment would be extremely challenging to current stellar evolution theories$^{9-11}$.

Gravitational-wave detections are now starting to probe the mass distribution of stellar-mass black holes (BHs). Robust predictions from stellar models are needed to interpret these. Theory predicts the existence of a gap in the BH mass distribution because of pair-instability supernova. The maximum BH mass below the gap is the result of pulsational mass loss. We evolve massive helium stars through their late hydrodynamical phases of evolution using the open-source MESA stellar evolution code. We find that the location of the lower edge of the mass gap at 45$M_\odot$ is remarkably robust against variations in the metallicity ($\approx 3M_\odot$), the treatment of internal mixing ($\approx 1M_\odot$), stellar wind mass loss ($\approx 4M_\odot$), making it the most robust predictions for the final stages of massive star evolution. The reason is that the onset of the instability is dictated by the near-final core mass, which in turn sets the resulting BH mass. However, varying $^{12}C\left(\alpha,\gamma\right)^{16}O$ reaction rate within its $1\sigma$ uncertainties shifts the location of the gap between $40M_\odot$ and $56M_\odot$. We provide updated analytic fits for population synthesis simulations. Our results imply that the detection of merging BHs can provide constraints on nuclear astrophysics. Furthermore, the robustness against metallicity suggests that there is a universal maximum for the location of the lower edge of the gap, which is insensitive to the formation environment and redshift for first-generation BHs. This is promising for the possibility to use the location of the gap as a "standard siren" across the Universe.

Gravitational-wave observations of double compact object (DCO) mergers are providing new insights into the physics of massive stars and the evolution of binary systems. Making the most of expected near-future observations for understanding stellar physics will rely on comparisons with binary population synthesis models. However, the vast majority of simulated binaries never produce DCOs, which makes calculating such populations computationally inefficient. We present an importance sampling algorithm, STROOPWAFEL, that improves the computational efficiency of population studies of rare events, by focusing the simulation around regions of the initial parameter space found to produce outputs of interest. We implement the algorithm in the binary population synthesis code COMPAS, and compare the efficiency of our implementation to the standard method of Monte Carlo sampling from the birth probability distributions. STROOPWAFEL finds $\sim$25-200 times more DCO mergers than the standard sampling method with the same simulation size, and so speeds up simulations by up to two orders of magnitude. Finding more DCO mergers automatically maps the parameter space with far higher resolution than when using the traditional sampling. This increase in efficiency also leads to a decrease of a factor $\sim$3-10 in statistical sampling uncertainty for the predictions from the simulations. This is particularly notable for the distribution functions of observable quantities such as the black hole and neutron star chirp mass distribution, including in the tails of the distribution functions where predictions using standard sampling can be dominated by sampling noise.

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 massive O4.5 V + O5.5 V binary VFTS 352 in the Tarantula nebula is one of the shortest-period and most massive overcontact binaries known. Recent theoretical studies indicate that some of these systems could ultimately lead to the formation of gravitational waves via black hole binary mergers through the chemically homogeneous evolution pathway. By analyzing ultraviolet-optical phase-resolved spectroscopic data, we aim to constrain atmospheric and wind properties that could be later used to confront theoretical predictions from binary evolution. In particular, surface abundances are powerful diagnostics of the evolutionary status, mass transfer and the internal mixing processes. From a set of 32 VLT/FLAMES visual and 8 HST/COS ultraviolet spectra, we used spectral disentangling to separate the primary and secondary components. Using a genetic algorithm wrapped around the NLTE model atmosphere and spectral synthesis code FASTWIND, we perform an 11-parameter optimization to derive the atmospheric and wind parameters of both components, including the surface abundances of He, C, N, O and Si. We find that both components are hotter than expected compared to single-star evolutionary models indicating that additional mixing processes may be at play. However the derived chemical abundances do not show significant indications of mixing when adopting baseline values typical for the system environment.

Interactions between massive stars in binaries are thought to be responsible for much of the observed diversity of supernovae. As surveys probe rarer populations of events, we should expect to see supernovae arising from increasingly uncommon progenitor channels. Here we examine a scenario in which massive stars merge after they have both formed a hydrogen-exhausted core. We suggest this could produce stars which explode as pair-instability supernovae (PISNe) with significantly more hydrogen, at a given metallicity, than in single-star models with the same pre-explosion oxygen-rich core mass. We investigate the subset of those stellar mergers which later produce pulsational PISNe, and estimate that the rate of such post-merger, hydrogen-rich pulsational PISNe could approach a few in a thousand of all core-collapse supernovae. The nature and predicted rate of such hydrogen-rich pulsational PISNe are reminiscent of the very unusual supernova iPTF14hls. For plausible assumptions, PISNe from similar mergers might dominate the rate of PISNe in the local Universe.

Next-generation observations will revolutionize our understanding of binary black holes and will detect new sources, such as intermediate-mass black holes. Primary science goals include: Discover binary black holes throughout the observable Universe; Reveal the fundamental properties of black holes; Uncover the seeds of supermassive black holes.

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.

Massive stars in binaries can give rise to extreme phenomena such as X-ray binaries and gravitational wave sources after one or both stars end their lives as core-collapse supernovae. Stars in close orbit around a stellar or compact companion are expected to explode as "stripped-envelope supernovae", showing no (Type Ib/c) or little (Type IIb) signs of hydrogen in the spectra, because hydrogen-rich progenitors are too large to fit. The physical processes responsible for the stripping process and the fate of the companion are still very poorly understood. Aiming to find new clues, we investigate Cas~A, which is a very young ($\sim$340 \,yr) and near ($\sim$3.4\,kpc) remnant of a core collapse supernova. Cas~A has been subject to several searches for possible companions, all unsuccessfully. We present new measurements of the proper motions and photometry of stars in the vicinity based on deep HST ACS/WFC and WFC3-IR data. We identify stellar sources that are close enough in projection, but using their proper motions we show that none are compatible with being at the location of center at the time of explosion, in agreement with earlier findings. Our photometric measurements allow us to place much deeper (order of magnitude) upper limits on the brightness of possible undetected companions. We systematically compare them with model predictions for a wide variety of scenarios. We can confidently rule out the presence of any stellar companion of any reasonable mass and age (main sequence, pre main sequence or stripped). Although this finding is not in direct conflict with model predictions, it does rule out what many considered to be the most likely evolutionary scenario for type IIb supernova. (abstract abbreviated)

How very massive stars form is still an open question in astrophysics. VFTS682 is among the most massive stars known, with an inferred initial mass of $\sim$150 $M_\odot$ . It is located in 30 Doradus at a projected distance of 29 pc from the central cluster R136. Its apparent isolation led to two hypotheses: either it formed in relative isolation or it was ejected dynamically from the cluster. We investigate the kinematics of VFTS682 as obtained by Gaia and Hubble Space Telescope astrometry. We derive a projected velocity relative to the cluster of $38 \pm 17 \mathrm{km \ s^{-1}}$ (1$\sigma$ confidence interval). Although the error bars are substantial, two independent measures suggest that VFTS682 is a runaway ejected from the central cluster. This hypothesis is further supported by a variety of circumstantial clues. The central cluster is known to harbor other stars more massive than 150 $M_\odot$ of similar spectral type and recent astrometric studies on VFTS16 and VFTS72 provide direct evidence that the cluster can eject some of its most massive members, in agreement with theoretical predictions. If future data confirm the runaway nature, this would make VFTS682 the most massive runaway star known to date.

Double neutron stars (DNSs) have been observed as Galactic radio pulsars, and the recent discovery of gravitational waves from the DNS merger GW170817 adds to the known DNS population. We perform rapid population synthesis of massive binary stars and discuss model predictions, including formation rates, mass distributions, and delay time distributions. We vary assumptions and parameters of physical processes such as mass transfer stability criteria, supernova kick distributions, remnant mass distributions and common-envelope energetics. We compute the likelihood of observing the orbital period-eccentricity distribution of the Galactic DNS population under each of our population synthesis models, allowing us to quantitatively compare the models. We find that mass transfer from a stripped post-helium-burning secondary (case BB) onto a neutron star is most likely dynamically stable. We also find that a natal kick distribution composed of both low (Maxwellian $\sigma=30\rm~km~s^{-1}$) and high ($\sigma=265\rm~km~s^{-1}$) components is preferred over a single high-kick component. We find that the observed DNS mass distribution can place strong constraints on model assumptions.

After the initial fast spiral-in phase experienced by a common-envelope binary, the system may enter a slow, self-regulated phase, possibly lasting 100s of years, in which all the energy released by orbital decay can be efficiently transported to the surface, where it is radiated away. If the remaining envelope is to be removed during this phase, this removal must occur through some as-yet-undetermined mechanism. We carried out 1-d hydrodynamic simulations of a low-mass red giant undergoing a synthetic common-envelope event in such a slow spiral-in phase, using the stellar evolutionary code MESA. We simulated the heating of the envelope due to frictional dissipation from a binary companion's orbit in multiple configurations and investigated the response of the giant's envelope. We find that our model envelopes become dynamically unstable and develop large-amplitude pulsations, with periods in the range 3-20 years and very short growth time-scales of similar order. The shocks and associated rebounds that emerge as these pulsations grow are in some cases strong enough to dynamically eject shells of matter of up to 0.1 $\mathrm{M}_{\odot}$, $\sim 10$ % of the mass of the envelope, from the stellar surface at above escape velocity. These ejections are seen to repeat within a few decades, leading to a time-averaged mass-loss rate of order $10^{-3}$ $\mathrm{M}_{\odot} \: \mathrm{yr}^{-1}$ which is sufficiently high to represent a candidate mechanism for removing the entire envelope over the duration of the slow spiral-in phase.

We present spectra of the extreme polar AR Ursae Majoris (AR UMa) which display a clear Al I absorption doublet, alongside spectra taken less than a year earlier in which that feature is not present. Re-examination of earlier SDSS spectra indicates that the Al I absorption doublet was also present $\approx$8 years before our first non-detection. We conclude that this absorbing material is unlikely to be on the surface of either the white dwarf (WD) or the donor star. We suggest that this Al I absorption feature arises in circumstellar material, perhaps produced by the evaporation of asteroids as they approach the hot WD. The presence of any remaining reservoir of rocky material in AR UMa might help to constrain the prior evolution of this unusual binary system. We also apply spectral decomposition to find the stellar parameters of the M dwarf companion, and attempt to dynamically measure the mass of the WD in AR UMa by considering both the radial velocity curves of the H$_\beta$ emission line and the Na I absorption line. Thereby we infer a mass range for the WD in AR UMa of 0.91 $M_{\odot}$ $<$ $M_{\mathrm{WD}}$ $<$ 1.24 $M_{\odot}$.

The formation of relativistic jets by an accreting compact object is one of the fundamental mysteries of astrophysics. While the theory is poorly understood, observations of relativistic jets from systems known as microquasars have led to a well-established phenomenology. Relativistic jets are not expected from sources with soft or supersoft X-ray spectra, although two such systems are known to produce relatively low-velocity bipolar outflows. Here we report optical spectra of an ultraluminous supersoft X-ray source (ULS) in the nearby galaxy M81 (M81 ULS-1) showing blueshifted broad H\alpha\ emission lines, characteristic of baryonic jets with relativistic speeds. The time variable jets have projected velocities ~17 per cent of the speed of light, and seem similar to those in the prototype microquasar SS 433. Such relativistic jets are not expected to be launched from white dwarfs, but an origin from a black hole or neutron star in M81 ULS-1 is hard to reconcile with its constant soft X-rays. The completely unexpected presence of relativistic jets in a ULS challenges the canonical theories for jet formation, but may possibly be explained by a long speculated super-critically accreting black hole with optically thick outflows

We investigate a scenario in which feedback from black-hole X-ray binaries (BHXBs) sometimes begins inside young star clusters before strong supernova feedback. Those BHXBs could reduce the gas fraction inside embedded young clusters whilst maintaining virial equilibrium, which may help globular clusters (GCs) to stay bound when supernova-driven gas ejection subsequently occurs. Adopting a simple toy model with parameters guided by BHXB population models, we produce GC formation efficiencies consistent with empirically-inferred values. The metallicity dependence of BHXB formation could naturally explain why GC formation efficiency is higher at lower metallicity. For reasonable assumptions about that metallicity dependence, our toy model can produce a GC metallicity bimodality in some galaxies without a bimodality in the field-star metallicity distribution.

We consider X-ray binaries (XBs) as potential sources of stellar feedback. XBs observationally appear able to deposit a high fraction of their power output into their local interstellar medium, which may make them a non-negligible source of energy input. The formation rate of the most luminous XBs rises with decreasing metallicity, which should increase their significance during galaxy formation in the early universe. We also argue that stochastic effects are important to XB feedback (XBF) and may dominate the systematic changes due to metallicity in many cases. Large stochastic variation in the magnitude of XBF at low absolute star formation rates provides a natural reason for diversity in the evolution of dwarf galaxies which were initially almost identical, with several percent of such halos experiencing energy input from XBs roughly two orders of magnitude above the most likely value. These probability distributions suggest that the effect of XBF is most commonly significant for total stellar masses between ~10^7 and 10^8 Msun, which might resolve a current problem with modelling populations of such galaxies. We explain how XBs might inject energy before luminous supernovae (SNe) contribute significantly to feedback and how XBs can assist in keeping gas hot long after the last core-collapse SN has exploded. [...] XBF could be especially important to some dwarf galaxies, potentially heating gas without expelling it; the properties of XBF also match those previously derived as allowing episodic star formation. We also argue that the efficiency of SN feedback (SNF) might be reduced when XBF has had the opportunity to act first. In addition, we note that the effect of SNF is unlikely to be scale-free; galaxies smaller than ~100 pc might well experience less effective SNF. (Slightly abbreviated to fit arXiv size limit.)

Evidence suggests that the direct progenitor stars of some core-collapse supernovae (CCSNe) are luminous blue variables (LBVs), perhaps including some `superluminous supernovae' (SLSNe). We examine models in which massive stars gain mass soon after the end of core hydrogen burning. These are mainly intended to represent mergers following a brief contact phase during early Case B mass transfer, but may also represent stars which gain mass in the Hertzsprung Gap or extremely late during the main-sequence phase for other reasons. The post-accretion stars spend their core helium-burning phase as blue supergiants (BSGs), and many examples are consistent with being LBVs at the time of core collapse. Other examples are yellow supergiants at explosion. We also investigate whether such post-accretion stars may explode successfully after core collapse. The final core properties of post-accretion models are broadly similar to those of single stars with the same initial mass as the pre-merger primary star. More surprisingly, when early Case B accretion does affect the final core properties, the effect appears likely to favour a successful SN explosion, i.e., to make the core properties more like those of a lower-mass single star. However, the detailed structures of these cores sometimes display qualitative differences to any single-star model we have calculated. The rate of appropriate binary mergers may match the rate of SNe with immediate LBV progenitors; for moderately optimistic assumptions we estimate that the progenitor birthrate is ~1% of the CCSN rate.

The energy budget in common-envelope events (CEEs) is not well understood, with substantial uncertainty even over to what extent the recombination energy stored in ionised hydrogen and helium might be used to help envelope ejection. We investigate the reaction of a red-giant envelope to heating which mimics limiting cases of energy input provided by the orbital decay of a binary during a CEE, specifically during the post-plunge-in phase during which the spiral-in has been argued to occur on a time-scale longer than dynamical. We show that the outcome of such a CEE depends less on the total amount of energy by which the envelope is heated than on how rapidly the energy was transferred to the envelope and on where the envelope was heated. The envelope always becomes dynamically unstable before receiving net heat energy equal to the envelope's initial binding energy. We find two types of outcome, both of which likely lead to at least partial envelope ejection: "runaway" solutions in which the expansion of the radius becomes undeniably dynamical, and superficially "self-regulated" solutions, in which the expansion of the stellar radius stops but a significant fraction of the envelope becomes formally dynamically unstable. Almost the entire reservoir of initial helium recombination energy is used for envelope expansion. Hydrogen recombination is less energetically useful, but is nonetheless important for the development of the dynamical instabilities. However, this result requires the companion to have already plunged deep into the envelope; therefore this release of recombination energy does not help to explain wide post-common-envelope orbits.

The potential importance of the angular momentum which is gained by accreting white dwarfs (WDs) has been increasingly recognized in the context of type Ia supernova (SN Ia) single-degenerate model. The expectation that the spin of the WD can delay the explosion should help the single-degenerate model to be consistent with the observed properties of most SNe Ia, in particular by avoiding hydrogen contamination. In this article, we attempt to study the most prominent single-degenerate supersoft (WD + MS) channel when the rotation of accreting WDs is considered. We present a detailed binary population synthesis study to examine the predicted population of SNe Ia for this channel. For our standard model, we find that 77% of these SNe Ia explode with WD masses which are low enough to be supported by solid-body rotation (1.378-1.5Msun); this is a substantially higher proportion than found by previous work. Only 2% have WD explosion masses >2.0Msun; these require the initial WD mass to be larger than 1.0 Msun. We further discuss the possible origin of the diversity of SNe Ia from the pre- and post- accretion properties of the WDs in this population. We also suggest that some SN Ia progenitors with substantial circumstellar hydrogen, including some apparent type IIn SNe, might be related to WDs which required support from differential rotation to avoid explosion, since these can still be accreting from hydrogen-rich donors with a relatively high mass-transfer rate at the time of the SN explosion.

We estimate the fraction of F,G,K stars with close binary companions by analysing multi-epoch stellar spectra from SDSS and LAMOST for radial velocity (RV) variations. We employ a Bayesian method to infer the maximum likelihood of the fraction of binary stars with orbital periods of 1000 days or shorter, assuming a simple model distribution for a binary population with circular orbits. The overall inferred fraction of stars with such a close binary companion is 43.0% \pm 2.0% for a sample of F, G, K stars from SDSS SEGUE, and 30% \pm 8.0% in a similar sample from LAMOST. The apparent close binary fraction decreases with the stellar effective temperature. We divide the SEGUE and LEGUE data into three subsamples with different metallicity ([Fe/H] < -1.1; -1.1 < [Fe/H] < -0.6; -0.6 < [Fe/H]), for which the inferred close binary fractions are 56% \pm 5.0%, 56.0% \pm 3%, and 30% \pm 5.7%. The metal-rich stars from our sample are therefore substantially less likely to possess a close binary companion than otherwise similar stars drawn from metal-poor populations. The different ages and formation environments of the Milky Way's thin disk, thick disk and halo may contribute to explaining these observations. Alternatively metallicity may have a significant effect on the formation and/or evolution of binary stars.

Secondary stars in cataclysmic variables (CVs) follow a well defined period-density relation. Thus, canonical donor stars in CVs are generally low-mass stars of spectral type M. However, several CVs have been observed containing secondary stars which are too hot for their inferred masses. This particular configuration can be explained if the donor stars in these systems underwent significant nuclear evolution before they reached contact. In this paper we present SDSSJ001153.08-064739.2 as an additional example belonging to this peculiar type of CV and discuss in detail its evolutionary history. We perform spectroscopic and photometric observations and make use of available Catalina Real-Time Transient Survey photometry to measure the orbital period of SDSSJ001153.08-064739.2 as 2.4 hours and estimate the white dwarf (Mwd>0.65Msun) and donor star (0.21Msun<Mdon<0.45Msun) masses, the mass ratio (q = 0.32 +- 0.08), the orbital inclination (47 degrees < i < 70 degrees), derive an accurate orbital ephemeris (T0 = 2453383.578 + E x 0.10028081), and report the detection of an outburst. We show that SDSSJ001153.08-064739.2 is one of the most extreme cases in which the donor star is clearly too hot for its mass. SDSSJ001153.08-064739.2 is therefore not only a peculiar CV containing an evolved donor star but also an accreting CV within the period gap. Intriguingly, approximately half of the total currently-observed sample of these peculiar CVs are located in the period gap with nearly the same orbital period.

There are two proposed explanations for ultraluminous X-ray sources (ULXs) with luminosities in excess of $10^{39}$ erg s$^{-1}$. They could be intermediate-mass black holes (more than 100-1,000, solar masses, $M_\odot$) radiating at sub-maximal (sub-Eddington) rates, as in Galactic black-hole X-ray binaries but with larger, cooler accretion disks. Alternatively, they could be stellar-mass black holes radiating at Eddington or super-Eddington rates. On its discovery, M101 ULX-1 had a luminosity of $3\times10^{39}$ erg s$^{-1}$ and a supersoft thermal disk spectrum with an exceptionally low temperature -- uncomplicated by photons energized by a corona of hot electrons -- more consistent with the expected appearance of an accreting intermediate-mass black hole. Here we report optical spectroscopic monitoring of M101 ULX-1. We confirm the previous suggestion that the system contains a Wolf-Rayet star, and reveal that the orbital period is 8.2 days. The black hole has a minimum mass of $5M_\odot$, and more probably a mass of $20-30M_\odot$, but we argue that it is very unlikely to be an intermediate-mass black hole. Therefore its exceptionally soft spectra at high Eddinton ratios violate the expectations for accretion onto stellar-mass black holes. Accretion must occur from captured stellar wind, which has hitherto been thought to be so inefficient that it could not power an ultraluminous source.

It has recently been proposed that one sub-class of type Ia supernovae (SNe Ia) is sufficiently both distinct and common to be classified separately from the bulk of SNe Ia, with a suggested class name of "type Iax supernovae" (SNe Iax), after SN 2002cx. However, their progenitors are still uncertain. We study whether the population properties of this class might be understood if the events originate from a subset of sub-Chandrasekhar mass explosions. In this potential progenitor population, a carbon--oxygen white dwarf (CO WD) accumulates a helium layer from a non-degenerate helium star; ignition of that helium layer then leads to ignition of the CO WD. We incorporated detailed binary evolution calculations for the progenitor systems into a binary population synthesis model to obtain rates and delay times for such events. The predicted Galactic event rate of these explosions is ~1.5\times10^{-3}{yr}^{-1} according to our standard model, in good agreement with the measured rates of SNe Iax. In addition, predicted delay times are ~70Myr-800Myr, consistent with the fact that most of SNe Iax have been discovered in late-type galaxies. If the explosions are assumed to be double-detonations -- following current model expectations -- then based on the CO WD masses at explosion we also estimate the distribution of resulting SN brightness (-13 \gtrsim M_{bol} \gtrsim -19mag), which can reproduce the empirical diversity of SNe Iax. We speculate on why binaries with non-degenerate donor stars might lead to SNe Iax if similar systems with degenerate donors do not. We suggest that the high mass of the helium layer necessary for ignition at the lower accretion rates typically delivered from non-degenerate donors might be necessary to produce SN 2002cx-like characteristics, perhaps even by changing the nature of the CO ignition.

It has recently been proposed that one sub-class of type Ia supernovae (SNe Ia) is sufficiently both distinct and common to be classified separately from the bulk of SNe Ia, with a suggested class name of "type Iax supernovae" (SNe Iax), after SN 2002cx. We show that the population properties of this class can be understood if the events originate from helium double-detonation sub-Chandrasekhar mass explosions, in which a carbon--oxygen white dwarf (CO WD) accumulates a helium layer from a non-degenerate helium star. We have incorporated detailed binary evolution calculations for the progenitor systems into a binary population synthesis model to obtain birthrates and delay times for such events. The predicted Galactic event rate is $\sim$$0.6-1.8\times 10^{-3}\,{\rm yr}^{-1}$, in good agreement with the measured rates of SNe Iax. In addition, predicted delay times are $\sim$70\,Myr$-$710\,Myr, consistent with the fact that SNe Iax have so far only been discovered in late-type galaxies. Based on the CO WD mass at explosion and previous detonation models, we also estimate the distribution of resulting SN brightness ($-13 \gtrsim M_{\rm bol} \gtrsim -19$ mag), which can reproduce the empirical diversity of SNe Iax. We speculate on why binaries with non-degenerate donor stars might lead to SNe Iax if similar systems with degenerate donors do not and suggest that the higher mass of the helium layer necessary for ignition at the lower accretion rates typically delivered from non-degenerate donors might provide the explanation.

Common-envelope events (CEEs), during which two stars temporarily orbit within a shared envelope, are believed to be vital for the formation of a wide range of close binaries. For decades, the only evidence that CEEs actually occur has been indirect, based on the existence of systems that could not be otherwise explained. Here we propose a direct observational signature of CEE arising from a physical model where emission from matter ejected in a CEE is controlled by a recombination front as the matter cools. The natural range of timescales and energies from this model, as well the expected colors, light-curve shapes, ejection velocities and event rate, match those of a recently-recognized class of red transient outbursts.

This work aims to present our current best physical understanding of common-envelope evolution (CEE). We highlight areas of consensus and disagreement, and stress ideas which should point the way forward for progress in this important but long-standing and largely unconquered problem. Unusually for CEE-related work, we mostly try to avoid relying on results from population synthesis or observations, in order to avoid potentially being misled by previous misunderstandings. As far as possible we debate all the relevant issues starting from physics alone, all the way from the evolution of the binary system immediately before CEE begins to the processes which might occur just after the ejection of the envelope. In particular, we include extensive discussion about the energy sources and sinks operating in CEE, and hence examine the foundations of the standard energy formalism. Special attention is also given to comparing the results of hydrodynamic simulations from different groups and to discussing the potential effect of initial conditions on the differences in the outcomes. We compare current numerical techniques for the problem of CEE and also whether more appropriate tools could and should be produced (including new formulations of computational hydrodynamics, and attempts to include 3D processes within 1D codes). Finally we explore new ways to link CEE with observations. We compare previous simulations of CEE to the recent outburst from V1309 Sco, and discuss to what extent post-common-envelope binaries and nebulae can provide information, e.g. from binary eccentricities, which is not currently being fully exploited.

The lack of hydrogen in spectra of type Ia supernovae (SN Ia) is often seen as troublesome for single-degenerate (SD) progenitor models. We argue that, since continued accretion of angular momentum can prevent explosion of the white dwarf, it may be natural for the donor stars in SD progenitors of SN Ia to exhaust their envelopes and shrink rapidly before the explosion. This outcome seems most likely for SD SN Ia progenitors where mass-transfer begins from a giant donor star, and might extend to other SD systems. Not only is the amount of hydrogen left in such a system below the current detection limit, but the donor star is typically orders of magnitude smaller than its Roche lobe by the point when a SD SN Ia occurs, in which case attempts to observe collisions between SN shocks and giant donor stars seem unlikely to succeed. We consider the constraints on this model from the circumstellar structures seen in spectra of SN 2006X and suggest a novel explanation for the origin of this material.

Type Ia supernovae are generally thought to be due to the thermonuclear explosions of carbon-oxygen white dwarfs with masses near the Chandrasekhar mass. This scenario, however, has two long-standing problems. First, the explosions do not naturally produce the correct mix of elements, but have to be finely tuned to proceed from sub-sonic deflagration to super-sonic detonation. Second, population models and observations give formation rates of near-Chandrasekhar white dwarfs that are far too small. Here, we suggest that type Ia supernovae instead result from mergers of roughly equal-mass carbon-oxygen white dwarfs, including those that produce sub-Chandrasekhar mass remnants. Numerical studies of such mergers have shown that the remnants consist of rapidly rotating cores that contain most of the mass and are hottest in the center, surrounded by dense, small disks. We argue that the disks accrete quickly, and that the resulting compressional heating likely leads to central carbon ignition. This ignition occurs at densities for which pure detonations lead to events similar to type Ia supernovae. With this merger scenario, we can understand the type Ia rates, and have plausible reasons for the observed range in luminosity and for the bias of more luminous supernovae towards younger populations. We speculate that explosions of white dwarfs slowly brought to the Chandrasekhar limit---which should also occur---are responsible for some of the "atypical" type Ia supernovae.

We propose a formation channel for the previously unexplained helium-rich subdwarf O (He-rich sdO) stars in which post-subdwarf B (sdB) stars (i.e. hybrid COHe white dwarfs) reignite helium burning in a shell after gaining matter from their helium white-dwarf (WD) companions. Such short-period binaries containing post-sdB WDs and helium WDs are predicted by one of the major binary formation channels for sdB stars. In the majority of cases, mass transfer is expected to lead to a dynamically unstable merger event, leaving a single-star remnant. Calculations of the evolution of these stars show that their properties are consistent with the observed He-rich sdO stars. The luminosity of these stars is about an order of magnitude higher than that of canonical sdB stars. We also suggest that binary systems such as PG 1544+488 (Ahmad et al. 2004) and HE 0301-3039 (Lisker et al. 2004), which each contain two hot subdwarfs, could be the outcome of a double-core common-envelope phase. Since this favours intermediate-mass progenitors, this may also explain why the subdwarfs in these systems are He-rich.

Among the initial results from Kepler were two striking lightcurves, for KOI 74 and KOI 81, in which the relative depths of the primary and secondary eclipses showed that the more compact, less luminous object was hotter than its stellar host. That result became particularly intriguing because a substellar mass had been derived for the secondary in KOI 74, which would make the high temperature challenging to explain; in KOI 81, the mass range for the companion was also reported to be consistent with a substellar object. We re-analyze the Kepler data and demonstrate that both companions are likely to be white dwarfs. We also find that the photometric data for KOI 74 show a modulation in brightness as the more luminous star orbits, due to Doppler boosting. The magnitude of the effect is sufficiently large that we can use it to infer a radial velocity amplitude accurate to 1 km/s. As far as we are aware, this is the first time a radial-velocity curve has been measured photometrically. Combining our velocity amplitude with the inclination and primary mass derived from the eclipses and primary spectral type, we infer a secondary mass of 0.22+/-0.03 Msun. We use our estimates to consider the likely evolutionary paths and mass-transfer episodes of these binary systems.

We present a new mechanism for the ejection of a common envelope in a massive binary, where the energy source is nuclear energy rather than orbital energy. This can occur during the slow merger of a massive primary with a secondary of 1-3 Msun when the primary has already completed helium core burning. We show that, in the final merging phase, hydrogen-rich material from the secondary can be injected into the helium-burning shell of the primary. This leads to a nuclear runaway and the explosive ejection of both the hydrogen and the helium layer, producing a close binary containing a CO star and a low-mass companion. We argue that this presents a viable scenario to produce short-period black-hole binaries and long-duration gamma-ray bursts (LGRBs). We estimate a LGRB rate of about 1.e-6 per year at solar metallicity, which implies that this may account for a significant fraction of all LGRBs, and that this rate should be higher at lower metallicity.

We study the properties of a sample of 211 heavily-obscured Active Galactic Nucleus (AGN) candidates in the Extended Chandra Deep Field-South selecting objects with f_24/f_R>1000 and R-K>4.5. Of these, 18 were detected in X-rays and found to be obscured AGN with neutral hydrogen column densities of ~10^23 cm^-2. In the X-ray undetected sample, the following evidence suggests a large fraction of heavily-obscured (Compton Thick) AGN: (i) The stacked X-ray signal of the sample is strong, with an observed ratio of soft to hard X-ray counts consistent with a population of ~90% heavily obscured AGN combined with 10% star-forming galaxies. (ii) The X-ray to mid-IR ratios for these sources are significantly larger than that of star-forming galaxies and ~2 orders of magnitude smaller than for the general AGN population, suggesting column densities of N_H>5x10^24 cm^-2. (iii) The Spitzer near- and mid-IR colors of these sources are consistent with those of the X-ray-detected sample if the effects of dust self-absorption are considered. Spectral fitting to the rest-frame UV/optical light (dominated by the host galaxy) returns stellar masses of ~10^11 M_sun and <E(B-V)> =0.5, and reveals evidence for a significant young stellar population, indicating that these sources are experiencing considerable star-formation. This sample of heavily-obscured AGN candidates implies a space density at z~2 of ~10^-5 Mpc^-3, finding a strong evolution in the number of L_X>10^44 erg/s sources from z=1.5 to 2.5, possibly consistent with a short-lived heavily-obscured phase before an unobscured quasar is visible.

There is still considerable debate over the progenitors of type Ia supernovae (SNe Ia). Likewise, it is not agreed how single white dwarfs with masses less than ~0.5 Msun can be formed in the field, even though they are known to exist. We consider whether single low-mass white dwarfs (LMWDs) could have been formed in binary systems where their companions have exploded as a SN Ia. In this model, the observed single LMWDs are the remnants of giant-branch donor stars whose envelopes have been stripped off by the supernova explosion. We investigate the likely remnants of SNe Ia, including the effects of the explosion on the envelope of the donor star. We also use evolutionary arguments to examine alternative formation channels for single LMWDs. In addition, we calculate the expected kinematics of the potential remnants of SNe Ia. SN Ia in systems with giant-branch donor stars can naturally explain the production of single LMWDs. It seems difficult for any other formation mechanism to account for the observations, especially for those single LMWDs with masses less than ~0.4 Msun. Independent of those results, we find that the kinematics of one potentially useful population containing single LMWDs is consistent with our model. Studying remnant white-dwarf kinematics seems to be a promising way to investigate SN Ia progenitors. The existence of single LMWDs appears to constitute evidence for the production of SNe Ia in binary systems with a red-giant donor star. Other single white dwarfs with higher space velocities support a second, probably dominant, population of SN Ia progenitors which contained main-sequence or subgiant donor stars at the time of explosion. The runaway stars LP400-22 and US 708 suggest the possibility of a third formation channnel for some SNe Ia in systems where the donor stars are hot subdwarfs.

Massive stars undergo a violent death when the supply of nuclear fuel in their cores is exhausted, resulting in a catastrophic "core-collapse" supernova. Such events are usually only detected at least a few days after the star has exploded. Observations of the supernova SNLS-04D2dc with the Galaxy Evolution Explorer space telescope reveal a radiative precursor from the supernova shock before the shock reached the surface of the star and show the initial expansion of the star at the beginning of the explosion. Theoretical models of the ultraviolet light curve confirm that the progenitor was a red supergiant, as expected for this type of supernova. These observations provide a way to probe the physics of core-collapse supernovae and the internal structures of their progenitor stars

The nature of the companion stars in Type Ia Supernova (SNe) progenitor systems remains unclear. One possible way to discriminate between different scenarios is the presence (or absence) of circumstellar material, the left overs from the progenitor evolution that may be revealed by their interaction with the SN. A new method to probe the circumstellar environment has been exploited for the normal Type Ia SN 2006X, leading for the first time to the direct detection of material which escaped the progenitor system. In this paper we apply the same analysis to the peculiar Type Ia SN 2000cx, with the aim of constraining the properties of its progenitor system. Using multi-epoch, high-resolution spectroscopy we have studied the spectral region where narrow, time-variable NaID absorption features are expected in case circumstellar material is present along the line of sight. No NaID absorption is detected in the rest-frame of the host galaxy to a level of a few mA, setting a stringent upper limit to the column density of the absorbing material (N(NaI)<=2 x 10^10 cm-2). In this respect the peculiar Type Ia SN 2000cx is different from the normal Ia SN 2006X. Whether this is to be attributed to a different progenitor system, to viewing-angle effects or to a low metallicity remains to be clarified.

Type Ia supernovae are thought to be thermonuclear explosions of accreting white dwarfs that reach a critical mass limit. Despite their importance as cosmological distance indicators, the nature of their progenitors has remained controversial. Here we report the detection of circumstellar material in a normal Type Ia supernova. The expansion velocities, densities and dimensions of the circumstellar envelope indicate that this material was ejected from the progenitor system. The relatively low expansion velocities appear to favor a progenitor system where a white dwarf accretes material from a companion star which is in the red-giant phase at the time of explosion.

Using a comprehensive binary population synthesis scheme, we investigate the statistical properties of a sample of eclipsing binaries that is detectable by an idealised extrasolar planet transit survey with specifications broadly similar to those of the SuperWASP (Wide Angle Search for Planets) project. In this idealised survey the total number of detectable single stars in the Galactic disc is of the order of 10^6-10^7, while, for a flat initial mass ratio distribution, the total number of detectable eclipsing binaries is of the order of 10^4-10^5. The majority of the population of detectable single stars is made up of main-sequence stars (60%), horizontal-branch stars (20%), and giant-branch stars (10%). The largest contributions to the population of detectable eclipsing binaries stem from detached double main-sequence star binaries (60%), detached giant-branch main-sequence star binaries (20%), and detached horizontal-branch main-sequence star binaries (10%). The ratio of the number of eclipsing binaries to the number of single stars detectable by the idealised SuperWASP survey varies by less than a factor of 2.5 across the sky, and decreases with increasing Galactic latitude. It is found to be largest in the direction of the Galactic longitude l=-7.5deg and the Galactic latitude b=-22.5deg. We also show that the fractions of systems in different subgroups of eclipsing binaries are sensitive to the adopted initial mass ratio distribution, which is one of the poorest constrained input parameters in present-day binary population synthesis calculations. This suggests that once statistically meaningful results from transit surveys are available, they will be able to significantly improve the predictive power of population synthesis studies of interacting binaries and related objects. (abridged)

We examine the proposal that the subset of neutron-star and black-hole X-ray binaries that form with Ap or Bp star companions will experience systemic angular-momentum losses due to magnetic braking, not otherwise operative with intermediate-mass companion stars. We suggest that for donor stars possessing the anomalously high magnetic fields associated with Ap and Bp stars, a magnetically coupled, irradiation-driven stellar wind can lead to substantial systemic loss of angular-momentum. In this paper we apply this mechanism to a specific astrophysics problem involving the formation of compact black-hole binaries with low-mass donor stars. At present, it is not understood how these systems form, given that low-mass companion stars are not likely to provide sufficient gravitational potential to unbind the envelope of the massive progenitor of the black hole during a prior `common-envelope' phase. However, in the absence of magnetic braking, such systems tend to evolve to long orbital periods. We show that, with the proposed magnetic braking properties afforded by Ap and Bp companions, such a scenario can lead to the formation of compact black-hole binaries with orbital periods, donor masses, lifetimes, and production rates that are in accord with the observations. In spite of these successes, our models reveal a significant discrepancy between the calculated effective temperatures and the observed spectral types of the donor stars. Finally, we show that this temperature discrepancy would still exist for other scenarios invoking initially intermediate-mass donor stars, and this presents a substantial unresolved mystery.

We have computed models for ultraluminous X-ray sources ("ULXs") consisting of a black-hole accretor of intermediate mass ("IMBH"; e.g., ~1000 Msun) and a captured donor star. For each of four different sets of initial donor masses and orbital separations, we computed 30,000 binary evolution models using a full Henyey stellar evolution code. To our knowledge this is the first time that a population of X-ray binaries this large has been carried out with other than approximation methods, and it serves to demonstrate the feasibility of this approach to large-scale population studies of mass-transfer binaries. In the present study, we find that in order to have a plausible efficiency for producing active ULX systems with IMBHs having luminosities > 10^{40} ergs/sec, there are two basic requirements for the capture of companion/donor stars. First, the donor stars should be massive, i.e., > 8 Msun. Second, the initial orbital separations, after circularization, should be close, i.e., < 6-30 times the radius of the donor star when on the main sequence. Even under these optimistic conditions, we show that the production rate of IMBH-ULX systems may fall short of the observed values by factors of 10-100.