Researcher profile

Floor S. Broekgaarden

Floor S. Broekgaarden contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
6works
0followers
5topics
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

6 published item(s)

preprint2026arXiv

The Stellar Winds Atlas I: Current uncertainties in mass-loss rates

Stellar winds are a major source of uncertainty in understanding the life and deaths of massive stars. Across studies in the field, prescriptions for stellar winds differ substantially in both their physical assumptions and implementation, making them a dominant contributor to model-to-model variation. In this work, we present a systematic analysis of the physical assumptions underlying commonly adopted wind prescriptions for optically thin and optically thick winds of hot stars, as well as the winds of cool supergiants. Our analysis reveals substantial discrepancies across all regimes: predicted mass-loss rates for optically thin winds differ by more than an order of magnitude, while rates for cool supergiants vary by several orders of magnitude, with even wider uncertainties arising in extrapolation regimes beyond the Humphreys-Davidson limit. These disparities introduce significant ambiguity into the predicted formation of Wolf-Rayet (WR) stars, a problem further compounded by the inconsistent application of transition criteria. A central issue is the "cool Wolf-Rayet problem", a temperature regime where the classical electron-scattering Eddington factor ($Γ_{\rm e}$) loses physical consistency. Because this factor is widely used to determine WR mass-loss rates, its failure forces models to rely on uncertain extrapolations and ad-hoc corrections. We conclude that the dominant stellar wind uncertainties arise from a mismatch between the physical assumptions in stellar wind models and the structure of the stars to which they are applied. Our framework clarifies the origins of current theoretical discrepancies and identifies the key physical bottlenecks that must be addressed to improve mass-loss modeling for massive stars.

preprint2025arXiv

Multi-messenger observations in the Einstein Telescope era: binary neutron star and black hole - neutron star mergers

The Einstein Telescope (ET), a proposed next-generation gravitational wave (GW) observatory, will expand the reach of GW astronomy of stellar-mass compact object binaries to unprecedented distances, enhancing opportunities for multi-messenger observations. Here we investigate multi-messenger emission properties of binary neutron star (NSNS) and black hole-neutron star (BHNS) mergers detectable by ET, providing projections to optimize observational strategies and maximize scientific insights from these sources. Using a synthetic population of compact binary mergers, we model each source's GW signal-to-noise ratio, sky localization uncertainty, kilonova (KN) light curves in optical and near-infrared bands, fluence of the relativistic jet gamma-ray burst (GRB) prompt emission and afterglow light curves across radio, optical, X-ray and very high energy wavelengths. We analyze multi-messenger detectability prospects for ET as a standalone observatory with two different configurations and within a network of next-generation GW detectors. ET will detect over $10^4$ NSNS mergers annually, enabling potential observation of tens to hundreds of electromagnetic (EM) counterparts. BHNS mergers have more limited multi-messenger prospects, but joint GW-EM rates will increase by an order of magnitude compared to current-generation instruments. We quantify uncertainties due to the NS equation of state (EoS) and mass distribution of NSNSs, as well as the NS EoS and BH spin for BHNSs. While a single ET will achieve an impressive GW detection rate, the fraction of well-localized events is orders of magnitude lower than in a network with additional detectors. This significantly limits efficient EM follow-up and science cases requiring well-characterized counterparts or early observations. The challenge is even greater for BHNS mergers due to their low EM rate.

preprint2022arXiv

Rates of Compact Object Coalescences

Gravitational-wave detections are enabling measurements of the rate of coalescences of binaries composed of two compact objects -- neutron stars and/or black holes. The coalescence rate of binaries containing neutron stars is further constrained by electromagnetic observations, including Galactic radio binary pulsars and short gamma-ray bursts. Meanwhile, increasingly sophisticated models of compact objects merging through a variety of evolutionary channels produce a range of theoretically predicted rates. Rapid improvements in instrument sensitivity, along with plans for new and improved surveys, make this an opportune time to summarise the existing observational and theoretical knowledge of compact-binary coalescence rates.

preprint2021arXiv

Evidence from Disrupted Halo Dwarfs that $r$-process Enrichment via Neutron Star Mergers is Delayed by $\gtrsim500$ Myrs

The astrophysical origins of $r$-process elements remain elusive. Neutron star mergers (NSMs) and special classes of core-collapse supernovae (rCCSNe) are leading candidates. Due to these channels' distinct characteristic timescales (rCCSNe: prompt, NSMs: delayed), measuring $r$-process enrichment in galaxies of similar mass, but differing star-formation durations might prove informative. Two recently discovered disrupted dwarfs in the Milky Way's stellar halo, Kraken and \textit{Gaia}-Sausage Enceladus (GSE), afford precisely this opportunity: both have $M_{\star}\approx10^{8}M_{\rm{\odot}}$, but differing star-formation durations of ${\approx}2$ Gyrs and ${\approx}3.6$ Gyrs. Here we present $R\approx50,000$ Magellan/MIKE spectroscopy for 31 stars from these systems, detecting the $r$-process element Eu in all stars. Stars from both systems have similar [Mg/H]$\approx-1$, but Kraken has a median [Eu/Mg]$\approx-0.1$ while GSE has an elevated [Eu/Mg]$\approx0.2$. With simple models we argue NSM enrichment must be delayed by $500-1000$ Myrs to produce this difference. rCCSNe must also contribute, especially at early epochs, otherwise stars formed during the delay period would be Eu-free. In this picture, rCCSNe account for $\approx50\%$ of the Eu in Kraken, $\approx25\%$ in GSE, and $\approx15\%$ in dwarfs with extended star-formation durations like Sagittarius. The inferred delay time for NSM enrichment is $10-100\times$ longer than merger delay times from stellar population synthesis -- this is not necessarily surprising because the enrichment delay includes time taken for NSM ejecta to be incorporated into subsequent generations of stars. For example, this may be due to natal kicks that result in $r$-enriched material deposited far from star-forming gas, which then takes $\approx10^{8}-10^{9}$ years to cool in these galaxies.

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.