Researcher profile

Simon Stevenson

Simon Stevenson contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

astro-ph.HE astro-ph.SR gr-qc astro-ph.IM astro-ph.GA

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author

13works

0followers

5topics

4close collaborators

Actions

Decide how to stay connected

Follow researcher0

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

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

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.

preprint2022arXiv

Biases in estimates of black hole kicks from the spin distribution of binary black holes

A population of more than 50 binary black hole mergers has now been observed by the LIGO and Virgo gravitational-wave observatories. While neutron stars are known to have large velocities associated with impulsive kicks imparted to them at birth in supernovae, whether black holes receive similar kicks, and of what magnitude, remains an open question. Recently, Callister et al. (2021) analysed the binary black hole population under the hypothesis that they were all formed through isolated binary evolution and claimed that large black hole kicks (greater than 260 km/s at 99% confidence) were required for the spin distribution of merging binary black holes to match observations. Here we highlight that a key assumption made by Callister et al. (2021) -- that all secondary black holes can be tidally spun up -- is not motivated by physical models, and may lead to a bias in their estimate of the magnitudes of black hole kicks. We make only minor changes to the Callister et al. (2021) model, accounting for a population of wider merging binaries where tidal synchronisation is ineffective. We show that this naturally produces a bimodal spin distribution for secondary black holes, and that the spin-orbit misalignments observed in the binary black hole population can be explained by more typical black hole kicks of order 100 km/s, consistent with kicks inferred from Galactic X-ray binaries containing black holes. We conclude that the majority of the binary black hole population is consistent with forming through isolated binary evolution.

preprint2022arXiv

Dynamical double black holes and their host cluster properties

We investigate the relationship between the global properties of star clusters and their double black hole (DBH) populations. We use the code {\tt NBODY6} to evolve a suite of star cluster models with an initial mass of $\mathcal{O}(10^4)$M$_\odot$ and varying initial parameters. We conclude that cluster metallicity plays the most significant role in determining the lifespan of a cluster, while the initial half-mass radius is dominant in setting the rate of BH exchange interactions in the central cluster regions. We find that the mass of interacting BHs, rather than how frequently their interactions with other BHs occur, is more crucial in the thermal expansion and eventual evaporation of the cluster. We formulate a novel approach to easily quantify the degree of BH-BH dynamical activity in each model. We report 12 in-cluster and three out-of-cluster (after ejection from the cluster) DBH mergers, of different types (inspiral, eccentric, hierarchical) across the ten $N$-body models presented. Our DBH merger efficiency is 3--4$\times10^{-5}$ mergers per M$_\odot$. We note the cluster initial density plays the most crucial role in determining the number of DBH mergers, with the potential presence of a transitional density point (between 1.2-3.8$\times10^3$M$_\odot$/pc$^3$) below which the number of in-cluster mergers increases with cluster density and above which the increased stellar density acts to prevent in-cluster BH mergers. The importance of the history of dynamical interactions within the cluster in setting up the pathways to ejected DBH mergers is also discussed. Our findings show a broad match with observed LIGO-Virgo DBH mergers.

preprint2022arXiv

Explaining the differences in massive star models from various simulations

The evolution of massive stars is the basis of several astrophysical investigations, from predicting gravitational-wave event rates to studying star-formation and stellar populations in clusters. However, uncertainties in massive star evolution present a significant challenge when accounting for these models' behaviour in stellar population studies. In this work, we present a comparison between five published sets of stellar models from the BPASS, BoOST, Geneva, MIST, and PARSEC simulations at near-solar metallicity. The different sets of stellar models have been computed using slightly different physical inputs in terms of mass-loss rates and internal mixing properties. Moreover, these models also employ various pragmatic methods to overcome the numerical difficulties that arise due to the presence of density inversions in the outer layers of stars more massive than 40 M$_\odot$. These density inversions result from the combination of inefficient convection in the low-density envelopes of massive stars and the excess of radiative luminosity to the Eddington luminosity. We find that the ionizing radiation released by the stellar populations can change by up to 18 percent, the maximum radial expansion of a star can differ between 100-1600 R$_\odot$, and the mass of the stellar remnant can vary up to 20 M$_\odot$ between the five sets of simulations. We conclude that any attempts to explain observations that rely on the use of models of stars more massive than 40 M$_\odot$ should be made with caution.

preprint2022arXiv

Modelling the formation of the first two neutron star-black hole mergers, GW200105 and GW200115: metallicity, chirp masses and merger remnant spins

The two neutron star-black hole mergers (GW200105 and GW200115) observed in gravitational waves by advanced LIGO and Virgo, mark the first ever discovery of such binaries in nature. We study these two neutron star-black hole systems through isolated binary evolution, using a grid of population synthesis models. Using both mass and spin observations (chirp mass, effective spin and remnant spin) of the binaries, we probe their different possible formation channels in different metallicity environments. Our models only support LIGO data when assuming the black hole is non spinning. Our results show a strong preference that GW200105 and GW200115 formed from stars with sub-solar metallicities $Z\lesssim 0.005$. Only two metal-rich ($Z=0.02$) models are in agreement with GW200115. We also find that chirp mass and remnant spins jointly aid in constraining the models, whilst the effective spin parameter does not add any further information. We also present the observable (i.e. post selection effects) median values of spin and mass distribution from all our models, which maybe used as a reference for future mergers. Further, we show that the remnant spin parameter distribution exhibits distinguishable features in different neutron star-black hole sub-populations. We find that non-spinning, first born black holes dominate significantly the merging neutron star-black hole population, ensuring electromagnetic counterparts to such mergers a rare affair.

preprint2022arXiv

Wide binary pulsars from electron-capture supernovae

Neutron stars receive velocity kicks at birth in supernovae. Those formed in electron-capture supernovae from super asymptotic giant branch stars -- the lowest mass stars to end their lives in supernovae -- may receive significantly lower kicks than typical neutron stars. Given that many massive stars are members of wide binaries, this suggests the existence of a population of low-mass ($1.25 < M_\mathrm{psr} / $M$_\odot < 1.3$), wide ($P_\mathrm{orb} \gtrsim 10^{4}$\,day), eccentric ($e \sim 0.7$), unrecycled ($P_\mathrm{spin} \sim 1$\,s) binary pulsars. The formation rate of such binaries is sensitive to the mass range of (effectively) single stars leading to electron capture supernovae, the amount of mass lost prior to the supernova, and the magnitude of any natal kick imparted on the neutron star. We estimate that one such binary pulsar should be observable in the Milky Way for every 10,000 isolated pulsars, assuming that the width of the mass range of single stars leading to electron-capture supernovae is $\lesssim 0.2$\,M$_\odot$, and that neutron stars formed in electron-capture supernovae receive typical kicks less than 10\,km s$^{-1}$. We have searched the catalog of observed binary pulsars, but find no convincing candidates that could be formed through this channel, consistent with this low predicted rate. Future observations with the Square Kilometre Array may detect this rare sub-class of binary pulsar and provide strong constraints on the properties of electron-capture supernovae and their progenitors.

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

Detecting Double Neutron Stars with LISA

We estimate the properties of the double neutron star (DNS) population that will be observable by the planned space-based interferometer LISA. By following the gravitational radiation driven evolution of DNSs generated from rapid population synthesis of massive binary stars, we estimate that around 35 DNSs will accumulate a signal-to-noise ratio above 8 over a four-year LISA mission. The observed population mainly comprises Galactic DNSs (94 per cent), but detections in the LMC (5 per cent) and SMC (1 per cent) may also be expected. The median orbital frequency of detected DNSs is expected to be 0.8 mHz, and many of them will be eccentric (median eccentricity of $0.11$). The orbital properties will provide insights into DNS progenitors and formation channels. LISA is expected to localise these DNSs to a typical angular resolution of $2^\circ$, with best-constrained sources localised to a few arcminutes. These localisations may allow neutron star natal kick magnitudes to be constrained through the Galactic distribution of DNSs, and make it possible to follow up the sources with radio pulsar searches. However, LISA is also expected to resolve $\sim 10^4$ Galactic double white dwarfs, many of which may have binary parameters that resemble DNSs; we discuss how the combined measurement of binary eccentricity, chirp mass, and sky location may aid the identification of a DNS. We expect the best-constrained DNSs to have eccentricities known to a few parts in a thousand, chirp masses measured to better than 1 per cent fractional uncertainty, and sky localisation at the level of a few arcminutes.

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

Modelling Double Neutron Stars: Radio and Gravitational Waves

We have implemented prescriptions for modelling pulsars in the rapid binary population synthesis code COMPAS. We perform a detailed analysis of the double neutron star (DNS) population, accounting for radio survey selection effects. The surface magnetic field decay timescale (${\sim}1000$\,Myr) and mass scale (${\sim}0.02$\,M$_\odot$) are the dominant uncertainties in our model. Mass accretion during common envelope evolution plays a non-trivial role in recycling pulsars. We find a best-fit model that is in broad agreement with the observed Galactic DNS population. Though the pulsar parameters (period and period derivative) are strongly biased by radio selection effects, the observed orbital parameters (orbital period and eccentricity) closely represent the intrinsic distributions. The number of radio observable DNSs in the Milky Way at present is about 2500 in our model, corresponding to approximately 10\% of the predicted total number of DNSs in the galaxy. Using our model calibrated to the Galactic DNS population, we make predictions for DNS mergers observed in gravitational waves. The DNS chirp mass distribution varies from 1.1M$_\odot$ to 2.1M$_\odot$ and the median is found to be 1.14\,M$_\mathrm{\odot}$. The expected effective spin $χ_\mathrm{eff}$ for isolated DNSs is $\lesssim$0.03 from our model. We predict that 34\% of the current Galactic isolated DNSs will merge within a Hubble time, and have a median total mass of 2.7\,M$_\mathrm{\odot}$. Finally, we discuss implications for fast radio bursts and post-merger remnant gravitational-waves.

preprint2020arXiv

On the origin of GW190425

The LIGO/Virgo collaborations recently announced the detection of a likely binary neutron star merger, GW190425. The total mass of GW190425 is significantly larger than the masses of Galactic double neutron stars known through radio astronomy. This suggests that GW190425 formed differently from Galactic double neutron stars. We hypothesize that GW190425 formed via unstable "case BB" mass transfer. According to this hypothesis, the progenitor of GW190425 was a binary consisting of a neutron star and a ${\sim} 4-5 {M_\odot}$ helium star, which underwent a common-envelope process. Following the supernovae of the helium star core, a tight, eccentric, double neutron star was formed, which merged in ${\lesssim}$ 10 Myr. The helium star progenitor may explain the unusually large mass of GW190425, while the short time to merger may explain why we do not see similar systems in radio. In order to test this hypothesis, we measure the eccentricity of GW190425 using publicly available LIGO/Virgo data. We constrain the eccentricity at 10 Hz to be $e \leq 0.007$ with $90\%$ confidence. This result provides no evidence for or against the unstable mass transfer scenario because the binary is likely to have circularized to $e\lesssim10^{-4}$ by the time it entered the LIGO/Virgo band. Future detectors operating in lower frequency bands will enable us to discern the formation channel of mergers similar to GW190425 using eccentricity measurements.

preprint2020arXiv

The fates of massive stars: exploring uncertainties in stellar evolution with METISSE

In the era of advanced electromagnetic and gravitational wave detectors, it has become increasingly important to effectively combine and study the impact of stellar evolution on binaries and dynamical systems of stars. Systematic studies dedicated to exploring uncertain parameters in stellar evolution are required to account for the recent observations of the stellar populations. We present a new approach to the commonly used Single-Star Evolution (SSE) fitting formulae, one that is more adaptable: Method of Interpolation for Single Star Evolution (METISSE). It makes use of interpolation between sets of pre-computed stellar tracks to approximate evolution parameters for a population of stars. We have used METISSE with detailed stellar tracks computed by the Modules for Experiments in Stellar Astrophysics (MESA), Bonn Evolutionary Code (BEC) and Cambridge STARS code. METISSE better reproduces stellar tracks computed using the STARS code compared to SSE, and is on average three times faster. Using stellar tracks computed with MESA and BEC, we apply METISSE to explore the differences in the remnant masses, the maximum radial expansion, and the main-sequence lifetime of massive stars. We find that different physical ingredients used in the evolution of stars, such as the treatment of radiation dominated envelopes, can impact their evolutionary outcome. For stars in the mass range 9 to 100 M$_\odot$, the predictions of remnant masses can vary by up to 20 M$_\odot$, while the maximum radial expansion achieved by a star can differ by an order of magnitude between different stellar models.

preprint2020arXiv

The origin of spin in binary black holes: Predicting the distributions of the main observables of Advanced LIGO

We study the formation of coalescing binary black holes via the evolution of isolated field binaries that go through the common envelope phase in order to obtain the combined distributions of the main observables of Advanced LIGO. We used a hybrid technique that combines the parametric binary population synthesis code COMPAS with detailed binary evolution simulations performed with the MESA code. We then convolved our binary evolution calculations with the redshift- and metallicity-dependent star-formation rate and the selection effects of gravitational-wave detectors to obtain predictions of observable properties. By assuming efficient angular momentum transport, we are able to present a model that is capable of simultaneously predicting the three main gravitational-wave observables: the effective inspiral spin parameter $χ_{eff}$, the chirp mass $M_{chirp}$, and the cosmological redshift of merger $z_{merger}$. We find an excellent agreement between our model and the ten events from the first two advanced detector observing runs. We make predictions for the third observing run O3 and for Advanced LIGO design sensitivity. We expect around 80% of events with $χ_{eff} < 0.1$, while the remaining 20% of events with $χ_{eff} \ge 0.1$ are split into ~10% with $M_{chirp} < 15$ M$_\odot$ and ~10% with $M_{chirp} \ge 15$ M$_\odot$. In conclusion, the favorable comparison of the existing LIGO/Virgo observations with our model predictions gives support to the idea that the majority, if not all of the observed mergers, originate from the evolution of isolated binaries. The first-born black hole has negligible spin because it lost its envelope after it expanded to become a giant star, while the spin of the second-born black hole is determined by the tidal spin up of its naked helium star progenitor by the first-born black hole companion after the binary finished the common-envelope phase.

Simon Stevenson

Quick read

Decide how to stay connected

How to connect with this researcher

Open a focused conversation when the fit is right

See the researcher in context

Building this graph slice

13 published item(s)

Biases in estimates of black hole kicks from the spin distribution of binary black holes

Dynamical double black holes and their host cluster properties

Explaining the differences in massive star models from various simulations

Modelling the formation of the first two neutron star-black hole mergers, GW200105 and GW200115: metallicity, chirp masses and merger remnant spins

Wide binary pulsars from electron-capture supernovae

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

Rapid stellar and binary population synthesis with COMPAS

Detecting Double Neutron Stars with LISA

Luminous Red Novae: population models and future prospects

Modelling Double Neutron Stars: Radio and Gravitational Waves

On the origin of GW190425

The fates of massive stars: exploring uncertainties in stellar evolution with METISSE

The origin of spin in binary black holes: Predicting the distributions of the main observables of Advanced LIGO