Source author record

Stephen Justham

Stephen Justham appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

25works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

25 published item(s)

preprint2025arXiv

Evolution of binaries containing a hot subdwarf and a white dwarf to double white dwarfs, and double detonation supernovae with hypervelocity runaway stars

Compact binaries containing hot subdwarfs and white dwarfs have the potential to evolve into a variety of explosive transients. These systems could also explain hypervelocity runaway stars such as US 708. We use the detailed binary evolution code MESA to evolve hot subdwarf and white dwarf stars interacting in binaries. We explore their evolution towards double detonation supernovae, helium novae, or double white dwarfs. Our grid of 3120 models maps from initial conditions such as orbital period and masses of hot subdwarf and white dwarf to these outcomes. The minimum amount of helium required to ignite the helium shell that leads to a double detonation supernova in our grid is $\approx 0.05 \, \mathrm{M_{\odot}}$, likely too large to produce spectra similar to normal type Ia supernovae, but compatible with inferred helium shell masses from some observed peculiar type I supernovae. We also provide the helium shell masses for our double white dwarf systems, with a maximum He shell mass of $\approx 0.18\,\mathrm{M_{\odot}}$. In our double detonation systems, the orbital velocity of the surviving donor star ranges from $\approx 450 \, \mathrm{km\,s^{-1}}$ to $\approx 1000 \, \mathrm{km\,s^{-1}}$. Among the surviving donors, we also estimate the runaway velocities of proto-white dwarfs, which have higher runaway velocities than hot subdwarf stars of the same mass. Our grid will provide a first-order estimate of the potential outcomes for the observation of binaries containing hot subdwarfs and white dwarfs from future missions like Gaia, LSST, and LISA.

preprint2025arXiv

The role of triple evolution in the formation of LISA double white dwarfs

Galactic double white dwarfs will be prominent gravitational-wave sources for the Laser Interferometer Space Antenna (LISA). While previous studies have primarily focused on formation scenarios in which binaries form and evolve in isolation, we present the first detailed study of the role of triple stellar evolution in forming the population of LISA double white dwarfs. In this work, we present the first detailed study of the role of triple stellar evolution in forming the population of LISA double white dwarfs. We use the multiple stellar evolution code (MSE) to model the stellar evolution, binary interactions, and the dynamics of triple star systems then use a Milky Way-like galaxy from the TNG50 simulations to construct a representative sample of LISA double white dwarfs. In our simulations about $7\times10^6$ Galactic double white dwarfs in the LISA frequency bandwidth originate from triple systems, whereas $\sim4\times10^6$ form from isolated binary stars. The properties of double white dwarfs formed in triples closely resemble those formed from isolated binaries, but we also find a small number of systems $\sim\mathcal{O}(10)$ that reach extreme eccentricities $(>0.9)$, a feature unique to the dynamical formation channels. Our population produces $\approx 10^{4} $ individually resolved double white dwarfs (from triple and binary channels) and an unresolved stochastic foreground below the level of the LISA instrumental noise. About $57\,\%$ of double white dwarfs from triple systems retain a bound third star when entering the LISA frequency bandwidth. However, we expect the tertiary stars to be too distant to have a detectable imprint in the gravitational-wave signal of the inner binary.

preprint2022arXiv

The circumstellar material around the Type IIP SN 2021yja

The majority of Type II-plateau supernovae (SNe IIP) have light curves that are not compatible with the explosions of stars in a vacuum; instead, the light curves require the progenitors to be embedded in circumstellar matter (CSM). We report on the successful fitting of the well-observed SN IIP 2021yja as a core-collapse explosion of a massive star with an initial mass of ~15 Msun and a pre-explosion radius of 631 Rsun. To explain the early-time behaviour of the broad-band light curves, the presence of 0.55 Msun CSM within ~2x10^14 cm is needed. Like many other SNe IIP, SN 2021yja exhibits an early-time flux excess including ultraviolet wavelengths. This, together with the short rise time (<2 days) in the gri bands, indicates the presence of a compact component in the CSM, essentially adjacent to the progenitor. We discuss the origin of the pre-existing CSM, which is most likely a common property of highly convective red supergiant envelopes. We argue that the difficulty in fitting the entire light curve with one spherical distribution indicates that the CSM around the SN 2021yja progenitor was asymmetric.

preprint2021arXiv

Impact of Massive Binary Star and Cosmic Evolution on Gravitational Wave Observations II: Double Compact Object Rates and Properties

Making the most of the rapidly increasing population of gravitational-wave detections of black hole (BH) and neutron star (NS) mergers requires comparing observations with population synthesis predictions. In this work we investigate the combined impact from the key uncertainties in population synthesis modelling of the isolated binary evolution channel: the physical processes in massive binary-star evolution and the star formation history as a function of metallicity, $Z$, and redshift $z, \mathcal{S}(Z,z)$. Considering these uncertainties we create 560 different publicly available model realizations and calculate the rate and distribution characteristics of detectable BHBH, BHNS, and NSNS mergers. We find that our stellar evolution and $\mathcal{S}(Z,z)$ variations can impact the predicted intrinsic and detectable merger rates by factors $10^2$-$10^4$. We find that BHBH rates are dominantly impacted by $\mathcal{S}(Z,z)$ variations, NSNS rates by stellar evolution variations and BHNS rates by both. We then consider the combined impact from all uncertainties considered in this work on the detectable mass distribution shapes (chirp mass, individual masses and mass ratio). We find that the BHNS mass distributions are predominantly impacted by massive binary-star evolution changes. For BHBH and NSNS we find that both uncertainties are important. We also find that the shape of the delay time and birth metallicity distributions are typically dominated by the choice of $\mathcal{S}(Z,z)$ for BHBH, BHNS and NSNS. We identify several examples of robust features in the mass distributions predicted by all 560 models, such that we expect more than 95% of BHBH detections to contain a BH $\gtrsim 8\,\rm{M}_{\odot}$ and have mass ratios $\lesssim 4$. Our work demonstrates that it is essential to consider a wide range of allowed models to study double compact object merger rates and properties.

preprint2021arXiv

Rapid stellar and binary population synthesis with COMPAS

Compact Object Mergers: Population Astrophysics and Statistics (COMPAS; https://compas.science) is a public rapid binary population synthesis code. COMPAS generates populations of isolated stellar binaries under a set of parametrized assumptions in order to allow comparisons against observational data sets, such as those coming from gravitational-wave observations of merging compact remnants. It includes a number of tools for population processing in addition to the core binary evolution components. COMPAS is publicly available via the github repository https://github.com/TeamCOMPAS/COMPAS/, and is designed to allow for flexible modifications as evolutionary models improve. This paper describes the methodology and implementation of COMPAS. It is a living document which will be updated as new features are added to COMPAS; the current document describes COMPAS v02.21.00.

preprint2020arXiv

Effect of binary evolution on the inferred initial and final core masses of hydrogen-rich, Type~II supernova progenitors

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

preprint2020arXiv

Looks can be deceiving: underestimating the age of single white dwarfs due to binary mergers

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

preprint2020arXiv

Luminous Red Novae: population models and future prospects

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

preprint2020arXiv

Non-interacting black hole binaries with Gaia and LAMOST

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

preprint2016arXiv

Time-Variable Aluminum Absorption in the Polar AR Ursae Majoris, and an Updated Estimate for the Mass of the White Dwarf

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

preprint2015arXiv

Another thread in the tapestry of stellar feedback: X-ray binaries

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

preprint2015arXiv

Globular cluster formation efficiencies from black-hole X-ray binary feedback

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

preprint2015arXiv

Relativistic baryonic jets from an ultraluminous supersoft X-ray source

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

preprint2014arXiv

Luminous Blue Variables and superluminous supernovae from binary mergers

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.

preprint2014arXiv

On the evolution of rotating accreting white dwarfs and type Ia supernovae

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

preprint2014arXiv

The binarity of Milky Way F,G,K stars as a function of effective temperature and metallicity

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

preprint2013arXiv

Double-detonation explosions as progenitors of type Iax supernovae

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

preprint2013arXiv

Producing type Iax supernovae from a specific class of helium-ignited WD explosions?

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

preprint2013arXiv

Puzzling accretion onto a black hole in the ultraluminous X-ray source M101 ULX-1

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

preprint2011arXiv

Single-degenerate type Ia supernovae without hydrogen contamination

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

preprint2010arXiv

Explosive Common-Envelope Ejection: Implications for Gamma-Ray Bursts and Low-Mass Black-Hole Binaries

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.

preprint2010arXiv

Observations of Doppler Boosting in Kepler Lightcurves

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.

preprint2010arXiv

On the formation of single and binary helium-rich sdO stars

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.

preprint2010arXiv

Sub-Chandrasekhar White Dwarf Mergers as the Progenitors of Type Ia Supernovae

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

preprint2009arXiv

Heavily Obscured AGN in Star-Forming Galaxies at z~2

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