Researcher profile

Silvia Toonen

Silvia Toonen contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

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

18 published item(s)

preprint2026arXiv

Gaia Barium Dwarfs and Their Ostensibly Ordinary Counterparts

The recently identified Gaia population of main-sequence--white dwarf (MS+WD) binaries at separations of ${\sim}\,1~{\rm AU}$, often with moderate eccentricities, is not readily reproduced by binary population synthesis models. Barium stars represent a closely related population whose enrichment in $s$-process elements both confirms the presence of a WD companion and attests to past binary interaction. It also indicates that mass transfer occurred at least during the late- and post-AGB phases of the WD progenitor, when $s$-process elements are dredged up. In this work, we further explore the connection between the astrometrically identified Gaia MS+WD binaries and the classical barium star population. To achieve this, we used high-resolution FEROS spectroscopy to measure abundances for 30 Gaia DR3 non-single-star binaries, identifying 10 as Ba-enriched. Together with our recent analysis of archival GALAH data, this yields a sample of 38 barium dwarfs with dynamically measured WD masses, compared to only 6 previously known systems with known WD masses at these separations. We find that, in cases where metallicity is sufficiently low to facilitate efficient $s$-process production, barium and yttrium enrichment is often detected. This enrichment is also identified in eccentric systems, suggesting that post-AGB mass transfer mechanisms are capable of pumping eccentricity into the orbit or occur without erasing it. Our results indicate that the Gaia MS+WD binaries trace the population from which barium stars emerge. Treating the large Gaia-discovered population as an extension of known $s$-process enriched dwarfs opens an avenue to empirically constrain their formation and evolution.

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

A gap in the double white dwarf separation distribution caused by the common-envelope evolution: astrometric evidence from Gaia

The trajectory of the center of light of an unresolved binary is different from that of its center of mass. Binary-induced stellar centroid wobbling can therefore be detected as an excess in the goodness-of-fit of the single-star astrometric model. We use reduced $χ^2$ of the astrometric fit in the {\it Gaia} Early Data Release 3 to detect the likely unresolved double white dwarfs (DWDs). Using parallax-based distances we convert the excess of reduced $χ^2$ into the amplitude of the centroid wobble $δa$, which is proportional to the binary separation $a$. The measured $δa$ distribution drops towards larger wobble amplitudes and shows a break around $δa \approx 0.2$ where it steepens. The integral of the distribution yields DWD fraction of $6.5 \pm 3.7$ per cent in the range $0.01 < a (\text{au}) < 2$. Using synthetic models of the Galactic DWDs we demonstrate that the break in the $δa$ distribution corresponds to one side of a deep gap in the DWD separation distribution at around $a\approx 1$ au. Model DWDs with separations less than several au shrink dramatically due to (al least one) common envelope phase, reshaping the original separation distribution, clearing a gap and creating a pile-up of systems with $a\approx 0.01$ au and $δa < 0.01$. Our models reproduce the overall shape of the observed $δa$ distribution and its normalisation, however the predicted drop in the numbers of DWDs beyond the break is steeper than in the data.

preprint2022arXiv

An X-ray quiet black hole born with a negligible kick in a massive binary within the Large Magellanic Cloud

Stellar-mass black holes are the final remnants of stars born with more than 15 solar masses. Billions are expected to reside in the Local Group, yet only few are known, mostly detected through X-rays emitted as they accrete material from a companion star. Here, we report on VFTS 243: a massive X-ray faint binary in the Large Magellanic Cloud. With an orbital period of 10.4-d, it comprises an O-type star of 25 solar masses and an unseen companion of at least nine solar masses. Our spectral analysis excludes a non-degenerate companion at a 5-sigma confidence level. The minimum companion mass implies that it is a black hole. No other X-ray quiet black hole is unambiguously known outside our Galaxy. The (near-)circular orbit and kinematics of VFTS 243 imply that the collapse of the progenitor into a black hole was associated with little or no ejected material or black-hole kick. Identifying such unique binaries substantially impacts the predicted rates of gravitational-wave detections and properties of core-collapse supernovae across the Cosmos.

preprint2022arXiv

New Clues to the Evolution of Dwarf Carbon Stars From Their Variability and X-ray Emission

As main-sequence stars with C$>$O, dwarf carbon (dC) stars are never born alone but inherit carbon-enriched material from a former asymptotic giant branch (AGB) companion. In contrast to M dwarfs in post-mass transfer binaries, C$_2$ and/or CN molecular bands allow dCs to be identified with modest-resolution optical spectroscopy, even after the AGB remnant has cooled beyond detectability. Accretion of substantial material from the AGB stars should spin up the dCs, potentially causing a rejuvenation of activity detectable in X-rays. Indeed, a few dozen dCs have recently been found to have photometric variability with periods under a day. However, most of those are likely post-common-envelope binaries (PCEBs), spin-orbit locked by tidal forces, rather than solely spun-up by accretion. Here, we study the X-ray properties of a sample of the five nearest known dCs with $Chandra$. Two are detected in X-rays, the only two for which we also detected short-period photometric variability. We suggest that the coronal activity detected so far in dCs is attributable to rapid rotation due to tidal locking in short binary orbits after a common-envelope phase, late in the thermally pulsing (TP) phase of the former C-AGB primary (TP-AGB).

preprint2022arXiv

Observationally driven Galactic double white dwarf population for LISA

Realistic models of the Galactic double white dwarf (DWD) population are crucial for testing and quantitatively defining the science objectives of the Laser Interferometer Space Antenna (LISA), a future European Space Agency&#39;s gravitational-wave observatory. In addition to numerous individually detectable DWDs, LISA will also detect an unresolved confusion foreground produced by the underlying Galactic population, which will affect the detectability of all LISA sources at frequencies below a few mHz. So far, the modelling of the DWD population for LISA has been based on binary population synthesis (BPS) techniques. The aim of this study is to construct an observationally driven population. To achieve this, we employ a model developed by Maoz, Hallakoun and Badenes (2018) for the statistical analysis of the local DWD population using two complementary large, multi-epoch, spectroscopic samples: the Sloan Digital Sky Survey (SDSS), and the Supernova Ia Progenitor surveY (SPY). We calculate the number of LISA-detectable DWDs and the Galactic confusion foreground, based on their assumptions and results. We find that the observationally driven estimates yield 1) 2 - 5 times more individually detectable DWDs than various BPS forecasts, and 2) a significantly different shape of the DWD confusion foreground. Both results have important implications for the LISA mission. A comparison between several variations to our underlying assumptions shows that our observationally driven model is robust, and that the uncertainty on the total number of LISA-detectable DWDs is in the order of 20 per cent.

preprint2022arXiv

Planet Hunters TESS IV: A massive, compact hierarchical triple star system TIC 470710327

We report the discovery and analysis of a massive, compact, hierarchical triple system (TIC 470710327) initially identified by citizen scientists in data obtained by NASA&#39;s Transiting Exoplanet Survey Satellite (TESS). Spectroscopic follow-up observations obtained with the HERMES spectrograph, combined with eclipse timing variations (ETVs), confirm that the system is comprised of three OB stars, with a compact 1.10 d eclipsing binary and a non-eclipsing tertiary on a 52.04 d orbit. Dynamical modelling of the system (from radial velocity and ETVs) reveal a rare configuration wherein the tertiary star (O9.5-B0.5V; 14-17 M$_{\odot}$) is more massive than the combined mass of the inner binary (10.9-13.2 M$_{\odot}$). Given the high mass of the tertiary, we predict that this system will undergo multiple phases of mass transfer in the future, and likely end up as a double neutron star gravitational wave progenitor or an exotic Thorne-Zytkow object. Further observational characterisation of this system promises constraints on both formation scenarios of massive stars as well as their exotic evolutionary end-products.

preprint2021arXiv

Near-Chandrasekhar-Mass Type Ia Supernovae from the Double-Degenerate Channel

Recent observational evidence has demonstrated that white dwarf (WD) mergers are a highly efficient mechanism for mass accretion onto WDs in the galaxy. In this paper, we show that WD mergers naturally produce highly-magnetized, uniformly-rotating WDs, including a substantial population within a narrow mass range close to the Chandrasekhar mass ($M_{\rm Ch}$). These near-$M_{\rm Ch}$ WD mergers subsequently undergo rapid spin up and compression on a $\sim 10^2$ yr timescale, either leading to central ignition and a normal SN Ia via the DDT mechanism, or alternatively to a failed detonation and SN Iax through pure deflagration. The resulting SNe Ia and SNe Iax will have spectra, light curves, polarimetry, and nucleosynthetic yields similar to those predicted to arise through the canonical near-$M_{\rm Ch}$ single degenerate (SD) channel, but with a $t^{-1}$ delay time distribution characteristic of the double-degenerate (DD) channel. Furthermore, in contrast to the SD channel, WD merger near-$M_{\rm Ch}$ SNe Ia and SNe Iax will not produce observable companion signatures. We discuss a range of implications of these findings, from SNe Ia explosion mechanisms, to galactic nucleosynthesis of iron peak elements including manganese.

preprint2021arXiv

Spatially Resolved X-ray Study of Supernova Remnant G306.3$-$0.9 with Unusually High Calcium Abundance

G306.3$-$0.9 is an asymmetric Galactic supernova remnant (SNR), whose progenitor has been thought to be a Type Ia supernova (SN), but its high Ca abundance appears inconsistent with the Type Ia origin. Hoping to uncover the reason for its asymmetry and the origin of this SNR, we performed a spatially resolved X-ray spectroscopic analysis of XMM-Newton and Chandra observation data. We divided the SNR into 13 regions and analyzed the spectra using two-temperature models (0.2 keV + 1 keV). Compared to the southwestern regions, the northeastern regions have higher metal abundances and a lower gas density. This suggests that the asymmetric morphology results from the non-uniform ambient environment. We found that neither Type Ia nor core-collapse SN models can account for the abnormally high abundance ratios of Ar/Si, Ca/Si, or the shape of the abundance curve. A comparison with the Ca-rich transient models shows that G306.3$-$0.9 is likely to be the first identified Galactic &#34;Ca-rich transient&#34; remnant, although the theoretical production of element S is lower. We also note that the conclusion for the SNR&#39;s origin relies on the measured abundance ratios and existing nucleosynthesis models. Between two groups of Ca-rich transient explosion models, we prefer the He shell detonation for an accreting WD, rather than the merger of a white dwarf and a neutron star.

preprint2021arXiv

Weighing Milky Way Satellites with LISA

White dwarf stars are a well-established tool for studying Galactic stellar populations. Two white dwarfs in a tight binary system offer us an additional messenger - gravitational waves - for exploring the Milky Way and its immediate surroundings. Gravitational waves produced by double white dwarf (DWD) binaries can be detected by the future Laser Interferometer Space Antenna (LISA). Numerous and widespread DWDs have the potential to probe shapes, masses and formation histories of the stellar populations in the Galactic neighbourhood. In this work we outline a method for estimating the total stellar mass of Milky Way satellite galaxies based on the number of DWDs detected by LISA. To constrain the mass we perform a Bayesian inference using binary population synthesis models and considering the number of detected DWDs associated with the satellite and the measured distance to the satellite as the only inputs. Using a fiducial binary population synthesis model we find that for large satellites the stellar masses can be recovered to within 1) a factor two if the star formation history is known and 2) an order of magnitude when marginalising over different star formation history models. For smaller satellites we can place upper limits on their stellar mass. Gravitational wave observations can provide mass measurements for large satellites that are comparable, and in some cases more precise, than standard electromagnetic observations.

preprint2020arXiv

$Gaia$ white dwarfs within 40 pc II: the volume-limited northern hemisphere sample

We present an overview of the sample of northern hemisphere white dwarfs within 40 pc of the Sun detected from $Gaia$ Data Release 2 (DR2). We find that 521 sources are spectroscopically confirmed degenerate stars, 111 of which were first identified as white dwarf candidates from $Gaia$ DR2 and followed-up recently with the William Herschel Telescope and Gran Telescopio Canarias. Three additional white dwarf candidates remain spectroscopically unobserved and six unresolved binaries are known to include a white dwarf but were not in our initial selection of white dwarfs in the $Gaia$ DR2 Hertzsprung-Russell diagram (HRD). Atmospheric parameters are calculated from $Gaia$ and Pan-STARRS photometry for all objects in the sample, confirming most of the trends previously observed in the much smaller 20 pc sample. Local white dwarfs are overwhelmingly consistent with Galactic disc kinematics, with only four halo candidates. We find that DAZ white dwarfs are significantly less massive than the overall DA population ($\overline{M}_\mathrm{DAZ} = 0.59\,\mathrm{M}_\odot$, $\overline{M}_\mathrm{DA} = 0.66\,\mathrm{M}_\odot$). It may suggest that planet formation is less efficient at higher mass stars, producing more massive white dwarfs. We detect a sequence of crystallised white dwarfs in the mass range from $0.6\ \lesssim M/\mathrm{M}_\odot \lesssim\ 1.0\,$ and find that the vast majority of objects on the sequence have standard kinematic properties that correspond to the average of the sample, suggesting that their nature can be explained by crystallisation alone. We also detect 56 wide binaries including a white dwarf and 26 double degenerates.

preprint2020arXiv

Double white dwarf merger products among high-mass white dwarfs

Double white dwarf (double-WD) binaries may merge within a Hubble time and produce high-mass WDs. Compared to other high-mass WDs, the double-WD merger products have higher velocity dispersion because they are older. With the power of Gaia data, we show strong evidence for double-WD merger products among high-mass WDs by analyzing the transverse-velocity distribution of more than a thousand high-mass WDs (0.8--1.3 $M_\odot$). We estimate that the fraction of double-WD merger products in our sample is about 20 %. We also obtain a precise double-WD merger rate and its mass dependence. Our merger rate estimates are close to binary population synthesis results and support the idea that double-WD mergers may contribute to a significant fraction of type Ia supernovae.

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

Massive Stellar Triples Leading to Sequential Binary Black-Hole Mergers in the Field

Stellar triples with massive stellar components are common, and can lead to sequential binary black-hole mergers. Here, we outline the evolution towards these sequential mergers, and explore these events in the context of gravitational-wave astronomy and the pair-instability mass gap. We find that binary black-hole mergers in the pair-instability mass gap can be of triple origin and therefore are not exclusively formed in dense dynamical environments. We discuss the sequential merger scenario in the context of the most massive gravitational-wave sources detected to date: GW170729 and GW190521. We propose that the progenitor of GW170729 is a low-metallicity field triple. We support the premise that GW190521 could not have been formed in the field. We conclude that triple stellar evolution is fundamental in the understanding of gravitational-wave sources, and likely, other energetic transientsas well.

preprint2020arXiv

Milky Way Satellites Shining Bright in Gravitational Waves

The population of Milky Way satellite galaxies is of great interest for cosmology, fundamental physics, and astrophysics. They represent the faint end of the galaxy luminosity function, are the most dark-matter dominated objects in the local Universe, and contain the oldest and most metal-poor stellar populations. Recent surveys have revealed around 60 satellites, but this could represent less than half of the total. Characterization of these systems remains a challenge due to their low luminosity. We consider the gravitational wave observatory LISA as a potential tool for studying these satellites through observations of their short-period double white dwarf populations. LISA will observe the entire sky without selection effects due to dust extinction, complementing optical surveys, and could potentially discover massive satellites hidden behind the disk of the galaxy.

preprint2020arXiv

The nearest discovered black hole is likely not in a triple configuration

HR 6819 was recently claimed to be a hierarchical triple system of a Be star in a wide orbit around an inner binary system of a black hole (BH) and a B III type star. We argue that this system is unlikely to be a hierarchical triple due to three reasons: (i) Given that this system is discovered in a magnitude limited Bright Star Catalog, the expected number of such systems in the Milky Way amounts to about $10^4$ while the estimate for the MW budget for such systems is between $10^2-10^3$ systems under generous assumptions. Such a large gap cannot be reconciled as it would otherwise likely overflow the MW budget for BHs; (ii) The dynamical stability of this system sets lower bounds on the orbital separation of the outer Be star, while it not being resolved by Gaia places an upper limit on its projected sky separation. We show that these two constraints would imply a narrow range for the outer orbit without resorting to geometrical fine-tuning; (iii) The triple system should have survived the stellar evolution prior to the formation of the BH in the inner binary. We perform numerical simulations starting with conservative initial conditions of this system and show that a small parameter space for BH progenitor star&#39;s mass loss, BH natal kicks, and initial orbital separation can reproduce HR 6819. Therefore, we propose this system is a chance superposition of a Be star with a binary.

preprint2019arXiv

An Empirical Fit for Viscoelastic Simulations of Tertiary Tides

Tertiary tides (TTs), or the continuous tidal distortion of the tertiary in a hierarchical triple system, can extract energy from the inner binary, inducing within it a proclivity to merge. Despite previous work on the subject, which established that it is significant for certain close triple systems, it is still not a well-understood process. A portion of our ignorance in this regard stems from our inability to integrate a simulation of this phenomenon into conventional stellar evolution codes, since full calculations of these tidal interactions are computationally expensive on stellar evolution timescales. Thus, to attain a better understanding of how these TTs act on longer timescales, an empirical expression of its effects as a function of parameters of the triple system involved is required. In our work, we evaluate the rate at which TTs extract energy from the inner binary within a series of constructed hierarchical triple systems under varying parameters, and study the rate at which the inner binary orbital separation shrinks as a function of those parameters. We find that this rate varies little with the absolute values of the masses of the three component objects, but is very sensitive to the mass ratio of the inner binary $q$, the tertiary radius $R_{\rm 3}$, the inner binary orbital separation $a_{\rm 1}$, the outer orbital separation $a_{\rm 2}$, and the viscoelastic relaxation time of the tertiary $τ$. More specifically, we find that the percentage by which $a_{\rm 1}$ shrinks per unit time can be reasonably approximated by (1/$a_{\rm 1}$)(d$a_{\rm 1}$/d$t$)=$\left(2.22{\times}10^{-8}{\rm yrs}^{-1}\right)4q\left(1+q\right)^{-2}(R_{\rm 3}/100{\rm R}_{\odot})^{5.2}(a_{\rm 1}/0.2{\rm AU})^{4.8}(a_{\rm 2}/2{\rm AU})^{-10.2}$ $(τ/0.534{\rm yrs})^{-1.0}$. We also provide tests of how precise this fitting function is.

preprint2008arXiv

The Tremaine-Weinberg method for pattern speeds using H-alpha emission from ionized gas

The Fabry-Perot interferometer FaNTOmM was used at the 3.6m Canada France Hawaii Telescope and the 1.6m Mont Megantic Telescope to obtain data cubes in H-alpha of 9 nearby spiral galaxies from which maps in integrated intensity, velocity, and velocity dispersion were derived. We then applied the Tremaine-Weinberg method, in which the pattern speed can be deduced from its velocity field, by finding the integrated value of the mean velocity along a slit parallel to the major axis weighted by the intensity and divided by the weighted mean distance of the velocity points from the tangent point measured along the slit. The measured variables can be used either to make separate calculations of the pattern speed and derive a mean, or in a plot of one against the other for all the points on all slits, from which a best fit value can be derived. Linear fits were found for all the galaxies in the sample. For two galaxies a clearly separate inner pattern speed with a higher value, was also identified and measured.