Evolution Of Massive Stars Research Articles

X-Shooting ULLYSES: Massive stars at low metallicity. VI. Atmosphere and mass-loss properties of O-type giants in the Small Magellanic Cloud

Mass loss through a stellar wind is an important physical process that steers the evolution of massive stars and controls the properties of their end-of-life products, such as the supernova type and the mass of compact remnants. To probe its role in stellar evolution over cosmic time, mass loss needs to be studied as function of metallicity. For mass loss to be accurately quantified, the wind structure needs to be established jointly with the characteristics of small-scale inhomogeneities in the outflow, which are known as wind clumping. We aim to improve empirical estimates of mass loss and wind clumping for hot main-sequence massive stars, study the dependence of both properties on the metal content, and compare the theoretical predictions of mass loss as a function of metallicity to our findings. Using the model atmosphere code Fastwind and the genetic algorithm fitting method Kiwi-GA we analyzed the optical and ultraviolet spectra of 13 O-type giant to supergiant stars in the Small Magellanic Cloud galaxy, which has a metallicity of approximately one-fifth of that of the Sun. We quantified the stellar global outflow properties, such as the mass-loss rate and terminal wind velocity, and the wind clumping properties. To probe the role of metallicity, our findings were compared to studies of Galactic and Large Magellanic Cloud samples that were analyzed with similar methods, including the description of clumping. We find significant variations in the wind clumping properties of the target stars, with clumping starting at flow velocities $0.01 - 0.23$ of the terminal wind velocity and reaching clumping factors $f_ cl = 2 - 30$. In the luminosity ($ L / odot = 5.0 - 6.0$) and metallicity ($Z/Z_ odot = 0.2 - 1$) range we considered, we find that the scaling of the mass loss $ M $ with metallicity $Z$ varies with luminosity. At $ odot = 5.75$, we find $ M Z^m$ with $m = 1.02 0.30$, in agreement with pioneering work in the field within the uncertainties. For lower luminosities, however, we obtain a significantly steeper scaling of $m > 2$. The monotonically decreasing $m(L)$ behavior adds a complexity to the functional description of the mass-loss rate of hot massive stars. Although the trend is present in the predictions, it is much weaker than we found here. However, the luminosity range for which $m$ is significantly larger than previously assumed (at $ odot 5.4$) is still poorly explored, and more studies are needed to thoroughly map the empirical behavior, in particular, at Galactic metallicity.

Searching for triple systems unbound by supernovae

Abstract A large fraction of massive stars are found in higher order systems where the presence of a tertiary may significantly modify the system’s evolution. In particular, it can lead to increased numbers of compact object binaries and accelerate their mergers with important implications for gravitational wave observations. Using Gaia, we constrain the number of Galactic supernovae that produce unbound triples. We do this by searching 8 supernova remnants for stars with consistent Gaia parallaxes and paths intersecting near the center of the supernova remnant at a time consistent with the age of the remnant. We find no candidates for unbound triple systems. Combined with prior work, less than 11.4% of supernovae leave behind unbound triples at a 90% confidence limit. The absence of such systems limits their role in the evolution of massive stars and the formation of merging compact objects.

Lower-mass-gap Black Holes in Dense Star Clusters

The existence of compact stellar remnants in the mass range 2–5 M ⊙ has long been debated. This so-called lower-mass gap (LMG) was initially suggested by the lack of low-mass X-ray binary observations with accretors about 2–5 M ⊙, but it has recently been called into question following newer observations, including an LMG candidate with a millisecond pulsar (MSP) companion in the dense globular cluster NGC 1851. Here, we model NGC 1851 with a grid of similar dense star clusters utilizing the state-of-the-art Monte Carlo N-body code Cluster Monte Carlo, and we specifically study the formation of LMG black holes (BHs). We demonstrate that both massive star evolution and dynamical interactions can contribute to forming LMG BHs. In general, the collapse of massive remnants formed through mergers of neutron stars (NSs) or massive white dwarfs produces the largest number of LMG BHs among all formation channels. However, in more massive clusters, supernova core collapse can contribute comparable numbers. Our NGC 1851-like models can reproduce MSP—LMG BH binaries similar to the observed system. Additionally, the LMG BHs can also become components of dynamically assembled binaries, and some will be in merging BH–NS systems similar to the recently detected gravitational wave source GW230529. However, the corresponding merger rate is probably ≲1 Gpc−3 yr−1.

Imaging the innermost circumstellar environment of the red supergiant WOH G64 in the Large Magellanic Cloud

Significant mass loss in the red supergiant (RSG) phase has great influence on the evolution of massive stars and their final fate as supernovae. We present near-infrared interferometric imaging of the circumstellar environment of the dust-enshrouded RSG WOH G64 in the Large Magellanic Cloud. WOH G64 was observed with the GRAVITY instrument at ESO's Very Large Telescope Interferometer (VLTI) at 2.0--2.45 mu m . We succeeded in imaging the innermost circumstellar environment of WOH G64 -- the first interferometric imaging of an RSG outside the Milky Way. The reconstructed image reveals elongated compact emission with a semimajor and semiminor axis of sim 2 and sim 1.5 mas (sim 13 and 9 $R_ star ), respectively. The GRAVITY data show that the stellar flux contribution at 2.2 mu m \ at the time of our observations in 2020 is much lower than predicted by the optically and geometrically thick dust torus model based on the VLTI/MIDI data taken in 2005 and 2007. We found a significant change in the near-infrared spectrum of WOH G64: while the (spectro)photometric data taken at 1--2.5 mu m \ before 2003 show the spectrum of the central RSG with H$_2$O absorption, the spectra and prime photometric data taken after 2016 are characterized by a monotonically rising continuum with very weak signatures of H$_2$O . This spectral change likely took place between December 2009 and 2016. On the other hand, the mid-infrared spectrum obtained in 2022 with VLT/VISIR agrees well with the spectra obtained before 2007. The compact emission imaged with GRAVITY and the near-infrared spectral change suggest the formation of hot new dust close to the star, which gives rise to the monotonically rising near-infrared continuum and the high obscuration of the central star. The elongation of the emission may be due to the presence of a bipolar outflow or effects of an unseen companion.

Constraining the overcontact phase in massive binary evolution. III. Period stability of known B+B and O+B overcontact systems

Binary systems play a crucial role in massive star evolution. Systems composed of B-type and O-type stars are of particular interest due to their potential to lead to very energetic phenomena or the merging of exotic compact objects. We aim to determine the orbital period variations of a sample of B+B and O+B massive overcontact binaries, with the primary objectives of characterizing the evolutionary timescales of these systems and addressing the existing discrepancy between observational data and theoretical predictions derived from population synthesis models. We used Period04 to analyze archival photometric data going back a century for a sample of seven binary systems to measure their orbital periods. We then determine the period variations using a linear fit. We find that the period variation timescales of five truly overcontact binary systems align with the nuclear timescale, in agreement with previous findings for more massive overcontact binaries. Additionally, we noticed a clear distinction between the five systems that had been unambiguously classified as overcontact systems and both SV Cen and VFTS 066, which seem to be evolving on thermal timescales and might be misclassified as overcontact systems. In the case of the five overcontact binaries, our results indicate a noticeable mismatch between the observational data and the theoretical predictions derived from population synthesis models. Furthermore, our results suggest that additional physical mechanisms must be investigated to compare the observed variations more thoroughly with theoretical predictions.

Boron depletion in Galactic early B-type stars reveals two different main sequence star populations

Context. The evolution and fate of massive stars are thought to be affected by rotationally induced internal mixing. The surface boron abundance is a sensitive tracer of this in early B-type main sequence stars. Aims. We test current stellar evolution models of massive main sequence stars which include rotational mixing through a systematic study of their predicted surface boron depletion. Methods. We construct a dense grid of rotating single star models using MESA, for which we employ a new nuclear network which follows all the stable isotopes up to silicon, including lithium, beryllium, boron, as well as the radioactive isotope aluminium-26. We also compile the measured physical parameters of the 90 Galactic early B-type stars with boron abundance information. We then compare each observed stars with our models through a Bayesian analysis, which yields the mixing efficiency parameter with which the star is reproduced the best, and the probability that it is represented by the stellar models. Results. We find that about two-thirds of the sample stars are well represented by the stellar models, with the best agreement achieved for a rotational mixing efficiency of ∼50% compared to the widely adopted value. The remaining one third of the stars, of which many are strongly boron depleted slow rotators, are largely incompatible with our models, for any rotational mixing efficiency. We investigate the observational incidence of binary companions and surface magnetic fields, and discuss their evolutionary implications. Conclusions. Our results confirm the concept of rotational mixing in radiative stellar envelopes. On the other hand, we find that a different boron depletion mechanism, and likely a different formation path, is required to explain about one-third of the sample stars. The large spread in the surface boron abundances of these stars may hold a clue to understanding their origin.

The Big Bang versus the Steady State: Gamow, Hoyle and Ryle, rivals in cosmology

Abstract This account of the controversy concerning the Big Bang and Steady State theories contrasts the different approaches to science of George Gamow, Fred Hoyle and Martin Ryle. Einstein’s general relativity offered theoretical physicists an indispensable tool kit for making mathematical models of the Universe. Georges Lemaître made an astonishing impact in 1927 with his proposal that the expanding Universe started from a decaying super-radioactive primeval atom. Gamow refined Lemaître’s model by proposing that primordial nucleosynthesis in the early expanding Universe had formed the chemical elements. Hoyle, who by contrast had discovered the conditions under which stellar nucleosynthesis could occur via the evolution and explosion of massive stars, publicly dismissed Gamow’s cosmology as an unnecessary speculation. Hoyle, Hermann Bondi and Thomas Gold proposed the alternative steady state theory in 1948. In 1950 Ryle constructed the first radio telescope with the sensitivity to discriminate between the rival models of the Universe. Steady state cosmology survived Ryle’s first surveys of cosmic radio sources because Ryle’s data was inadequate. However results from the much improved Third and Fourth Cambridge Surveys demonstrated that our Universe had evolved. The Ryle-Hoyle debates became intensely passionate and personal. In 1965 the discovery of the cosmic microwave background radiation settled the argument in favour of the Lemaître Big Bang evolving Universe.

The YMDB catalog: Young massive detached binaries for the determination of high-precision absolute stellar parameters

Context. Massive stars play a crucial role in the cosmic dynamics and chemical evolution of galaxies. Despite their significance, our understanding of their evolution and properties remains limited. An accurate determination of stellar parameters, such as the mass and radius, is essential for advancing our knowledge. Detached eclipsing binaries (DEBs) are particularly valuable for these determinations due to the minimal interaction between their stellar components, allowing for precise measurements. Aims. This study aims to introduce the Young Massive Detached Binary (YMDB) catalog, designed to address the gap in the high-precision absolute parameter determination for young massive stars. By focusing on DEBs within the spectral range O9-B1, this catalog seeks to provide a reliable database for future astronomical studies and improve our understanding of massive star evolution. Methods. We conducted a photometric analysis of 87 young massive stars in detached eclipsing systems using TESS light curves (LCs) that were processed through a custom pipeline. This analysis involved determining the amplitude of magnitude variations, orbital periods, times of minima, eccentricities, and the presence of apsidal motion and heartbeat phenomena. A thorough literature review was performed to obtain MK spectral classifications. We performed our own spectral classification of 19 systems to support the sample where a new classification was lacking or inconclusive. Results. The analysis identified 20 previously unreported binary systems, with 13 newly recognized as variable stars. Among the 87 stars examined, 30 are confirmed as YMDB members, and 25 are candidates pending spectral classification. The exclusion of the remaining 32 stars is attributed to unsuitable spectral types or their nondetached binary nature. Notable findings include the identification of new LC classifications, eccentricities in 13 systems, and heartbeat phenomena in several targets. Conclusions. The YMDB catalog offers a resource of high-quality LCs and reliable stellar classifications, serving as a valuable tool for the astronomical community.

EWOCS-II: X-ray properties of the Wolf–Rayet stars in the young Galactic super star cluster Westerlund 1

Context. Wolf–Rayet (WR) stars are massive evolved stars that exhibit particularly fast and dense stellar winds. Although they constitute a very short phase near the end of a massive star’s life, they play a crucial role in the evolution of massive stars and have a substantial impact on their surrounding environment. Aims. We present the most comprehensive and deepest X-ray study to date of the properties of the richest Wolf–Rayet population observed in a single stellar cluster, Westerlund 1 (Wd1). By examining the X-ray signatures of WR stars, we aim to shed light on the hottest plasma in their stellar winds and gain insights into whether they exist as single stars or within binary systems. Methods. This work is based on 36 Chandra observations obtained from the “Extended Westerlund 1 and 2 Open Clusters Survey” (EWOCS) project, plus 8 archival Chandra observations. The overall exposure depth Ms) and baseline of the EWOCS observations extending over more than one year enable us to perform a detailed photometric, colour, and spectral analysis, as well as to search for short- and long-term periodicity. Results. In X-rays, we detect 20 out of the 24 known Wolf–Rayet stars in Wd1 down to an observed luminosity of ~7 × 1029 erg s−1 (assuming a distance of 4.23 kpc to Wd1), with 8 WR stars being detected in X-rays for the first time. Nine stars show clear evidence of variability over the year-long baseline, with clear signs of periodicity. The X-ray colours and spectral analysis reveal that the vast majority of the WR stars are hard X-ray sources (kT≥2.0 keV). The Fe XXV emission line at ~6.7 keV, which commonly originates from the wind–wind collision zone in binary systems, is detected for the first time in the spectra of 17 WR stars in Wd1. In addition the ~6.4 keV fluorescent line is observed in the spectra of three stars, which are among the very few massive stars exhibiting this line, indicating that dense cold material coexists with the hot gas in these systems. Overall, our X-ray results alone suggest a very high binary fraction (≥80%) for the WR star population in Wd1. When combining our results with properties of the WR population from other wavelengths, we estimate a binary fraction of ≥92%, which could even reach unity. This suggests that either all the most massive stars are found in binary systems within Wd1, or that binarity is essential for the formation of such a rich population of WR stars.

The Scale of Stellar Yields: Implications of the Measured Mean Iron Yield of Core Collapse Supernovae

The scale of α-element yields is difficult to predict from theory because of uncertainties in massive star evolution, supernova physics, and black hole formation, and it is difficult to constrain empirically because the impact of higher yields can be compensated by greater metal loss in galactic winds. We use a recent measurement of the mean iron yield of core collapse supernovae (CCSN) by Rodriguez et al., y¯Fecc=0.058±0.007M⊙ , to infer the scale of α-element yields by assuming that the plateau of [α/Fe] abundance ratios observed in low-metallicity stars represents the yield ratio of CCSN. For a plateau at [α/Fe]cc = 0.45, we find that the population-averaged yields of O and Mg are about equal to the solar abundance of these elements, logyOcc/ZO,⊙=logyMgcc/ZMg,⊙=−0.01±0.1 , where yXcc is the mass of element X produced by massive stars per unit mass of star formation. The inferred O and Fe yields agree with predictions of the Sukhbold et al. CCSN models assuming their Z9.6+N20 neutrino-driven engine, a scenario in which many progenitors with M < 40M ⊙ implode to black holes rather than exploding. The yields are lower than assumed in many models of the galaxy mass–metallicity relation, reducing the level of outflows needed to match observed abundances. Our one-zone chemical evolution models with η=Ṁout/Ṁ*≈0.6 evolve to solar metallicity at late times. By further requiring that models reach [α/Fe] ≈ 0 at late times, we infer a Hubble-time integrated Type Ia supernova rate of 1.1×10−3M⊙−1 , compatible with estimates from supernova surveys.

The Discovery of Three Galactic Wolf–Rayet Stars

Wolf–Rayet stars (WRs) are evolved massive stars in the brief stage before they undergo core collapse. Not only are they rare, but they also can be particularly difficult to find due to the high extinction in the Galactic plane. This paper discusses the discovery of three new Galactic WRs previously classified as Hα emission stars, but thanks to Gaia spectra, we were able to identify the broad, strong emission lines that characterize WRs. Using the Lowell Discovery Telescope and the DeVeny spectrograph, we obtained spectra for each star. Two are WC9s, and the third is a WN6 + O6.5 V binary. The latter is a known eclipsing system with a 4.4 day period from ASAS-SN data. We calculate absolute visual magnitudes for all three stars to be between −7 and −6, which is consistent with our expectations of these subtypes. These discoveries highlight the incompleteness of the WR census in our local volume of the Milky Way and suggest the potential for future Galactic WR discoveries from Gaia low-dispersion spectra. Furthermore, radial velocity studies of the newly found binary will provide direct mass estimates and orbital parameters, adding to our knowledge of the role that binarity plays in massive star evolution.

The Evolution of Massive Binary Stars

Massive stars play a major role in the evolution of their host galaxies and serve as important probes of the distant Universe. It has been established that the majority of massive stars reside in close binaries and interact with their companion stars during their lifetimes. Such interactions drastically alter their life cycles and complicate our understanding of their evolution, but are also responsible for the production of interesting and exotic interaction products. ▪Extensive observation campaigns with well-understood detection sensitivities have enabled the conversion of observed properties into intrinsic characteristics, facilitating a direct comparison to theory.▪Studies of large samples of massive stars in our Galaxy and the Magellanic Clouds have unveiled new types of interaction products, providing critical constraints on the mass transfer phase and the formation of compact objects.▪The direct detection of gravitational waves has revolutionized the study of stellar mass compact objects, providing a new window to study massive star evolution. Their formation processes are, however, still unclear. The known sample of compact object mergers will increase by orders of magnitude in the coming decade, which is vastly outgrowing the number of stellar-mass compact objects detected through electromagnetic radiation.

The evolution of accreting population III stars at 10−6–103 M⊙ yr−1

Context. The first stars formed over five orders of magnitude in mass by accretion in primordial dark matter halos. Aims. We study the evolution of massive, very massive and supermassive primordial (Pop III) stars over nine orders of magnitude in accretion rate. Methods. We use the stellar evolution code GENEC to evolve accreting Pop III stars from 10−6–103 M⊙ yr−1 and study how these rates determine final masses. The stars are evolved until either the end central Si burning or they encounter the general relativistic instability (GRI). We also examine how metallicity affects the evolution of the star at one accretion rate. Results. At rates below ∼2.5 × 10−5 M⊙ yr−1 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⊙ is ∼5 × 10−5 M⊙ yr−1, well within the range of infall rates found in numerical simulations of halos that cool via H2, ≲10−3 M⊙ yr−1. At rates of 5 × 10−5 M⊙ yr−1 to 4 × 10−2 M⊙ yr−1, like those expected for halos cooling by both H2 and Lyα, 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(Teff) = 3.8 and luminosities that can reach 1010.5 L⊙. At accretion rates of 100–1000 M⊙ yr−1 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 ∼106 M⊙ without ever having become a star. Conclusions. Our models corroborate previous studies of Pop III stellar evolution with and without hydrodynamics over separate, smaller ranges in accretion rate. They also 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.

Exploring the stellar rotation of early-type stars in the LAMOST medium-resolution survey

Context. Stellar rotation significantly shapes the evolution of massive stars, yet the interplay of mass and metallicity remains elusive, limiting our capacity to construct accurate stellar evolution models and to better estimate the impact of rotation on the chemical evolution of galaxies. Aims. Our goal is to investigate how mass and metallicity influence the rotational evolution of A-type stars on the main sequence (MS). We seek to identify deviations in rotational behaviors that could serve as new constraints for existing stellar models. Methods. Using the LAMOST Median-Resolution Survey Data Release 9, we derived stellar parameters for a population of 104 752 A-type stars. Our study focused on the evolution of surface rotational velocities and their dependence on mass and metallicity in 84 683 “normal” stars. Results. Normalizing surface rotational velocities to zero age main sequence (ZAMS) values revealed a prevailing evolutionary profile from 1.7 to 4.0 M⊙. This profile features an initial rapid acceleration until t/tMS = 0.25 ± 0.1 and potentially a second acceleration peak near t/tMS = 0.55 ± 0.1 for stars heavier than 2.5 M⊙, followed by a steady decline and a “hook” feature at the end. Surpassing theoretical expectations, the initial acceleration likely stems from a concentrated distribution of angular momentum at the ZAMS, resulting in a prolonged increase in speed. A transition phase for stars with 2.0 &lt; M/M⊙ &lt; 2.3 emerged, a region where evolutionary tracks remain uncertain. Stellar expansion primarily drives the spin down in the latter half of the MS, accompanied by significant influence from inverse meridional circulation. The inverse circulation becomes more efficient at lower metallicities, explaining the correlation of the slope of this deceleration phase with metallicity from –0.3 dex up to 0.1 dex. The metal-poor subsample (−0.3 dex &lt; [M/H]&lt; − 0.1 dex) starts with lower velocities at the ZAMS, suggesting that there is a metallicity-dependent mechanism that removes angular momentum during star formation. The proportion of fast rotators decreases with an increase in metallicity, up to log(Z/Z⊙)∼ − 0.2, a trend consistent with observations of OB-type stars found in the Small and Large Magellanic Clouds.

Investigation of nuclear structure and [formula omitted]-decay properties of As isotopes

The nuclear ground state properties of 67−80As nuclei have been investigated within the framework of relativistic mean field (RMF) approach. The RMF model with density-dependent (DD-ME2) interaction is utilized for the calculation of potential energy curves and the nuclear ground-state deformation parameters (β2) of selected As isotopes. Later, the β-decay properties of As isotopes were studied using the proton–neutron quasi particle random phase approximation (pn-QRPA) model. These include Gamow Tellar (GT) strength distributions, log ft values, β-decay half-lives, stellar β± decays and stellar electron/positron capture rates. The β2 values computed from RMF model were employed in the pn-QRPA model as an input parameter for the calculations of β-decay properties for 67−80As. The calculated log ft values were in decent agreement with the measured data. The predicted β-decay half-lives matched the experimental values within a factor of 10. The stellar rates were compared with the shell model results. Only at high temperature and density values, the sum of β+ and electron capture rates had a finite contribution. On the other hand, the sum of β− and positron capture rates were sizeable only at low density and high temperature values. For all such cases, the pn-QRPA rates were found to be bigger than the shell model rates up to a factor of 33 or more. The findings reported in the current investigation could prove valuable for simulating the late-stage stellar evolution of massive stars.

Impact of main sequence mass loss on the appearance, structure, and evolution of Wolf-Rayet stars

Context. Stellar winds are one of the most important drivers of massive star evolution and are a vital source of chemical, mechanical, and radiative feedback on the galactic scale. Despite its significance, mass loss remains a major uncertainty in stellar evolution models. In particular, the interdependencies between the different approaches and the subsequent evolutionary stages and predicted observable phenomena are far from being systematically understood. Aims. In this study, we examine the impact of main sequence mass loss on the structure of massive stars throughout their entire evolution. Particular focus is placed on the consequences in terms of entrance into the Wolf-Rayet (WR) regime and the subsequent evolution. Methods. Using the Geneva stellar evolution code (GENEC), we computed grids of single, nonrotating stellar models at solar and Large Magellanic Cloud (LMC) metallicities of initial masses between 20 and 120 solar masses, with two representative prescriptions for high and low main sequence mass loss. Results. We obtain detailed numerical predictions regarding the structure and evolution of massive stars, and infer the role of main sequence mass loss by comparison of the mass-loss rate prescriptions. We present implications for the overall evolutionary trajectory, including the evolution of WR stars, as well as the effect on stellar yields and stellar populations. Conclusions. Mass loss during the main sequence plays an important role because of its ability to affect the sequence and duration of all subsequent phases. We identify several distinct evolutionary paths for massive stars, which are significantly influenced by the chosen main sequence mass-loss description. We also discuss the impact of uncertainties – other than that regarding mass loss – on the evolution, in particular those relating to convection. We further demonstrate that not only the total mass loss but also the specific mass-loss history throughout a star’s life is a crucial determinant of many aspects, such as the resulting stellar yields.

The effects of surface fossil magnetic fields on massive star evolution: V. Models at low metallicity

Abstract At metallicities lower than that of the Small Magellanic Cloud, it remains essentially unexplored how fossil magnetic fields, forming large-scale magnetospheres, could affect the evolution of massive stars, thereby impacting the fundamental building blocks of the early Universe. We extend our stellar evolution model grid with representative calculations of main-sequence, single-star models with initial masses of 20 and 60 M⊙, including appropriate changes for low-metallicity environments (Z = 10−3–10−6). We scrutinise the magnetic, rotational, and chemical properties of the models. When lowering the metallicity, the rotational velocities can become higher and the tendency towards quasi-chemically homogeneous evolution increases. While magnetic fields aim to prevent the development of this evolutionary channel, the weakening stellar winds lead to less efficient magnetic braking in our models. Since the stellar radius is almost constant during a blueward evolution caused by efficient chemical mixing, the surface magnetic field strength remains unchanged in some models. We find core masses at the terminal-age main sequence between 22 and 52 M⊙ for initially 60 M⊙ models. This large difference is due to the vastly different chemical and rotational evolution. We conclude that in order to explain chemical species and, in particular, high nitrogen abundances in the early Universe, the adopted stellar models need to be under scrutiny. The assumptions regarding wind physics, chemical mixing, and magnetic fields will strongly impact the model predictions.

Exploring the Origin of Ultralong Gamma-Ray Bursts: Lessons from GRB 221009A

The brightest gamma-ray burst (GRB) ever, GRB 221009A, displays ultralong GRB (ULGRB) characteristics, with a prompt emission duration exceeding 1000 s. To constrain the origin and central engine of this unique burst, we analyze its prompt and afterglow characteristics and compare them to the established set of similar GRBs. To achieve this, we statistically examine a nearly complete sample of Swift-detected GRBs with measured redshifts. We categorize the sample to bronze, silver, and gold by fitting a Gaussian function to the log-normal of T 90 duration distribution and considering three subsamples respectively to 1, 2, and 3 times of the standard deviation to the mean value. GRB 221009A falls into the gold subsample. Our analysis of prompt emission and afterglow characteristics aims to identify trends between the three burst groups. Notably, the gold subsample (a higher likelihood of being ULGRB candidates) suggests a collapsar scenario with a hyperaccreting black hole as a potential central engine, while a few GRBs (GRB 060218, GRB 091024A, and GRB 100316D) in our gold subsample favor a magnetar. Late-time near-IR observations from 3.6 m Devasthal Optical Telescope rule out the presence of any bright supernova associated with GRB 221009A in the gold subsample. To further constrain the physical properties of ULGRB progenitors, we employ the tool MESA to simulate the evolution of low-metallicity massive stars with different initial rotations. The outcomes suggest that rotating (Ω ≥ 0.2 Ωc) massive stars could potentially be the progenitors of ULGRBs within the considered parameters and initial inputs to MESA.

Modeling the Progenitor Stars of Observed Type IIP Supernovae

Type IIP supernovae (SNe IIP) are thought to originate from the explosion of massive stars >10 M ⊙. Their luminosity is primarily powered by the explosion energy and the radioactive decay energy of 56Co, with the photosphere location regulated by hydrogen recombination. However, the physical connections between SNe IIP and their progenitor stars remain unclear. This paper presents a comprehensive study of SNe IIP and their progenitor stars by using the one-dimensional stellar evolution code, MESA. Our model grids consider the effects of stellar metallicity, mass, and rotation in the evolution of massive stars, as well as the explosion energy and 56Ni production in modeling supernovae. To elucidate the observed SNe IIP and their origins, we compare their light curves (LCs) with our models. Furthermore, we investigate the impact of stellar parameters on LCs by considering stellar mass, metallicity, rotation, explosion energy, and 56Ni production. We find that more massive stars exhibit longer plateaus due to increased photon diffusion time caused by massive ejecta. Higher metallicity leads to increased opacity and mass loss of progenitor stars. Rapid rotation affects internal stellar structures, enhancing convective mixing and mass loss, potentially affecting the plateau’s brightness and duration. Higher explosion energy results in brighter but shorter plateaus due to faster-moving ejecta. 56Ni mass affects late-time luminosity and plateau duration, with larger masses leading to slower declines.

Shell mergers in the late stages of massive star evolution: new insight from 3D hydrodynamic simulations

ABSTRACT One-dimensional (1D) stellar evolution models are widely used across various astrophysical fields, however they are still dominated by important uncertainties that deeply affect their predictive power. Among those, the merging of independent convective regions is a poorly understood phenomenon predicted by some 1D models but whose occurrence and impact in real stars remain very uncertain. Being an intrinsically multi-D phenomenon, it is challenging to predict the exact behaviour of shell mergers with 1D models. In this work, we conduct a detailed investigation of a multiple shell merging event in a 20 M$_\odot$ star using 3D hydrodynamic simulations. Making use of the active tracers for composition and the nuclear network included in the 3D model, we study the merging not only from a dynamical standpoint but also considering its nucleosynthesis and energy generation. Our simulations confirm the occurrence of the merging also in 3D, but reveal significant differences from the 1D case. Specifically, we identify entrainment and the erosion of stable regions as the main mechanisms that drive the merging, we predict much faster convective velocities compared to the mixing-length theory velocities, and observe multiple burning phases within the same merged shell, with important effects for the chemical composition of the star, which presents a strongly asymmetric (dipolar) distribution. We expect that these differences will have important effects on the final structure of massive stars and thus their final collapse dynamics and possible supernova explosion, subsequently affecting the resulting nucleosynthesis and remnant.