Source author record

Selma E. de Mink

Selma E. de Mink 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
6topics
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

Carbon from massive binary-stripped stars over cosmic time: effect of metallicity

The origin of carbon in the Universe remains uncertain. At solar metallicity, binary-stripped massive stars -- stars that lost their envelope through stable interaction with a companion -- have been suggested to produce twice as much carbon as their single-star counterparts. However, understanding the chemical evolution of galaxies over cosmic time requires examining stellar yields across a range of metallicities. Using the stellar evolution code MESA, we compute the carbon yields from wind mass loss and supernova explosions of single and binary-stripped stars across a wide range of initial masses ($10$-$46,M_\odot$), metallicities ($Z = 0.0021$, $0.0047$, $0.0142$), and initial orbital periods ($10$-$5000$ days). We find that metallicity is the dominant factor influencing the carbon yields of massive stars, outweighing the effects of binarity and detailed orbital parameters. Since the chemical yields from binary massive stars are highly sensitive to metallicity, we caution that yields predicted at solar metallicity should not be directly extrapolated to lower metallicities. At sub-solar metallicities, particularly below $1/7$ solar, weak stellar winds and inefficient binary stripping result in carbon yields from binary-stripped stars that closely resemble those of single stars. This suggests that binary-stripped massive stars are unlikely to explain the presence of carbon-enhanced metal-poor stars or the carbon enrichment observed in high-redshift galaxies as probed by the James Webb Space Telescope. Our findings only concern the stripped stars in massive binaries. The impact of other outcomes of binary star evolution, in particular stellar mergers and accretors, remains largely unexplored but will be necessary for a full understanding of the role of massive binaries in nucleosynthesis.

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.

preprint2023arXiv

Close Encounters of Tight Binary Stars with Stellar-mass Black Holes

Strong dynamical interactions among stars and compact objects are expected in a variety of astrophysical settings, such as star clusters and the disks of active galactic nuclei. Via a suite of 3D hydrodynamics simulations using the moving-mesh code AREPO, we investigate the formation of transient phenomena and their properties in close encounters between an $2M_{\odot}$ or $20M_{\odot}$ equal-mass circular binary star and single $20M_{\odot}$ black hole (BH). Stars can be disrupted by the BH during dynamical interactions, naturally producing electromagnetic transient phenomena. Encounters with impact parameters smaller than the semimajor axis of the initial binary frequently lead to a variety of transients whose electromagnetic signatures are qualitatively different from those of ordinary disruption events involving just two bodies. These include the simultaneous or successive disruptions of both stars and one full disruption of one star accompanied by successive partial disruptions of the other star. On the other hand, when the impact parameter is larger than the semimajor axis of the initial binary, the binary is either simply tidally perturbed or dissociated into bound and unbound single stars ("micro-Hills" mechanism). The dissociation of $20M_{\odot}$ binaries can produce a runaway star and an active BH moving away from one another. Also, the binary dissociation can either produce an interacting binary with the BH, or a non-interacting, hard binary; both could be candidates of BH high- and low-mass X-ray binaries. Hence our simulations especially confirm that strong encounters can lead to the formation of the (generally difficult to form) BH low-mass X-ray binaries.

preprint2022arXiv

Formation and fate of low metallicity stars in IllustrisTNG50

Low metallicity stars give rise to unique spectacular transients and are of immense interest for understanding stellar evolution. Their importance has only grown further with the recent detections of mergers of stellar mass black holes that likely originate mainly from low metallicity progenitor systems. Moreover, the formation of low metallicity stars is intricately linked to galaxy evolution, in particular to early enrichment and to later accretion and mixing of lower metallicity gas. Because low metallicity stars are difficult to observe directly, cosmological simulations are crucial for understanding their formation. Here we quantify the rates and locations of low metallicity star formation using the high-resolution TNG50 magnetohydrodynamical cosmological simulation, and we examine where low metallicity stars end up at $z=0$. We find that $20\%$ of stars with $Z_*<0.1\,\mathrm{Z_\odot}$ form after $z=2$, and that such stars are still forming in galaxies of all masses at $z=0$ today. Moreover, most low-metallicity stars at $z=0$ reside in massive galaxies. We analyse the radial distribution of low metallicity star formation, and discuss the curious case of seven galaxies in TNG50 that form stars from primordial gas even at $z=0$.

preprint2022arXiv

LEGWORK: A python package for computing the evolution and detectability of stellar-origin gravitational-wave sources with space-based detectors

We present LEGWORK (LISA Evolution and Gravitational Wave Orbit Kit), an open-source Python package for making predictions about stellar-origin gravitational wave sources and their detectability in LISA or other space-based gravitational wave detectors. LEGWORK can be used to evolve the orbits of sources due to gravitational wave emission, calculate gravitational wave strains (using post-Newtonian approximations), compute signal-to-noise ratios and visualise the results. It can be applied to a variety of potential sources, including binaries consisting of white dwarfs, neutron stars and black holes. Although we focus on double compact objects, in principle LEGWORK can be used for any system with a user-specified orbital evolution, such as those affected by a third object or gas drag. We optimised the package to make it efficient for use in population studies which can contain tens-of-millions of sources. This paper describes the package and presents several potential use cases. We explain in detail the derivations of the expressions behind the package as well as identify and clarify some discrepancies currently present in the literature. We hope that LEGWORK will enable and accelerate future studies triggered by the rapidly growing interest in gravitational wave sources.

preprint2022arXiv

Stellar mergers as the origin of the blue main-sequence band in young star clusters

Recent high-quality Hubble Space Telescope (HST) photometry shows that the main sequences (MS) stars of young star clusters form two discrete components in the color-magnitude diagram (CMD). Based on their distribution in the CMD, we show that stars of the blue MS component can be understood as slow rotators originating from stellar mergers. We derive the masses of the blue MS stars, and find that they follow a nearly flat mass function, which supports their unusual formation path. Our results imply that the cluster stars gain their mass in two different ways, by disk accretion leading to rapid rotation, contributing to the red MS, or by binary merger leading to slow rotation and populating the blue MS. We also derive the approximate merger time of the individual stars of the blue MS component, and find a strong early peak in the merger rate, with a lower level merger activity prevailing for tens of Myr. This supports recent binary formation models, and explains new velocity dispersion measurements for members of young star clusters. Our findings shed new light on the origin of the bi-modal mass, spin, and magnetic field distributions of main-sequence stars.

preprint2022arXiv

Successful Common Envelope Ejection and Binary Neutron Star Formation in 3D Hydrodynamics

A binary neutron star merger has been observed in a multi-messenger detection of gravitational wave (GW) and electromagnetic (EM) radiation. Binary neutron stars that merge within a Hubble time, as well as many other compact binaries, are expected to form via common envelope evolution. Yet five decades of research on common envelope evolution have not yet resulted in a satisfactory understanding of the multi-spatial multi-timescale evolution for the systems that lead to compact binaries. In this paper, we report on the first successful simulations of common envelope ejection leading to binary neutron star formation in 3D hydrodynamics. We simulate the dynamical inspiral phase of the interaction between a 12$M_\odot$ red supergiant and a 1.4$M_\odot$ neutron star for different initial separations and initial conditions. For all of our simulations, we find complete envelope ejection and final orbital separations of $a_{\rm f} \approx 1.3$-$5.1 R_\odot$ depending on the simulation and criterion, leading to binary neutron stars that can merge within a Hubble time. We find $α_{\rm CE}$-equivalent efficiencies of $\approx 0.1$-$2.7$ depending on the simulation and criterion, but this may be specific for these extended progenitors. We fully resolve the core of the star to $\lesssim 0.005 R_\odot$ and our 3D hydrodynamics simulations are informed by an adjusted 1D analytic energy formalism and a 2D kinematics study in order to overcome the prohibitive computational cost of simulating these systems. The framework we develop in this paper can be used to simulate a wide variety of interactions between stars, from stellar mergers to common envelope episodes leading to GW sources.

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

Candidate LBV stars in galaxy NGC 7793 found via HST photometry + MUSE spectroscopy

Only about 19 Galactic and 25 extra-galactic bona-fide Luminous Blue Variables (LBVs) are known to date. This incomplete census prevents our understanding of this crucial phase of massive star evolution which leads to the formation of heavy binary black holes via the classical channel. With large samples of LBVs one could better determine the duration and maximum stellar luminosity which characterize this phase. We search for candidate LBVs (cLBVs) in a new galaxy, NGC 7793. For this purpose, we combine high spatial resolution images from two Hubble Space Telescope (HST) programs with optical spectroscopy from the Multi Unit Spectroscopic Explorer (MUSE). By combining PSF-fitting photometry measured on F547M, F657N, and F814W images, with restrictions on point-like appearance (at HST resolution) and Hα luminosity, we find 100 potential cLBVs, 36 of which fall in the MUSE fields. Five of the latter 36 sources are promising cLBVs which have M(V) less than or equal to -7 and a combination of: Hα with a P-Cygni profile; no [O I] 6300 emission; weak or no [O III] 5007 emission; large [N II]/Hα relative to H II regions; and [S II] 6716 / [S II] 6731 similar to 1. It is not clear if these five cLBVs are isolated from O-type stars, which would favor the binary formation scenario of LBVs. Our study, which approximately covers one fourth of the optical disc of NGC 7793, demonstrates how by combining the above HST surveys with multi-object spectroscopy from 8-m class telescopes, one can efficiently find large samples of cLBVs in nearby galaxies.

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

How stellar rotation shapes the colour magnitude diagram of the massive intermediate-age star cluster NGC 1846

We present a detailed study of stellar rotation in the massive 1.5 Gyr old cluster NGC 1846 in the Large Magellanic Cloud. Similar to other clusters at this age, NGC 1846 shows an extended main sequence turn-off (eMSTO), and previous photometric studies have suggested it could be bimodal. In this study, we use MUSE integral-field spectroscopy to measure the projected rotational velocities (vsini) of around 1400 stars across the eMSTO and along the upper main sequence of NGC 1846. We measure vsini values up to ~250 km/s and find a clear relation between the vsini of a star and its location across the eMSTO. Closer inspection of the distribution of rotation rates reveals evidence for a bimodal distribution, with the fast rotators centred around vsini = 140 km/s and the slow rotators centred around vsini = 60 km/s. We further observe a lack of fast rotating stars along the photometric binary sequence of NGC 1846, confirming results from the field that suggest that tidal interactions in binary systems can spin down stars. However, we do not detect a significant difference in the binary fractions of the fast and slowly rotating sub-populations. Finally, we report on the serendipitous discovery of a planetary nebula associated with NGC 1846.

preprint2020arXiv

HST astrometry in the Orion Nebula Cluster: census of low-mass runaways

We present a catalog of high-precision proper motions in the Orion Nebula Cluster (ONC), based on Treasury Program observations with the Hubble Space Telescope's (HST) ACS/WFC camera. Our catalog contains 2,454 objects in the magnitude range of $14.2<m_{\rm F775W}<24.7$, thus probing the stellar masses of the ONC from $\sim$0.4 $M_\odot$ down to $\sim$0.02 $M_\odot$ over an area of $\sim$550 arcmin$^2$. We provide a number of internal velocity dispersion estimates for the ONC that indicate a weak dependence on the stellar location and mass. There is good agreement with the published velocity dispersion estimates, although nearly all of them (including ours at $σ_{v,x}=0.94$ and $σ_{v,y}=1.25$ mas yr$^{-1}$) might be biased by the overlapping young stellar populations of Orion A. We identified 4 new ONC candidate runaways based on HST and the Gaia DR2 data, all with masses less than $\sim$1 $M_\odot$. The total census of known candidate runaway sources is 10 -- one of the largest samples ever found in any Milky Way open star cluster. Surprisingly, none of them has the tangential velocity exceeding 20 km s$^{-1}$. If most of them indeed originated in the ONC, it may compel re-examination of dynamical processes in very young star clusters. It appears that the mass function of the ONC is not significantly affected by the lost runaways.

preprint2020arXiv

HST survey of the Orion Nebula Cluster in the H$_2$O 1.4 $μ$m absorption band: I. A census of substellar and planetary mass objects

In order to obtain a complete census of the stellar and sub-stellar population, down to a few M$_{Jup}$ in the $\sim1$ Myr old Orion Nebula Cluster, we used the infrared channel of the Wide Field Camera 3 of the Hubble Space Telescope with the F139M and F130N filters. These bandpasses correspond to the $1.4 μ$m H$_2$O absorption feature and an adjacent line-free continuum region. Out of $4,504$ detected sources, $3,352$ (about $75\%$) appear fainter than m$_{130}=14$ (Vega mag) in the F130N filter, a brightness corresponding to the hydrogen-burning limit mass (M$\simeq 0.072 M_\odot$) at $\sim 1$ Myr. Of these, however, only $742$ sources have a negative F130M-139N color index, indicative of the presence of H$_2$O vapor in absorption, and can therefore be classified as bona-fide M and L dwarfs, with effective temperatures T$\lesssim 2850$ K at an assumed $1$ Myr cluster age. On our color-magnitude diagram, this population of sources with H$_2$O absorption appears clearly distinct from the larger background population of highly reddened stars and galaxies with positive F130M-F139N color index, and can be traced down to the sensitivity limit of our survey, m$_{130}\simeq 21.5$, corresponding to a $1$ Myr old $\simeq 3 $M$_{Jup}$, planetary mass object under about 2 magnitudes of visual extinction. Theoretical models of the BT-Settl family predicting substellar isochrones of $1, 2$ and $3$ Myr (down to $\sim 1 $M$_{Jup}$) fail to reproduce the observed H$_2$O color index at M$\lesssim 20 $M$_{Jup}$. We perform a Bayesian analysis to determine extinction, mass and effective temperature of each sub-stellar member of our sample, together with its membership probability.

preprint2016arXiv

Large proper motion of the Thorne-Żytkow object candidate HV 2112 reveals its likely nature as foreground Galactic S-star

Using the Southern Proper Motion (SPM) catalog, we show that the candidate Thorne-Żytkow object HV~2112 has a proper motion implying a space velocity of about 3000\kms if the object is located at the distance of the Small Magellanic Cloud. The proper motion is statistically different from that of the SMC at approximately $4σ$ in SPM, although the result can drop to about $3σ$ significance by including the UCAC4 data and considering systematic uncertainties in addition to the statistical ones. Assuming the measurement is robust, this proper motion is sufficient to exclude its proposed membership of the Small Magellanic Cloud and to argue instead that it is likely to be a foreground star in the Milky Way halo. The smaller distance and therefore lower brightness argue against its proposed nature as a Thorne-Żytkow object (the hypothesized star-like object formed when a normal star and a neutron star merge) or a super-Asymptotic Giant Branch (AGB) star. Instead we propose a binary scenario where this star is the companion of a former massive AGB star, which polluted the object with via its stellar wind, i.e. a special case of an extrinsic S star. Our new scenario solves two additional problems with the two existing scenarios for its nature as Thorne-Żytkow object or present-day super-AGB star. The puzzling high ratio of the strength of calcium to iron absorption lines is unexpected for SMC supergiants, but is fully consistent with the expectations for halo abundances. Secondly, its strong variability can now be explained naturally as a manifestation of the Mira phenomenon.

preprint2016arXiv

Merging binary black holes formed through chemically homogeneous evolution in short-period stellar binaries

We explore a newly proposed channel to create binary black holes of stellar origin. This scenario applies to massive, tight binaries where mixing induced by rotation and tides transports the products of hydrogen burning throughout the stellar envelopes. This slowly enriches the entire star with helium, preventing the build-up of an internal chemical gradient. The stars remain compact as they evolve nearly chemically homogeneously, eventually forming two black holes, which, we estimate, typically merge 4--11 Gyr after formation. Like other proposed channels, this evolutionary pathway suffers from significant theoretical uncertainties, but could be constrained in the near future by data from advanced ground-based gravitational-wave detectors. We perform Monte Carlo simulations of the expected merger rate over cosmic time to explore the implications and uncertainties. Our default model for this channel yields a local binary black hole merger rate of about $10$ Gpc$^{-3}$ yr$^{-1}$ at redshift $z=0$, peaking at twice this rate at $z=0.5$. This means that this channel is competitive, in terms of expected rates, with the conventional formation scenarios that involve a common-envelope phase during isolated binary evolution or dynamical interaction in a dense cluster. The events from this channel may be distinguished by the preference for nearly equal-mass components and high masses, with typical total masses between 50 and 110 $\textrm{M}_\odot$. Unlike the conventional isolated binary evolution scenario that involves shrinkage of the orbit during a common-envelope phase, short time delays are unlikely for this channel, implying that we do not expect mergers at high redshift.

preprint2015arXiv

A High-Resolution Multiband Survey of Westerlund 2 With the Hubble Space Telescope I: Is the Massive Star Cluster Double?

We present first results from a high resolution multi-band survey of the Westerlund 2 region with the Hubble Space Telescope. Specifically, we imaged Westerlund 2 with the Advanced Camera for Surveys through the $F555W$, $F814W$, and $F658N$ filters and with the Wide Field Camera 3 in the $F125W$, $F160W$, and $F128N$ filters. We derive the first high resolution pixel-to-pixel map of the color excess $E(B-V)_g$ of the gas associated with the cluster, combining the H$α$ ($F658N$) and Pa$β$ ($F128N$) line observations. We demonstrate that, as expected, the region is affected by significant differential reddening with a median of $E(B-V)_g=1.87$~mag. After separating the populations of cluster members and foreground contaminants using a $(F814W-F160W)$ vs. $F814W$ color-magnitude diagram, we identify a pronounced pre-main-sequence population in Westerlund 2 showing a distinct turn-on. After dereddening each star of Westerlund 2 individually in the color-magnitude diagram we find via over-plotting PARSEC isochrones that the distance is in good agreement with the literature value of $\sim4.16 \pm 0.33$~kpc. With zero-age-main-sequence fitting to two-color-diagrams, we derive a value of total to selective extinction of $R_V=3.95 \pm 0.135$. A spatial density map of the stellar content reveals that the cluster might be composed of two clumps. We estimate the same age of 0.5-2.0 Myr for both clumps. While the two clumps appear to be coeval, the northern clump shows a $\sim 20 \%$ lower stellar surface density.

preprint2015arXiv

Constraints on the Binary Companion to the SN Ic 1994I Progenitor

Core-collapse supernovae (SNe), marking the deaths of massive stars, are among the most powerful explosions in the Universe, responsible, e.g., for a predominant synthesis of chemical elements in their host galaxies. The majority of massive stars are thought to be born in close binary systems. To date, putative binary companions to the progenitors of SNe may have been detected in only two cases, SNe 1993J and 2011dh. We report on the search for a companion of the progenitor of the Type Ic SN 1994I, long considered to have been the result of binary interaction. Twenty years after explosion, we used the Hubble Space Telescope to observe the SN site in the ultraviolet (F275W and F336W bands), resulting in deep upper limits on the expected companion: F275W > 26.1 mag and F336W > 24.7 mag. These allows us to exclude the presence of a main sequence companion with a mass >~ 10 Msun. Through comparison with theoretical simulations of possible progenitor populations, we show that the upper limits to a companion detection exclude interacting binaries with semi-conservative (late Case A or early Case B) mass transfer. The limits tend to favor systems with non-conservative, late Case B mass transfer with intermediate initial orbital periods and mass ratios. The most likely mass range for a putative main sequence companion would be ~5--12 Msun, the upper end of which corresponds to the inferred upper detection limit.

preprint2015arXiv

LEGUS Discovery of a Light Echo Around Supernova 2012aw

We have discovered a luminous light echo around the normal Type II-Plateau Supernova (SN) 2012aw in Messier 95 (M95; NGC 3351), detected in images obtained approximately two years after explosion with the Wide Field Channel 3 on-board the Hubble Space Telescope (HST) by the Legacy ExtraGalactic Ultraviolet Survey (LEGUS). The multi-band observations span from the near-ultraviolet through the optical (F275W, F336W, F438W, F555W, and F814W). The apparent brightness of the echo at the time was ~21--22 mag in all of these bands. The echo appears circular, although less obviously as a ring, with an inhomogeneous surface brightness, in particular, a prominent enhanced brightness to the southeast. The SN itself was still detectable, particularly in the redder bands. We are able to model the light echo as the time-integrated SN light scattered off of diffuse interstellar dust in the SN environment. We have assumed that this dust is analogous to that in the Milky Way with R_V=3.1. The SN light curves that we consider also include models of the unobserved early burst of light from the SN shock breakout. Our analysis of the echo suggests that the distance from the SN to the scattering dust elements along the echo is ~45 pc. The implied visual extinction for the echo-producing dust is consistent with estimates made previously from the SN itself. Finally, our estimate of the SN brightness in F814W is fainter than that measured for the red supergiant star at the precise SN location in pre-SN images, possibly indicating that the star has vanished and confirming it as the likely SN progenitor.

preprint2015arXiv

New Constraints on the Star Formation History of the Star Cluster NGC 1856

We use the Wide Field Camera 3 onboard the Hubble Space Telescope to obtain deep, high-resolution photometry of the young (age ~ 300 Myr) star cluster NGC1856 in the Large Magellanic Cloud. We compare the observed colour-magnitude diagram (CMD), after having applied a correction for differential reddening, with Monte Carlo simulations of simple stellar populations (SSPs) of various ages. We find that the main sequence turn-off (MSTO) region is wider than that derived from the simulation of a single SSP. Using constraints based on the distribution of stars in the MSTO region and the red clump, we find that the CMD is best reproduced using a combination of two different SSPs with ages separated by 80 Myr (0.30 and 0.38 Gyr, respectively). However, we can not formally exclude that the width of the MSTO could be due to a range of stellar rotation velocities if the efficiency of rotational mixing is higher than typically assumed. Using a King-model fit to the surface number density profile in conjunction with dynamical evolution models, we determine the evolution of cluster mass and escape velocity from an age of 10 Myr to the present age, taking into account the possible effects of primordial mass segregation. We find that the cluster has an escape velocity Vesc ~ 17 km/s at an age of 10 Myr, and it remains high enough during a period of ~ 100 Myr to retain material ejected by slow winds of first-generation stars. Our results are consistent with the presence of an age spread in NGC1856, in contradiction to the results of Bastian & Silva-Villa (2013).

preprint2013arXiv

Massive Binary Stars and Self-Enrichment of Globular Clusters

Globular clusters contain many stars with surface abundance patterns indicating contributions from hydrogen burning products, as seen in the anti-correlated elemental abundances of e.g. sodium and oxygen, and magnesium and aluminium. Multiple generations of stars can explain this phenomenon, with the second generation forming from a mixture of pristine gas and ejecta from the first generation. We show that massive binary stars may be a source of much of the material that makes this second generation of stars. Mass transfer in binaries is often non-conservative and the ejected matter moves slowly enough that it can remain inside a globular cluster and remain available for subsequent star formation. Recent studies show that there are more short-period massive binaries than previously thought, hence also more stars that interact and eject nuclear-processed material.

preprint2013arXiv

The ages of young star clusters, massive blue stragglers and the upper mass limit of stars: analysing age dependent stellar mass functions

Massive stars rapidly change their masses through strong stellar winds and mass transfer in binary systems. We show that such mass changes leave characteristic signatures in stellar mass functions of young star clusters which can be used to infer their ages and to identify products of binary evolution. We model the observed present day mass functions of the young Galactic Arches and Quintuplet star clusters using our rapid binary evolution code. We find that shaping of the mass function by stellar wind mass loss allows us to determine the cluster ages to 3.5$\pm$0.7 Myr and 4.8$\pm$1.1 Myr, respectively. Exploiting the effects of binary mass exchange on the cluster mass function, we find that the most massive stars in both clusters are rejuvenated products of binary mass transfer, i.e. the massive counterpart of classical blue straggler stars. This resolves the problem of an apparent age spread among the most luminous stars exceeding the expected duration of star formation in these clusters. We perform Monte Carlo simulations to probe stochastic sampling, which support the idea of the most massive stars being rejuvenated binary products. We find that the most massive star is expected to be a binary product after 1.0$\pm$0.7 Myr in Arches and after 1.7$\pm$1.0 Myr in Quintuplet. Today, the most massive 9$\pm$3 stars in Arches and 8$\pm$3 in Quintuplet are expected to be such objects. Our findings have strong implications for the stellar upper mass limit and solve the discrepancy between the claimed 150 $\mathrm{M}_\odot$ limit and observations of fours stars with initial masses of 165-320 $\mathrm{M}_\odot$ in R136 and of SN 2007bi, which is thought to be a pair-instability supernova from an initial 250 $\mathrm{M}_\odot$ star. Using the stellar population of R136, we revise the upper mass limit to values in the range 200-500 $\mathrm{M}_\odot$.

preprint2011arXiv

Rotating Massive Main-Sequence Stars I: Grids of Evolutionary Models and Isochrones

We present a dense grid of evolutionary tracks and isochrones of rotating massive main-sequence stars. We provide three grids with different initial compositions tailored to compare with early OB stars in the Small and Large Magellanic Clouds and in the Galaxy. Each grid covers masses ranging from 5 to 60 Msun and initial rotation rates between 0 and about 600 km/s. To calibrate our models we used the results of the VLT-FLAMES Survey of Massive Stars. We determine the amount of convective overshooting by using the observed drop in rotation rates for stars with surface gravities log g < 3.2 to determine the width of the main sequence. We calibrate the efficiency of rotationally induced mixing using the nitrogen abundance determinations for B stars in the Large Magellanic cloud. We describe and provide evolutionary tracks and the evolution of the central and surface abundances. In particular, we discuss the occurrence of quasi-chemically homogeneous evolution, i.e. the severe effects of efficient mixing of the stellar interior found for the most massive fast rotators. We provide a detailed set of isochrones for rotating stars. Rotation as an initial parameter leads to a degeneracy between the age and the mass of massive main sequence stars if determined from its observed location in the Hertzsprung-Russell diagram. We show that the consideration of surface abundances can resolve this degeneracy.

preprint2011arXiv

Rotating Massive Main-Sequence Stars II: Simulating a Population of LMC early B-type Stars as a Test of Rotational Mixing

Rotational mixing in massive stars is a widely applied concept, with far reaching consequences for stellar evolution. Nitrogen surface abundances for a large and homogeneous sample of massive B-type stars in the LMC were obtained by the VLT-FLAMES Survey of Massive Stars. This sample is the first covering a broad range of projected stellar rotational velocities, with a large enough sample of high quality data to allow for a statistically significant analysis. We use the sample to provide the first rigorous test of the theory of rotational mixing in massive stars. We calculated a grid of stellar evolution models, using the FLAMES sample to calibrate some of the uncertain mixing processes. We developed a new population-synthesis code, which uses this grid to simulate a large population of stars with masses, ages and rotational velocity distributions consistent with those from the FLAMES sample. The synthesized population is then filtered by the selection effects in the observed sample, to enable a direct comparison between the empirical results and theoretical predictions. Our simulations reproduce the fraction of stars without significant nitrogen enrichment. The predicted number of rapid rotators with enhanced nitrogen is about twice as large as found observationally. Furthermore, a group of stars consisting of slowly rotating, nitrogen-enriched objects and another consisting of rapidly rotating un-enriched objects can not be reproduced by our single-star population synthesis. Additional physical processes appear to be required to understand the population of massive main-sequence stars from the FLAMES sample.We discuss the possible role of binary stars and magnetic fields in the interpretation of our results. We find that the population of slowly rotating nitrogen-enriched stars is unlikely produced via mass transfer and subsequent tidal spin-down in close binary systems