Most massive stars exchange mass with a companion, leading to evolution which is altered drastically from that expected of stars in isolation. Such systems are the result of unusual binary evolution pathways and, as such, may be used to place stringent constraints on the physics of these interactions. We use the R4 system's B[e] supergiant, which has been postulated to be the product of a binary stellar merger, to guide our understanding of such outcomes by comparing observations of R4 to the results of simulations of mergers performed with the 3d hydrodynamics code FLASH. Our approach tailors the simulation initial conditions to the observed properties of R4 and implements realistic stellar profiles generated by the 1d stellar evolution code MESA onto the 3d grid, resolving the merger inspiral to within $0.02\, R_{\odot}$. We then map the merger remnant into MESA to track its evolution on the HR diagram over a period of $10^4$ years. This generates models for a B[e] supergiant with stellar properties, age, and nebula structure in qualitative agreement with that of the R4 system. Our calculations provide concrete evidence to support the idea that R4 was originally a member of a triple system in which the inner binary merged after its most massive member evolved off the main sequence, producing a new object that is of similar mass yet significantly more luminous than the A supergiant companion. The potential applications of the code framework presented in this paper are wide ranging and can be used to generate models of a variety of merger stellar remnants.

About 10% of stars more massive than ${\approx}\,1.5\,\mathrm{M}_\odot$ have strong, large-scale surface magnetic fields and are being discussed as progenitors of highly-magnetic white dwarfs and magnetars. The origin of these fields remains uncertain. Recent 3D magnetohydrodynamical simulations have shown that strong magnetic fields can be generated in the merger of two massive stars. Here, we follow the long-term evolution of such a 3D merger product in a 1D stellar evolution code. During a thermal relaxation phase after the coalescence, the merger product reaches critical surface rotation, sheds mass and then spins down primarily because of internal mass readjustments. The spin of the merger product after thermal relaxation is mainly set by the co-evolution of the star-torus structure left after coalescence. This evolution is still uncertain, so we also consider magnetic braking and other angular-momentum-gain and -loss mechanisms that may influence the final spin of the merged star. Because of core compression and mixing of carbon and nitrogen in the merger, enhanced nuclear burning drives a transient convective core that greatly contributes to the rejuvenation of the star. Once the merger product relaxed back to the main sequence, it continues its evolution similar to that of a genuine single star of comparable mass. It is a slow rotator that matches the magnetic blue straggler $\tau$ Sco. Our results show that merging is a promising mechanism to explain some magnetic massive stars and it may also be key to understand the origin of the strong magnetic fields of highly-magnetic white dwarfs and magnetars.

Be X-ray binaries (BeXRBs) consist of rapidly rotating Be stars with neutron star companions accreting from the circumstellar emission disk. We compare the observed population of BeXRBs in the Small Magellanic Cloud with simulated populations of BeXRB-like systems produced with the COMPAS population synthesis code. We focus on the apparently higher minimal mass of Be stars in BeXRBs than in the Be population at large. Assuming that BeXRBs experienced only dynamically stable mass transfer, their mass distribution suggests that at least $\sim$ 30% of the mass donated by the progenitor of the neutron star is typically accreted by the B-star companion. We expect these results to affect predictions for the population of double compact object mergers. A convolution of the simulated BeXRB population with the star formation history of the Small Magellanic Cloud shows that the excess of BeXRBs is most likely explained by this galaxy's burst of star formation $\sim$ 40 Myr ago, rather than by its low metallicity.

During the first three observing runs of the Advanced gravitational-wave detector network, the LIGO/Virgo collaboration detected several black hole binary (BHBH) mergers. As the population of detected BHBH mergers grows, it will become possible to constrain different channels for their formation. Here we consider the chemically homogeneous evolution (CHE) channel in close binaries, by performing population synthesis simulations that combine realistic binary models with detailed cosmological calculations of the chemical and star-formation history of the Universe. This allows us to constrain population properties, as well as cosmological and aLIGO detection rates of BHBH mergers formed through this pathway. We predict a BHBH merger rate at redshift zero of $5.8 \hspace{1mm} \textrm{Gpc}^{-3} \textrm{yr}^{-1}$ through the CHE channel, to be compared with aLIGO's measured rate of ${53.2}_{-28.2}^{+55.8} \hspace{1mm} \text{Gpc}^{-3}\text{yr}^{-1}$, and find that eventual merger systems have BH masses in the range $17 - 43 \hspace{1mm} \textrm{M}_{\odot}$ below the pair-instability supernova (PISN) gap, and $>124 \hspace{1mm} \textrm{M}_{\odot}$ above the PISN gap. We further investigate the effects of momentum kicks during black hole formation, calculate cosmological and magnitude limited PISN rates and investigate the effects of high-redshift deviations in the star formation rate. We find that momentum kicks tend to increase delay times of BHBH systems, and our magnitude limited PISN rate estimates indicate that current deep surveys should be able to detect such events. Lastly, we find that our cosmological merger rate estimates change by at most $\sim 8\%$ for mild deviations of the star formation rate in the early Universe, and by up to $\sim 40\%$ for extreme deviations.

Type Ibn supernovae (SNe Ibn) are intriguing stellar explosions whose spectra exhibit narrow helium lines with little hydrogen. They trace the presence of circumstellar material (CSM) formed via pre-SN eruptions of their stripped-envelope progenitors. Early work has generally assumed that SNe Ibn come from massive Wolf-Rayet (WR) stars via single star evolution. In this paper, we report ultraviolet (UV) and optical observations of two nearby Type Ibn SNe 2006jc and 2015G conducted with the Hubble Space Telescope (HST) at late times. A point source is detected at the position of SN 2006jc, and we confirm the conclusion of Maund et al. that it is the progenitor's binary companion. Its position on the Hertzsprung-Russell (HR) diagram corresponds to a star that has evolved off the main sequence (MS); further analysis implies a low initial mass for the companion star ($M_2$ $\le$ 11.9$^{+1.2}_{-0.8}$ $M_\odot$) and a secondary-to-primary initial mass ratio very close to unity ($q$ = $M_2/M_1$ $\sim$ 1); the SN progenitor's hydrogen envelope had been stripped through binary interaction. We do not detect the binary companion of SN 2015G. For both SNe, the surrounding stellar populations have relatively old ages and argue against any massive WR stars as their progenitors. These results suggest that SNe Ibn may have lower-mass origins in interacting binaries. As a result, they also provide evidence that the giant eruptions commonly seen in massive luminous blue variables (LBVs) can also occur in much lower-mass, stripped-envelope stars just before core collapse.

Magnetic fields are ubiquitous in the Universe. The Sun's magnetic field drives the solar wind and causes solar flares and other energetic surface phenomena that profoundly affect space weather here on Earth. The first magnetic field in a star other than the Sun was detected in 1947 in the peculiar A-type star 78 Vir. It is now known that the magnetic fields of the Sun and other low-mass stars (<1.5 solar masses) are generated in-situ by a dynamo process in their turbulent, convective envelopes. Unlike such stars, intermediate-mass and high-mass stars (>1.5 solar masses; referred to as "massive" stars here) have relatively quiet, radiative envelopes where a solar-like dynamo cannot operate. However, about 10% of them, including 78 Vir, have strong, large-scale surface magnetic fields whose origin has remained a major mystery. The massive star $\tau$ Sco is a prominent member of this group and appears to be surprisingly young compared to other presumably coeval members of the Upper Scorpius association. Here, we present the first 3D magneto-hydrodynamical simulations of the coalescence of two massive main-sequence stars and 1D stellar evolution computations of the subsequent evolution of the merger product that can explain $\tau$ Sco's magnetic field, apparent youth and other observed characteristics. We argue that field amplification in stellar mergers is a general mechanism to form strongly-magnetised massive stars. These stars are promising progenitors of those neutron stars that host the strongest magnetic fields in the Universe, so-called magnetars, and that may give rise to some of the enigmatic fast radio bursts. Strong magnetic fields affect the explosions of core-collapse supernovae and, moreover, those magnetic stars that have rapidly-rotating cores at the end of their lives might provide the right conditions to power long-duration gamma-ray bursts and super-luminous supernovae.

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.

Recently, the black hole X-ray binary (BHXB) Nova Muscae 1991 has been reported to be experiencing an extremely rapid orbital decay. So far, three BHXBs have anomalously high orbital period derivatives, which can not be interpreted by the standard stellar evolution theory. In this work, we investigate whether the resonant interaction between the binary and a surrounding circumbinary (CB) disk could produce the observed orbital period derivatives. Analytical calculations indicate that the observed orbital period derivatives of XTE J1118+480 and A0620-00 can originate from the tidal torque between the binary and a CB disk with a mass of $10^{-9}~\rm M_{\odot}$, which is approximately in agreement with the dust disk mass detected in these two sources. However, Nova Muscae 1991 was probably surrounded by a heavy CB disk with a mass of $10^{-7}~\rm M_{\odot}$. Based on the CB disk model and the anomalous magnetic braking theory, we simulate the evolution of the three BHXBs with intermediate-mass donor stars by using the MESA code. Our simulated results are approximately consistent with the observed donor star masses, orbital periods, and orbital-period derivatives. However, the calculated effective temperatures of the donor stars are higher than indicated by the observed spectral types of two sources.

Motivated by the eruptive mass loss inferred from Luminous Blue Variable (LBV) stars, we present 1D hydrodynamical simulations of the response from sudden energy injection into the interior of a very massive ($100 \, M_\odot$) star. For a fiducial case with total energy addition set to a factor $f=0.5$ of the net stellar binding energy, and applied within the stellar envelope, we detail the dynamical response that leads to ejection of the outermost $7.2 \, M_\odot$. We find that the ejecta's variations in time $t$ and radius $r$ for the velocity $v$, density $\rho$, and temperature $T$ are quite well fit by similarity forms in the variable $r/t \approx v$. Specifically the scaled density follows a simple exponential decline $\rho t^{3} \sim \exp (-r/v_{\rm o} t)$. This `exponential similarity' leads to analytic scaling relations for total ejecta mass $\Delta M$ and kinetic energy $\Delta K$ that agree well with the hydrodynamical simulations, with the specific-energy-averaged speed related to the exponential scale speed $v_{\rm o}$ through ${\bar v} \equiv \sqrt{2 \Delta K/\Delta M} = \sqrt{12} \, v_{\rm o}$, and a value comparable to the star's surface escape speed, $v_{\rm esc}$. Models with energy added in the core develop a surface shock breakout that propels an initial, higher-speed ejecta ($>$5000km s$^{-1}$), but the bulk of the ejected material still follows the same exponential similarity scalings with ${\bar v} \approx v_{\rm esc}$. A broader parameter study examines how the ejected mass and energy depends on the energy-addition factor $f$, for three distinct model series that locate the added energy in either the core, envelope, or near-surface. We conclude by discussing the relevance of these results for understanding LBV outbursts and other eruptive phenomena, such as failed supernovae and pulsational pair instability events.

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)

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.

We carry out a comprehensive study of supernova ejecta-companion interaction in massive binary systems. We aim to physically understand the kinematics of the interaction and predict observational signatures. To do this we perform simulations over a vast parameter space of binary configurations, varying the masses of the progenitor and companion, structure of the companion, explosion energy, and orbital separation. Our results were not so consistent with classical models by Wheeler et al. 1975, sometimes deviating by an order of magnitude. We construct an alternative simple model which explains the simulated results reasonably well and can be used to estimate impact velocities for arbitrary explosion profiles and companion star structures. We then investigate the long term evolution after the supernova, where the companion can be inflated by the energy injected into the star. We find that the companion can become more than an order of magnitude overluminous straight after the supernova, but quickly fades away after $\sim10$ years and returns to its original luminosity in about a thermal timescale of the star. Finally, we also discuss the possible surface contamination of heavy elements from the slower ejecta.

Farr and Mandel reanalyse our data, finding initial-mass-function slopes for high mass stars in 30 Doradus that agree with our results. However, their reanalysis appears to underpredict the observed number of massive stars. Their technique results in more precise slopes than in our work, strengthening our conclusion that there is an excess of massive stars above $30\,\mathrm{M}_\odot$ in 30 Doradus.

SN 2002cx-like and SN Ia-CSM objects show similar early spectra and both belong to a young stellar population, suggesting that they could share the same progenitor origin. Adopting the framework of the common-envelope-wind (CEW) model developed in \citet{MENGXC17}, we here propose that both subclasses of SNe Ia are caused by the explosion of hybrid carbon-oxygen-neon white dwarfs (CONe WDs) in single-degenerate systems, where SNe Ia-CSM explode in systems with a massive common envelope (CE) of $\sim1~M_{\odot}$, while SN 2002cx-like events correspond to those events where most of the CE has been lost in a wind. Using binary-population-synthesis (BPS) calculations, we estimate a number ratio of SNe Ia-CSM to SN 2002cx-like objects between $\frac{1}{3}$ and $\frac{2}{3}$, consistent with observational constraints, and an overall contribution from hybrid CONe WDs to the total SN Ia population that also matches the observed number from these peculiar objects. Our model predicts a statistical sequence of CSM density from SN Ia-CSM to SN 2002cx-like events and normal SNe Ia, consistent with existing radio constraints. We also find a new subclass of hybrid SNe which share the properties of Type II and Type Ia SNe, consistent with some observed SNe, which do not have a surviving companion. In some cases these could even produce SNe Ia from apparently single WDs.

The 30 Doradus star-forming region in the Large Magellanic Cloud is a nearby analogue of large star-formation events in the distant Universe. We determine the recent formation history and the initial mass function (IMF) of massive stars in 30 Doradus based on spectroscopic observations of 247 stars more massive than 15 solar masses ($\mathrm{M}_\odot$). The main episode of massive star formation started about $8\,\mathrm{Myr}$ ago and the star-formation rate seems to have declined in the last $1\,\mathrm{Myr}$. The IMF is densely sampled up to $200\,\mathrm{M}_\odot$ and contains $32\pm12\%$ more stars above $30\,\mathrm{M}_\odot$ than predicted by a standard Salpeter IMF. In the mass range $15-200\,\mathrm{M}_\odot$, the IMF power-law exponent is $1.90^{+0.37}_{-0.26}$, shallower than the Salpeter value of 2.35.

Gravitational waves were discovered with the detection of binary black hole mergers and they should also be detectable from lower mass neutron star mergers. These are predicted to eject material rich in heavy radioactive isotopes that can power an electromagnetic signal called a kilonova. The gravitational wave source GW170817 arose from a binary neutron star merger in the nearby Universe with a relatively well confined sky position and distance estimate. Here we report observations and physical modelling of a rapidly fading electromagnetic transient in the galaxy NGC4993, which is spatially coincident with GW170817 and a weak short gamma-ray burst. The transient has physical parameters broadly matching the theoretical predictions of blue kilonovae from neutron star mergers. The emitted electromagnetic radiation can be explained with an ejected mass of 0.04 +/- 0.01 Msol, with an opacity of kappa <= 0.5 cm2/gm at a velocity of 0.2 +/- 0.1c. The power source is constrained to have a power law slope of beta = -1.2 +/- 0.3, consistent with radioactive powering from r-process nuclides. We identify line features in the spectra that are consistent with light r-process elements (90 < A < 140). As it fades, the transient rapidly becomes red, and emission may have contribution by a higher opacity, lanthanide-rich ejecta component. This indicates that neutron star mergers produce gravitational waves, radioactively powered kilonovae, and are a nucleosynthetic source of the r-process elements.

Double neutron star (DNS) systems represent extreme physical objects and the endpoint of an exotic journey of stellar evolution and binary interactions. Large numbers of DNS systems and their mergers are anticipated to be discovered using the Square-Kilometre-Array searching for radio pulsars and high-frequency gravitational wave detectors (LIGO/VIRGO), respectively. Here we discuss all key properties of DNS systems, as well as selection effects, and combine the latest observational data with new theoretical progress on various physical processes with the aim of advancing our knowledge on their formation. We examine key interactions of their progenitor systems and evaluate their accretion history during the high-mass X-ray binary stage, the common envelope phase and the subsequent Case BB mass transfer, and argue that the first-formed NSs have accreted at most $\sim 0.02\;M_{\odot}$. We investigate DNS masses, spins and velocities, and in particular correlations between spin period, orbital period and eccentricity. Numerous Monte Carlo simulations of the second supernova (SN) events are performed to extrapolate pre-SN stellar properties and probe the explosions. All known close-orbit DNS systems are consistent with ultra-stripped exploding stars. Although their resulting NS kicks are often small, we demonstrate a large spread in kick magnitudes which may, in general, depend on the past interaction history of the exploding star and thus correlate with the NS mass. We analyze and discuss NS kick directions based on our SN simulations. Finally, we discuss the terminal evolution of close-orbit DNS systems until they merge and possibly produce a short $\gamma$-ray burst.

He accretion onto carbon-oxygen white dwarfs (CO WDs) plays a fundamental role when studying the formation of type Ia supernovae (SNe Ia). Employing the MESA stellar evolution code, we calculated the long-term evolution of He-accreting CO WDs. Previous studies usually supposed that a WD can grow in mass to the Chandrasekhar limit in the stable He burning region and finally produce a SN Ia. However, in this study we find that off-centre carbon ignition occurs in the stable He burning region if the accretion rate is above a critical value (~2.05*10^{-6}{M$_\odot$}{yr}^-1), resulting in accretion-induced collapse rather than a SN Ia. If the accretion rate is below the critical value, explosive carbon ignition will eventually happen in the centre producing a SN Ia. Taking into account the possibility of off-centre carbon ignition, we have re-determined the initial parameter space that produces SNe Ia in the He star donor channel, one of the promising channels to produce SNe Ia in young populations. Since this parameter space is smaller than was found in the previous study of Wang et al. (2009), the SN Ia rates are also correspondingly smaller. We also determined the chemical abundance profile of the He-accreting WDs at the moment of explosive carbon ignition, which can be used as initial input for SN Ia explosion models.

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.

The single-degenerate (SD) model is one of the principal models for the progenitors of type Ia supernovae (SNe Ia), but some of the predictions in the most widely studied version of the SD model, i.e. the optically thick wind (OTW) model, have not been confirmed by observations. Here, we propose a new version of the SD model in which a common envelope (CE) is assumed to form when the mass-transfer rate between a carbon-oxygen white dwarf (CO WD) and its companion exceeds a critical accretion rate. The WD may gradually increase its mass at the base of the CE. Due to the large nuclear luminosity for stable hydrogen burning, the CE may expand to giant dimensions and will lose mass from the surface of the CE by a CE wind (CEW). Because of the low CE density, the binary system will avoid a fast spiral-in phase and finally re-emerge from the CE phase. Our model may share the virtues of the OTW model but avoid some of its shortcomings. We performed binary stellar evolution calculations for more than 1100 close WD + MS binaries. Compared with the OTW model, the parameter space for SNe Ia from our CEW model extends to more massive companions and less massive WDs. Correspondingly, the Galactic birth rate from the CEW model is higher than that from the OTW model by $\sim$30\%. Finally, we discuss the uncertainties of the CEW model and the differences between our CEW model and the OTW model.

Gravitational waves from the binary black hole (BH) merger GW150914 may enlighten our understanding of ultra-luminous X-ray sources (ULXs), as BHs>30Msun can reach luminosities>4x10^39 erg s^-1 without exceeding their Eddington limit. It is then important to study variations of evolutionary channels for merging BHs, which might instead form accreting BHs and become ULXs. It was recently shown that massive binaries with mass ratios close to unity and tight orbits can undergo efficient rotational mixing and evolve chemically homogeneously, resulting in a compact BH binary. We study similar systems by computing ~120000 detailed binary models with the MESA code covering a wide range of initial parameters. For initial mass ratios M2/M1~0.1-0.4, primaries >40Msun can evolve chemically homogeneously, remaining compact and forming a BH without undergoing Roche-lobe overflow. The secondary then expands and transfers mass to the BH, initiating a ULX phase. We predict that ~1 out of 10^4 massive stars evolves this way, and that in the local universe 0.13 ULXs per Msun yr^-1 of star-formation rate are observable, with a strong preference for low-metallicities. At metallicities log Z>-3, BH masses in ULXs are limited to 60Msun due to the occurrence of pair-instability supernovae which leave no remnant, resulting in an X-ray luminosity cut-off. At lower metallicities, very massive stars can avoid exploding as pair-instability supernovae and instead form BHs with masses above 130Msun, producing a gap in the ULX luminosity distribution. After the ULX phase, neutron-star-BH binaries that merge in less than a Hubble time are produced with a formation rate <0.2 Gpc^-3 yr^-1. We expect that upcoming X-ray observatories will test these predictions, which together with additional gravitational wave detections will provide strict constraints on the origin of the most massive BHs that can be produced by stars.

Strong winds from massive stars are a topic of interest to a wide range of astrophysical fields. In High-Mass X-ray Binaries the presence of an accreting compact object on the one side allows to infer wind parameters from studies of the varying properties of the emitted X-rays; but on the other side the accretor's gravity and ionizing radiation can strongly influence the wind flow. Based on a collaborative effort of astronomers both from the stellar wind and the X-ray community, this presentation attempts to review our current state of knowledge and indicate avenues for future progress.

Some short-period, detached binary systems have recently been reported to experience very rapid orbital decay, much faster than is expected from the angular-momentum loss caused by gravitational radiation alone. As these systems contain fully convective stars, magnetic braking is not believed to be operative, making the large orbital-period derivative puzzling. Here we explore whether a resonant interaction between the binary and a surrounding circumbinary (CB) disk could account for the observed orbital decay. Our calculations indicate that the observed orbital-period derivatives in seven detached binaries can be produced by the resonant interaction between the binary and a CB disk if the latter has a mass in the range of $10^{-4}-10^{-2}~ M_{\odot}$, which is of the same order as the inferred disk mass ($\sim2.4\times 10^{-4}~ M_{\odot}$) in the post-common-envelope binary NN Ser.

Since the majority of massive stars are members of binary systems, an understanding of the intricacies of binary interactions is essential for understanding the large variety of supernova types and sub-types. I therefore briefly review the basic elements of binary evolution theory and discuss how binary interactions affect the presupernova structure of massive stars and the resulting supernovae. SN 1987A was a highly anomalous supernova, almost certainly because of a previous binary interaction. The most likely scenario at present is that the progenitor was a member of a massive close binary that experienced dynamical mass transfer during its second red-supergiant phase and merged completely with its companion as a consequence. This can naturally explain the three main anomalies of SN 1987A: the blue color of the progenitor, the chemical anomalies and the complex triple-ring nebula.

From the canonical binary scenario, the majority of sdBs are produced from low-mass stars with degenerate cores where helium is ignited in a way of flashes. Due to numerical difficulties, the models of produced sdBs are generally constructed from more massive stars with non-degenerate cores, leaving several uncertainties on the exact characteristics of sdB stars. Employing MESA, we systematically studied the characteristics of sdBs produced from the common envelope (CE) ejection channel, and found that the sdB stars produced from the CE ejection channel appear to form two distinct groups on the effective temperature-gravity diagram. One group (the flash-mixing model) almost has no H-rich envelope and crows at the hottest temperature end of the extremely horizontal branch (EHB), while the other group has significant H-rich envelope and spreads over the whole canonical EHB region. The key factor for the dichotomy of the sdB properties is the development of convection during the first helium flash, which is determined by the interior structure of the star after the CE ejection. For a given initial stellar mass and a given core mass at the onset of the CE, if the CE ejection stops early, the star has a relatively massive H-rich envelope, resulting in a canonical sdB generally. The fact of only a few short-orbital-period sdB binaries being in the flash-mixing sdB region and the lack of He-rich sdBs in short-orbital-period binaries indicate that the flash mixing is not very often in the products of the CE ejection. A falling back process after the CE ejection, similar to that happened in nova, is an appropriate way of increasing the envelope mass, then prevents the flash mixing.

We investigate light-curve and spectral properties of ultra-stripped core-collapse supernovae. Ultra-stripped supernovae are the explosions of heavily stripped massive stars which lost their envelopes via binary interactions with a compact companion star. They eject only ~ 0.1 Msun and may be the main way to form double neutron-star systems which eventually merge emitting strong gravitational waves. We follow the evolution of an ultra-stripped supernova progenitor until iron core collapse and perform explosive nucleosynthesis calculations. We then synthesize light curves and spectra of ultra-stripped supernovae using the nucleosynthesis results and present their expected properties. Ultra-stripped supernovae synthesize ~ 0.01 Msun of radioactive 56Ni, and their typical peak luminosity is around 1e42 erg/s or -16 mag. Their typical rise time is 5 - 10 days. Comparing synthesized and observed spectra, we find that SN 2005ek, some of the so-called calcium-rich gap transients, and SN 2010X may be related to ultra-stripped supernovae. If these supernovae are actually ultra-stripped supernovae, their event rate is expected to be about 1 per cent of core-collapse supernovae. Comparing the double neutron-star merger rate obtained by future gravitational-wave observations and the ultra-stripped supernova rate obtained by optical transient surveys identified with our synthesized light-curve and spectral models, we will be able to judge whether ultra-stripped supernovae are actually a major contributor to the binary neutron star population and provide constraints on binary stellar evolution.

The recently detected gravitational wave signals (GW150914 and GW151226) of the merger event of a pair of relatively massive stellar-mass black holes (BHs) calls for an investigation of the formation of such progenitor systems in general. We analyse the common envelope (CE) stage of the "traditional" formation channel in binaries where the first-formed compact object undergoes an in-spiral inside the envelope of its evolved companion star and ejects the envelope in that process. We calculate envelope binding energies of donor stars with initial masses between 4 and 115 Msun for metallicities of Z=Zsun/2 and Z=Zsun/50, and derive minimum masses of in-spiralling objects needed to eject these envelopes. We find that CE evolution, besides from producing WD-WD and NS-NS binaries, may, in principle, also produce massive BH-BH systems with individual BH component masses up to ~50-60 Msun, in particular for donor stars evolved to giants. However, the physics of envelope ejection of massive stars remains uncertain. We discuss the applicability of the energy-budget formalism, the location of the bifurcation point, the recombination energy and the accretion energy during in-spiral as possible energy sources, and also comment on the effect of inflated helium cores. Massive stars in a wide range of metallicities and with initial masses up to at least 115 Msun may possibly shed their envelopes and survive CE evolution, depending on their initial orbital parameters, similarly to the situation for intermediate mass and low-mass stars with degenerate cores. We conclude that based on stellar structure calculations, and in the view of the usual simple energy budget analysis, events like GW150914 and GW151226 could possibly be produced from the CE channel. Calculations of post-CE orbital separations, however, and thus the estimated LIGO detection rates, remain highly uncertain. [Abridged]

It is generally believed that Ultracompact X-ray binaries (UCXBs) evolved from binaries consisting of a neutron star accreting from a low-mass white dwarf or helium star where mass transfer is driven by gravitational radiation. However, the standard white-dwarf evolutionary channel cannot produce the relatively long-period ($40 - 60$\,min) UCXBs with high time-averaged mass-transfer rate. In this work, we explore an alternative evolutionary route toward UCXBs where the companions evolve from intermediate-mass Ap/Bp stars with an anomalously strong magnetic field ($100 - 10000$\,G). Including the magnetic braking caused by the coupling between the magnetic field and an irradiation-driven wind induced by the X-ray flux from the accreting component, we show that intermediate-mass X-ray binaries (IMXBs) can evolve into UCXBs. Using the \emph{MESA} code, we have calculated evolutionary sequences for a large number of IMXBs. The simulated results indicate that, for a small wind-driving efficiency $f=10^{-5}$, the anomalous magnetic braking can drive IMXBs to an ultra-short period of 11 min. Comparing our simulated results with the observed parameters of fifteen identified UCXBs, the anomalous magnetic braking evolutionary channel can account for the formation of seven and eight sources with $f=10^{-3}$, and $10^{-5}$, respectively. In particular, a relatively large value of $f$ can fit three of the long-period, persistent sources with high mass-transfer rate. Though the proportion of Ap/Bp stars in intermediate-mass stars is only 5\%, the lifetime of the UCXB phase is $\gtrsim$ 2 Gyr, producing a relatively high number of observable systems, making this an alternative evolutionary channel for the formation of UCXBs.

Recent observations suggest that some type Ia supernovae (SNe Ia) originate from the merging of two carbon-oxygen white dwarfs (CO WDs). Meanwhile, recent hydrodynamical simulations have indicated that the accretion-induced collapse may be avoided under certain conditions when double WDs merge violently. However, the properties of SNe Ia from this violent merger scenario are highly dependent on a particular mass-accretion stage, the so-called WD + He subgiant channel, during which the primary WD is able to increase its mass by accreting He-rich material from a He subgiant before the systems evolves into a double WD system. In this article, we aim to study this particular evolutionary stage systematically and give the properties of violent WD mergers. By employing the Eggleton stellar evolution code, we followed a large number of binary calculations and obtained the regions in parameter space for producing violent mergers based on the WD + He subgiant channel. According to these simulations, we found that the primary WDs can increase their mass by ~0.10-0.45 Msun during the mass-accretion stage. We then conducted a series of binary population synthesis calculations and found that the Galactic SN Ia birthrate from this channel is about 0.01-0.4*10-3 yr-1. This suggests that the violent WD mergers from this channel may only contribute to ~0.3%-10% of all SNe Ia in our Galaxy. The delay times of violent WD mergers from this channel are >1.7 Gyr, contributing to the SNe Ia in old populations. We also found that the WD + He subgiant channel is the dominent way for producing violent WD mergers that may be able to eventually explode as SNe Ia.

Recent advances in gravitational-wave astronomy make the direct detection of gravitational waves from the merger of two stellar-mass compact objects a realistic prospect. Evolutionary scenarios towards mergers of double compact objects generally invoke common-envelope evolution which is poorly understood, leading to large uncertainties in merger rates. We explore the alternative scenario of massive overcontact binary (MOB) evolution, which involves two very massive stars in a very tight binary which remain fully mixed due to their tidally induced high spin. We use the public stellar-evolution code MESA to systematically study this channel by means of detailed simulations. We find that, at low metallicity, MOBs produce double-black-hole (BH+BH) systems that will merge within a Hubble time with mass ratios close to one, in two mass ranges, ~25...60msun and >~ 130msun, with pair instability supernovae (PISNe) being produced in-between. Our models are also able to reproduce counterparts of various stages in the MOB scenario in the local Universe, providing direct support for it. We map the initial parameter space that produces BH+BH mergers, determine the expected chirp mass distribution, merger times, Kerr parameters and predict event rates. We typically find that for Z~<Z_sun/10, there is one BH+BH merger for ~1000 core-collapse supernovae. The advanced LIGO (aLIGO) detection rate is more uncertain and depends on the metallicity evolution. Deriving upper and lower limits from a local and a global approximation for the metallicity distribution of massive stars, we estimate aLIGO detection rates (at design limit) of ~19-550 yr^(-1) for BH+BH mergers below the PISN gap and of ~2.1-370 yr^(-1) above the PISN gap. Even with conservative assumptions, we find that aLIGO should soon detect BH+BH mergers from the MOB scenario and that these could be the dominant source for aLIGO detections.

Approximately 10% of massive OBA main-sequence (MS) and pre-MS stars harbour strong, large-scale magnetic fields. At the same time there is a dearth of magnetic stars in close binaries. A process generating strong magnetic fields only in some stars must be responsible with the merging of pre-MS and MS stars being suggested as one such channel. Stars emerging from the coalescence of two MS stars are rejuvenated, appearing younger than they are. They can therefore be identified by comparison with reference clocks. Here we predict the rejuvenation of MS merger products over a wide range of masses and binary configurations calibrated to smoothed-particle-hydrodynamical merger models. We find that the rejuvenation is of the order of the nuclear timescale and is strongest in the lowest-mass mergers and the most evolved binary progenitors with the largest mass ratios. These predictions allow us to put constraints on the binary progenitors of merger products. We show that the magnetic stars HR 2949 and $\tau$ Sco are younger than the potential binary companion HR 2948 and the Upper-Sco association, respectively, making them promising merger candidates. We find that the age discrepancies and the potential binary progenitors of both are consistent with them being rejuvenated merger products, implying that their magnetic fields may originate from this channel. Searching for age discrepancies in magnetic stars is therefore a powerful way to explore which fraction of magnetic stars may have obtained their strong magnetic field in MS mergers and to improve our understanding of magnetism in massive stars and their remnants.

Recent interpretations of narrow, variable absorption lines detected in some Type Ia supernovae suggest that their progenitors are surrounded by dense, circumstellar material. Similar variations detected in the symbiotic recurrent nova system RS Oph, which undergoes thermonuclear outbursts every ~20 years, making it an ideal candidate to investigate the origin of these lines. To this end, we present simulations of multiple mass transfer-nova cycles in RS Oph. We find that the quiescent mass transfer produces a dense, equatorial outflow, i.e., concentrated towards the binary orbital plane, and an accretion disc forms around the white dwarf. The interaction of a spherical nova outburst with these aspherical circumstellar structures produces a bipolar outflow, similar to that seen in HST imaging of the 2006 outburst. In order to produce an ionization structure that is consistent with observations, a mass-loss rate of $5 \times 10^{-7}\,\mathrm{M}_{\odot}\,\mathrm{yr}^{-1}$ from the red giant is required. The simulations also produce a polar accretion flow, which may explain the broad wings of the quiescent H {\alpha} line and hard X-rays. By comparing simulated absorption line profiles to observations of the 2006 outburst, we are able to determine which components arise in the wind and which are due to the novae. We explore the possible behaviour of absorption line profiles as they may appear should a supernova occur in a system like RS Oph. Our models show similarities to supernovae like SN 2006X, but require a high mass-loss rate, $\dot{M} \sim 10^{-6}$ to $10^{-5}\,\mathrm{M_\odot}\,\mathrm{yr}^{-1}$, to explain the variability in SN 2006X.

The explosion of ultra-stripped stars in close binaries may explain new discoveries of weak and fast optical transients. We have demonstrated that helium star companions to neutron stars (NSs) may evolve into naked metal cores as low as ~1.5 Msun, barely above the Chandrasekhar mass limit, by the time they explode. Here we present a new systematic investigation of the progenitor evolution leading to such ultra-stripped supernovae (SNe), in some cases yielding pre-SN envelopes of less than 0.01 Msun. We discuss the nature of these SNe (electron-capture vs iron core-collapse) and their observational light-curve properties. Ultra-stripped SNe are highly relevant for binary pulsars, as well as gravitational wave detection of merging NSs by LIGO/VIRGO, since these events are expected to produce mainly low-kick NSs in the mass range 1.10-1.80 Msun.

The explosion of ultra-stripped stars in close binaries can lead to ejecta masses < 0.1 M_sun and may explain some of the recent discoveries of weak and fast optical transients. In Tauris et al. (2013), it was demonstrated that helium star companions to neutron stars (NSs) may experience mass transfer and evolve into naked ~1.5 M_sun metal cores, barely above the Chandrasekhar mass limit. Here we present a systematic investigation of the progenitor evolution leading to ultra-stripped supernovae (SNe). In particular, we examine the binary parameter space leading to electron-capture (EC SNe) and iron core-collapse SNe (Fe CCSNe), respectively, and determine the amount of helium ejected with applications to their observational classification as Type Ib or Type Ic. We mainly evolve systems where the SN progenitors are helium star donors of initial mass M_He = 2.5 - 3.5 M_sun in tight binaries with orbital periods of P_orb = 0.06 - 2.0 days, and hosting an accreting NS, but we also discuss the evolution of wider systems and of both more massive and lighter - as well as single - helium stars. In some cases we are able to follow the evolution until the onset of silicon burning, just a few days prior to the SN explosion. We find that ultra-stripped SNe are possible for both EC SNe and Fe CCSNe, and that the amount of helium ejected is correlated with P_orb - the tightest systems even having donors being stripped down to envelopes of less than 0.01 M_sun. We estimate the rise time of ultra-stripped SNe to be in the range 12 hr - 8 days, and light curve decay times between 1 and 50 days. Ultra-stripped SNe may produce NSs in the mass range 1.10 - 1.80 M_sun and are highly relevant for LIGO/VIRGO since most (possibly all) merging double NS systems have evolved through this phase. Finally, we discuss the low-velocity kicks which might be imparted on these resulting NSs at birth. [Abridged]

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.

Considering the uncertainties of the C-burning rate (CBR) and the treatment of convective boundaries, Chen et al. (2014) found that there is a regime where it is possible to form hybrid CONe white dwarfs (WDs), i.e. ONe WDs with carbon-rich cores. As these hybrid WDs can be as massive as 1.30 $M_{\odot}$, not much mass needs to be accreted for these objects to reach the Chandrasekhar limit and to explode as Type Ia supernovae (SNe Ia). We have investigated their contribution to the overall SN Ia birth rate and found that such SNe Ia tend to be relatively young with typical time delays between 0.1 and 1 Gyr, where some may be as young as 30 Myr. SNe Ia from hybrid CONe WDs may contribute several percent to all SNe Ia, depending on the common-envelope ejection efficiency and the CBR. We suggest that these SNe Ia may produce part of the 2002cx-like SN Ia class.

Recent discoveries of weak and fast optical transients raise the question of their origin. We investigate the minimum ejecta mass associated with core-collapse supernovae (SNe) of Type Ic. We show that mass transfer from a helium star to a compact companion can produce an ultra-stripped core which undergoes iron core collapse and leads to an extremely fast and faint SN Ic. In this Letter, a detailed example is presented in which the pre-SN stellar mass is barely above the Chandrasekhar limit, resulting in the ejection of only ~0.05-0.20 M_sun of material and the formation of a low-mass neutron star. We compute synthetic light curves of this case and demonstrate that SN 2005ek could be explained by our model. We estimate that the fraction of such ultra-stripped to all SNe could be as high as 0.001-0.01. Finally, we argue that the second explosion in some double neutron star systems (for example, the double pulsar PSR J0737-3039B) was likely associated with an ultra-stripped SN Ic.

Justham (2011) and DiStefano et al.\ (2011) proposed that the white-dwarf progenitor of a Type Ia supernova (SN Ia) may have to spin down before it can explode. As the white dwarf spin-down timescale is not well known theoretically, we here try to constrain it empirically (within the framework of this spin-down model) for progenitor systems that contain a giant donor and for which circumbinary material has been detected after the explosion: we obtain an upper limit of a few $10^{\rm 7} {\rm yr}$. Based on the study of Di Stefano & Kilic (2012) this means that it is too early to rule out the existence of a surviving companion in SNR 0509-67.5.

A key tracer of the elusive progenitor systems of Type Ia supernovae (SNe Ia) is the detection of narrow blueshifted time-varying Na I D absorption lines, interpreted as evidence of circumstellar material (CSM) surrounding the progenitor system. The origin of this material is controversial, but the simplest explanation is that it results from previous mass loss in a system containing a white dwarf and a non-degenerate companion star. We present new single-epoch intermediate-resolution spectra of 17 low-redshift SNe Ia taken with XShooter on the ESO Very Large Telescope. Combining this sample with events from the literature, we confirm an excess (~20 per cent) of SNe Ia displaying blueshifted narrow Na I D absorption features compared to non-blueshifted Na I D features. The host galaxies of SNe Ia displaying blueshifted absorption profiles are skewed towards later-type galaxies, compared to SNe Ia that show no Na I D absorption, and SNe Ia displaying blueshifted narrow Na I D absorption features have broader light curves. The strength of the Na I D absorption is stronger in SNe Ia displaying blueshifted Na I D absorption features than those without blueshifted features, and the strength of the blueshifted Na I D is correlated with the B-V colour of the SN at maximum light. This strongly suggests the absorbing material is local to the SN. In the context of the progenitor systems of SNe Ia, we discuss the significance of these findings and other recent observational evidence on the nature of SN Ia progenitors. We present a summary that suggests there are at least two distinct populations of normal, cosmologically useful SNe Ia.

We present a large parameter study where we investigate the structure of white dwarf (WD) merger remnants after the dynamical phase. A wide range of WD masses and compositions are explored and we also probe the effect of different initial conditions. We investigated the degree of mixing between the WDs, the conditions for detonations as well as the amount of gas ejected. We find that systems with lower mass ratios have more total angular momentum and as a result more mass is flung out in a tidal tail. Nuclear burning can affect the amount of mass ejected. Many WD binaries that contain a helium-rich WD achieve the conditions to trigger a detonation. In contrast, for carbon-oxygen transferring systems only the most massive mergers with a total mass above ~2.1 solar masses detonate. Even systems with lower mass may detonate long after the merger if the remnant remains above the Chandrasekhar mass and carbon is ignited at the centre. Finally, our findings are discussed in the context of several possible observed astrophysical events and stellar systems, such as hot subdwarfs, R Coronae Borealis stars, single massive white dwarfs, supernovae of type Ia and other transient events. A large database containing 225 white dwarf merger remnants is made available via a dedicated web page.

Long-orbital-period subdwarf B (sdB) stars with main-sequence companions are believed to be the product of stable Roche Lobe overflow (RLOF), a scenario challenged by recent observations. Here we represent the results of a systematic study of the orbital-period distribution of sdB binaries in this channel using detailed binary evolution calculations. We show that the observed orbital-period distribution of long-period sdB binaries can be well explained by this scenario. Furthermore, we find that, if the progenitors of the sdB stars have initial masses below the helium flash mass, the sdB binaries produced from stable RLOF follow a unique mass -- orbital period relation for a given metallicity $Z$; increasing the orbital period from $\sim 400$ to $\sim 1100$\,d corresponds to increasing the mass of the sdB star from $\sim 0.40$ to $\sim 0.49\,M_\odot$ for $Z=0.02$. We suggest that the longest sdB binaries (with orbital period $> 1100$\,d) could be the result of atmospheric RLOF. The mass -- orbital period relation can be tested observationally if the mass of the sdB star can be determined precisely, e.g.\ from asteroseismology. Using this relation, we revise the orbital period distribution of sdB binaries produced by the first stable RLOF channel for the best fitting model of Han et al (2003), and show that the orbital period has a peak around 830\,d.

We put constraints on the properties of the progenitors of peculiar calcium-rich transients using the distribution of locations within their host galaxies. We confirm that this class of transients do not follow the galaxy stellar mass profile and are more likely to be found in remote locations of their apparent hosts. We test the hypothesis that these transients are from low metallicity progenitors by comparing their spatial distributions with the predictions of self-consistent cosmological simulations that include star formation and chemical enrichment. We find that while metal-poor stars and our transient sample show a consistent preference for large offsets, metallicity alone cannot explain the extreme cases. Invoking a lower age limit on the progenitor helps to improve the match, indicating these events may result from a very old metal-poor population. We also investigate the radial distribution of globular cluster systems, and show that they too are consistent with the class of calcium-rich transients. Because photometric upper limits exist for globular clusters for some members of the class, a production mechanism related to the dense environment of globular clusters is not favoured for the calcium-rich events. However the methods developed in this paper may be used in the future to constrain the effects of low metallicity on radially distant core-collapse events or help establish a correlation with globular clusters for other classes of peculiar explosions.

The material that is ejected in a common-envelope (CE) phase in a close binary system provides an ideal environment for dust formation. By constructing a simple toy model to describe the evolution of the density and the temperature of CE ejecta and using the \emph{AGBDUST} code to model dust formation, we show that dust can form efficiently in this environment. The actual dust masses produced in the CE ejecta depend strongly on their temperature and density evolution. We estimate the total dust masses produced by CE evolution by means of a population synthesis code and show that, compared to dust production in AGB stars, the dust produced in CE ejecta may be quite significant and could even dominate under certain circumstances.

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.

In this paper, we report on our analysis using Hubble Space Telescope astrometry and Keck-I HIRES spectroscopy of the central six stars of Tycho's supernova remnant (SN 1572). With these data, we measured the proper motions, radial velocities, rotational velocities, and chemical abundances of these objects. Regarding the chemical abundances, we do not confirm the unusu- ally high [Ni/Fe] ratio previously reported for Tycho-G. Rather, we find that for all metrics in all stars, none exhibit the characteristics expected from traditional SN Ia single-degenerate-scenario calculations. The only possible exception is Tycho-B, a rare, metal-poor A-type star; however, we are unable to find a suitable scenario for it. Thus, we suggest that SN 1572 cannot be explained by the standard single-degenerate model.

We present an analysis of the maximum light, near ultraviolet (NUV; 2900-5500 A) spectra of 32 low redshift (0.001<z<0.08) Type Ia supernovae (SNe Ia), obtained with the Hubble Space Telescope (HST). We combine this spectroscopic sample with high-quality gri light curves obtained with robotic telescopes to measure photometric parameters, such as stretch, optical colour, and brightness. By comparing our data to a comparable sample of SNe Ia at intermediate-z (0.4<z<0.9), we detect modest spectral evolution (3-sigma), in the sense that our mean low-z NUV spectrum has a depressed flux compared to its intermediate-z counterpart. We also see a strongly increased dispersion about the mean with decreasing wavelength, confirming the results of earlier surveys. These trends are consistent with changes in metallicity as predicted by contemporary SN Ia spectral models. We also examine the properties of various NUV spectral diagnostics in the individual spectra. We find a general correlation between stretch and the velocity (or position) of many NUV spectral features. In particular, we observe that higher stretch SNe have larger Ca II H&K velocities, that also correlate with host galaxy stellar mass. This latter trend is probably driven by the well-established correlation between stretch and stellar mass. We find no trends between UV spectral features and optical colour. Mean spectra constructed according to whether the SN has a positive or negative Hubble residual show very little difference at NUV wavelengths, indicating that the NUV evolution and variation we identify do not directly correlate with Hubble residuals. Our work confirms and strengthens earlier conclusions regarding the complex behaviour of SNe Ia in the NUV spectral region, but suggests the correlations we find are more useful in constraining progenitor models than improving the use of SNe Ia as cosmological probes.

Type Ia supernovae play a significant role in the evolution of the Universe and have a wide range of applications. It is widely believed that these events are the thermonuclear explosions of carbon-oxygen white dwarfs close to the Chandrasekhar mass (1.38 M\odot). However, CO white dwarfs are born with masses much below the Chandrasekhar limit and thus require mass accretion to become Type Ia supernovae. There are two main scenarios for accretion. First, the merger of two white dwarfs and, second, a stable mass accretion from a companion star. According to predictions, this companion star (also referred to as donor star) survives the explosion and thus should be visible in the center of Type Ia remnants. In this paper we scrutinize the central stars (79 in total) of the SN 1006 remnant to search for the surviving donor star as predicted by this scenario. We find no star consistent with the traditional accretion scenario in SN1006.

There is a consensus that Type-Ia supernovae (SNe Ia) arise from the thermonuclear explosion of white dwarf stars that accrete matter from a binary companion. However, direct observation of SN Ia progenitors is lacking, and the precise nature of the binary companion remains uncertain. A temporal series of high-resolution optical spectra of the SN Ia PTF 11kx reveals a complex circumstellar environment that provides an unprecedentedly detailed view of the progenitor system. Multiple shells of circumsteller are detected and the SN ejecta are seen to interact with circumstellar material (CSM) starting 59 days after the explosion. These features are best described by a symbiotic nova progenitor, similar to RS Ophiuchi.

We present broad-band observations of the afterglow and environment of the short GRB 111020A. An extensive X-ray light curve from Swift/XRT, XMM-Newton and Chandra, spanning ~100 seconds to 10 days after the burst, reveals a significant break at t~2 days with pre- and post-break decline rates of alphaX,1 ~ -0.78 and alphaX,2<-1.7, respectively. Interpreted as a jet break, we infer a collimated outflow with an opening angle of ~3-8 degrees. The resulting beaming-corrected gamma-ray (10-1000 keV band) and blastwave kinetic energies are (2-3)e48 erg and (0.3-2)e49 erg, respectively, with the range depending on the unknown redshift of the burst. We report a radio afterglow limit of <39 microJy (3-sigma) from EVLA observations which, along with our finding that vc<vX, constrains the circumburst density to n~0.01-0.1 cm^(-3). Optical observations provide an afterglow limit of i>24.4 mag at 18 hours after the burst, and reveal a potential host galaxy with i~24.3 mag. The sub-arcsecond localization from Chandra provides a precise offset of 0.80"+/-0.11" (1-sigma) from this galaxy corresponding to an offset of 5-7 kpc for z=0.5-1.5. We find a high excess neutral Hydrogen column density of (7.5+/-2.0)e21 cm^(-2) (z=0). Our observations demonstrate that a growing fraction of short GRBs are collimated which may lead to a true event rate of >100-1000 Gpc^(-3) yr^(-1), in good agreement with the NS-NS merger rate of ~200-3000 Gpc^(-3) yr^(-1). This consistency is promising for coincident short GRB-gravitational wave searches in the forthcoming era of Advanced LIGO/VIRGO.

On August 24 (UT) the Palomar Transient Factory (PTF) discovered PTF11kly (SN 2011fe), the youngest and most nearby type Ia supernova (SN Ia) in decades. We followed this event up in the radio (centimeter and millimeter bands) and X-ray bands, starting about a day after the estimated explosion time. We present our analysis of the radio and X-ray observations, yielding the tightest constraints yet placed on the pre-explosion mass-loss rate from the progenitor system of this supernova. We find a robust limit of dM/dt<10^-8 (w/100 km/s) [M_solar/yr] from sensitive X-ray non-detections, as well as a similar limit from radio data, which depends, however, on assumptions about microphysical parameters. We discuss our results in the context of single-degenerate models for SNe Ia and find that our observations modestly disfavor symbiotic progenitor models involving a red giant donor, but cannot constrain systems accreting from main-sequence or sub-giant stars, including the popular supersoft channel. In view of the proximity of PTF11kly and the sensitivity of our prompt observations we would have to wait for a long time (decade or longer) in order to more meaningfully probe the circumstellar matter of Ia supernovae.

Two types of supernova are thought to produce the overwhelming majority of neutron stars in the Universe. The first type, iron-core collapse supernovae, occurs when a high-mass star develops a degenerate iron core that exceeds the Chandrasekhar limit. The second type, electron-capture supernovae, is associated with the collapse of a lower-mass oxygen-neon-magnesium core as it loses pressure support owing to the sudden capture of electrons by neon and/or magnesium nuclei. It has hitherto been impossible to identify the two distinct families of neutron stars produced in these formation channels. Here we report that a large, well-known class of neutron-star-hosting X-ray pulsars is actually composed of two distinct sub-populations with different characteristic spin periods, orbital periods and orbital eccentricities. This class, the Be/X-ray binaries, contains neutron stars that accrete material from a more massive companion star. The two sub-populations are most probably associated with the two distinct types of neutron-star-forming supernovae, with electron-capture supernovae preferentially producing system with short spin period, short orbital periods and low eccentricity. Intriguingly, the split between the two sub-populations is clearest in the distribution of the logarithm of spin period, a result that had not been predicted and which still remains to be explained.

Type Ia supernovae (SNe Ia) have been used empirically as standardized candles to reveal the accelerating universe even though fundamental details, such as the nature of the progenitor system and how the star explodes, remained a mystery. There is consensus that a white dwarf star explodes after accreting matter in a binary system, but the secondary could be anything from a main sequence star to a red giant, or even another white dwarf. The uncertainty stems from the fact that no recent SN Ia has been discovered close enough to detect the stars before explosion. Here we report early observations of SN 2011fe (PTF11kly) in M101 at a distance of 6.4 Mpc, the closest SN Ia in the past 25 years. We find that the exploding star was likely a carbon-oxygen white dwarf, and from the lack of an early shock we conclude that the companion was most likely a main sequence star. Early spectroscopy shows high-velocity oxygen that varies on a time scale of hours and extensive mixing of newly synthesized intermediate mass elements in the outermost layers of the supernova. A companion paper uses pre-explosion images to rule out luminous red giants and most helium stars as companions.

Type Ia supernovae (SNe) serve as a fundamental pillar of modern cosmology, owing to their large luminosity and a well-defined relationship between light-curve shape and peak brightness. The precision distance measurements enabled by SNe Ia first revealed the accelerating expansion of the universe, now widely believed (though hardly understood) to require the presence of a mysterious "dark" energy. General consensus holds that Type Ia SNe result from thermonuclear explosions of a white dwarf (WD) in a binary system; however, little is known of the precise nature of the companion star and the physical properties of the progenitor system. Here we make use of extensive historical imaging obtained at the location of SN 2011fe/PTF11kly, the closest SN Ia discovered in the digital imaging era, to constrain the visible-light luminosity of the progenitor to be 10-100 times fainter than previous limits on other SN Ia progenitors. This directly rules out luminous red giants and the vast majority of helium stars as the mass-donating companion to the exploding white dwarf. Any evolved red companion must have been born with mass less than 3.5 times the mass of the Sun. These observations favour a scenario where the exploding WD of SN 2011fe/PTF11kly, accreted matter either from another WD, or by Roche-lobe overflow from a subgiant or main-sequence companion star.

PTF10ops is a Type Ia supernova (SN Ia), whose lightcurve and spectral properties place it outside the current SN Ia subtype classifications. Its spectra display the characteristic lines of subluminous SNe Ia, but it has a normal-width lightcurve with a long rise-time, typical of normal luminosity SNe Ia. The early-time optical spectra of PTF10ops were modelled using a spectral fitting code and found to have all the lines typically seen in subluminous SNe Ia, without the need to invoke more uncommon elements. The host galaxy environment of PTF10ops is also unusual with no galaxy detected at the position of the SN down to an absolute limiting magnitude of r \geq -12.0 mag, but a very massive galaxy is present at a separation of ~148 kpc and at the same redshift as suggested by the SN spectral features. The progenitor of PTF10ops is most likely a very old star, possibly in a low metallicity environment, which affects its explosion mechanism and observational characteristics. PTF10ops does not easily fit into any of the current models of either subluminous or normal SN Ia progenitor channels.

We present spectroscopic observations of GRB091127 (z=0.490) at the peak of the putative associated supernova, SN2009nz. Subtracting a late-time spectrum of the host galaxy, we isolate the contribution of SN2009nz and uncover broad features typical of nearby GRB-SNe. This establishes unambiguously that GRB091127 was accompanied by a broad-lined Type Ic SN, and links a cosmological long burst with a standard energy release (E_gamma,iso ~ 1.1e52 erg) to a massive star progenitor. The spectrum of SN2009nz closely resembles that of SN2006aj, with SN2003dh also providing an acceptable match, but has significantly narrower features than SNe 1998bw and 2010bh, indicative of a lower expansion velocity. The photospheric velocity inferred from the SiII 6355 absorption feature, v_ph ~ 17,000 km/s, is indeed closer to that of SNe 2006aj and 2003dh than to the other GRB-SNe. Combining the measured velocity with the light curve peak brightness and width, we estimate the following explosion parameters: M_Ni ~ 0.35 M_sun, E_K ~ 2.3e51 erg, and M_ej ~ 1.4 M_sun, similar to those of SN2006aj. These properties indicate that SN2009nz follows a trend of lower M_Ni for GRB-SNe with lower E_K and M_ej. Equally important, since GRB091127 is a typical cosmological burst, the similarity of SN2009nz to SN2006aj either casts doubt on the claim that XRF060218/SN2006aj was powered by a neutron star, or indicates that the nature of the central engine is encoded in the SN properties but not in the prompt emission. Future spectra of GRB-SNe at z > 0.3, including proper subtraction of the host galaxy contribution, will shed light on the full dispersion of SN properties for standard long GRBs, on the relation between SNe associated with sub-energetic and standard GRBs, and on a potential dispersion in the associated SN types.

Gamma-ray bursts (GRBs) serve as powerful probes of the early Universe, with their luminous afterglows revealing the locations and physical properties of star forming galaxies at the highest redshifts, and potentially locating first generation (Population III) stars. Since GRB afterglows have intrinsically very simple spectra, they allow robust redshifts from low signal to noise spectroscopy, or photometry. Here we present a photometric redshift of z~9.4 for the Swift-detected GRB 090429B based on deep observations with Gemini-North, the Very Large Telescope, and the GRB Optical and Near-infrared Detector. Assuming a Small Magellanic Cloud dust law (which has been found in a majority of GRB sight-lines), the 90% likelihood range for the redshift is 9.06 < z < 9.52, although there is a low-probability tail to somewhat lower redshifts. Adopting Milky Way or Large Magellanic Cloud dust laws leads to very similar conclusions, while a Maiolino law does allow somewhat lower redshift solutions, but in all cases the most likely redshift is found to be z>7. The non-detection of the host galaxy to deep limits (Y_AB >~ 28 mag, which would correspond roughly to 0.001 L* at z=1) in our late time optical and infrared observations with the Hubble Space Telescope strongly supports the extreme redshift origin of GRB 090429B, since we would expect to have detected any low-z galaxy, even if it were highly dusty. Finally, the energetics of GRB 090429B are comparable to those of other GRBs, and suggest that its progenitor is not greatly different to those of lower redshift bursts.

Aims: Recurrent nova systems like RS Oph have been proposed as a possible channel to Type Ia Supernova explosions, based on the high mass of the accreting white dwarf. Additional support to this hypothesis has been recently provided by the detection of circumstellar material around SN2006X and SN2007le, showing a structure compatible with that expected for recurrent nova outbursts.In this paper we investigate the circumstellar environment of RS Oph and its structure, with the aim of establishing a firmer and independent link between this class of objects and Type Ia SN progenitors. Methods: We study the time evolution of CaII, NaI and KI absorption features in RS Oph, before, during, and after the last outburst, using multi-epoch, high-resolution spectroscopy, and applying the same method adopted for SN2006X and SN2007le. Results: A number of components, blue-shifted with respect to the systemic velocity of RS Oph, are detected. In particular, one feature strongly weakens in the first two weeks after the outburst, simultaneously with the disappearance of very narrow P-Cyg profiles overimposed on the much wider nova emission lines of H, He, FeII and other elements. Conclusions: We interpret these facts as the signature of density enhancements in the circumstellar material, suggesting that the recurrent eruptions might indeed create complex structures within the material lost by the donor star. This establishes a strong link between RS Oph and the progenitor system of the Type Ia SN2006X, for which similar features have been detected.

PTF09dav is a peculiar subluminous type Ia supernova (SN) discovered by the Palomar Transient Factory (PTF). Spectroscopically, it appears superficially similar to the class of subluminous SN1991bg-like SNe, but it has several unusual features which make it stand out from this population. Its peak luminosity is fainter than any previously discovered SN1991bg-like SN Ia (M_B -15.5), but without the unusually red optical colors expected if the faint luminosity were due to extinction. The photospheric optical spectra have very unusual strong lines of Sc II and Mg I, with possible Sr II, together with stronger than average Ti II and low velocities of ~6000 km/s. The host galaxy of PTF09dav is ambiguous. The SN lies either on the extreme outskirts (~41kpc) of a spiral galaxy, or in an very faint (M_R>-12.8) dwarf galaxy, unlike other 1991bg-like SNe which are invariably associated with massive, old stellar populations. PTF09dav is also an outlier on the light-curve-width--luminosity and color--luminosity relations derived for other sub-luminous SNe Ia. The inferred 56Ni mass is small (0.019+/-0.003Msun), as is the estimated ejecta mass of 0.36Msun. Taken together, these properties make PTF09dav a remarkable event. We discuss various physical models that could explain PTF09dav. Helium shell detonation or deflagration on the surface of a CO white-dwarf can explain some of the features of PTF09dav, including the presence of Sc and the low photospheric velocities, but the observed Si and Mg are not predicted to be very abundant in these models. We conclude that no single model is currently capable of explaining all of the observed signatures of PTF09dav.

We have computed an extensive grid of binary evolution tracks to represent low- and intermediate mass X-ray binaries (LMXBs and IMXBs). The grid includes 42,000 models which covers 60 initial donor masses over the range of 1-4 solar masses and, for each of these, 700 initial orbital periods over the range of 10-250 hours. These results can be applied to understanding LMXBs and IMXBs: those that evolve analogously to CVs; that form ultracompact binaries with orbital periods in the range of 6-50 minutes; and that lead to wide orbits with giant donors. We also investigate the relic binary recycled radio pulsars into which these systems evolve. To evolve the donor stars in this study, we utilized a newly developed stellar evolution code called "MESA" that was designed, among other things, to be able to handle very low-mass and degenerate donors. This first application of the results is aimed at an understanding of the newly discovered pulsar PSR J1614-2230 which has a 1.97 solar masses neutron star, orbital period = 8.7 days, and a companion star of 0.5 solar mass. We show that (i) this system is a cousin to the LMXB Cyg X-2; (ii) for neutron stars of canonical birth mass 1.4 solar masses, the initial donor stars which produce the closest relatives to PSR J1614-2230 have a mass between 3.4-3.8 solar masses; (iii) neutron stars as massive as 1.97 solar masses are not easy to produce in spite of the initially high mass of the donor star, unless they were already born as relatively massive neutron stars; (iv) to successfully produce a system like PSR J1614-2230 requires a minimum initial neutron star mass of at least 1.6+-0.1 solar masses, as well as initial donor masses and orbital period of ~ 4.25+-0.10 solar masses and ~49+-2 hrs, respectively; and (v) the current companion star is largely composed of CO, but should have a surface H abundance of ~10-15%.

(abridged) PSR J1903+0327, a millisecond pulsar in an eccentric (e = 0.44) 95-day orbit with a (~ 1Msun) companion poses a challenge to our understanding of stellar evolution in binary and multiple-star systems. Here we describe optical and radio observations which rule out most of the scenarios proposed to explain formation of this system. Radio timing measurements of three post-Keplerian effects yield the most precise measurement of the mass of a millisecond pulsar to date: 1.667 +/- 0.021 solar masses (99.7% confidence limit) (...). Optical spectroscopy of a proposed main sequence counterpart show that its orbital motion mirrors the pulsar's 95-day orbit; being therefore its binary companion (...) The optical detection also provides a measurement of the systemic radial velocity of the binary; this and the proper motion measured from pulsar timing allow the determination of the systemic 3-D velocity in the Galaxy. We find that the system is always within 270 pc of the plane of the Galaxy, but always more than 3 kpc away from the Galactic centre. Thus an exchange interaction in a dense stellar environment (like a globular cluster or the Galactic centre) is not likely to be the origin of this system. We suggest that after the supernova that formed it, the neutron star was in a tight orbit with a main-sequence star, the present companion was a tertiary farther out. The neutron star then accreted matter from its evolving inner MS companion, forming a millisecond pulsar. The former donor star then disappears, either due to a chaotic 3-body interaction with the outer star (caused by the expansion of the inner orbit that necessarily results from mass transfer), or in the case of a very compact inner system, due to ablation/accretion by the newly formed millisecond pulsar.

This paper presents the first results from a new citizen science project: Galaxy Zoo Supernovae. This proof of concept project uses members of the public to identify supernova candidates from the latest generation of wide-field imaging transient surveys. We describe the Galaxy Zoo Supernovae operations and scoring model, and demonstrate the effectiveness of this novel method using imaging data and transients from the Palomar Transient Factory (PTF). We examine the results collected over the period April-July 2010, during which nearly 14,000 supernova candidates from PTF were classified by more than 2,500 individuals within a few hours of data collection. We compare the transients selected by the citizen scientists to those identified by experienced PTF scanners, and find the agreement to be remarkable - Galaxy Zoo Supernovae performs comparably to the PTF scanners, and identified as transients 93% of the ~130 spectroscopically confirmed SNe that PTF located during the trial period (with no false positive identifications). Further analysis shows that only a small fraction of the lowest signal-to-noise SN detections (r > 19.5) are given low scores: Galaxy Zoo Supernovae correctly identifies all SNe with > 8{\sigma} detections in the PTF imaging data. The Galaxy Zoo Supernovae project has direct applicability to future transient searches such as the Large Synoptic Survey Telescope, by both rapidly identifying candidate transient events, and via the training and improvement of existing machine classifier algorithms.

We present 3D numerical simulations of the early evolution of long-duration gamma-ray bursts in the collapsar scenario. Starting from the core-collapse of a realistic progenitor model, we follow the formation and evolution of a central black hole and centrifugally balanced disk. The dense, hot accretion disk produces freely-escaping neutrinos and is hydrodynamically unstable to clumping and to forming non-axisymmetric (m=1, 2) modes. We show that these spiral structures, which form on dynamical timescales, can efficiently transfer angular momentum outward and can drive the high required accretion rates (>=0.1-1 M_sun) for producing a jet. We utilise the smoothed particle hydrodynamics code, Gadget-2, modified to implement relevant microphysics, such as cooling by neutrinos, a plausible treatment approximating the central object and relativistic effects. Finally, we discuss implications of this scenario as a source of energy to produce relativistically beamed gamma-ray jets.

We introduce two simplified nuclear networks that can be used in hydrostatic carbon burning reactions occurring in white dwarf interiors. They model the relevant nuclear reactions in carbon-oxygen white dwarfs (COWDs) approaching ignition in Type Ia supernova (SN Ia) progenitors, including the effects of the main e-captures and \beta-decays that drive the convective Urca process. They are based on studies of a detailed nuclear network compiled by the authors and are defined by approximate sets of differential equations whose derivations are included in the text. The first network, N1, provides a good first order estimation of the distribution of ashes and it also provides a simple picture of the main reactions occurring during this phase of evolution. The second network, N2, is a more refined version of N1 and can reproduce the evolution of the main physical properties of the full network to the 5% level. We compare the evolution of the mole fraction of the relevant nuclei, the neutron excess, the photon energy generation and the neutrino losses between both simplified networks and the detailed reaction network in a fixed temperature and density parcel of gas.

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.

The planet formation process and subsequent planet migration may lead to configurations resulting in strong dynamical interactions among the various planets. Well-studied possible outcomes include collisions between planets, scattering events that eject one or more of the planets, and a collision of one or more of the planets with the parent star. In this work we consider one other possibility that has seemingly been overlooked in the various scattering calculations presented in the literature: the tidal capture of two planets which leads to the formation of a binary planet (or binary brown dwarf) in orbit about the parent star. We carry out extensive numerical simulations of such dynamical and tidal interactions to explore the parameter space for the formation of such binary planets. We show that tidal formation of binary planets is possible for typical planet masses and distances from the host star. The detection (or lack thereof) of planet-planet binaries can thus be used to constrain the properties of planetary systems, including their mutual spacing during formation, and the fraction of close planets in very eccentric orbits which are believed to form by a closely related process.

We use a collection of 14 well-measured neutron star masses to strengthen the case that a substantial fraction of these neutron stars was formed via electron-capture supernovae (SNe) as opposed to Fe-core collapse SNe. The e-capture SNe are characterized by lower resultant gravitational masses and smaller natal kicks, leading to lower orbital eccentricities when the e-capture SN has led to the formation of the second neutron star in a binary system. Based on the measured masses and eccentricities, we identify four neutron stars, which have a mean post-collapse gravitational mass of ~1.25 solar masses, as the product of e-capture SNe. We associate the remaining ten neutron stars, which have a mean mass of 1.35 solar masses, with Fe-core collapse SNe. If the e-capture supernova occurs during the formation of the first neutron star, then this should substantially increase the formation probability for double neutron stars, given that more systems will remain bound with the smaller kicks. However, this does not appear to be the case for any of the observed systems, and we discuss possible reasons for this.

The origin of hot subdwarf B stars (sdBs) is still unclear. About half of the known sdBs are in close binary systems for which common envelope ejection is the most likely formation channel. Little is known about this dynamic phase of binary evolution. Due to the tidal influence of the companion in close binary systems, the rotation of the primary becomes synchronised to its orbital motion. In this case it is possible to constrain the mass of the companion, if the primary mass, its projected rotational velocity as well as its surface gravity are known. For the first time we measured the projected rotational velocities of a large sdB binary sample from high resolution spectra. We analysed a sample of 51 sdB stars in close binaries, 40 of which have known orbital parameters comprising half of all such systems known today. Synchronisation in sdB binaries is discussed both from the theoretical and the observational point of view. The masses and the nature of the unseen companions could be constrained in 31 cases. The companions to seven sdBs could be clearly identified as late M stars. One binary may have a brown dwarf companion. The unseen companions of nine sdBs are white dwarfs with typical masses. The mass of one white dwarf companion is very low. In eight cases the companion mass exceeds 0.9 solar masses, four of which even exceed the Chandrasekhar limit indicating that they may be neutron stars. Even stellar mass black holes are possible for the most massive companions. The distribution of the inclinations of the systems with low mass companions appears to be consistent with expectations, whereas a lack of high inclinations becomes obvious for the massive systems. We show that the formation of such systems can be explained with common envelope evolution and present an appropriate formation channel including two phases of unstable mass transfer and one supernova explosion.

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 the spectroscopic discovery of a broad-lined Type Ic supernova (SN 2010bh) associated with the nearby long-duration gamma-ray burst (GRB) 100316D. At z = 0.0593, this is the third-nearest GRB-SN. Nightly optical spectra obtained with the Magellan telescopes during the first week after explosion reveal the gradual emergence of very broad spectral features superposed on a blue continuum. The supernova features are typical of broad-lined SNe Ic and are generally consistent with previous supernovae associated with low-redshift GRBs. However, the inferred velocities of SN 2010bh at 21 days after explosion are a factor of ~2 times larger than those of the prototypical SN 1998bw at similar epochs, with v ~ 26,000 km/s, indicating a larger explosion energy or a different ejecta structure. A near-infrared spectrum taken 13.8 days after explosion shows no strong evidence for He I at 1.083 microns, implying that the progenitor was largely stripped of its helium envelope. The host galaxy is of low luminosity (M_R ~ -18.5 mag) and low metallicity (Z < 0.4 Z_solar), similar to the hosts of other low-redshift GRB-SNe.

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.

Most stars are members of binaries, and the evolution of a star in a close binary system differs from that of an ioslated star due to the proximity of its companion star. The components in a binary system interact in many ways and binary evolution leads to the formation of many peculiar stars, including blue stragglers and hot subdwarfs. We will discuss binary evolution and the formation of blue stragglers and hot subdwarfs, and show that those hot objects are important in the study of evolutionary population synthesis (EPS), and conclude that binary interactions should be included in the study of EPS. Indeed, binary interactions make a stellar population younger (hotter), and the far-ultraviolet (UV) excess in elliptical galaxies is shown to be most likely resulted from binary interactions. This has major implications for understanding the evolution of the far-UV excess and elliptical galaxies in general. In particular, it implies that the far-UV excess is not a sign of age, as had been postulated prviously and predicts that it should not be strongly dependent on the metallicity of the population, but exists universally from dwarf ellipticals to giant ellipticals.

It is thought that the first generations of massive stars in the Universe were an important, and quite possibly dominant, source of the ultra-violet radiation that reionized the hydrogen gas in the intergalactic medium (IGM); a state in which it has remained to the present day. Measurements of cosmic microwave background anisotropies suggest that this phase-change largely took place in the redshift range z=10.8 +/- 1.4, while observations of quasars and Lyman-alpha galaxies have shown that the process was essentially completed by z=6. However, the detailed history of reionization, and characteristics of the stars and proto-galaxies that drove it, remain unknown. Further progress in understanding requires direct observations of the sources of ultra-violet radiation in the era of reionization, and mapping the evolution of the neutral hydrogen fraction through time. The detection of galaxies at such redshifts is highly challenging, due to their intrinsic faintness and high luminosity distance, whilst bright quasars appear to be rare beyond z~7. Here we report the discovery of a gamma-ray burst, GRB 090423, at redshift z=8.26 -0.08 +0.07. This is well beyond the redshift of the most distant spectroscopically confirmed galaxy (z=6.96) and quasar (z=6.43). It establishes that massive stars were being produced, and dying as GRBs, ~625 million years after the Big Bang. In addition, the accurate position of the burst pinpoints the location of the most distant galaxy known to date. Larger samples of GRBs beyond z~7 will constrain the evolving rate of star formation in the early universe, while rapid spectroscopy of their afterglows will allow direct exploration of the progress of reionization with cosmic time.

It is widely believed that Type Ia supernovae (SN Ia) originate in binary systems where a white dwarf accretes material from a companion star until its mass approaches the Chandrasekhar mass and carbon is ignited in the white dwarf's core. This scenario predicts that the donor star should survive the supernova explosion, providing an opportunity to understand the progenitors of Type Ia supernovae.In this paper we argue that rotation is a generic signature expected of most non-giant donor stars that is easily measurable. \citep{2004Natur.431.1069R} examined stars in the center of the remnant of SN 1572 (Tycho's SN) and showed evidence that a subgiant star (Star G by their naming convention) near the remnant's centre was the system's donor star. We present high-resolution (R \simeq 40000) spectra taken with the High Dispersion Spectrograph on Subaru of this candidate donor star and measure the star's radial velocity as $79\pm 2$ \kms with respect to the LSR and put an upper limit on the star's rotation of 7.5 \kms. In addition, by comparing images that were taken in 1970 and 2004, we measure the proper motion of Star G to be $\mu_l = -1.6 \pm 2.1$ \masyr and $\mu_b = -2.7 \pm 1.6$ \masyr. We demonstrate that all of the measured properties of Star G presented in this paper are consistent with those of a star in the direction of Tycho's SN that is not associated with the supernova event. However, we discuss an unlikely, but still viable scenario for Star G to be the donor star, and suggest further observations that might be able to confirm or refute it.

We show how the recent discovery of a likely close white dwarf companion to the well known star Regulus, one of the brightest stars in the sky, leads to considerable insight into the prior evolutionary history of this star, including the cause of its current rapid rotation. We infer a relatively narrow range for the initial masses of the progenitor system: M_{10} = 2.3 +/- 0.2 M_sun and M_{20} = 1.7 +/- 0.2 M_sun, where M_{10} and M_{20} are the initial masses of the progenitors of the white dwarf and Regulus, respectively. In this scenario, the age of the Regulus system would exceed 1 Gyr. We also show that Regulus, with a current orbital period of 40 days, has an interesting future ahead of it. This includes (i) a common envelope phase, and, quite possibly, (ii) an sdB phase, followed by (iii) an AM CVn phase with orbital periods < 1 hr. Binary evolution calculations are presented in support of this scenario. We also discuss alternative possibilities, emphasizing the present uncertainties in binary evolution theory. Thus, this one particular star system illustrates many different aspects of binary stellar evolution.

Neutron stars and stellar-mass black holes are the remnants of massive stars, which ended their lives in supernova explosions. These exotic objects can only be studied in relatively rare cases. If they are interacting with close companions they become bright X-ray sources. If they are neutron stars, they may be detected as pulsars. Only a few hundred such systems are presently known in the Galaxy. However, there should be many more binaries with basically invisible compact objects in non-interacting binaries. Here we report the discovery of unseen compact companions to hot subdwarfs in close binary systems. Hot subdwarfs are evolved helium-core-burning stars that have lost most of their hydrogen envelopes, often due to binary interactions. Using high-resolution spectra and assuming tidal synchronisation of the subdwarfs, we were able to constrain the companion masses of 32 binaries. While most hot subdwarf binaries have white-dwarf or late-type main sequence companions, as predicted by binary evolution models, at least 5% of the observed subdwarfs must have very massive companions: unusually heavy white dwarfs, neutron stars and, in some cases, even black holes. We present evolutionary models which show that such binaries can indeed form if the system has evolved through two common-envelope phases. This new connection between hot subdwarfs, which are numerous in the Galaxy, and massive compact objects may lead to a tremendous increase in the number of known neutron stars and black holes and shed some light on this dark population and its evolutionary link to the X-ray binary population.

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.

The excess of far-ultraviolet (far-UV) radiation in elliptical galaxies has remained one of their most enduring puzzles. In contrast, the origin of old blue stars in the Milky Way, hot subdwarfs, is now reasonably well understood: they are hot stars that have lost their hydrogen envelopes by various binary interactions. Here, we review the main evolutionary channels that produce hot subdwarfs in the Galaxy and present the results of binary population synthesis simulations that reproduce the main properties of the Galactic hot-subdwarf population. Applying the same model to elliptical galaxies, we show how this model can explain the main observational properties of the far-UV excess, including the far-UV spectrum, without the need to invoke ad hoc physical processes. The model implies that the UV excess is not a sign of age, as has been postulated previously, and predicts that it should not be strongly dependent on the metallicity of the population.

I describe the IR and X-ray observational campaign we have undertaken for the purpose of determining the nature of the faint discrete X-ray source population discovered by Chandra in the Galactic Center (GC). Data obtained for this project includes a deep Chandra survey of the Galactic Bulge; deep, high resolution IR imaging from VLT/ISAAC, CTIO/ISPI, and the UKIDSS Galactic Plane Survey (GPS); and IR spectroscopy from VLT/ISAAC and IRTF/SpeX. By cross-correlating the GC X-ray imaging from Chandra with our IR surveys, we identify candidate counterparts to the X-ray sources via astrometry. Using a detailed IR extinction map, we are deriving magnitudes and colors for all the candidates. Having thus established a target list, we will use the multi-object IR spectrograph FLAMINGOS-2 on Gemini-South to carry out a spectroscopic survey of the candidate counterparts, to search for emission line signatures which are a hallmark of accreting binaries. By determining the nature of these X-ray sources, this FLAMINGOS-2 Galactic Center Survey will have a dramatic impact on our knowledge of the Galactic accreting binary population.

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

In this talk, we present the general principles of binary evolution and give two examples. The first example is the formation of subdwarf B stars (sdBs) and their application to the long-standing problem of ultraviolet excess (also known as UV-upturn) in elliptical galaxies. The second is for the progenitors of type Ia supernovae (SNe Ia). We discuss the main binary interactions, i.e., stable Roche lobe overflow (RLOF) and common envelope (CE) evolution, and show evolutionary channels leading to the formation of various binary-related objects. In the first example, we show that the binary model of sdB stars of Han et al. (2002, 2003) can reproduce field sdB stars and their counterparts, extreme horizontal branch (EHB) stars, in globular clusters. By applying the binary model to the study of evolutionary population synthesis, we have obtained an ``a priori'' model for the UV-upturn of elliptical galaxies and showed that the UV-upturn is most likely resulted from binary interactions. This has major implications for understanding the evolution of the UV excess and elliptical galaxies in general. In the second example, we introduce the single degenerate channel and the double degenerate channel for the progenitors of SNe Ia. We give the birth rates and delay time distributions for each channel and the distributions of companion stars at the moment of SN explosion for the single degenerate channel, which would help to search for the remnant companion stars observationally.

The collapsar model requires rapidly rotating Wolf-Rayet stars as progenitors of long gamma-ray bursts. However, Galactic Wolf-Rayet stars rapidly lose angular momentum due to their intense stellar winds. We investigate whether the tidal interaction of a Wolf-Rayet star with a compact object in a binary system can spin up the Wolf-Rayet star enough to produce a collapsar. We compute the evolution of close Wolf-Rayet binaries, including tidal angular momentum exchange, differential rotation of the Wolf-Rayet star, internal magnetic fields, stellar wind mass loss, and mass transfer. The Wolf-Rayet companion is approximated as a point mass. We then employ a population synthesis code to infer the occurrence rates of the various relevant binary evolution channels. We find that the simple scenario -- i.e., the Wolf-Rayet star being tidally spun up and producing a collapsar -- does not occur at solar metallicity and may only occur with low probability at low metallicity. It is limited by the widening of the binary orbit induced by the strong Wolf-Rayet wind or by the radius evolution of the Wolf-Rayet star that most often leads to a binary merger. The tidal effects enhance the merger rate of Wolf-Rayet stars with black holes such that it becomes comparable to the occurrence rate of long gamma-ray bursts.

We investigate the evolution of the properties of model populations of ultraluminous X-ray sources (ULXs) consisting of a black-hole accretor in a binary with a donor star. We have computed models corresponding to three different populations of black-hole binaries; two invoke stellar-mass (~10 Msun) black hole accretors, and the third utilizes intermediate-mass (~1000 Msun) black holes (IMBHs). For each of the three populations, we computed 30,000 binary evolution sequences using a full Henyey stellar evolution code. The optical flux from the model ULXs includes contributions from the accretion disk, due to x-ray irradiation as well as intrinsic viscous heating, and that due to the donor star. We present "probability images" for the ULX systems in planes of color-magnitude, orbital period vs. X-ray luminosity, and luminosity vs. evolution time. Estimates of the numbers of ULXs in a typical galaxy as functions of time and of X-ray luminosity are also presented. Our model CMDs are compared with six ULX counterparts that have been discussed in the literature. Overall, the observed systems seem more closely related to model systems with very high-mass donors (> ~25 Msun) in binaries with IMBH accretors. However, significant difficulties remain with both the IMBH and stellar-mass black hole models.

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.

The discovery of an excess of light in the far-ultraviolet (UV) spectrum in elliptical galaxies was a major surprise in 1969. While it is now clear that this UV excess is caused by an old population of hot helium-burning stars without large hydrogen-rich envelopes rather than young stars, their origin has remained a mystery. Here we show that these stars most likely lost their envelopes because of binary interactions, similar to the hot subdwarf population in our own Galaxy. This has major implications for understanding the evolution of the UV excess and of elliptical galaxies in general. In particular, it implies that the UV excess is not a sign of age, as had been postulated previously, and predicts that it should not be strongly dependent on the metallicity of the population.

The discovery of a flux excess in the far-ultraviolet (UV) spectrum of elliptical galaxies was a major surprise in 1969. While it is now clear that this UV excess is caused by an old population of hot helium-burning stars without large hydrogen-rich envelopes, rather than young stars, their origin has remained a mystery. Here we show that these stars most likely lost their envelopes because of binary interactions, similar to the hot subdwarf population in our own Galaxy. We have developed an evolutionary population synthesis model for the far-UV excess of elliptical galaxies based on the binary model developed by Han et al (2002, 2003) for the formation of hot subdwarfs in our Galaxy. Despite its simplicity, it successfully reproduces most of the properties of elliptical galaxies with a UV excess: the range of observed UV excesses, both in $(1550-V)$ and $(2000-V)$, and their evolution with redshift. We also present colour-colour diagrams for use as diagnostic tools in the study of elliptical galaxies. The model has major implications for understanding the evolution of the UV excess and of elliptical galaxies in general. In particular, it implies that the UV excess is not a sign of age, as had been postulated previously, and predicts that it should not be strongly dependent on the metallicity of the population, but exists universally from dwarf ellipticals to giant ellipticals.

We systematically explore the evolution of the merger of two carbon-oxygen (CO) white dwarfs. The dynamical evolution of a 0.9 Msun + 0.6 Msun CO white dwarf merger is followed by a three-dimensional SPH simulation. We use an elaborate prescription in which artificial viscosity is essentially absent, unless a shock is detected, and a much larger number of SPH particles than earlier calculations. Based on this simulation, we suggest that the central region of the merger remnant can, once it has reached quasi-static equilibrium, be approximated as a differentially rotating CO star, which consists of a slowly rotating cold core and a rapidly rotating hot envelope surrounded by a centrifugally supported disc. We construct a model of the CO remnant that mimics the results of the SPH simulation using a one-dimensional hydrodynamic stellar evolution code and then follow its secular evolution. The stellar evolution models indicate that the growth of the cold core is controlled by neutrino cooling at the interface between the core and the hot envelope, and that carbon ignition in the envelope can be avoided despite high effective accretion rates. This result suggests that the assumption of forced accretion of cold matter that was adopted in previous studies of the evolution of double CO white dwarf merger remnants may not be appropriate. Our results imply that at least some products of double CO white dwarfs merger may be considered good candidates for the progenitors of Type Ia supernovae. In this case, the characteristic time delay between the initial dynamical merger and the eventual explosion would be ~10^5 yr. (Abridged).

The redshift distribution of the short-duration GRBs is a crucial, but currently fragmentary, clue to the nature of their progenitors. Here we present optical observations of nine short GRBs obtained with Gemini, Magellan, and the Hubble Space Telescope. We detect the afterglows and host galaxies of two short bursts, and host galaxies for two additional bursts with known optical afterglow positions, and five with X-ray positions (<6'' radius). In eight of the nine cases we find that the most probable host galaxies are faint, R~23-26.5 mag, and are therefore starkly different from the first few short GRB hosts with R~17-22 mag and z<0.5. Indeed, we measure spectroscopic redshifts of z~0.4-1.1 for the four brightest hosts. A comparison to large field galaxy samples, as well as the hosts of long GRBs and previous short GRBs, indicates that the fainter hosts likely reside at z>1. Our most conservative limit is that at least half of the five hosts without a known redshift reside at z>0.7 (97% confidence level), suggesting that about 1/3-2/3 of all short GRBs originate at higher redshifts than previously determined. This has two important implications: (i) We constrain the acceptable age distributions to a wide lognormal (sigma>1) with tau~4-8 Gyr, or to a power law, P(tau)~tau^n, with -1<n<0; and (ii) the inferred isotropic energies, E_{gamma,iso}~10^50-10^52 erg, are significantly larger than ~10^48-10^49 erg for the low redshift short GRBs, indicating a large spread in energy release or jet opening angles. Finally, we re-iterate the importance of short GRBs as potential gravitational wave sources and find a conservative Advanced LIGO detection rate of ~2-6 yr^-1.

The discovery of an excess of light in the far-ultraviolet (UV) in 1969 in elliptical galaxies was a major surprise. While it is now clear that this UV excess (UV-upturn) is probably caused by an old population of helium-burning stars. Han et al (2002, 2003) proposed a binary model for the formation of hot subdwarfs (helium burning stars) and the model can reproduce the observations in our Galaxy. By applying the binary model to the study of evolutionary population synthesis, we have obtained an {\it a priori} model for the UV-upturn of elliptical galaxies. The model shows that the UV-upturn is most likely resulted from binary interactions and it is universal (not very much metallicity-dependant) in ellipticals. This has major implications for understanding the evolution of the UV-upturn and elliptical galaxies in general; contrary to previous postulates, it implies that the UV-upturn is not a sign of age, but could be a potentially powerful indicator for a recent minor burst of star-forming activity.

Supernova 1987A, the first naked-eye supernova observed since Kepler's supernova in 1604, defies a number of theoretical expectations. Its anomalies have long been attributed to a merger between two massive stars that occurred some 20,000 years before the explosion, but so far there has been no conclusive proof that this merger took place. Here, we present three-dimensional hydrodynamical simulations of the mass ejection associated with such a merger and the subsequent evolution of the ejecta, and we show that this accurately reproduces the properties of the triple-ring nebula surrounding the supernova.

Although there is strong support for the collapsar engine as the power source of long-duration gamma-ray bursts (GRBs), we still do not definitively know the progenitor of these explosions. Here we review the current set of progenitor scenarios for long-duration GRBs and the observational constraints on these scenarios. Examining these, we find that single-star models cannot be the only progenitor for long-duration GRBs. Several binary progenitors can match the solid observational constraints and also have the potential to match the trends we are currently seeing in the observations. Type Ib/c supernovae are also likely to be produced primarily in binaries; we discuss the relationship between the progenitors of these explosions and those of the long-duration GRBs.

We present results of a multi-wavelength program to study the faint discrete X-ray source population discovered by Chandra in the Galactic Centre (GC). From IR imaging obtained with the VLT we identify candidate K-band counterparts to 75% of the X-ray sources in our sample. By combining follow-up VLT K-band spectroscopy of a subset of these candidate counterparts with the magnitude limits of our photometric survey, we suggest that only a small percentage of the sources are HMXBs, while the majority are likely to be canonical LMXBs and CVs at the distance of the GC. In addition, we present our discovery of highly structured small-scale (5-15") extinction towards the Galactic Centre. This is the finest-scale extinction study of the Galactic Centre to date. Finally, from these VLT observations we are able to place constraints on the stellar counterpart to the ``bursting pulsar'' GRO J1744-28.

We investigate the difference between the host galaxy properties of core-collapse supernovae and long-duration gamma-ray bursts (LGRBs), and quantify a possible metallicity dependence of the efficiency of producing LGRBs. We use a sample of 16 CC SNe and 16 LGRBs from Fruchter et al. (2006) which have similar redshift distributions to eliminate galaxy evolution biases. We make a forward prediction of their host galaxy luminosity distributions from the overall cosmic metallicity distribution of star formation. This appoach is supported by the finding that LGRB hosts follow the L-Z relations of star-forming galaxies. We then compare predictions for metallicity-dependent event efficiencies with the observed host data. We find that UV-based SFR estimates predict the hosts distribution of CC SNe perfectly well in a metallicity-independent form. In contrast, LGRB hosts are fainter on average by one magnitude, almost as faint as the Large Magellanic Cloud. Assuming this is a metallicity effect, the present data are insufficient to discriminate between a sharp metallicity cutoff and a soft decrease in efficiency towards higher metallicity. For a sharp cut-off, however, we find a best value for the cutoff metallicity, as reflected in the oxygen abundance, 12+log (O/H)_lim ~ 8.7+/-0.3 at 95% confidence including systematic uncertainties on the calibration of Kobulnicky & Kewley (2004). This value is somewhat lower than the traditionally quoted value for the Sun, but is comparable to the revised solar oxygen abundance (Asplund, Grevesse & Sauval 2005). LGRB models that require sharp metallicity cutoffs well below ~1/2 the revised solar metallicity appear to be effectively ruled out (abridged).

Theoretical models predict that the initial metallicity of the progenitor of a Type Ia supernova (SN Ia) affects the peak of the supernova light curve. This can cause a deviation from the standard light curve calibration employed when using SNe Ia as standardizable distance candles and, if there is a systematic evolution of the metallicity of SN Ia progenitors, could affect the determination of cosmological parameters. Here we show that this metallicity effect can be substantially larger than has been estimated previously, when the neutronisation in the immediate pre-explosion phase in the CO white dwarf is taken into account, and quantitatively assess the importance of metallicity evolution for determining cosmological parameters. We show that, in principle, a moderate and plausible amount of metallicity evolution could mimic a lambda-dominated, flat Universe in an open, lambda-free Universe. However, the effect of metallicity evolution appears not large enough to explain the high-z SN Ia data in a flat Universe, for which there is strong independent evidence, without a cosmological constant. We also estimate the systematic uncertainties introduced by metallicity evolution in a lambda-dominated, flat Universe. We find that metallicity evolution may limit the precision with which Omega_m and w can be measured and that it will be difficult to distinguish evolution of the equation of state of dark energy from metallicity evolution, at least from SN Ia data alone.

A variation of the relative content of (54Fe+58Ni) versus 56Ni may be responsible for the observed scatter of Type Ia Supernovae (SNe Ia) about a mean relation between their intrinsic brightness and the shape of their light curve. Synthetic light curves are computed of parametrised Chandrasekhar-mass explosion models of constant kinetic energy, where the ejecta are divided into an inner NSE zone, composed of (54Fe+58Ni) inside and 56Ni outside, an outer zone with Intermediate Mass Elements and a CO zone. Both the size of the NSE zone and the fraction of (54Fe+58Ni) v. 56Ni are varied systematically. Models with the same original NSE content but different (54Fe+58Ni)/56Ni ratios reach different peak brightness but have similar light curve shapes. Synthetic spectra indicate that the V-band decline rate is not affected by the (54Fe+58Ni)/56Ni ratio. While the 56Ni mass and the total NSE mass are the dominant parameters determining the peak luminosity and the shape of the light curve, respectively, a variation in the (54Fe+58Ni)/56Ni ratio, which may depend on the metallicity of the progenitor (Timmes, Brown & Truran 2003) is likely to account for a significant part of the observed scatter of local SNe Ia about the mean brightness--decline rate relation.

We re-assess the question of a systematic time delay between the formation of the progenitor and its explosion in a type Ia supernova (SN Ia) using the Hubble Higher-z Supernova Search sample (Strolger et al. 2004). While the previous analysis indicated a significant time delay, with a most likely value of 3.4 Gyr, effectively ruling out all previously proposed progenitor models, our analysis shows that the time-delay estimate is dominated by systematic errors, in particular due to uncertainties in the star-formation history. We find that none of the popular progenitor models under consideration can be ruled out with any significant degree of confidence. The inferred time delay is mainly determined by the peak in the assumed star-formation history. We show that, even with a much larger Supernova sample, the time delay distribution cannot be reliably reconstructed without better constraints on the star-formation history.

We present the results of a systematic exploration of an alternative evolutionary scenario to form double neutron-star binaries, first proposed by Brown (1995), which does not involve a neutron star passing through a common envelope. In this scenario, the initial binary components have very similar masses, and both components have left the main sequence before they evolve into contact; preferably the primary has already developed a CO core. We have performed population synthesis simulations to study the formation of double neutron star binaries via this channel and to predict the orbital properties and system velocities of such systems. We obtain a merger rate for DNSs in this channel in the range of 0.1 - 12/Myr. These rates are still subject to substantial uncertainties such as the modelling of the contact phase.

Supernova 2002ic was an atypical Type Ia supernova (SN Ia) with evidence for substantial amounts of hydrogen associated with the system. Contrary to previous claims, we show that its unusual properties can be understood within the framework of one of the most favoured progenitor models, the so-called supersoft channel. This requires that the donor star was initially relatively massive ($\sim 3\Ms$) and that the system experienced a delayed dynamical instability, leading to a large amount of mass loss from the system in the last few $10^4 $yr before the explosion. This can produce the inferred hydrogen-rich circumstellar environment, most likely with a disc-like geometry. However, to apply these models requires a larger accretion efficiency onto the white dwarf than is assumed in present parameterizations. If this is confirmed, it would most likely increase estimates for the frequency of the single-degenerate channel. Based on population synthesis simulations we estimate that not more than 1 in 100 SNe Ia should belong to this subgroup of SNe Ia.

Thanks to a stellar evolution code able to compute through the C-flash we link the binary population synthesis of single degenerate progenitors of type Ia supernovae (SNe Ia) to their physical condition at the time of ignition. We show that there is a large range of possible ignition densities and we detail how their probability distribution depends on the accretion properties. The low density peak of this distribution qualitatively reminds of the clustering of the luminosities of Branch-normal SNe Ia. We tighten the possible range of initial physical conditions for explosion models: they form a one-parameter family, independent of the metallicity. We discuss how these results may be modified if we were to relax our hypothesis of a permanent Hachisu wind or if we were to include electron captures.

We investigate the mass loss from a rotationally distorted envelope following the early, rapid in-spiral of a companion star inside a common envelope. For initially wide, massive binaries (M_1+M_2=20M_{\odot}, P\sim 10 yr), the primary has a convective envelope at the onset of mass transfer and is able to store much of the available orbital angular momentum in its expanded envelope. Three-dimensional smoothed particle hydrodynamics calculations show that mass loss is enhanced at mid-latitudes due to shock reflection from a torus-shaped outer envelope. Mass ejection in the equatorial plane is completely suppressed if the shock wave is too weak to penetrate the outer envelope in the equatorial direction (typically when the energy deposited in the star is less than about one-third of the binding energy of the envelope). We present a parameter study to show how the geometry of the ejecta depends on the angular momentum and the energy deposited in the envelope during a merging event. Applications to the nearly axisymmetric, but very non-spherical nebulae around SN1987A and Sheridan 25 are discussed, as well as possible links to RY Scuti and the Small Magellanic Cloud object R4.

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.

We present near-IR imaging of a sample of the faint, hard X-ray sources discovered in the 2001 Chandra ACIS-I survey towards the Galactic Centre (GC) (Wang et al. 2002). These ~800 discrete sources represent an important and previously undetected population within the Galaxy. From our VLT observations of 77 X-ray sources, we identify candidate K-band counterparts to 75% of the Chandra sources in our sample. The near-IR magnitudes and colours of the majority of candidate counterparts are consistent with highly reddened stars, indicating that most of the Chandra sources are likely to be accreting binaries at or near the GC.

Double neutron stars (DNSs), binary systems consisting of a radio pulsar and a generally undetected second neutron star (NS), have proven to be excellent laboratories for testing the theory of general relativity. The seven systems discovered in our Galaxy exhibit a remarkably well-defined relation between the pulsar spin period and the orbital eccentricity. Here we show, using a simple model where the pulsar is spun up by mass transfer from a helium-star companion, that this relation can only be produced if the second neutron star received a kick that is substantially smaller (with a velocity dispersion of less than 50 km/s) than the standard kick received by a single radio pulsar. This demonstrates that the kick mechanism depends on the evolutionary history of the NS progenitor and that the orbital parameters of DNSs are completely determined by the evolution in the preceding helium star - neutron star phase. This has important implications for estimating the rates of NS-NS mergers, one of the major potential astrophysical sources for the direct detection of gravitational waves, and for short-period gamma-ray bursts.

The double pulsar J0737--3039 has become an important astrophysical laboratory for testing fundamental physics. Here we demonstrate that the low measured mass of Pulsar B can be used to constrain the equation of state of neutron star matter {\em under the assumption} that it formed in an electron-capture supernova. We show that the observed orbital parameters as well as the likely evolutionary history of the system support such a hypothesis and discuss future refinements that will improve the constraints this test may provide.

Binaries containing a stellar-mass black hole and a recycled radio pulsar have so far eluded detection. We present a focused investigation of the formation and evolution of these systems in the Galactic disk, highlighting the factors that limit their numbers and the reasons why they may be extremely rare. We surmise that the birthrate of black-hole/recycled-pulsar binaries in the Galactic disk is probably no higher than ~10^{-7}/yr, and may be much less. Simple arguments regarding common-envelope evolution suggest that these binaries should have orbital periods less than 10 hr and an average lifetime of <~10^8 yr before coalescence due to the emission of gravitational radiation. We expect that fewer than ~10 of these compact, relativistic binaries currently reside in the Galactic disk, less than 0.1--1% of the number of double neutron stars. The discovery of two or more black-hole/recycled-pulsar binaries using current radio telescopes would tightly constrain certain ideas regarding the evolution of massive stars, dynamical mass transfer, and black-hole formation.

Evidence of high-velocity features such as those seen in the near-maximum spectra of some Type Ia Supernovae (eg SN 2000cx) has been searched for in the available SNIa spectra observed earlier than one week before B maximum. Recent observational efforts have doubled the number of SNeIa with very early spectra. Remarkably, all SNeIa with early data (7 in our RTN sample and 10 from other programmes) show signs of such features, to a greater or lesser degree, in CaII IR, and some also in SiII 6255A line. High-velocity features may be interpreted as abundance or density enhancements. Abundance enhancements would imply an outer region dominated by Si and Ca. Density enhancements may result from the sweeping up of circumstellar material by the highest velocity SN ejecta. In this scenario, the high incidence of HVFs suggests that a thick disc and/or a high-density companion wind surrounds the exploding white dwarf, as may be the case in Single Degenerate systems. Large-scale angular fluctuations in the radial density and abundance distribution may also be responsible: this could originate in the explosion, and would suggest a deflagration as the more likely explosion mechanism. CSM-interaction and surface fluctuations may coexist, possibly leaving different signatures on the spectrum. In some SNe the HVFs are narrowly confined in velocity, suggesting the ejection of blobs of burned material.

We present results from the first near-IR imaging of the weak X-ray sources discovered in the Chandra/ACIS-I survey (Wang et al. 2002) towards the Galactic Centre (GC). These ~800 discrete sources, which contribute significantly to the GC X-ray emission, represent an important and previously unknown population within the Galaxy. From our VLT observations we will identify likely IR counterparts to a sample of the hardest sources, which are most likely X-ray binaries. With these data we can place constraints on the nature of the discrete weak X-ray source population of the GC.

One possible fate of an accreting white dwarf is explosion in a type Ia supernova. However, the route to the thermonuclear runaway has always been uncertain owing to the lack of a convective model consistent with the Urca process. We derive a formalism for convective motions involving two radial flows. This formalism provides a framework for convective models that guarantees self-consistency for chemistry and energy budget, allows time-dependence and describes the interaction of convective motions with the global contraction or expansion of the star. In the one-stream limit, we reproduce several already existing convective models and allow them to treat chemistry. We also suggest as a model easy to implement in a stellar evolution code. We apply this formalism to convective Urca cores in Chandrasekhar mass white dwarfs. We stress that in degenerate matter, nuclear reactions that change the number of electrons strongly influence the convective velocities. We point out the sensitivity of the energy budget on the mixing. We illustrate our model by computing {\it stationary} convective cores with Urca nuclei. We show that even a very small mass fraction of Urca nuclei ($10^{-8}$) strongly influences the convective velocities. Finally, we present preliminary computations of the late evolution of a close to Chandrasekhar mass C+O white dwarf including the convective Urca process.

We derive a new formalism for convective motions involving two radial flows. This formalism provides a framework for convective models that guarantees consistency for the chemistry and the energy budget in the flows, allows time-dependence and accounts for the interaction of the convective motions with the global contraction or expansion of the star. In the one-stream limit the formalism reproduces several existing convective models and allows them to treat the chemistry in the flows. We suggest a version of the formalism that can be implemented easily in a stellar evolution code. We then apply the formalism to convective Urca cores in Chandrasekhar mass white dwarfs and compare it to previous studies. We demonstrate that, in degenerate matter, nuclear reactions that change the number of electrons strongly influence the convective velocities and we show that the net energy budget is sensitive to the mixing. We illustrate our model by computing stationary convective cores with Urca nuclei. Even a very small mass fraction of Urca nuclei (as little as $10^{-8}$) strongly influences the convective velocities. We conclude that the proper modelling of the Urca process is essential for determining the ignition conditions for the thermonuclear runaway in Chandrasekhar-mass white dwarfs.

Ultraluminous X-ray sources (ULXs) with Lx > 10^{39} ergs/s have been discovered in great numbers in external galaxies with ROSAT, Chandra, and XMM. The central question regarding this important class of sources is whether they represent an extension in the luminosity function of binary X-ray sources containing neutron stars and stellar-mass black holes (BHs), or a new class of objects, e.g., systems containing intermediate-mass black holes (100-1000 Msun). We have carried out a theoretical study to test whether a large fraction of the ULXs, especially those in galaxies with recent star formation activity, can be explained with binary systems containing stellar-mass black holes. To this end, we have applied a unique set of binary evolution models for black-hole X-ray binaries, coupled to a binary population synthesis code, to model the ULXs observed in external galaxies. We find that for donor stars with initial masses >10 Msun the mass transfer driven by the normal nuclear evolution of the donor star is sufficient to potentially power most ULXs. This is the case during core hydrogen burning and, to an even more pronounced degree, while the donor star ascends the giant branch, though the latter phases lasts only ~5% of the main sequence phase. We show that with only a modest violation of the Eddington limit, e.g., a factor of ~10, both the numbers and properties of the majority of the ULXs can be reproduced. One of our conclusions is that if stellar-mass black-hole binaries account for a significant fraction of ULXs in star-forming galaxies, then the rate of formation of such systems is ~3 x 10^{-7} per year normalized to a core-collapse supernova rate of 0.01 per year.

A critical comparison of estimates for the rates of hypernovae (HNe) and gamma-ray bursts (GRBs) is presented. Within the substantial uncertainties, the estimates are shown to be quite comparable and give a Galactic rate of $10^{-6}$ -- $10^{-5} $yr$^{-1}$ for both events. These rates are several orders of magnitude lower than the rate of core-collapse supernovae, suggesting that the evolution leading to a HN/GRB requires special circumstances, very likely due to binary interactions. Various possible binary channels are discussed, and it is shown that these are generally compatible with the inferred rates.

We systematically examine how the presence in a binary affects the final core structure of a massive star and its consequences for the subsequent supernova explosion. Interactions with a companion star may change the final rate of rotation, the size of the helium core, the strength of carbon burning and the final iron core mass. Stars with initial masses larger than \sim 11\Ms that experiece core collapse will generally have smaller iron cores at the time of the explosion if they lost their envelopes due to a previous binary interaction. Stars below \sim 11\Ms, on the other hand, can end up with larger helium and metal cores if they have a close companion, since the second dredge-up phase which reduces the helium core mass dramatically in single stars does not occur once the hydrogen envelope is lost. We find that the initially more massive stars in binary systems with masses in the range 8 - 11\Ms are likely to undergo an electron-capture supernova, while single stars in the same mass range would end as ONeMg white dwarfs. We suggest that the core collapse in an electron-capture supernova (and possibly in the case of relatively small iron cores) leads to a prompt explosion rather than a delayed neutrino-driven explosion and that this naturally produces neutron stars with low-velocity kicks. This leads to a dichotomous distribution of neutron star kicks, as inferred previously, where neutron stars in relatively close binaries attain low kick velocities. We illustrate the consequences of such a dichotomous kick scenario using binary population synthesis simulations and discuss its implications. This scenario has also important consequences for the minimum initial mass of a massive star that becomes a neutron star. (Abbreviated.)

We report on recent progress in our theoretical understanding of X-ray binaries, which has largely been driven by new observations, and illustrate the interplay between theory and observations considering as examples intermediate-mass X-ray binaries, irradiation-driven evolution, ultraluminous X-ray sources and neutron stars with low-velocity kicks.

The massive star which underwent core-collapse to produce SN1993J was identified as a non-variable red supergiant star in images of the galaxy M81 taken before explosion. However the stellar source showed an excess in UV and B-band colours that suggested it had either a hot, massive companion star or was embedded in an unresolved young stellar association. The spectra of SN1993J underwent a remarkable transformation between a hydrogen-rich Type II supernova and a helium-rich (hydrogen-deficient) Type Ib. The spectral and photometric peculiarities were explained by models in which the 13-20 solar mass supergiant had lost almost its entire hydrogen envelope to a close binary companion. The binary scenario is currently the best fitting model for the production of such type IIb supernovae, however the hypothetical massive companion stars have so far eluded discovery. Here we report the results of new photometric and spectroscopic observations of SN1993J, 10 years after explosion. At the position of the fading SN we detect the unambiguous signature of a massive star, the binary companion to the progenitor. This is evidence that this type of SN originate in interacting binary systems.

We present early results from the first IR imaging of the weak X-ray sources discovered in a recent Chandra survey towards the Galactic Centre. From our VLT observations we will identify likely counterparts to a sample of the hardest sources in order to place constraints on the nature of this previously unknown population.

We have carried out a detailed study of one of the most popular evolutionary channels for the production of Type Ia supernova (SN Ia) progenitors, the semi-degenerate channel (CO+MS), where a carbon/oxygen (CO) white dwarf (WD) accretes matter from an unevolved or slightly evolved non-degenerate star until it reaches the Chandrasekhar mass limit. Employing Eggleton's stellar evolution code and adopting the prescription of Hachisu et al. (1999) for the accretion efficiency, we have carried out full binary evolution calculations for about 2300 close WD binary systems and mapped out the initial parameters in the orbital period -- secondary mass ($P$--$M_2$) plane (for a range of WD masses) which lead to a successful Type Ia supernova. We obtained accurate, analytical fitting formulae to describe this parameter range which can be used for binary population synthesis (BPS) studies. The contours in the $P$--$M_2$ plane differ from those obtained by Hachisu et al. (1999) for low-mass CO WDs, which are more common than massive CO WDs. We show that white dwarfs with a mass as low as $0.67M_\odot$ can accrete efficiently and reach the Chandrasekhar limit. We have implemented these results in a BPS study to obtain the birthrates for SNe Ia and the evolution of birthrates with time of SNe Ia for both a constant star formation rate and a single star burst. The birthrates are somewhat lower than (but comparable to) those inferred observationally.

One promising method to test progenitor models for Type Ia supernovae is to identify surviving companion stars in historical supernova remnants. A surviving companion will have been strongly affected by its interaction with the supernova ejecta. Here we systematically explore the evolution and appearance of a typical companion star that has been stripped and heated by the supernova interaction during the post-impact re-equilibration phase. We show that, depending on the amount of heating and the amount of mass stripped by the impact or the previous binary mass transfer, such a star may be significantly overluminous or underluminous 10^3 - 10^4 yr after the supernova relative to its pre-supernova luminosity (by up to two orders of magnitude) and discuss the implications of these results for the strategies to be employed in searches for such companions.

(abridged) We present the first study that combines binary population synthesis in the Galactic disk and detailed evolutionary calculations of low- and intermediate-mass X-ray binaries (L/IMXBs). We show that the formation probability of IMXBs with initial donor masses of 1.5--4 Msun is typically >~5 times higher than that of standard LMXBs, and suggest that the majority of the observed systems may have descended from IMXBs. Distributions at the current epoch of the orbital periods, donor masses, and mass accretion rates have been computed, as have orbital-period distributions of BMPs. Several significant discrepancies between the theoretical and observed distributions are discussed. The orbital-period distribution of observed BMPs strongly favors cases where the envelope of the neutron-star progenitor is more easily ejected during the common-envelope phase. However, this leads to a >~100-fold overproduction of the theoretical number of luminous X-ray sources relative to the total observed number of LMXBs. X-ray irradiation of the donor star may result in a dramatic reduction in the X-ray active lifetime of L/IMXBs, thus possibly resolving the overproduction problem, as well as the long-standing BMP/LMXB birthrate problem.

We present the results of a systematic study of the formation and evolution of binaries containing black holes and normal-star companions with a wide range of masses. We first reexamine the standard formation scenario for close black-hole binaries, where the spiral-in of the companion in the envelope of a massive star causes the ejection of the envelope. We estimate the formation rates for different companion masses and different assumptions about the common-envelope structure and other model parameters. We find that black-hole binaries with intermediate- and high-mass secondaries can form for a wide range of assumptions, while black-hole binaries with low-mass secondaries can only form with apparently unrealistic assumptions (in agreement with previous studies). We then present detailed binary evolution sequences for black-hole binaries with secondaries of 2 to 17 Msun and demonstrate that in these systems the black hole can accrete appreciably even if accretion is Eddington limited (up to 7 Msun for an initial black-hole mass of 10 Msun) and that the black holes can be spun up significantly in the process. We discuss the implications of these calculations for well-studied black-hole binaries (in particular GRS 1915+105), ultra-luminous X-ray sources and Cygnus X-1. Finally, we discuss how some of the assumptions in the standard model could be relaxed to allow the formation of low-mass, short-period black-hole binaries which appear to be very abundant in Nature. (Abstract abridged)

We have carried out a detailed binary populations synthesis (BPS) study of the formation of subdwarf B (sdB) stars and related objects (sdO, sdOB stars) using the latest version of the BPS code developed by Han et al.(1994, 1995a, 1995b, 1998, 2001). We systematically investigate the importance of the five main evolutionary channels in which the sdB stars form after one or two common-envelope (CE) phases, one or two phases of stable Roche-lobe overflow (RLOF) or as the result of the merger of two helium white dwarfs (WD) (see Han et al. 2002, Paper I). Our best BPS model can satisfactorily explain the main observational characteristics of sdB stars, in particular their distributions in the orbital period - minimum companion mass diagram and in the effective temperature - surface gravity diagram, their distributions of orbital period, log (g theta^4), and mass function, their binary fraction and the fraction of sdB binaries with WD companions, their birthrates and their space density. We obtain a Galactic formation rate, a total number in the Galaxy, the intrinsic binary fraction for sdB stars. We also predict a distribution of masses for sdB stars that is wider than is commonly assumed and that some sdB stars have companions of spectral type as early as B. The percentage of A type stars with sdB companions can in principle be used to constrain some of the important parameters in the binary evolution model. We conclude that (a) the first RLOF phase needs to be more stable than is commonly assumed; (b) mass transfer in the first stable RLOF phase is non-conservative, and the mass lost from the system takes away a specific angular momentum similar to that of the system; (c) common-envelope ejection is very efficient.

We study the complete merger of two massive stars inside a common envelope and the subsequent evolution of the merger product, a rapidly rotating massive supergiant. Three qualitatively different types of mergers have been identified and investigated in detail, and the post-merger evolution has been followed to the immediate presupernova stage. The ``quiet merger'' case does not lead to significant changes in composition, and the star remains a red supergiant. In the case of a ``moderate merger'', the star may become a blue supergiant and end its evolution as a blue supergiant, depending on the core to total mass ratio (as may be appropriate for the progenitor of SN 1987A). In the case of the most effective ``explosive merger'', the merger product stays a red giant. In last two cases, the He abundance in the envelope is increased drastically, but significant s-processing is mainly expected in the ``explosive merger'' case.

We present unique radio observations of SS433, using MERLIN, the VLBA, and the VLA, which allow us to, for the first time, properly image and derive a meaningful spectral index for the `ruff' of equatorial emission which surrounds SS433's jet. We interpret this smooth ruff as a wind-like outflow from the binary.

Subdwarf B (sdB) stars (and related sdO/sdOB stars) are believed to be helium core-burning objects with very thin hydrogen-rich envelopes. In recent years it has become increasingly clear from observational surveys that a large fraction of these objects are members of binary systems. To better understand their formation, we here present the results of a detailed investigation of the three main binary evolution channels that can lead to the formation of sdB stars: the common envelope (CE) ejection channel, the stable Roche lobe overflow (RLOF) channel and the double helium white dwarfs (WDs) merger channel. We obtained the conditions for the formation of sdB stars from each of these channels using detailed stellar and binary evolution calculations where we modelled the detailed evolution of sdB stars and carried out simplified binary population synthesis simulations. The observed period distribution of sdB stars in compact binaries strongly constrains the CE ejection parameters. We also present the distribution of sdB stars in the $T_{\rm eff}$ - $\log g$ diagram, the Hertzsprung-Russell diagram and the distribution of mass functions.

We describe the theory and implementation of a three-dimensional fluid dynamics code which we have developed for calculating the surface geometry and circulation currents in the secondaries of interacting binary systems. The main method is based on an Eulerian-Lagrangian scheme to solve the advective and force terms in Euler's equation. Surface normalised spherical polar coordinates are used to allow the accurate modelling of the surface of the star, as is necessary when free surfaces and irradiation effects are to be considered. The potential and its gradient are expressed as sums of Legendre polynomials, which allows a very efficient solution of Poisson's equation. The basic solution scheme, based on operator splitting, is outlined, and standard numerical tests are presented.

We have searched for evidence of emission at 15 microns with ISOCAM and at 90 microns with ISOPHOT from dust orbiting six nearby pulsars, both in binaries and in isolation, located at distances between about 100 to 1000 pc. No emission was detected at any of the pulsar positions, and for the nearest pulsar J0108-1431 the 3 sigma upper limits on the flux density is about 66 mJy at 15 microns and 22.5 mJy at 90 microns. Upper limits on the masses of circumpulsar dust are inferred at a given temperature using a simple modelling of the radiated flux; they are compared to upper limits of orbiting mass obtained with the dust heating model of Foster & Fisher (1996). These results suggest that it is unlikely that any of them have sufficiently massive, circumpulsar discs out of which planets may form in the future.

We explore the possibility that neutron stars accreting from the winds of main-sequence stellar companions account for a significant fraction of low-luminosity, hard X-ray sources (L_X <~ 10^35 ergs/s; 1-10 keV) in the Galaxy. This work was motivated by recent Chandra observations of the Galactic center by Wang et al. (2002). Our calculations indicate that many of the discrete X-ray sources detected in this survey may be wind-accreting neutron stars, and that many more may be discovered with deeper X-ray observations. We propose that an infrared observing campaign be undertaken to search for the stellar counterparts of these X-ray sources.

We study the nucleosynthesis and the induced mixing during the merging of massive stars inside a common envelope. The systems of interest are close binaries, initially consisting of a massive red supergiant and a main-sequence companion of a few solar masses. We apply parameterized results based on hydrodynamical simulations to model the stream-core interaction and the response of the star in a standard stellar-evolution code. Preliminary results are presented illustrating the possibility of unusual nucleosynthesis and post-merging dredge-up which can cause composition anomalies in the supergiant's envelope.

We have imaged the X-ray binary SS433 with unprecedented Fourier-plane coverage at 6cm using simultaneously the VLBA, MERLIN, and the VLA, and also at 20cm with the VLBA. At both wavelengths we have securely detected smooth, low-surface brightness emission having the appearance of a `ruff' or collar attached perpendicularly to the well-studied knotty jets in this system, extending over at least a few hundred AU. We interpret this smooth emission as a wind-like outflow from the binary, and discuss its implications for the present evolutionary stage of this system.

We present detailed simulations of the interaction of a stream emanating from a mass-losing secondary with the core of a massive supergiant in the slow merger of the two stars inside a common envelope. The dynamics of the stream can be divided into a ballistic phase, starting at the L_1 point, and a hydrodynamical phase where the stream interacts strongly with the core. Considering the merger of a 1 and 5Msun star with a 20Msun evolved supergiant, we present two-dimensional hydrodynamical simulations using the PROMETHEUS code to demonstrate how the penetration depth and post-impact conditions depend on the initial properties of stream material (e.g. entropy, angular momentum, stream width) and the properties of the core (e.g. density structure and rotation rate). Using these results, we present a fitting formula for the entropy generated in the stream--core interaction and a recipe for the determination of the penetration depth based on a modified Bernoulli integral.

We investigate an interesting new class of high-mass X-ray binaries (HMXBs) with long orbital periods (P_orb > 30 days) and low eccentricities (e <~ 0.2). The orbital parameters suggest that the neutron stars in these systems did not receive a large impulse, or ``kick,'' at the time of formation. We develop a self-consistent phenomenological picture wherein the neutron stars born in the observed wide HMXBs receive only a small kick (<~ 50 km/s), while neutron stars born in isolation, in the majority of low-mass X-ray binaries, or in many of the well-known HMXBs with P_orb <~ 30 days receive the conventional large kicks, with a mean speed of ~ 300 km/s. We propose that the magnitude of the natal kick to a neutron star born in a binary system depends on the rotation rate of the pre-collapse core. We further suggest that the rotation rate of the core is a strong, well-defined function of the evolutionary path of the progenitor star.

As has only recently been recognized, X-ray binaries with intermediate-mass secondaries are much more important than previously believed. To assess the relative importance of low- and intermediate-mass X-ray binaries (LMXBs and IMXBs), we have initiated a systematic study of these systems consisting of two parts: an exploration of the evolution of LMXBs and IMXBs for a wide range of initial masses and orbital periods using detailed binary stellar evolution calculations, and an integration of these results into a Monte-Carlo binary population synthesis code. Here we present some of the main results of our binary calculations and some preliminary results of the population synthesis study for a ``standard'' reference model. While the inclusion of IMXBs improves the agreement with the observed properties of ``LMXBs'', several significant discrepancies remain, which suggests that additional physical processes need to be included in the model.

We present a method for the calculation of the effects of external irradiation on the geometrical shape of the secondary in a close binary containing a compact star, the source of the radiation, and a normal companion star, where we include the possibility of shadowing by an accretion disc. The model is based on a simple modification of the standard Roche binary potential in which the radiation-pressure force is parameterised using the ratio of the radiation to the gravitational force. We have constructed numerical solutions of approximate 3-dimensional irradiated equipotential surfaces to demonstrate the main geometrical effects of external radiation pressure. For systems in which the inner Lagrangian point is irradiated directly and for sufficiently high irradiation fluxes, the critical condition for which the secondary fills its tidal lobe changes from an inner to an outer critical configuration, where the critical equipotential surface connects to one of the outer rather than the inner Lagrangian point. Such a situation may apply to evaporating binary pulsar systems (e.g. PSR 1957+20), stars orbiting supermassive black holes in AGN and some high-mass X-ray binaries (including Centaurus X-3). For systems containing an accretion disc, which shadows the inner Lagrangian point from the external irradiation, the presence of significant radiation pressure causes a non-axisymmetric deformation of the stellar surface. We have calculated modified Roche potentials to determine the main effects on the optical lightcurves and radial velocity curves for typical binary parameters. Compared to previous studies, the inclusion of irradiation-pressure effects results in changes in the derived system parameters (e.g. component masses, radial velocities) that may be as high as 30 per cent.

Israelian et al. (1999) showed that the stellar companion of the black-hole binary Nova Sco is polluted with material ejected in the supernova that accompanied the formation of the black-hole primary. Here we systematically investigate the implications of these observations for the black-hole formation process. Using a variety of supernova models, including both standard as well as hypernova models (for different helium-star masses, explosion energies, and explosion geometries) and a simple model for the evolution of the binary and the pollution of the secondary, we show that most of the observed abundance anomalies can be explained for a large range of model parameters (apart from the abundance of Ti). The best models are obtained for He star masses of 10 to 16 Msun, where spherical hypernova models are generally favoured over standard supernova ones. Aspherical hypernova models also produce acceptable fits, provided there is extensive lateral mixing. All models require substantial fallback and that the fallback material either reached the orbit of the secondary or was mixed efficiently with material that escaped. The black hole therefore formed in a two-step process, where the initial mass of the collapsed remnant was increased substantially by matter that fell back after the initial collapse. This may help to explain the high observed space velocity of Nova Sco either because of a neutrino-induced kick (if a neutron star was formed first) or by asymmetric mass ejection in an asymmetric supernova explosion.

We present the results of a systematic study of cataclysmic variables (CVs) and related systems, combining detailed binary-population synthesis (BPS) models with a grid of 120 binary evolution sequences calculated with a Henyey-type stellar evolution code. In these sequences, we used 3 masses for the white dwarf (0.6, 0.8, 1.0 Msun) and seven masses for the donor star in the range of 0.6-1.4 Msun. The shortest orbital periods were chosen to have initially unevolved secondaries, and the longest orbital period for each secondary mass was taken to be just longer than the bifurcation period (16 - 22 hr), beyond which systems evolve towards long orbital periods. These calculations show that systems which start with evolved secondaries near the end or just after their main-sequence phase become ultra-compact systems with periods as short as 7 min. These systems are excellent candidates for AM CVn stars. Using a standard BPS code, we show how the properties of CVs at the beginning of mass transfer depend on the efficiency for common-envelope (CE) ejection and the efficiency of magnetic braking. In our standard model, where CE ejection is efficient, some 10 per cent of all CVs have initially evolved secondaries (with a central hydrogen abundance X_c < 0.4) and ultimately become ultra-compact systems (implying a Galactic birthrate for AM CVn-like stars of 10^{-3} yr^{-1}). Almost all CVs with orbital periods longer than 5 hr are found to have initially evolved or relatively massive secondaries. We show that their distribution of effective temperatures is in good agreement with the distribution of spectral types obtained by Beuermann et al. (1998).

We present ellipsoidal light-curve fits to the quiescent B, V, R and I light curves of GRO J1655-40 (Nova Scorpii 1994). The fits are based on a simple model consisting of a Roche-lobe filling secondary and an accretion disc around the black-hole primary. Unlike previous studies, no assumptions are made about the interstellar extinction or the distance to the source; instead these are determined self-consistently from the observed light curves. In order to obtain tighter limits on the model parameters, we used the distance determination from the kinematics of the radio jet as an additional constraint. We obtain a value for the extinction that is lower than was assumed previously; this leads to lower masses for both the black hole and the secondary star of 5.4 +/- 0.3 Msun and 1.45 +/- 0.35 Msun, respectively. The errors in the determination of the model parameters are dominated by systematic errors, in particular due to uncertainties in the modeling of the disk structure and uncertainties in the atmosphere model for the chemically anomalous secondary in the system. A lower mass of the secondary naturally explains the transient nature of the system if it is either in a late case A or early case B mass-transfer phase.

We present the results of a systematic study of the evolution of low- and intermediate-mass X-ray binaries (LMXBs and IMXBs). Using a standard Henyey-type stellar-evolution code and a standard model for binary interactions, we have calculated 100 binary evolution sequences containing a neutron star and a normal-type companion star, where the initial mass of the secondary ranges from 0.6 to 7\Ms and the initial orbital period from $\sim 4 $hr to $\sim 100 $d. This grid of models samples the entire range of parameters one is likely to encounter for LMXBs and IMXBs. The sequences show an enormous variety of evolutionary histories and outcomes, where different mass-transfer mechanisms dominate in different phases. Very few sequences resemble the classical evolution of cataclysmic variables, where the evolution is driven by magnetic braking and gravitational radiation alone. Many systems experience a phase of mass transfer on a thermal timescale and may briefly become detached immediately after this phase (for the more massive secondaries). In agreement with previous results (\markcite{Tauris1}Tauris & Savonije 1999), we find that all sequences with (sub-)giant donors up to $\sim 2\Ms$ are stable against dynamical mass transfer. Sequences where the secondary has a radiative envelope are stable against dynamical mass transfer for initial masses up to $\sim 4\Ms$. (abridged)

(abridged) Observations of very high speeds among pulsars in the Galactic disk present a puzzle regarding neutron stars in globular clusters. The inferred characteristic speed of single pulsars in the Galaxy is $\sim 5-10$ times as large as the central escape speed from the most massive globular clusters. It is then reasonable to ask why any pulsars are seen in globular clusters, whereas, in fact, quite a large number have been detected and as many as $\sim 1000$ are thought to be present in some of the richest clusters. We use a Monte Carlo approach to generate a population of massive primordial binaries. If we utilize the convention assumptions regarding mass transfer and neutron star kicks, we find that < 5% of the neutron stars initially formed in a massive cluster can be retained. We suggest that this fraction is too low to account for what is observed, and we speculate on possible alternative solutions to the retention problem.

We study the merging of massive stars inside a common envelope for binary systems consisting of a red supergiant with a mass of 15-20 Msun and a main-sequence companion of 1-5 Msun. We are particularly interested in the stage when the secondary, having overfilled its Roche lobe inside the common envelope, starts to transfer mass to the core of the primary at a very high mass-transfer rate and the subsequent nucleo-synthesis in the core-impact region. Using a parametrized model for the structure of the envelope at this stage, we perform 2-dimensional hydrodynamical calculations with the Munich Prometheus code to calculate the dynamics of the stream emanating from the secondary and its impact on the core of the primary. We find that, for the lower end of the estimated mass-transfer rate, low-entropy, hydrogen-rich material can penetrate deep into the primary core where nucleosynthesis through the hot CNO cycle can take place and that the associated neutron exposure may be sufficiently high for significant s-processing. For mass-transfer rates at the high end of our estimated range and higher densities in the stream, the stream impact can lead to the dredge-up of helium, but the neutron production is too low for significant s-processing.

We report here on two complementary population synthesis studies which relate directly to the formation and evolution of neutron star binaries in globular clusters. In the first, we compute the probability of retaining neutron stars in globular clusters, and quantitatively confirm the idea that the retention fraction for neutron stars born in binary systems is greatly enhanced over those born in isolated stars. However, the retention fraction may well be insufficient to explain the current population of neutron star binaries. In the second study, we follow a large population of primordial binaries and neutron stars throughout the lifetime of a globular cluster whose properties may be similar to 47 Tuc. We directly compute all 3-body interactions among binary systems, neutron stars, and isolated field stars throughout the history of the cluster. The evolution of certain types of neutron star binaries is followed up to the current epoch. The numbers of close, recycled, binary radio pulsars are evaluated and compared with the results of radio observations.

The X-ray binary Cygnus X-2 (Cyg X-2) has recently been shown to contain a secondary that is much more luminous and hotter than is appropriate for a low-mass subgiant. We present detailed binary-evolution calculations which demonstrate that the present evolutionary state of Cyg X-2 can be understood if the secondary had an initial mass of around 3.5 M_sun and started to transfer mass near the end of its main-sequence phase (or, somewhat less likely, just after leaving the main sequence). Most of the mass of the secondary must have been ejected from the system during an earlier rapid mass-transfer phase. In the present phase, the secondary has a mass of around 0.5 M_sun with a non-degenerate helium core. It is burning hydrogen in a shell, and mass transfer is driven by the advancement of the burning shell. Cyg X-2 therefore is related to a previously little studied class of intermediate-mass X-ray binaries (IMXBs). We suggest that perhaps a significant fraction of X-ray binaries presently classified as low-mass X-ray binaries may be descendants of IMXBs and discuss some of the implications.

We present a reanalysis of the outburst radial velocity data for X-Ray Nova Scorpii 1994. Using a model based on X-ray heating of the secondary star we suggest a more realistic treatment of the radial velocity data. Solutions are obtained in the (K_2,q) plane which, when combined with the published value for the binary mass ratio and inclination, constrain the mass of the black hole to within the region 4.1<M_1<6.6 Msun (90 per cent confidence), which is significantly lower than the value obtained by Orosz & Bailyn (1997). This reduced lower bound for the black hole mass together with the high space velocity of the system is consistent with the idea that it was formed by the post-supernova collapse of a neutron star.

We systematically investigate the effects of high supernova kick velocities on the orbital parameters of post-supernova neutron-star binaries. Using Monte- Carlo simulations, we determine the post-supernova distributions of orbital parameters for progeneitors of HMXBs and LMXBs. With the recent distribution of pulsar birth velocities by Lyne & Lorimer (1994), only about 27% of massive systems remain bound after the supernova, of which about 26% immediately experience dynamical mass transfer and possibly merge to become Thorne-Zytkow objects. The correlations between various orbital parameters can be compared with observational samples to yield information about supernova kick velocities and pre-supernova orbital-period distributions. After the supernova, the spins of most stars in massive systems have large inclinations with respect to their orbital axes, and a significant fraction of systems (about 20%) contain stars with retrograde spins. In the case of low-mass companions, we find that about 19% of systems remain bound after the supernova, of which about 57% experience immediate dynamical mass transfer. There is a tight correlation between the eccentricity and the post-supernova orbital period in LMXB progeneitors. All LMXBs with post-supernova periods longer than a few days initially have large eccentricities (e.g. Cir X-1 is discussed in this context). We use the results of these calculations to simulate the sky distributions of HMXBs and LMXBs. The simulated distributions agree with observed samples. Most importantly, the distribution of Galactic LMXBs is consistent with an ordinary Galactic disk population which has been widened because of large supernova kicks and does not require a special population of progenitors.