Researcher profile

Jorick S. Vink

Jorick S. Vink contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
14works
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

14 published item(s)

preprint2025arXiv

Metal-poor single Wolf-Rayet stars: The interplay of optically thick winds and rotation

The Small Magellanic Cloud (SMC) hosts 12 known Wolf-Rayet (WR) stars, seven of which are apparently single. Their formation is a challenge for current stellar evolution models because line-driven winds are generally assumed to be quenched at a metallicity of Z < 0.004. Here, we present a set of mesa models of single stars with zero-age main sequence masses of 20 - 80 Msun considering different initial rotation speeds (Ω = 0 - 0.7 Ω_c), metallicities (Z = 0.002 - 0.0045), and wind mass-loss models (optically thin and thick winds). We show that if we account for optically thick winds, fast rotating (Ω = 0.6 Ω_c) single metal-poor O-type stars (with M > 20 Msun) shed their envelope and become WR stars even at the low metallicity of the SMC. The luminosity, effective temperature, evolutionary timescale, surface abundance, and rotational velocity of our simulated WR stars are compatible to the WRs observed in the SMC. We speculate that this scenario can also alleviate the excess of giant stars across the Humphreys-Davidson limit. Our results have key implications for black hole masses, (pair instability) supernova explosions, and other observable signatures.

preprint2023arXiv

Potential Subpopulations and Assembling Tendency of the Merging Black Holes

The origins of coalescing binary black holes (BBHs) detected by the advanced LIGO/Virgo are still under debate, and clues may be present in the joint mass-spin distribution of these merger events. Here we construct phenomenological models containing two sub-populations to investigate the BBH population detected in gravitational wave observations. We find that our models can explain the GWTC-3 data rather well, and several constraints to our model are required by the data: first, the maximum mass for the component with a stellar-origin, $m_{\rm max}$, is $39.1^{+2.4}_{-2.7}M_{\odot}$ at 90\% credibility; second, about $15\%$ of the mergers happen in dynamical environments, in which $7-16\%$ of events are hierarchical mergers, and these BHs have an average spin magnitude significantly larger than the first-generation mergers, with ${\rm d}μ_{\rm a} > 0.4 $ at $99\%$ credibility; third, the dynamical component BHs tend to pair with each other with larger total mass and higher mass ratio. An independent analysis focusing on spins is also carried out, and we find that the spin amplitude of component BHs can be divided into two groups according to a division mass $m_{\rm d} = 46.1^{+5.6}_{-5.1}M_{\odot}$. These constraints can be naturally explained by current formation channels, and our results suggest that some of the observed events were likely from AGN disks.

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

Mass-loss implementation and temperature evolution of very massive stars

Very massive stars (VMS) dominate the physics of young clusters due to their ionising radiation and extreme stellar winds. It is these winds that determine their lifepaths until expiration. Observations in the Arches cluster show that VMS all have similar temperatures. The VLT-Flames Tarantula survey analysed VMS in the 30 Dor region of the LMC also finding a narrow range of temperatures, albeit at higher values - likely a metallicity effect. Using MESA, we study the main-sequence evolution of VMS with a new mass-loss recipe that switches from optically-thin O-star winds to optically-thick Wolf-Rayet type winds through the model-independent transition mass-loss rate of Vink & Gräfener. We examine the temperature evolution of VMS with mass loss that scales with the luminosity-over-mass (L/M) ratio and the Eddington parameter ($Γ_{\rm e}$), assessing the relevance of the surface hydrogen (H) abundance which sets the number of free electrons. We present grids of VMS models at Galactic and LMC metallicity and compare our temperature predictions with empirical results. Models with a steep $Γ_{\rm e}$-dependence evolve horizontally in the Hertzsprung-Russel (HR) diagram at nearly constant luminosities, requiring a delicate and unlikely balance between envelope inflation and enhanced mass loss over the entire VMS mass range. By contrast, models with a steep L/M-dependent mass loss are shown to evolve vertically in the HR-diagram at nearly constant Teff, naturally reproducing the narrow range of observed temperatures, as well as the correct trend with metallicity. This distinct behavior of a steeply dropping luminosity is a self-regulatory mechanism that keeps temperatures constant during evolution in the HR-diagram.

preprint2022arXiv

Resolving the core of R136 in the optical

The sharpest optical images of the R136 cluster in the Large Magellanic Cloud are presented, allowing for the first time to resolve members of the central core, including R136a1, the most massive star known. These data were taken using the Gemini speckle imager Zorro in medium-band filters with effective wavelengths similar to BVRI achieving angular resolutions between 30-40 mas. All stars previously known in the literature, having $V<$16 mag within the central $2&#39;&#39;\times2&#39;&#39;$ were recovered. Visual companions ($\geq$40 mas; 2000 au) were detected for the WN5h stars R136 a1, and a3. Photometry of the visual companion of a1 suggests it is of mid O spectral type. Based on new photometric luminosities using the resolved Zorro imaging, the masses of the individual WN5h stars are estimated to be between 150-200 $M_{\odot}$, lowering significantly the present-day masses of some of the most massive stars known. These mass estimates are critical anchor points for establishing the stellar upper-mass function.

preprint2022arXiv

Rubin Observatory LSST Transients and Variable Stars Roadmap

The Vera C. Rubin Legacy Survey of Space and Time holds the potential to revolutionize time domain astrophysics, reaching completely unexplored areas of the Universe and mapping variability time scales from minutes to a decade. To prepare to maximize the potential of the Rubin LSST data for the exploration of the transient and variable Universe, one of the four pillars of Rubin LSST science, the Transient and Variable Stars Science Collaboration, one of the eight Rubin LSST Science Collaborations, has identified research areas of interest and requirements, and paths to enable them. While our roadmap is ever-evolving, this document represents a snapshot of our plans and preparatory work in the final years and months leading up to the survey&#39;s first light.

preprint2022arXiv

The R136 star cluster dissected with Hubble Space Telescope/STIS. III. The most massive stars and their clumped winds

Context: The star cluster R136 inside the LMC hosts a rich population of massive stars, including the most massive stars known. The strong stellar winds of these very luminous stars impact their evolution and the surrounding environment. We currently lack detailed knowledge of the wind structure that is needed to quantify this impact. Aims: To observationally constrain the stellar and wind properties of the massive stars in R136, in particular the parameters related to wind clumping. Methods: We simultaneously analyse optical and UV spectroscopy of 53 O-type and 3 WNh-stars using the FASTWIND model atmosphere code and a genetic algorithm. The models account for optically thick clumps and effects related to porosity and velocity-porosity, as well as a non-void interclump medium. Results: We obtain stellar parameters, surface abundances, mass-loss rates, terminal velocities and clumping characteristics and compare these to theoretical predictions and evolutionary models. The clumping properties include the density of the interclump medium and the velocity-porosity of the wind. For the first time, these characteristics are systematically measured for a wide range of effective temperatures and luminosities. Conclusions: We confirm a cluster age of 1.0-2.5 Myr and derive an initial stellar mass of $\geq 250 {\rm M}_\odot$ for the most massive star in our sample, R136a1. The winds of our sample stars are highly clumped, with an average clumping factor of $f_{\rm cl}=29\pm15$. We find tentative trends in the wind-structure parameters as a function of mass-loss rate, suggesting that the winds of stars with higher mass-loss rates are less clumped. We compare several theoretical predictions to the observed mass-loss rates and terminal velocities and find that none satisfactorily reproduces both quantities. The prescription of Krtička & Kubát (2018) matches best the observed mass-loss rates.

preprint2022arXiv

Theory and Diagnostics of Hot Star Mass Loss

Massive stars have strong stellar winds that direct their evolution through the upper Hertzsprung-Russell diagram and determine the black hole mass function. Secondly, wind strength dictates the atmospheric structure that sets the ionising flux. Thirdly, the wind directly intervenes with the stellar envelope structure, which is decisive for both single star and binary evolution, affecting predictions for gravitational wave events. Key findings of current hot-star research include: * The traditional line-driven wind theory is being updated with Monte Carlo and co-moving frame computations, revealing a rich multi-variate behaviour of the mass-loss rate dM/dt in terms of M, L, Eddington Gamma, Teff , and chemical composition Z. Concerning the latter, dM/dt is shown to depend on the iron (Fe) opacity, making Wolf-Rayet populations, and gravitational wave events dependent on host galaxy Z. * On top of smooth mass-loss behaviour, there are several transitions in the Hertzsprung-Russell diagram, involving bi-stability jumps around Fe recombination temperatures, leading to quasi-stationary episodic, and not necessarily eruptive, Luminous Blue Variable and pre-SN mass loss. * Moreover, there are kinks. At 100 Solar Masses, a high Gamma mass-loss transition implies that hydrogen-rich very massive stars have higher mass-loss rates than commonly considered. At the other end of the mass spectrum, low-mass stripped helium stars no longer appear as Wolf-Rayet stars, but as optically-thin stars. These stripped stars, in addition to very massive stars, are two newly identified sources of ionising radiation that could play a key role in local star formation as well as at high-redshift.

preprint2021arXiv

Ultraviolet Spectropolarimetry with Polstar: Clumping and Mass-loss Rate Corrections

The most massive stars are thought to lose a significant fraction of their mass in a steady wind during the main-sequence and blue supergiant phases. This in turn sets the stage for their further evolution and eventual supernova, with consequences for ISM energization and chemical enrichment. Understanding these processes requires accurate observational constraints on the mass-loss rates of the most luminous stars, which can also be used to test theories of stellar wind generation. In the past, mass-loss rates have been characterized via collisional emission processes such as H$α$ and free-free radio emission, but these so-called &#34;density squared&#34; diagnostics require correction in the presence of widespread clumping. Recent observational and theoretical evidence points to the likelihood of a ubiquitously high level of such clumping in hot-star winds, but quantifying its effects requires a deeper understanding of the complex dynamics of radiatively driven winds. Furthermore, large-scale structures arising from surface anisotropies and propagating throughout the wind can further complicate the picture by introducing further density enhancements, affecting mass-loss diagnostics. Time series spectroscopy of UV resonance lines with high resolution and high signal-to-noise are required to better understand this complex dynamics, and help correct &#34;density squared&#34; diagnostics of mass-loss rates. The proposed Polstar mission easily provides the necessary resolution at the sound-speed scale of 20 km s$^{-1}$, on three dozen bright targets with signal-to noise an order of magnitude above that of the celebrated IUE MEGA campaign, via continuous observations that track structures advecting through the wind in real time.

preprint2020arXiv

A theoretical investigation of the Humphreys-Davidson limit at high and low metallicity

Current massive star evolution grids are not able to simultaneously reproduce the empirical upper luminosity limit of red supergiants, the Humphrey-Davidson (HD) limit at high and low metallicity. In this study, we provide a better understanding of what drives massive star evolution to blue and red supergiant phases, with the ultimate aim of reproducing the HD limit at varied metallicities (Z). For solar, LMC, and SMC Z, we develop eight grids of MESA models for the mass range 20-60M to probe the effect of semiconvection and overshooting. We compare rotating and non-rotating models with efficient (alpha_semi = 100) and inefficient semi-convection (alpha_semi = 0.1), with high and low core overshooting (alpha_ov of 0.1 or 0.5). The red and blue supergiant evolutionary phases are investigated by comparing the fraction of core He-burning lifetimes spent in each phase. We find that the extension of the convective core by overshooting alpha_ov = 0.5 has an effect on the post-MS evolution which can disable semiconvection leading to more RSGs, but a lack of BSGs. We therefore implement alpha_ov = 0.1 which switches on semiconvective mixing, though for standard alpha_semi = 1, would result in an HD limit which is higher than observed at low Z. Therefore, we need to implement very efficient semiconvection of alpha_semi = 100 which reproduces the HD limit at log L ~ 5.5 for the Magellanic Clouds while simultaneously reproducing the Galactic HD limit of log L ~ 5.8 naturally. The effect of semiconvection is not active at high Z due to the depletion of the envelope structure by strong mass loss such that semiconvective regions could not form. Z-dependent mass loss plays an indirect, yet decisive role in setting the HD limit as a function of Z. For a combination of efficient semiconvection and low overshooting with standard Z-dependent mass loss, we find a natural HD limit at all metallicities.

preprint2020arXiv

Properties of OB star-black hole systems derived from detailed binary evolution models

The recent gravitational wave measurements have demonstrated the existence of stellar mass black hole binaries. It is essential for our understanding of massive star evolution to identify the contribution of binary evolution to the formation of double black holes. A promising way to progress is investigating the progenitors of double black hole systems and comparing predictions with local massive star samples such as the population in 30 Doradus in the Large Magellanic Cloud (LMC). Methods. To this purpose, we analyse a large grid of detailed binary evolution models at LMC metallicity with initial primary masses between 10 and 40 Msun, and identify which model systems potentially evolve into a binary consisting of a black hole and a massive main sequence star. We then derive the observable properties of such systems, as well as peculiarities of the OB star component. We find that about 3% of the LMC late O and early B stars in binaries are expected to possess a black hole companion, when assuming stars with a final helium core mass above 6.6 M to form black holes. While the vast majority of them may be X-ray quiet, our models suggest that these may be identified in spectroscopic binaries, either by large amplitude radial velocity variations ( > 50 km s ) and simultaneous nitrogen surface enrichment, or through a moderate radial velocity ( > 10 km/s ) and simultaneously rapid rotation of the OB star. The predicted mass ratios are such that main sequence companions could be excluded in most cases. A comparison to the observed OB+WR binaries in the LMC, Be/X-ray binaries, and known massive BH binaries supports our conclusion. We expect spectroscopic observations to be able to test key assumptions in our models, with important implications for massive star evolution in general, and for the formation of double-black hole mergers in particular.

preprint2020arXiv

Star Formation in the Ultraviolet

With the launch of JWST and the upcoming installation of extremely large telescopes, the first galaxies in our Universe will finally be revealed. Their light will be dominated by massive stars, which peak in in the ultra-violet (UV) part of the electromagnetic spectrum. Star formation is the key driver of the evolution of our Universe. At young ages, within 10 Million years, both high and low mass stars generate complex UV emission processes which are poorly understood yet are vital for interpreting high red-shift line emission. For these reasons, the Hubble Space Telescope (HST) will devote 1000 orbits to obtaining a UV Legacy Library of Young Stars as Essential Standards (ULLYSES). The purpose of this Overview is to outline the basic physical principles driving UV emission processes from local (within 100 parsecs of) star formation, ranging from huge star-forming complexes containing hundreds of massive and very-massive stars (VMS), such as 30 Doradus (the Tarantula Nebula) in the neighbouring Magellanic Clouds (only 50 kpc away), to galaxies near and far, out to the epoch of Cosmic Reionization.

preprint2018arXiv

Constraining the progenitor evolution of GW 150914

One of the largest surprises from the LIGO results regarding the first gravitational wave detection (GW 150914) was the fact the black holes (BHs) were &#34;heavy&#34;, of order 30 - 40 solar masses. The most promising explanation for this obesity is that the BH-BH merger occurred at low metallicity (Z): when the iron (Fe) contents is lower this is expected to result in weaker mass loss during the Wolf-Rayet (WR) phase. We therefore critically evaluate the claims for the reasons of heavy BHs as a function of Z in the literature. Furthermore, weaker stellar winds might lead to more rapid stellar rotation, allowing WR and BH progenitor evolution in a chemically homogeneous manner. However, there is as yet no empirical evidence for more rapid rotation amongst WR stars in the low Z environment of the Magellanic Clouds. Due to the intrinsic challenge of determining WR rotation rates from emission lines, the most promising avenue to constrain rotation-rate distributions amongst various WR subgroups is through the utilisation of their emission lines in polarised light. We thus provide an overview of linear spectro-polarimetry observations of both single and binary WRs in the Galaxy, as well as the Large and Small Magellanic Clouds, at 50% and 20% of solar Z, respectively. Initial results suggest that the route of chemically homogeneous evolution (CHE) through stellar rotation is challenging, whilst the alternative of a post-LBV or common envelope evolution is more likely.

preprint2018arXiv

Massive star evolution revealed in the Mass-Luminosity plane

Massive star evolution is dominated by key physical processes such as mass loss, convection and rotation, yet these effects are poorly constrained, even on the main sequence. We utilise a detached, eclipsing binary HD166734 as a testbed for single star evolution to calibrate new MESA stellar evolution grids. We introduce a novel method of comparing theoretical models with observations in the &#39;Mass-Luminosity Plane&#39;, as an equivalent to the HRD (see Higgins & Vink, 2018). We reproduce stellar parameters and abundances of HD166734 with enhanced overshooting (alpha=0.5), mass loss and rotational mixing. When comparing the constraints of our testbed to the systematic grid of models we find that a higher value of alpha = 0.5 (rather than 0.1) results in a solution which is more likely to evolve to a neutron star than a black hole, due to a lower value of the compactness parameter.