Search

The James Webb Space Telescope (JWST) has revealed an unexpectedly high abundance of UV luminous galaxies at redshifts $z\gtrsim 10$, challenging `standard' galaxy formation models. This study investigates the role of rapidly rotating (massive) stars undergoing chemically homogeneous evolution (CHE) in reconciling this potential tension. These stars are more compact, hotter, and exhibit enhanced UV emission. We find that the rest-frame UV luminosity of star-forming galaxies can be significantly enhanced by a factor of $\sim 3-6$ when CHE stars above a minimum initial mass of $m_{\star,\min}^{\rm CHE}\sim 2-10\ \rm M_\odot$ account for more than half of the total stellar mass following a Salpeter initial mass function. As a result, the UV luminosity functions observed at $z\sim 12-16$ can be reproduced with less extreme values of star formation efficiency and UV luminosity stochastic variability. Our results highlight the potential of CHE in explaining the UV-bright galaxy populations detected by JWST and call for future work to explore the broader astrophysical implications of CHE and its associated phenomena in the early universe, such as gamma-ray bursts, compact object binaries, and metal enrichment., Comment: 7+3 pages, 6 figures, to be submitted to ApJL, comments are welcome

Context. Carbon, nitrogen, and oxygen are the most abundant elements throughout the universe, after hydrogen and helium. Studying these elements in low-metallicity stars can provide crucial information on the chemical composition in the early Galaxy and possible internal mixing processes that can alter the surface composition of the stars. Aims. This work aims to investigate the chemical abundance patterns for CNO elements and Li in a homogeneously analyzed sample of 52 metal-poor halo giant stars. Methods. We used high-resolution spectra with a high signal-to-noise ratio (S/N) to carry out a spectral synthesis to derive detailed C, N, O, and Li abundances for a sample of stars with metallicities in the range of -3.58 <= [Fe/H] <= -1.79 dex. Our study was based on the assumption of one-dimensional (1D) local thermodynamic equilibrium (LTE) atmospheres. Results. Based on carbon and nitrogen abundances, we investigated the deep mixing taking place within stars along the red giant branch (RGB). The individual abundances of carbon decrease towards the upper RGB while nitrogen shows an increasing trend, indicating that carbon has been converted into nitrogen. No signatures of ON-cycle processed material were found for the stars in our sample. We computed a set of galactic chemical evolution (GCE) models, implementing different sets of massive star yields, both with and without including the effects of stellar rotation on nucleosynthesis. We confirm that stellar rotation is necessary to explain the highest [N/Fe] and [N/O] ratios observed in unmixed halo stars. The predicted level of N enhancement varies sensibly in dependence of the specific set of yields that are adopted. For stars with stellar parameters similar to those of our sample, heavy elements such as Sr, Y, and Zr appear to have unchanged abundances despite the stellar evolution mixing processes., Comment: 19 pages, 17 figures. Accepted for publication in A&A on Semptember 04, 2024. Table with derived abundances will be available at the CDS after the paper will be published

During the fourth observing run, the LIGO-Virgo-KAGRA Collaboration reported the detection of a coalescing compact binary (GW230529$_{-}$181500) with component masses estimated at $2.5-4.5\, M_\odot$ and $1.2-2.0\, M_\odot$ with 90\% credibility. Given the current constraints on the maximum neutron star (NS) mass, this event is most likely a lower-mass-gap (LMG) black hole-neutron star (BHNS) binary. The spin magnitude of the BH, especially when aligned with the orbital angular momentum, is critical in determining whether the NS is tidally disrupted. An LMG BHNS merger with a rapidly spinning BH is an ideal candidate for producing electromagnetic counterparts. However, no such signals have been detected. In this study, we employ a detailed binary evolution model, incorporating new dynamical tide implementations, to explore the origin of BH spin in an LMG BHNS binary. If the NS forms first, the BH progenitor (He-rich star) must begin in orbit shorter than 0.35 days to spin up efficiently, potentially achieving a spin magnitude of $\chi_{\rm BH} > 0.3$. Alternatively, if a non-spinning BH (e.g., $M_{\rm BH} = 3.6\, M_\odot$) forms first, it can accrete up to $\sim 0.2\, M_\odot$ via Case BA mass transfer (MT), reaching a spin magnitude of $\chi_{\rm BH} \sim 0.18$ under Eddington-limited accretion. With a higher Eddington accretion limit (i.e., 10.0 $\Dot{M}_{\rm Edd}$), the BH can attain a significantly higher spin magnitude of $\chi_{\rm BH} \sim\,0.65$ by accreting approximately $1.0\, M_\odot$ during Case BA MT phase., Comment: Accepted by A&A Letters (09/11/2024)

On April 25th, 2019, the LIGO-Virgo Collaboration discovered a Gravitational-wave (GW) signal from a binary neutron star (BNS) merger, i.e., GW190425. Due to the inferred large total mass, the origin of GW190425 remains unclear. We perform detailed stellar structure and binary evolution calculations that take into account mass-loss, internal differential rotation, and tidal interactions between a He-rich star and a NS companion. We explore the parameter space of the initial binary properties, including initial NS and He-rich masses and initial orbital period. We find that the immediate post-common-envelope progenitor system, consisting of a primary $\sim2.0\,M_\odot$ ($\sim1.7\,M_\odot$) NS and a secondary He-rich star with an initial mass of $\sim3.0-5.5\,M_\odot$ ($\sim5.5-6.0\,M_\odot$) in a close binary with an initial period of $\sim0.08-0.5\,{\rm{days}}$ ($\sim 0.08-0.4\,{\rm{days}}$), that experiences stable Case BB/BC mass transfer (MT) during binary evolution, can reproduce the formation of GW190425-like BNS events. Our studies reveal that the secondary He-rich star of the GW190425's progenitor before its core collapse can be efficiently spun up through tidal interaction, finally remaining as a NS with rotational energy even reaching $\sim10^{52}\,{\rm{erg}}$, which is always much higher than the neutrino-driven energy of the supernova (SN) explosion. If the newborn secondary NS is a magnetar, we expect that GW190425 can be the remnant of a magnetar-driven SN, e.g., a magnetar-driven ultra-stripped SN, a superluminous SN, or a broad-line Type Ic SN. Our results show that GW190425 could be formed through the isolated binary evolution, which involves a stable Case BB/BC MT just after the common envelope phase. On top of that, we show the He-rich star can be tidally spun up, potentially forming a spinning magnetized NS (magnetar) during the second SN explosion., Comment: 7 pages, 8 figures, accepted for publication in A&A

Stars with initial mass above roughly 8 solar masses will evolve to form a core made of iron group elements at which point no further exothermic nuclear reactions between charged nuclei may prevent the core collapse. Electron captures, neutrino losses, and the photo-disintegration of heavy nuclei trigger the collapse of these stars. Models at the brink of core collapse are produced using stellar evolution codes and these pre-collapse models may be used in the study of the subsequent dynamical evolution (including their explosion as supernovae and the formation of compact remnants such as neutron stars or black holes). We upgrade the physical ingredients employed by the GENeva stellar Evolution Code, GENEC, so that it may cover the regime of high temperatures and high densities required to produce progenitors of core-collapse. We have improved GENEC in three directions, equation of state, the nuclear reaction network and the radiative and conductive opacities adapted for the computation of the advanced phases of evolution. We produce a small grid of pre-supernova models of stars with zero-age main sequence masses of 15, 20 and 25 solar masses at solar and less than half solar metallicities. The results are compared with analogous models produced with the MESA code. The global properties of our new models, particularly of their inner cores, are comparable to models computed with MESA and pre-existing progenitors in the literature. Between codes the exact shell structure varies impacting explosion predictions. Using GENEC with state-of-the-art physics, we have produced massive stellar progenitors prior to collapse. These progenitors are suitable for follow-up studies, including the dynamical collapse and supernova phases. Larger grids of supernova progenitors are now feasible, with potential for further dynamical evolution., Comment: 22 pages, 17 figures

Grids of stellar evolution models with rotation using the Geneva stellar evolution code (Genec) have been published for a wide range of metallicities. We introduce the last remaining grid of Genec models, with a metallicity of $Z=10^{-5}$. We study the impact of this extremely metal-poor initial composition on various aspects of stellar evolution, and compare it to the results from previous grids at other metallicities. We provide electronic tables that can be used to interpolate between stellar evolution tracks and for population synthesis. Using the same physics as in the previous papers of this series, we computed a grid of stellar evolution models with Genec spanning masses between 1.7 and 500 $M_\odot$, with and without rotation, at a metallicity of $Z=10^{-5}$. Due to the extremely low metallicity of the models, mass-loss processes are negligible for all except the most massive stars. For most properties (such as evolutionary tracks in the Hertzsprung-Russell diagram, lifetimes, and final fates), the present models fit neatly between those previously computed at surrounding metallicities. However, specific to this metallicity is the very large production of primary nitrogen in moderately rotating stars, which is linked to the interplay between the hydrogen- and helium-burning regions. The stars in the present grid are interesting candidates as sources of nitrogen-enrichment in the early Universe. Indeed, they may have formed very early on from material previously enriched by the massive short-lived Population III stars, and as such constitute a very important piece in the puzzle that is the history of the Universe., Comment: Accepted in A&A. 20 pages, 14 figures, 3 tables

The first stars formed over five orders of magnitude in mass by accretion in primordial dark matter halos. We study the evolution of massive, very massive and supermassive primordial (Pop III) stars over nine orders of magnitude in accretion rate. We use the stellar evolution code GENEC to evolve accreting Pop III stars from 10$^{-6}$ - 10$^3$ M$_\odot$/yr and study how these rates determine final masses. The stars are evolved until either the end of central Si burning or until they encounter the general relativistic instability (GRI). We also examine how metallicity affects the evolution of the stars. At rates below $2.5 x 10^{-5}$ M$_\odot$/yr the final mass of the star falls below that required for pair-instability supernovae. The minimum rate required to produce black holes with masses above 250 M$_\odot$ is $5 x 10^{-5}$ M$_\odot$/yr, well within the range of infall rates found in numerical simulations of halos that cool via H$_2$, $10^{-3}$ M$_\odot$/yr. At rates of $5 x 10^{-5}$ M$_\odot$/yr to $4 x 10^{-2}$ \Ms\ yr$^{-1}$, like those expected for halos cooling by both H$_2$ and Ly-alpha, the star collapses after Si burning. At higher accretion rates the GRI triggers the collapse of the star during central H burning. Stars that grow at above these rates are cool red hypergiants with effective temperatures $log(T_{\text{eff}}) = 3.8$ and luminosities that can reach 10$^{10.5}$ L$_\odot$. At accretion rates of 100 - 1000 M$_\odot$/yr the gas encounters the general relativistic instability prior to the onset of central hydrogen burning and collapses to a black hole with a mass of 10$^6$ M$_\odot$ without ever having become a star. We reveal for the first time the critical transition rate in accretion above which catastrophic baryon collapse, like that which can occur during galaxy collisions in the high-redshift Universe, produces supermassive black holes via dark collapse., Comment: 11 pages, 8 figures

With stunning clarity, JWST has revealed the Universe's first billion years. The scientific community is analyzing a wealth of JWST imaging and spectroscopic data from that era, and is in the process of rewriting the astronomy textbooks. Here, 1.5 years into the JWST science mission, we provide a snapshot of the great progress made towards understanding the initial chapters of our cosmic history. We highlight discoveries and breakthroughs, topics and issues that are not yet understood, and questions that will be addressed in the coming years, as JWST continues its revolutionary observations of the Early Universe. While this compendium is written by a small number of authors, invited to ISSI Bern in March 2024 as part of the 2024 ISSI Breakthrough Workshop, we acknowledge the work of a large community that is advancing our collective understanding of the evolution of the Early Universe., Comment: review article written by the attendees of the 2024 ISSI breakthrough workshop "The first billion year of the Universe", submitted. Comments welcome

The first stars might have been fast rotators. This would have important consequences for their radiative, mechanical and chemical feedback. We discuss the impact of fast initial rotation on the evolution of massive Population III models and on their nitrogen and oxygen stellar yields. We explore the evolution of Population III stars with initial masses in the range of 9Msol < Mini < 120Msol starting with an initial rotation on the Zero Age Main Sequence equal to 70% of the critical one. We find that with the physics of rotation considered here, our rapidly-rotating Population III stellar models do not follow a homogeneous evolution. They lose very little mass in case mechanical winds are switched on when the surface rotation becomes equal or larger than the critical velocity. Impact on the ionising flux appears modest when compared to moderately-rotating models. Fast rotation favours, in models with initial masses above ~20Msol, the appearance of a very extended intermediate convective zone around the H-burning shell during the core He-burning phase. This shell has important consequences on the sizes of the He- and CO-cores and thus impacts the final fate of stars. Moreover, it has a strong impact on nucleosynthesis boosting the production of primary 14N. Fast initial rotation impacts significantly the chemical feedback of Population III stars. Observations of extremely metal-poor stars and/or starbursting regions are essential to provide constraints on the properties of the first stars., Comment: Accepted in A&A. Pages 17. Figures 15. Tables 2

We present a method for a non-linear asteroseismic inversion suitable for gravity-mode pulsators and apply it to slowly pulsating B-type (SPB) stars. Our inversion method is based on the iterative improvement of a parameterised static stellar structure model, which in turn is based on constraints from the observed oscillation periods. We present tests to demonstrate that the method is successful in recovering the properties of artificial targets both inside and outside the parameter space. We also present a test of our method on the well-studied SPB star KIC 7760680. We believe that this method is promising for carrying out detailed analyses of observations of SPB and $\gamma$ Dor stars and will provide complementary information to evolutionary models., Comment: Accepted in A&A

We use Geneva-evolution-code to run evolutionary tracks for stellar masses ranging from $20$ to $85$ $M_\odot$ at SMC metallicity ($Z=0.002$). We upgrade the recipe for stellar winds by adopting our self-consistent m-CAK prescription, which reduces the value of mass-loss rate by a factor between 2 and 6 depending on the mass range. The impact of our new winds is wide, and it can be divided between direct and indirect impact. For the most massive models ($60$ and $85$ $M_\odot$) with $\dot M\gtrsim2\times10^{-7}$ $M_\odot$ yr$^{-1}$, the impact is direct because lower mass loss make stars remove less envelope and therefore remain more massive and less chemically enriched at their surface at the end of their MS phase. For the less massive models ($20$ and $25$ $M_\odot$) with $\dot M\lesssim2\times10^{-8}$ $M_\odot$ yr$^{-1}$, the impact is indirect because lower mass loss make the stars keep high rotational velocities for a longer period of time, then extending the H-core burning lifetime and reaching the end of the MS with higher surface enrichment. Given that the conditions at the H-depletion change, the stars will lose more mass during their He-core burning stages anyways. For $M_\text{zams}=20$ to $40$ $M_\odot$, our models predict stars will evolve through the Hertzsprung gap, from O-type supergiants to BSG and finally RSG, with larger mass fractions of helium compared to old evolution models. New models also set down to $M_\text{zams}=85\,M_\odot$ the minimal initial mass required for a single star to become WR at metallicity $Z=0.002$. New values for $\dot M$ need to be complemented with upgrades in additional features such as convective core overshooting and distribution of rotational velocities, besides more detailed observations from projects such as XShootU, in order to provide a robust framework for the study of massive stars at low metallicity environments., Comment: Accepted for publication in Astronomy & Astrophysics

The chemical enrichment of the early Universe is a crucial element in the formation and evolution of galaxies, and Population III (PopIII) stars must play a vital role in this process. In this study, we examine metal enrichment from massive stars in the early Universe's embryonic galaxies. Using radiation hydrodynamic simulations and stellar evolution modelling, we calculated the expected metal yield from these stars. Specifically, we applied accretion rates from a previous radiation-hydrodynamic simulation to inform our stellar evolution modelling, executed with the Geneva code, across 11 selected datasets, with final stellar masses between 500 and 9000 Msol. Our results demonstrate that the first generation of Pop III stars within a mass range of 2000 to 9000 Msol result in N/O, C/O and O/H ratios compatible with the values observed in very high-z galaxies GN-z11 and CEERS 1019. The ejecta of these Pop III stars are predominantly composed of He, H, and N. Our Pop III chemical enrichment model of the halo can accurately reproduce the observed N/O and C/O ratios, and, by incorporating a hundred times more zero-metallicity interstellar material with the stellar ejecta, it accurately attains the observed O/H ratio. Thus, a sub-population of extremely massive PopIII stars, with masses surpassing approximately 2000 Msol, effectively reproduces the CNO elemental abundances observed in high-z JWST galaxies to date., Comment: Accepted in A&A, 13 pages, 4 figures

Stars formed with initial mass over 50 Msun are very rare today, but they are thought to be more common in the early universe. The fates of those early, metal-poor, massive stars are highly uncertain. Most are expected to directly collapse to black holes, while some may explode as a result of rotationally powered engines or the pair-creation instability. We present the chemical abundances of J0931+0038, a nearby low-mass star identified in early followup of SDSS-V Milky Way Mapper, which preserves the signature of unusual nucleosynthesis from a massive star in the early universe. J0931+0038 has relatively high metallicity ([Fe/H] = -1.76 +/- 0.13) but an extreme odd-even abundance pattern, with some of the lowest known abundance ratios of [N/Fe], [Na/Fe], [K/Fe], [Sc/Fe], and [Ba/Fe]. The implication is that a majority of its metals originated in a single extremely metal-poor nucleosynthetic source. An extensive search through nucleosynthesis predictions finds a clear preference for progenitors with initial mass > 50 Msun, making J0931+0038 one of the first observational constraints on nucleosynthesis in this mass range. However the full abundance pattern is not matched by any models in the literature. J0931+0038 thus presents a challenge for the next generation of nucleosynthesis models and motivates study of high-mass progenitor stars impacted by convection, rotation, jets, and/or binary companions. Though rare, more examples of unusual early nucleosynthesis in metal-poor stars should be found in upcoming large spectroscopic surveys., Comment: 11 pages + 22 page appendix, accepted to ApJL

The physics of stellar rotation plays a crucial role in the evolution of stars, their final fate and the properties of compact remnants. Diverse approaches have been adopted to incorporate the effects of rotation in stellar evolution models. This study seeks to explore the consequences of these various prescriptions for rotation on essential outputs of massive star models. We compute a grid of 15 and 60 M$_{\odot}$ stellar evolution models with the Geneva Stellar Evolution Code (GENEC), accounting for both hydrodynamical and magnetic instabilities induced by rotation. In both the 15 and 60 M$_{\odot}$ models, the choice of the vertical and horizontal diffusion coefficients for the non magnetic models strongly impacts the evolution of the chemical structure, but has a weak impact on the angular momentum transport and the rotational velocity of the core. In the 15 M$_{\odot}$ models, the choice of diffusion coefficient impacts the convective core size during the core H-burning phase, whether the model begins core He-burning as a blue or red supergiant and the core mass at the end of He-burning. In the 60 M$_{\odot}$ models, the evolution is dominated by mass loss and is less affected by the choice of diffusion coefficient. In the magnetic models, magnetic instability dominates the angular momentum transport and such models are found to be less mixed when compared to their rotating non-magnetic counterparts. Stellar models with the same initial mass, chemical composition, and rotation may exhibit diverse characteristics depending on the physics applied. By conducting thorough comparisons with observational features, we can ascertain which method(s) produce the most accurate results in different cases., Comment: Accepted in A&A, 18 pages

Binary evolution codes are essential tools to help in understanding the evolution of binary systems. They contain a great deal of physics, for example stellar evolution, stellar interactions, mass transfer, tides, orbital evolution. Since many of these processes are difficult to account for in detail, we often rely on prescriptions obtained in earlier studies. We highlight that the impact of the dynamical tides with radiative damping has been implemented inconsistently with respect to its original theoretical formulation in many studies. We derive a new analytical solution for the evolution toward synchronization in the case of circular orbits and propose turnkey equations for the case of eccentric orbits that can be used in population synthesis studies. We compare the strength of the tidal torque obtained with this new formula with respect to that obtained with the formula generally used in literature by studying how the evolution toward synchronization of main sequence stellar models is affected. We conclude that by using an incorrect formula for the tidal torque, as has been done in many binary codes, the strength of the dynamical tides with radiative damping is over- or underestimated depending on whether the star is close to or far from synchronization.

*Context*. In the framework of the PLATO mission, to be launched in late 2026, seismic inversion techniques will play a key role in the mission precision requirements of the stellar mass, radius, and age. It is therefore relevant to discuss the challenges of the automation of seismic inversions, which were originally developed for individual modelling.\\ *Aims*. We tested the performance of our newly developed quality assessment procedure of seismic inversions, which was designed in the perspective of a pipeline implementation.\\ *Methods*. We applied our assessment procedure on a testing set composed of 26 reference models. We divided our testing set into two categories, calibrator targets whose inversion behaviour is well known from the literature and targets for which we assessed manually the quality of the inversion. We then compared the results of our assessment procedure with our expectations as a human modeller for three types of inversions, the mean density inversion, the acoustic radius inversion, and the central entropy inversion.\\ *Results*. We found that our quality assessment procedure performs as well as a human modeller. The mean density inversion and the acoustic radius inversion are suited for a large-scale application, but not the central entropy inversion, at least in its current form.\\ *Conclusions*. Our assessment procedure showed promising results for a pipeline implementation. It is based on by-products of the inversion and therefore requires few numerical resources to assess quickly the quality of an inversion result., Comment: Accepted for publication in Astronomy & Astrophysics

We study the differences between models computed with Ledoux and Schwarzschild criteria on the internal structure, evolutionary track in the Hertzsprung-Russell diagram (HRD), lifetimes, evolution of the surface abundances and velocities, and masses of the He and CO cores. We investigate the consequences on the nature of the supernova (SN) progenitors and the type of SN events, as well as on the yields of light elements. We also study the impact on the outputs of population synthesis models. Models with initial masses between 7 and 120 M$_\odot$ at solar metallicity ($Z$=0.014) and with an initial rotation equal to 0 or 0.4 times the critical velocity at the zero-age main sequence were computed with either the Schwarzschild or the Ledoux criterion until the end of the C-burning phase. Models with initial masses between 15 and 32 M$_\odot$ computed with the Schwarzschild criterion show larger intermediate convective zones attached to the H-burning shell than models computed with the Ledoux criterion. Their CO cores and outer convective zones in the red supergiant (RSG) phase are also smaller. This impacts many outputs of stars during the core He-burning phase. Schwarzschild models have smaller CO cores and outer convective zones in the RSG phase, and their blue-to-red supergiant ratio is much higher than for Ledoux models. They also produce longer crossings of the Hertzsprung gap and favour blue loops. The upper luminosity of RSGs is little affected by the change in the convective criterion. The maximum luminosity of RSG progenitors for type II-P SN events is lowered from 5.2 to 4.95 when the Ledoux criterion is used instead of the Schwarzschild criterion in non-rotating models. The Schwarzschild criterion predicts longer-lasting, less nitrogen-enriched, and faster-rotating Cepheids. Rotational mixing decreases the differences between Schwarzschild and Ledoux models., Comment: Accepted for publication in A&A (17 October 2023); 20 pages, 14 figures

Efforts to understand the origin and growth of massive black holes observed in the early Universe have spurred a strong interest in the evolution and fate of rapidly-accreting primordial (metal-free) stars. Here, we investigate the evolution of such Population III stars under variable accretion rates, focusing on the thermal response and stellar structure, the impact of the luminosity wave encountered early in the pre-main sequence phase, and the influence of accretion on their subsequent evolution. We employ the Geneva stellar evolution code and simulate ten models with varying accretion histories, covering a final mass range from 491 M$_{\odot}$ to 6127 M$_{\odot}$. Our findings indicate that the critical accretion rate delineating the red and blue supergiant regimes during the pre-main sequence evolution is approximately $2.5\times10^{-2} $M$_{\odot}$/yr. Once core hydrogen burning commences, the value of this critical accretion rate drops to $7.0\times10^{-3}$M$_{\odot}$/yr. Moreover, we also confirm that the Kelvin-Helmholtz timescale in the outer surface layers is the more relevant timescale for determining the transition between red and blue phases. Regarding the luminosity wave, we find that it affects only the early pre-main sequence phase of evolution and does not directly influence the transition between red and blue phases, which primarily depends on the accretion rate. Finally, we demonstrate that variable accretion rates significantly impact the lifetimes, surface enrichment, final mass and time spent in the red phase. Our study provides a comprehensive understanding of the intricate evolutionary patterns of Population III stars subjected to variable accretion rates., Comment: Accepted in A&A

Asteroseismic modelling will be part of the pipeline of the PLATO mission and will play a key role in the mission precision requirements on stellar mass, radius and age. It is therefore crucial to compare how current modelling strategies perform, and discuss the limitations and remaining challenges for PLATO, such as the so-called surface effects, the choice of physical ingredients, and stellar activity. In this context, we carried out a systematic study of the impact of surface effects on the estimation of stellar parameters. In this work, we demonstrated how combining a mean density inversion with a fit of frequencies separation ratios can efficiently damp the surface effects and achieve precise and accurate stellar parameters for ten Kepler LEGACY targets, well within the PLATO mission requirements. We applied and compared two modelling approaches, directly fitting the individual frequencies, or coupling a mean density inversion with a fit of the ratios, to six synthetic targets with a patched 3D atmosphere from Sonoi et al. (2015) and ten actual targets from the LEGACY sample. The fit of the individual frequencies is unsurprisingly very sensitive to surface effects and the stellar parameters tend to be biased, which constitutes a fundamental limit to both accuracy and precision. In contrast, coupling a mean density inversion and a fit of the ratios efficiently damps the surface effects, and allows us to get both precise and accurate stellar parameters. The average statistical precision of our selection of LEGACY targets with this second strategy is 1.9% for the mass, 0.7% for the radius, and 4.1% for the age, well within the PLATO requirements. Using the mean density in the constraints significantly improves the precision of the mass, radius and age determinations, on average by 20%, 33%, and 16%, respectively., Comment: Accepted for publication in Astronomy & Astrophysics

Betelgeuse is a well known bright red supergiant that shows semi-regular variations with four approximate periods of 2200, 420, 230, and 185 days. While the longest period was customarily regarded as LSP (long secondary period) of unknown origin, we identify the ~2200-d period as the radial fundamental mode, and the three shorter periods as the radial first, second, and third overtones. From a linear nonadiabatic pulsation analysis including the pulsation/convection coupling, we have found that these radial pulsation modes are all excited in the envelope of a model in a late stage of the core-carbon burning. Models with similar pulsation property have masses of 11~12 M_\odot (19M_\odot at ZAMS) with luminosities (log L/L_\odot =5.27~5.28) and effective temperatures (log T_{eff}\approx 3.53) that are consistent with the range of the observational determinations. We also find that a synthetic light curve obtained by adding the fundamental and the first-overtone mode is comparable with the light curve of Betelgeuse up to the Great Dimming. We conclude that Betelgeuse is likely in the late stage of core carbon burning, and a good candidate for the next Galactic Type II supernova., Comment: 12 pages, 6 figures, version accepted for publication in MNRAS

Mass loss due to line-driven winds is central to our understanding of the evolution of massive stars. We extend the evolution models introduced in Paper I, where the mass loss recipe is based on the simultaneous calculation of the wind hydrodynamics and the line-acceleration, by incorporating the effects of stellar rotation. We introduce a grid of self-consistent line-force parameters for a set of standard evolutionary tracks. With that, we generate a new set of evolutionary tracks with rotation for $M_\text{ZAMS}=25,40,70,$ and $120\,M_\odot$, and metallicities $Z=0.014$ and $0.006$. The self-consistent approach gives lower mass loss rates than the standard values adopted in previous evolution models. This decrease impacts strongly on the tracks of the most massive models. Weaker winds allow the star to retain more mass, but also more angular momentum. As a consequence, weaker wind models rotate faster and show a less efficient mixing in their inner stellar structure. The new tracks predict an evolution of the rotational velocities through the MS in close agreement with the range of $\varv\sin i$ values found by recent surveys of Galactic O-type stars. As subsequent implications, the weaker winds from self-consistent models suggest a reduction of the contribution of the isotope $^{26}$Al to the ISM due to stellar winds of massive stars during the MS phase. Moreover, the higher luminosities found for the self-consistent evolutionary models suggest that some populations of massive stars might be less massive than previously thought, as in the case of Ofpe stars at the Galactic Centre. Therefore, this study opens a wide range of consequences for further research based on the evolution of massive stars., Comment: Accepted for publication in Astronomy & Astrophysics journal

Stars strongly impact their environment, and shape structures on all scales throughout the universe, in a process known as ``feedback''. Due to the complexity of both stellar evolution and the physics of larger astrophysical structures, there remain many unanswered questions about how feedback operates, and what we can learn about stars by studying their imprint on the wider universe. In this white paper, we summarize discussions from the Lorentz Center meeting `Bringing Stellar Evolution and Feedback Together' in April 2022, and identify key areas where further dialogue can bring about radical changes in how we view the relationship between stars and the universe they live in., Comment: Accepted to the Publications of the Astronomical Society of the Pacific

Mass-loss in red supergiants (RSGs) is generally recognized to be episodic, but mass-loss prescriptions fail to reflect this. Evolutionary models show that the total amount of mass lost during this phase determines if these stars evolve to warmer temperatures before undergoing core collapse. The current Geneva evolutionary models mimic episodic mass loss by enhancing the quiescent prescription rates whenever the star's outer layers exceed the Eddington luminosity by a large factor. This results in a 20 solar-mass model undergoing significantly more mass loss during the RSG phase than it would have otherwise, but has little effect on models of lower masses. We can test the validity of this approach observationally by measuring the proportion of high-luminosity RSGs to that predicted by the models. To do this, we use our recent luminosity-limited census of RSGs in M31 and M33, making modest improvements to membership, and adopting extinctions based on the recent panchromatic M31 and M33 Hubble surveys. We then compare the proportions of the highest luminosity RSGs found to that predicted by published Geneva models, as well as to a special set of models computed without the enhanced rates. We find good agreement with the models which include the supra-Eddington enhanced mass loss. The models with lower mass-loss rates predict a larger fraction of high-luminosity RSGs than observed, and thus can be ruled out. We also use these improved data to confirm that the upper luminosity limit of RSGs is log L/Lo~5.4, regardless of metallicity, using our improved data on M31 and M33 plus previous results on the Magellanic Clouds., Comment: Accepted for publication in the Astrophysical Journal. Revised version includes an additional reference

The LIGO, Virgo, and KAGRA (LVK) collaboration has announced 90 coalescing binary black holes (BBHs) with $p_{\rm astro} > 50\%$ to date, however, the origin of their formation channels is still an open scientific question. Given various properties of BBHs (BH component masses and individual spins) inferred using the default priors by the LVK, independent groups have been trying to explain the formation of the BBHs with different formation channels. Of all formation scenarios, the chemically homogeneous evolution (CHE) channel has stood out with distinguishing features, namely, nearly-equal component masses and preferentially high individual spins aligned with the orbital angular momentum. We perform Bayesian inference on the BBH events officially reported in GWTC-3 with astrophysically-predicted priors representing different formation channels of the isolated binary evolution (CEE: common-envelope evolution channel; CHE; SMT: stable mass transfer). Given assumed models, we report strong evidence for GW190517\_055101 being most likely to have formed through the CHE channel. Assuming the BBH events in the subsample are all formed through one of the isolated binary evolution channels, we obtain the lower limits on the local merger rate density of these channels at $11.45 ~\mathrm{Gpc^{-3}~yr^{-1}}$ (CEE), $0.18 ~\mathrm{Gpc^{-3}~yr^{-1}}$ (CHE), and $0.63 ~\mathrm{Gpc^{-3}~yr^{-1}}$ (SMT) at $90\%$ credible level., Comment: 13 pages, 4 figures, 1 table

A vital condition for life on Earth is the steady supply of radiative heat by the Sun. Like all other stars, the Sun generates its emitted energy in its central regions where densities and temperatures are high enough for nuclear fusion processes to take place. Because stellar cores are usually covered by an opaque envelope, most of our knowledge about them and their life-giving nuclear processes comes from theoretical modelling or from indirect observations such as the detection of solar neutrinos and the study of stellar pulsations, respectively. Only in very rare cases, stars may expose their cores, e.g., when a tiny fraction of them evolves into Wolf-Rayet or helium hot subdwarf stars. However, for the vast majority of stars, namely unevolved stars that burn hydrogen to helium in their centres, direct observational clues on the cores are still missing. Based on a comprehensive spectroscopic and asteroseismic analysis, we show here that the bright B-type star $\gamma$ Columbae is the stripped pulsating core (with a mass of $4$-$5\,M_\odot$, where $M_\odot$ is the mass of the Sun) of a previously much more massive star of roughly $12\,M_\odot$ that just finished central hydrogen fusion. The star's inferred parameters indicate that it is still in a short-lived post-stripping structural readjustment phase, making it an extremely rare object. The discovery of this unique star paves the way to obtain invaluable insights into the physics of both single and binary stars with respect to nuclear astrophysics and common-envelope evolution. In particular, it provides first observational constraints on the structure and evolution of stripped envelope stars., Comment: This preprint has not undergone peer review or any post-submission improvements or corrections. The Version of Record of this article is published in Nature Astronomy, and is available online at https://doi.org/10.1038/s41550-022-01809-6

Type IIb supernovae are believed to originate from core-collapse progenitors having kept only a very thin hydrogen envelope. We aim to explore how some physical factors, such as rotation, metallicity, overshooting, and the initial orbital period in binaries, significantly affect the Roche lobe overflow and the formation of type IIb supernovae. It is found that binaries are the main channel that capable of producing type typeIIb supernovae progenitors in the mass range for initial masses below 20 $M_{\odot}$. The formation of type IIb supernova progenitors is extremely sensitive to the initial orbital period. A less massive hydrogen indicates smaller radius and a higher effective temperatures, and vice versa. Binary systems with initial periods between 300 and 720 days produce type IIb progenitors that are a red supergiant. Those with an initial period between 50 and 300 days produce yellow supergiant progenitors and those with initial periods shorter than 50 days, blue supergiant progenitors. Both rapid rotation and larger overshooting can enlarge the carbon-oxygen core mass and lead to higher core temperature and lower central density at the pre-collapse phase. They are also beneficial to surface nitrogen enrichment but restrict the efficiency of the first dredge-up. SN IIb progenitors with low metallicity have smaller hydrogen envelope masses and radii than the high metallicity counterparts. Ultra-stripped binary models have systematically higher core mass fraction $\rm ^{12}C$ left, which has important influence on the compactness of type IIb progenitors., Comment: Accepted in ApJS

The center of the Milky Way contains stellar populations spanning a range in age and metallicity, with a recent star formation burst producing young and massive stars. Chemical abundances in the most luminous stellar member of the Nuclear Star Cluster (NSC), IRS 7, are presented for $^{19}$F, $^{12}$C, $^{13}$C, $^{14}$N, $^{16}$O, $^{17}$O, and Fe from an LTE analysis based on spherical modeling and radiative transfer with a 25M$_{\odot}$ model atmosphere, whose chemistry was tailored to the derived photospheric abundances. We find IRS 7 to be depleted heavily in both $^{12}$C (~-0.8 dex) and $^{16}$O (~-0.4 dex), while exhibiting an extremely enhanced $^{14}$N abundance (~+1.1 dex), which are isotopic signatures of the deep mixing of CNO-cycled material to the stellar surface. The $^{19}$F abundance is also heavily depleted by ~1 dex relative to the baseline fluorine of the Nuclear Star Cluster, providing evidence that fluorine along with carbon constrain the nature of the deep mixing in this very luminous supergiant. The abundances of the minor isotopes $^{13}$C and $^{17}$O are also derived, with ratios of $^{12}$C/$^{13}$C~5.3 and $^{16}$O/$^{17}$O~525. The derived abundances for IRS 7, in conjunction with previous abundance results for massive stars in the NSC, are compared with rotating and non-rotating models of massive stars and it is found that the IRS 7 abundances overall follow the behavior predicted by stellar models. The depleted fluorine abundance in IRS 7 illustrates, for the first time, the potential of using the $^{19}$F abundance as a mixing probe in luminous red giants., Comment: 11 pages, 6 figures, 2 tables

Context: Population III (Pop III) stars may be fast rotating. An expected consequence of fast rotation is strong internal mixing that deeply affects their evolutionary tracks in the Hertzsprung-Russell diagram and hence their ionising power. aims: We investigate the impact on the ionising power of Pop III stars in an extreme case of internal mixing, the one leading to chemically homogeneous evolution (CHE). In that situation, during the main sequence phase, the star keeps the same chemical composition from its center to its surface. Homogeneous stars have larger effective temperatures and luminosities than stars evolving non-homogeneously and thus are stronger ionising sources. Methods: The stellar evolution models are based on $n=3$ polytropes with a time varying hydrogen mass fraction. The ionisation model employs the self-similar champagne flow solution from Shu et al. (2002), as well as numerical simulations for the stochastic treatment of star clusters over a grid of redshifts and halo masses. Results: We find that haloes containing chemically homogeneous stars have an escape fraction of ionising photons about twice that of haloes containing classical Pop III stars. By extrapolating the high-$z$ ionisation history powered by Pop III stars (at $z\gtrsim 15$) to the post-reionisation epoch, we derive the Thomson scattering optical depth $\tau$, which is compared with the value measured by $\textit{Planck}$. We find that $\tau$ is overproduced by $\sim1.5- 5\sigma$, when all Pop III stars evolve homogeneously. This indicates that CHE is unlikely to be realised in the majority of Pop III stars, although the present study cannot exclude that a fraction of them undergoes CHE. Conclusions: Fast rotation might have a significant impact on the ionising budget of Pop III stars, and thus on early cosmic reionisation., Comment: Accepted for publication in A&A (12 August 2022); 13 pages, 9 figures

Context: Grids of stellar models, computed with the same physical ingredients, allow one to study the impact of a given physics on a broad range of initial conditions and are a key ingredient for modeling the evolution of galaxies. Aims: We present a grid of single star models for masses between 0.8 and 120 $M_\odot$, with and without rotation for a mass fraction of heavy element Z=0.006, representative of the Large Magellanic Cloud (LMC). Methods: We used the Geneva stellar evolution code. The evolution was computed until the end of the central carbon-burning phase, the early asymptotic giant branch phase, or the core helium-flash for massive, intermediate, and low mass stars, respectively. Results: The outputs of the present stellar models are well framed by the outputs of the two grids obtained by our group for metallicities above and below the one considered here. The models of the present work provide a good fit to the nitrogen surface enrichments observed during the main sequence for stars in the LMC with initial masses around 15 $M_\odot$. They also reproduce the slope of the luminosity function of red supergiants of the LMC well, which is a feature that is sensitive to the time-averaged mass loss rate over the red supergiant phase. The most massive black hole that can be formed from the present models at Z=0.006 is around 55 $M_\odot$. No model in the range of mass considered will enter into the pair-instability supernova regime, while the minimal mass to enter the region of pair pulsation instability is around 60 $M_\odot$ for the rotating models and 85 $M_\odot$ for the nonrotating ones. Conclusions: The present models are of particular interest for comparisons with observations in the LMC and also in the outer regions of the Milky Way. We provide public access to numerical tables that can be used for computing interpolated tracks and for population synthesis studies., Comment: 11 pages, 11 figures, published in A&A

We present a grid of stellar models at super-solar metallicity (Z = 0.020) extending the previous grids of Geneva models at solar and sub-solar metallicities. A metallicity of Z = 0.020 was chosen to match that of the inner Galactic disk. A modest increase of 43% (=0.02/0.014) in metallicity compared to solar models means that the models evolve similarly to solar models but with slightly larger mass loss. Mass loss limits the final total masses of the super-solar models to 35 M$_\odot$ even for stars with initial masses much larger than 100 M$_\odot$. Mass loss is strong enough in stars above 20 M$_\odot$ for rotating stars (25 M$_\odot$ for non-rotating stars) to remove the entire hydrogen-rich envelope. Our models thus predict SNII below 20 M$_\odot$ for rotating stars (25 M$_\odot$ for non-rotating stars) and SNIb (possibly SNIc) above that. We computed both isochrones and synthetic clusters to compare our super-solar models to the Westerlund 1 (Wd1) massive young cluster. A synthetic cluster combining rotating and non-rotating models with an age spread between log10 (age/yr) = 6.7 and 7.0 is able to reproduce qualitatively the observed populations of WR, RSG and YSG stars in Wd1, in particular their simultaneous presence at log10(L/L$_\odot$) = 5-5.5. The quantitative agreement is imperfect and we discuss the likely causes: synthetic cluster parameters, binary interactions, mass loss and their related uncertainties. In particular, mass loss in the cool part of the HRD plays a key role., Comment: 15 pages, 12 figures. Accepted for publication in MNRAS

The $\beta$ Cephei pulsators are massive main-sequence stars, presenting low radial-order modes. These modes probe in particular the chemical gradient at the edge of the convective core. They hence give constraints on macroscopic processes, such as hydrodynamic or magnetic instabilities. Yet, it is not clear to what extent the seismic inferences depend on the physics employed for the stellar modelling or on the observational dataset. We investigate the observational constraints which are necessary to provide accurate constraints on the mixing processes in $\beta$ Cephei stars. We explore the importance of the identification of the angular degree of modes. Depending on the quality of the seismic dataset and the classical constraints, we estimate the precision achievable with asteroseismology. We propose a method extending the forward approach classically used to model $\beta$ Cephei stars The probability distributions of the asteroseismic-derived stellar parameters are obtained. With these distributions, we provide a systemic way to estimate the errors of the modelling. A particular effort is made to also include the theoretical uncertainties of the models. We then estimate the accuracy and precision of asteroseismology for $\beta$ Cephei stars in a series of hare and hound exercises. The exercises show that a set of four to five oscillation frequencies with an identified angular degree already leads to accurate inferences on the stellar parameters. Without the identification of the modes, the addition of other classical observational constraints allow to succeed the seismic modelling. When the micro-physics of the star and stellar models used for the modelling differ, the constraints derived on the internal structure remain valid if expressed in terms of acoustic variables. The characterisation of the mixing processes at the boundary of the convective core remain model-dependent., Comment: 17 pages, 14 figures, accepted for publication in A&A

The observed late-type WC Wolf-Rayet stars (WC7-9) with low luminosity below $\rm \log L/L_{\odot} < 5.4$ in the HR diagram cannot be reproduced satisfactorily by the evolutionary track of single stars. The mass transfer due to Roche lobe overflow drastically modifies the internal structure and surface compositions of two components. Therefore, binaries provide a very promising evolutionary channel to produce these WC stars., Comment: Accept for publication in A&A

In the search of a sample of metal-poor bright giants using Str{\"o}mgren photometry, we serendipitously found a sample of 26 young (ages younger than 1 Gyr) metal-rich giants, some of which have high rotational velocities.We determined the chemical composition and rotational velocities of these stars in order to compare them with predictions from stellar evolution models. These stars where of spectral type A to B when on the main sequence, and we therefore wished to compare their abundance pattern to that of main-sequence A and B stars.Stellar masses were derived by comparison of the position of the stars in the colour-magnitude diagram with theoretical evolutionary tracks. These masses, together with Gaia photometry and parallaxes, were used to derive the stellar parameters. We used spectrum synthesis and model atmospheres to determine chemical abundances for 16 elements (C, N, O, Mg, Al, Ca, Fe, Sr, Y, Ba, La, Ce, Pr, Nd, Sm, and Eu) and rotational velocities.The age-metallicity degeneracy can affect photometric metallicity calibrations. We identify 15 stars as likely binary stars. All stars are in prograde motion around the Galactic centre and belong to the thin-disc population. All but one of the sample stars present low [C/Fe] and high [N/Fe] ratios together with constant [(C+N+O)/Fe], suggesting that they have undergone CNO processing and first dredge-up. The observed rotational velocities are in line with theoretical predictions of the evolution of rotating stars., Comment: Astronomy and Astrophysics - A\&A, EDP Sciences, 2021

The evolution of massive stars is affected by a variety of physical processes including convection, rotation, mass loss and binary interaction. Because these processes modify the internal chemical abundance profiles in multiple ways simultaneously, it can be challenging to determine which properties of the stellar interior are primarily driving the overall evolution. Building on previous work, we develop a new modelling approach called SNAPSHOT that allows us to isolate the key features of the internal abundance profile that drive the evolution of massive stars. Using our approach, we compute numerical stellar structure models in thermal equilibrium covering key phases of stellar evolution. For the main sequence, we demonstrate that models with the same mass and very similar surface properties can have different internal distributions of hydrogen and convective core masses. We discuss why massive stars expand after the main sequence and the fundamental reasons for why they become red, blue or yellow supergiants. For the post-main sequence, we demonstrate that small changes in the abundance profile can cause very large effects on the surface properties. We also discuss the effects that produce blue supergiants and the cause of blue loops. Our models show that massive stars with lower metallicity tend to be more compact due to the combined effect of lower CNO abundances in the burning regions and lower opacity in the envelope., Comment: Accepted in MNRAS

The LIGO Scientific Collaboration and Virgo Collaboration (LVC) have recently reported in GWTC-2.1 eight additional candidate events with a probability of astrophysical origin greater than 0.5 in the LVC deeper search on O3a running. In GWTC-2.1, the majority of the effective inspiral spins ($\chi_{\rm eff}$) show magnitudes consistent with zero, while two (GW190403$_{-}$051519 and GW190805$_{-}$211137) of the eight new events have $\chi_{\rm eff}$ $> 0$ (at 90% credibility). We note that GW190403$_{-}$051519 was reported with $\chi_{\rm eff}$ = $0.70^{+0.15}_{-0.27}$ and mass ratio $q$ = $0.25^{+0.54}_{-0.11}$, respectively. Assuming a uniform prior probability between 0 and 1 for each black hole's dimensionless spin magnitude, GW190403$_{-}$051519 was reported with the dimensionless spin of the more massive black hole, $\chi_1$ = $0.92^{+0.07}_{-0.22}$. This is the fastest first-born black hole ever measured in all current gravitational-wave events. If GW190403$_{-}$051519 is formed through isolated binary evolution channel, this extremely high spin challenges, at least in that case, the existence of efficient angular momentum transport mechanism between the stellar core and the radiative envelope of massive stars, as for instance predicted by the Tayler-Spruit dynamo (Spruit 2002) or its revised version by Fuller et al. 2019., Comment: Accepted for publication in ApJ on 11th November 2021

The first stars are thought to be one of the dominant sources of hydrogen reionization in the early Universe, with their high luminosities and surface temperatures expected to drive high ionizing photon production rates. In this work, we take our Geneva stellar evolution models of zero-metallicity stars and predict their production rates of photons capable to ionize H, He I and He II, based on a blackbody approximation. We present analytical fits in the range 1.7-500 solar masses. We then explore the impact of stellar initial mass, rotation, and convective overshooting for individual stars. We have found that ionizing photon production rates increase with increasing initial mass. For the rotational velocities considered we see changes of up to 25% to ionizing photons produced. This varies with initial mass and ionizing photon species and reflects changes to surface properties due to rotation. We have also found that higher convective overshooting increases ionizing photon production by approximately 20% for the change in overshooting considered here. For stellar populations, we explore how the production of ionizing photons varies as a function of the initial mass function (IMF) slope, and minimum and maximum initial masses. For a fixed population mass we have found changes of the order of 20-30% through varying the nature of the IMF. This work presents ionizing photon production predictions for the most up to date Geneva stellar evolution models of Population III stars, and provides insight into how key evolutionary parameters impact the contribution of the first stars to reionization., Comment: 20 pages, 18 figures. Accepted for publication in MNRAS

Context: The Sun is by far a privileged target for testing stellar models with unique precision. A recent concern appeared with the progress in the solar surface abundances derivation that has led to a decrease of the solar metallicity. While the ancient high-metallicity models were in fair agreement with other solar observational indicators, it is no longer the case for low-metallicity models. Recent collection of data are however promising to shed a new light on it. For instance, the Borexino collaboration released in 2020 the first-ever complete estimate of neutrinos emitted in the CNO cycle. It has reaffirmed the role of the neutrino constraints in the solar modelling process and its associated issues. In parallel, newly claimed detection of solar gravity modes of oscillations offers another opportunity of probing the stratification in the Sun's central layers. Aims: We propose to combine the diagnoses from neutrinos and helioseismology, both from pressure and gravity modes, for assessing the predictions of solar models. We compare in detail the different physical prescriptions currently at disposal for stellar model computations. Results: The CNO neutrino flux confirms a preference for high-metallicity models. Nevertheless, we found that mild modification of the nuclear screening factors can re-match low-metallicity model predictions to observed fluxes, although it does not restore the agreement with the helioseismic frequency ratios. Neither the high-metallicity or low-metallicity models are able to reproduce the gravity-mode period spacing. The disagreement is huge, more than 100$\sigma$ to the reported value. Reversely, the family of standard models narrows the expected range of the Sun's period spacing: between $\sim$2150 to $\sim$2190~s. Moreover, we show this indicator can constrain the chemical mixture, opacity, and to a lower extent nuclear reactions in solar models., Comment: 18 pages, 10 figures, accepted for publications in Astronomy and Astrophysics

Stellar winds from fast-rotating Population III (Pop III) stars have long been suspected to make important contributions to early metal enrichment, as features in the nucleosynthesis of such `spinstars' are consistent with the chemical abundance patterns of some metal-poor stars in the local Universe. Particularly, stellar winds rich in light elements can provide another pathway towards explaining the carbon enhancement in carbon-enhanced metal-poor (CEMP) stars. In this work, we focus on the feedback of Pop III stellar winds combined with supernovae (SNe), and derive the resulting chemical signatures in the enriched medium. We explore a large parameter space of Pop III star formation, feedback, yields from winds and SNe with a semi-analytical model. The predicted pattern of carbon and iron abundances of second-generation stars agrees well with observations of CEMP-no stars ($[\rm Ba/Fe]<0$) at $[\rm Fe/H]\lesssim -3$ and $A(\mathrm{C})\lesssim 7$, under the optimistic assumption of significant mass loss by winds from massive ($\gtrsim 25\ \rm M_{\odot}$) stars that collapse into BHs without SNe. In this scenario, carbon-rich but iron-free second-generation stars can form in systems dominated by enrichment from winds, gaining trace amounts of iron by accretion from the interstellar medium, to become the most iron-poor and carbon-enhanced stars seen in observations ($[\rm Fe/H]\lesssim -4$, $[\rm C/Fe]\gtrsim 2$). We conclude that the observed CEMP-no stars can be explained by both our winds + ISM accretion channel as well as the well-studied faint SN scenario. Wind feedback from Pop~III spinstars deserves more detailed modelling in early cosmic structure formation., Comment: 23 pages, 19 figures, accepted for publication in MNRAS

Tidal interactions and planet evaporation processes impact the evolution of close-in star-planet systems. We study the impact of stellar rotation on these processes. We compute the time evolution of star-planet systems consisting of a planet with initial mass between 0.02 and 2.5 M$_{ Jup}$ (6 and 800 M$_{ Earth}$), in a quasi-circular orbit with an initial orbital distance between 0.01 and 0.10 au, around a solar-type star evolving from the Pre-Main-Sequence (PMS) until the end of the Main-Sequence (MS) phase. We account for the evolution of: the stellar structure, the stellar angular momentum due to tides and magnetic braking, the tidal interactions (equilibrium and dynamical tides in stellar convective zones), the mass-evaporation of the planet, and the secular evolution of the planetary orbit. We consider that at the beginning of the evolution, the proto-planetary disk has fully dissipated and planet formation is complete. Both a rapid initial stellar rotation, and a more efficient angular momentum transport inside the star, in general, contribute toward the enlargement of the domain which is devoid of planets after the PMS phase, in the plane of planet mass vs. orbital distance. Comparisons with the observed distribution of exoplanets orbiting solar mass stars, in the plane of planet mass vs. orbital distance (addressing the "Neptunian desert" feature), show an encouraging agreement with the present simulations, especially since no attempts have been made to fine-tune initial parameters of the models to fit the observations. We also obtain an upper limit for the orbital period of bare-core planets, that agrees with observations of the "radius valley" feature in the plane of planet radius vs. orbital period. The two effects, tides and planet evaporation, should be accounted for simultaneously and in a consistent way, with a detailed model for the evolution of the star., Comment: 15 pages, 13 figures. This manuscript has been accepted for publication in Astronomy & Astrophysics section 10. Planets and planetary systems

Massive star evolution is still poorly understood, and observational tests are required to discriminate between different implementations of physical phenomenon in stellar evolution codes. By confronting stellar evolution models with observed properties of blue supergiants, such as pulsations, chemical composition and position in the Hertzsprung-Russell diagram, we aim at determining which of the criterion used for convection (Schwarzschild or Ledoux) is best able to explain the observations. We compute state-of-the-art stellar evolution models with either the Schwarzschild or the Ledoux criterion for convection. Models are for $14$ to $35\,M_\odot$ at solar or Large Magellanic Cloud metallicity. For each model, we compute the pulsation properties to know when radial modes are excited. We then compare our results with the position of blue supergiants in the Hertzsprung-Russell diagram, with their surface chemical composition, and with their variability. Our results at Large Magellanic Cloud metallicity shows only a slight preference for the Ledoux criterion over the Schwarzschild one in reproducing at the same time the observed properties of blue supergiants, even if the Schwarzschild criterion cannot be excluded at this metallicity. We check that changing the overshoot parameter at solar metallicity does not improve the situation. We also check that our models are able to reproduce the position of Galactic blue supergiant in the flux-weighted-gravity -- luminosity relation. We confirm that overall, models computed with the Ledoux criterion are slightly better in matching observations. Our results also support the idea that most Galactic $\alpha$ Cyg variables are blue supergiants of the group 2, i.e. stars that have been through a previous red supergiant phase where they have lost large amount of mass., Comment: 14 pages, 10 figures, accepted for publication in A&A

Understanding the nature of the first stars is key to understanding the early universe. With new facilities such as JWST we may soon have the first observations of the earliest stellar populations, but to understand these observations we require detailed theoretical models. Here we compute a grid of stellar evolution models using the Geneva code with the aim to improve our understanding of the evolution of zero-metallicity stars, with particular interest in how rotation affects surface properties, interior structure, and metal enrichment. We produce a range of models of initial masses (Mini) from 1.7 Msun to 120 Msun, focusing on massive models of 9 Msun < Mini < 120 Msun. Our grid includes models with and without rotation, with rotating models having an initial velocity of 40% of the critical velocity. We find that rotation strongly impacts the evolution of the first stars, mainly through increased core size and stronger H-burning shells during core He-burning. Without radiative mass loss, angular momentum builds at the surface in rotating models, thus models of initial masses Mini > 60 Msun reach critical rotation on the main sequence and experience mass loss. We find that rotational mixing strongly affects metal enrichment, but does not always increase metal production as we see at higher metallicities. This is because rotation leads to an earlier CNO boost to the H shell during He-burning, which may hinder metal enrichment depending on initial mass and rotational velocity. Electronic tables of this new grid of Population III models are publicly available., Comment: 20 pages, 15 figures. Accepted for publication in MNRAS

We use N-body simulations to study the dynamical evolution of Population III (Pop III) stellar systems and the resulting binary statistics. We design a physically-motivated framework for the initial conditions of Pop III star clusters, based on small-scale hydrodynamic simulations and the scale-free nature of disk evolution during Pop III star formation. Our novel approach enables us to explore the dependence of binary statistics on initial conditions and arrive at more robust predictions for the signals of Pop III X-ray binaries (XRBs) and binary black hole (BBH) mergers, compared to simple extrapolations of Pop III protostar systems. We find that binary properties are highly sensitive to the initial cluster size and distribution of binary separation, while the effect of initial mass function is relatively minor. Our simulations predict less close binaries, and thus, significantly lower efficiencies (by a factor of $\sim 10-10^{4}$) for the formation and accretion of Pop III XRBs, than found in previous studies, implying that the contribution of Pop III XRBs to the cosmic X-ray background is negligible and their feedback effects are unimportant. We estimate the efficiency of Pop III BBH mergers as $\sim 10^{-5}-10^{-4}\ \rm M_{\odot}^{-1}$, for which 3-body hardening by surrounding stars in dense star clusters or close binary interactions is required to facilitate in-spirals of BBHs. All simulation data, including catalogs of Pop III binaries and multiple systems, are publicly available., Comment: 19 + 4 pages, 15 + 3 figures, accepted for publication in MNRAS, data available at https://drive.google.com/drive/folders/1JHBjhSBDPT3jdipdQtmr2IJ5TV6Xe6Id?usp=sharing

GW190521 challenges our understanding of the late-stage evolution of massive stars and the effects of the pair-instability in particular. We discuss the possibility that stars at low or zero metallicity could retain most of their hydrogen envelope until the pre-supernova stage, avoid the pulsational pair-instability regime and produce a black hole with a mass in the mass gap by fallback. We present a series of new stellar evolution models at zero and low metallicity computed with the Geneva and MESA stellar evolution codes and compare to existing grids of models. Models with a metallicity in the range 0-0.0004 have three properties which favour higher BH masses as compared to higher metallicity models. These are (i) lower mass-loss rates during the post-MS phase, (ii) a more compact star disfavouring binary interaction and (iii) possible H-He shell interactions which lower the CO core mass. We conclude that it is possible that GW190521 may be the merger of black holes produced directly by massive stars from the first stellar generations. Our models indicate BH masses up to 70-75 Msun. Uncertainties related to convective mixing, mass loss, H-He shell interactions and pair-instability pulsations may increase this limit to ~85 Msun., Comment: 6 pages, 3 figures, accepted for publication in MNRAS Letters

We introduce SNAPSHOT, a technique to systematically compute stellar structure models in hydrostatic and thermal equilibrium based on 3 structural properties - core mass $M_{\rm core}$, envelope mass $M_{\rm env}$ and core composition. This approach allows us to connect these properties of stellar interiors to the luminosity and effective temperature $T_{\rm eff}$ in a more systematic way than with stellar evolution models. For MS models, we derive an analytical relationship between $M_{\rm core}$, $M_{\rm total}$ and central H abundance that can be used in rapid stellar evolution algorithms. Core-He burning models with $M_{\rm core}/M_{\rm total}$ from 0.2 to 0.8 have convective envelopes, low $T_{\rm eff}$ and will appear as red supergiants. For a given $M_{\rm core}$, they exhibit a small variation in luminosity (0.02 dex) and $T_{\rm eff}$ ($\sim 400\,\mathrm{K}$) over a wide range of $M_{\rm env}$ ($\sim 2 - 20\,\mathrm{M}_{\odot}$). This means that it is not possible to derive red supergiant masses from luminosities and $T_{\rm eff}$ alone. We derive the following relationship between $M_{\rm core}$ and the total luminosity of a red supergiant during core He burning: $\log M_{\rm core} \simeq 0.44 \log L/L_{\odot} - 1.38$. At $M_{\rm core}$/$M_{\rm total} \approx 0.2$, our models exhibit a bi-stability and jump from a RSG to a BSG structure. Our models with $M_{\rm core}/M_{\rm total} > 0.8$, which correspond to stripped stars produced by mass loss or binary interaction, show that $T_{\rm eff}$ has a strong dependence on $M_{\rm env}$, $M_{\rm core}$ and the core composition. We find the mass of one of these stripped stars in a binary system, HD 45166, to be less than its dynamical mass. When a large observational sample of stripped stars becomes available, our results can be used to constrain their masses and the physics of binary interaction., Comment: Accepted for publication in MNRAS

We show that it is not possible to determine the final mass $M_{\rm fin}$ of a red supergiant (RSG) at the pre-supernova (SN) stage from its luminosity $L$ and effective temperature $T_{\rm eff}$ alone. Using a grid of stellar models, we demonstrate that for a given value of $L$ and $T_{\rm eff}$, a RSG can have a range of $M_{\rm fin}$ as wide as 3 to $45~\mathrm{M}_{\odot}$. While the probability distribution within these limits is not flat, any individual determination of $M_{\rm fin}$ for a RSG will be degenerate. This makes it difficult to determine its evolutionary history and to map $M_{\rm fin}$ to an initial mass. Single stars produce a narrower range that is difficult to accurately determine without making strong assumptions about mass loss, convection, and rotation. Binaries would produce a wider range of RSG $M_{\rm fin}$. However, the final Helium core mass M$_{\rm He-core}$ is well determined by the final luminosity and we find $\log (\mathrm{M}_{\rm He-core}/M_{\odot}) = 0.659 \log (L/\mathrm{L}_{\odot}) -2.630$ Using this relationship, we derive M$_{\rm He-core}$ for directly imaged SN progenitors and one failed SN candidate. The value of $M_{\rm fin}$ for stripped star progenitors of SNe IIb is better constrained by $L$ and $T_{\rm eff}$ due to the dependence of $T_{\rm eff}$ on the envelope mass $M_{\rm env}$ for $M_{\rm env} \lesssim 1~$M$_{\odot}$. Given the initial mass function, our results apply to the majority of progenitors of core collapse SNe, failed SNe and direct collapse black holes., Comment: Re-submitted to MNRAS Letters after referee comments

We investigate the effects of mass loss during the main-sequence (MS) and post-MS phases of massive star evolution on black hole (BH) birth masses. We compute solar metallicity Geneva stellar evolution models of an 85 $M_{\odot}$ star with mass-loss rate ($\dot{M}$) prescriptions for MS and post-MS phases and analyze under which conditions such models could lead to very massive BHs. Based on the observational constraints for $\dot{M}$ of luminous stars, we discuss two possible scenarios that could produce massive BHs at high metallicity. First, if a massive BH progenitor evolves from the observed population of massive MS stars known as WNh stars, we show that its average post-MS mass-loss rate has to be less than $1\,\times10^{-5}\,M_{\odot}$/yr. However, this is lower than the typical observed mass-loss rates of luminous blue variables (LBV). Second, a massive BH progenitor could evolve from a yet undetected population of $80-85$ $M_{\odot}$ stars with strong surface magnetic fields, which could quench mass loss during the evolution. In this case, the average mass-loss rate during the post-MS LBV phase has to be less than $5\,\times10^{-5}\,M_{\odot}$/yr to produce 70 $M_{\odot}$ BHs. We suggest that LBVs that explode as SNe have large envelopes and small cores that could be prone to explosion, possibly evolving from binary interaction (either mergers or mass gainers that do not fully mix). Conversely, LBVs that directly collapse to BHs could have evolve from massive single stars or binary-star mergers that fully mix, possessing large cores that would favor BH formation., Comment: Accepted in ApJ

The mass-loss rates of red supergiant stars (RSGs) are poorly constrained by direct measurements, and yet the subsequent evolution of these stars depends critically on how much mass is lost during the RSG phase. In 2012 the Geneva evolutionary group updated their mass-loss prescription for RSGs with the result that a 20 solar mass star now loses 10x more mass during the RSG phase than in the older models. Thus, higher mass RSGs evolve back through a second yellow supergiant phase rather than exploding as Type II-P supernovae, in accord with recent observations (the so-called "RSG Problem"). Still, even much larger mass-loss rates during the RSG phase cannot be ruled out by direct measurements of their current dust-production rates, as such mass-loss is episodic. Here we test the models by deriving a luminosity function for RSGs in the nearby spiral galaxy M31 which is sensitive to the total mass loss during the RSG phase. We carefully separate RSGs from asymptotic giant branch stars in the color-magnitude diagram following the recent method exploited by Yang and collaborators in their Small Magellanic Cloud studies. Comparing our resulting luminosity function to that predicted by the evolutionary models shows that the new prescription for RSG mass-loss does an excellent job of matching the observations, and we can readily rule out significantly larger values., Comment: ApJ, accepted

Massive stars are the drivers of the chemical evolution of dwarf galaxies. We review here the basics of massive star evolution and the specificities of stellar evolution in low-Z environment. We discuss nucleosynthetic aspects and what observations could constrain our view on the first generations of stars., Comment: 8 pages, 4 figures. Proceedings of IAUS 344 "Dwarf Galaxies: From the Deep Universe to the Present" held in Vienna (Austria) August 20-24 2018

In the last decade, the Kepler and CoRoT space-photometry missions have demonstrated the potential of asteroseismology as a novel, versatile and powerful tool to perform exquisite tests of stellar physics, and to enable precise and accurate characterisations of stellar properties, with impact on both exoplanetary and Galactic astrophysics. Based on our improved understanding of the strengths and limitations of such a tool, we argue for a new small/medium space mission dedicated to gathering high-precision, high-cadence, long photometric series in dense stellar fields. Such a mission will lead to breakthroughs in stellar astrophysics, especially in the metal poor regime, will elucidate the evolution and formation of open and globular clusters, and aid our understanding of the assembly history and chemodynamics of the Milky Way's bulge and few nearby dwarf galaxies., Comment: White paper submitted in response to Voyage 2050 long-term plan in the ESA Science Programme, now published in Experimental Astronomy, see https://link.springer.com/article/10.1007/s10686-021-09711-1

Current burning issues in stellar physics, for both hot and cool stars, concern their magnetism. In hot stars, stable magnetic fields of fossil origin impact their stellar structure and circumstellar environment, with a likely major role in stellar evolution. However, this role is complex and thus poorly understood as of today. It needs to be quantified with high-resolution UV spectropolarimetric measurements. In cool stars, UV spectropolarimetry would provide access to the structure and magnetic field of the very dynamic upper stellar atmosphere, providing key data for new progress to be made on the role of magnetic fields in heating the upper atmospheres, launching stellar winds, and more generally in the interaction of cool stars with their environment (circumstellar disk, planets) along their whole evolution. UV spectropolarimetry is proposed on missions of various sizes and scopes, from POLLUX on the 15-m telescope LUVOIR to the Arago M-size mission dedicated to UV spectropolarimetry., Comment: 25 pages, 11 figures, white paper submitted to ESA Voyage 2050 call