BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Stephen Justham

Stephen Justham contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.SRastro-ph.HEastro-ph.GAastro-ph.IM

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

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

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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 &#34;age-dating&#34; 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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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&#39;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.

0 citations0 reviewsNo code linkedOpen access