Stellar Models Research Articles

The JWST Resolved Stellar Populations Early Release Science Program. VI. Identifying Evolved Stars in Nearby Galaxies

We present an investigation of evolved stars in the nearby star-forming galaxy Wolf–Lundmark–Melotte (WLM), using Near-Infrared Camera (NIRCam) imaging from the JWST Resolved Stellar Populations Early Release Science program. We find that various combinations of the F090W, F150W, F250M, and F430M filters can effectively isolate red supergiants and thermally pulsing asymptotic giant branch (TP-AGB) stars from one another, while also providing a reasonable separation of the primary TP-AGB subtypes: carbon-rich C-type stars and oxygen-rich M-type stars. The classification scheme we present here agrees very well with the well-established Hubble Space Telescope (HST) medium-band filter technique. The ratio of C to M-type stars is 0.8 ± 0.1 for both the new JWST and the HST classifications, which is within 1σ of empirical predictions from optical narrowband CN and TiO filters. The evolved star colors show good agreement with the predictions from the PARSEC + COLIBRI stellar evolutionary models, and the models indicate a strong metallicity dependence that makes stellar identification even more effective at higher metallicity. However, the models also indicate that evolved star identification with NIRCam may be more difficult at lower metallicities. We test every combination of NIRCam filters using the models and present additional filters that are also useful for evolved star studies. We also find that ≈90% of the dusty evolved stars are carbon rich, suggesting that carbonaceous dust dominates the present-day dust production in WLM, similar to the findings in the Magellanic Clouds. These results demonstrate the usefulness of NIRCam in identifying and classifying dust-producing stars without the need for mid-infrared data.

Constraining Wind-driven Accretion onto Gaia BH3 with Chandra

Gaia BH3 is the most massive known stellar-origin black hole in the Milky Way, with a mass M • ≈ 33 M ⊙. Detected from Gaia’s astrometry, this black hole is in the mass range of those observed via gravitational waves, whose nature is still highly debated. Hosted in a binary system with a companion giant star that is too far away for Roche-lobe mass transfer, this black hole could nonetheless accrete at low levels due to wind-driven mass loss from its companion star, thus accreting in advection-dominated accretion flow, or ADAF, mode. Using stellar wind models, we constrain its Eddington ratio in the range 10−9 < f Edd < 10−7, corresponding to radiative efficiencies 5 × 10−5 < ϵ < 10−3, compatible with radiatively inefficient accretion modes. Chandra ACIS-S observed this object and obtained the most sensitive upper bound of its [2–10] keV flux: F X < 3.25 × 10−15 erg s−1 cm−2 at 90% confidence level corresponding to L [2–10] < 2.10 × 1029 erg s−1. Using ADAF emission models, we constrained its accretion rate to f Edd < 4.91 × 10−7 at the apastron, in agreement with our theoretical estimate. At the periastron, we expect fluxes ∼50 times larger. Because of the inferred low rates, accretion did not significantly contribute to black hole growth over the system’s lifetime. Detecting the electromagnetic emission from Gaia BH3 will be fundamental to informing stellar wind and accretion disk models.

Anisotropic Durgapal-Fuloria compact stars in f(R) gravity

This study presented a new exact solution for anisotropic compact stellar objects within the framework of f(R)=R+αR2 gravity. In this context, the Durgapal-Fuloria metric potential has been employed to solve the field equation derived for f(R) theory. Furthermore, we have derived the generalized Darmois-Israel junction condition necessary for seamlessly connecting the interior region to the Schwarzschild exterior metric across the boundary hypersurface of the star in the context of f(R) gravity, and the interior solution is matched with the Schwarzschild exterior metric over the bounding surface of a compact star. These junction conditions stipulate that the pressure must not be zero at the boundary and should be proportional to the non-linear terms of f(R) gravity, a crucial aspect often overlooked by many researchers when investigating compact stellar models. Additionally, we derived the values of these parameters by using observational data of various compact stars (CSs), namely Her X-1, SAX J1808.4-3658, SMC X-1, LMC X-4, Cen X-3, 4U 1820-30, PSR J1903+327, 4U 1608-52, Vela X-1, and PSR J1416-2230. This approach enables us to investigate the comprehensive analysis of solutions numerically and graphically. We conducted various physical tests, including gradient of energy density and pressures, anisotropy, stability, equilibrium conditions, energy-density constraints, mass function, compactness, redshift, and adiabatic index, to assess the feasibility of our models. Our findings demonstrate the consistent behavior of our models provides a satisfactory physical situation as far as the observational results are confirmed.

Validating full-spectrum fitting with a synthetic integral-field spectroscopic observation of the milky way

Abstract Ongoing deep IFS observations of disk galaxies provide opportunities for comparison with the Milky Way (MW) to understand galaxy evolution. However, such comparisons are marred by many challenges such as selection effects, differences in observations and methodology, and proper validation of full-spectrum fitting methods. In this study, we present a novel code GalCraft to address these challenges by generating mock IFS data cubes of the MW using simple stellar population models and a mock MW stellar catalog derived from E-Galaxia. We use the widely adopted full-spectrum fitting code pPXF to investigate the ability to recover kinematics and stellar populations for an edge-on mock MW IFS observation. We confirm that differences in kinematics, mean age, [M/H], and [α/Fe] between thin and thick disks can be distinguished. However, the age distribution is overestimated in the ranges between 2 − 4 and 12 − 14 Gyr compared to the expected values. This is likely due to the age spacing and degeneracy of SSP templates. We find systematic offsets in the recovered kinematics due to insufficient spectral resolution and the variation of line-of-sight velocity distribution with age and [M/H]. With future higher resolution and multi-[α/Fe] SSP templates, GalCraft will be useful to validate key signatures such as [α/Fe]-[M/H] distribution at different R and |z| and potentially infer radial migration and kinematic heating efficiency to study detailed chemodynamical evolution of MW-like galaxies.

Towards a universal analytical model for Population III star formation: interplay between feedback and fragmentation

ABSTRACT JWST has brought us new insights into Cosmic Dawn with tentative detection of the unique signatures of metal-free Population III (Pop III) stars, such as strong He II emission, extremely blue ultraviolet spectrum, and enhanced nitrogen abundance. Self-consistent theoretical predictions of the formation rates, sites, and masses of Pop III stars are crucial for interpreting the observations, but are challenging due to complex physical processes operating over the large range of length-scales involved. One solution is to combine analytical models for the small-scale star formation process with cosmological simulations that capture the large-scale physics such as structure formation, radiation backgrounds, and baryon-dark matter streaming motion that regulate the conditions of Pop III star formation. We build an analytical model to predict the final masses of Pop III stars/clusters from the properties of star-forming clouds, based on the key results of small-scale star formation simulations and stellar evolution models. Our model for the first time considers the interplay between feedback and fragmentation and covers different modes of Pop III star formation ranging from ordinary small ($\sim\!{10{-}2000}\ \rm M_\odot$) clusters in molecular-cooling clouds to massive ($\gtrsim\!{10^{4}}\ \rm M_\odot$) clusters containing supermassive ($\sim\!{10^{4}{-}3}\times 10^{5}\ \rm M_\odot$) stars under violent collapse of atomic-cooling clouds with large gas accretion rates of $\gtrsim\!{0.1}\ \rm M_\odot \ yr^{-1}$. As an example, the model is applied to the Pop III star-forming clouds in the progenitors of typical haloes hosting high-z luminous quasars ($M_{\rm h}\sim 10^{12}\ \rm M_\odot$ at $z\sim 6$), which shows that formation of Pop III massive clusters is common ($\sim\!{20{-}70}{{\ \rm per\ cent}}$) in such biased ($\sim\!{4}\sigma$) regions, and the resulting heavy black hole seeds from supermassive stars can account for a significant fraction of observed luminous ($\gtrsim\!{10^{46}}\ \rm erg\ s^{-1}$) quasars at $z\sim 6$.

The rates and host galaxies of pair-instability supernovae through cosmic time: predictions from BPASS and IllustrisTNG

ABSTRACT Pair-instability supernovae (PISNe) have long been predicted to be the final fates of near-zero-metallicity very massive stars ($Z \lt Z_\odot /3$, M$_\mathrm{ZAMS} \gtrsim 140\, \text{M}_\odot$). However, no definite PISN has been observed to date, leaving theoretical modelling validation open. To investigate the observability of these explosive transients, we combine detailed stellar evolution models for PISNe formation, computed from the binary population and spectral synthesis code suite, bpass, with the star formation history of all individual computational elements in the Illustris-TNG simulation. This allows us to compute comic PISN rates and predict their host galaxy properties. Of particular importance is that IllustrisTNG galaxies do not have uniform metallicities throughout, with metal-enriched galaxies often harbouring metal-poor pockets of gas where PISN progenitors may form. Accounting for the chemical inhomogeneities within these galaxies, we find that the peak redshift of PISNe formation is $z=3.5$ instead of the value of $z=6$ when ignoring chemical inhomogeneities within galaxies. Furthermore, the rate increases by an order of magnitude from 1.9 to 29 PISN Gpc$^{-3}$ yr$^{-1}$ at $z=0$, if the chemical inhomogeneities are considered. Using state-of-the-art theoretical PISN light curves, we find an observed rate of 13.8 (1.2) visible PISNe per year for the Euclid-Deep survey, or 83 (7.3) over the 6-yr lifetime of the mission when considering chemically inhomogeneous (homogenous) systems. Interestingly, only 12 per cent of helium PISN progenitors are sufficiently massive to power a superluminous supernova event, which can potentially explain why PISN identification in time-domain surveys remains elusive and progress requires dedicated strategies.

The impact of asteroseismically calibrated internal mixing on nucleosynthetic wind yields of massive stars

Context. Asteroseismology gives us the opportunity to look inside stars and determine their internal properties, such as the radius and mass of the convective core. Based on these observations, estimations can be made for the amount of the convective boundary mixing and envelope mixing of such stars and for the shape of the mixing profile in the envelope. However, these results are not typically included in stellar evolution models. Aims. We aim to investigate the impact of varying convective boundary mixing and envelope mixing in a range based on asteroseismic modelling in stellar models up to the core collapse, both for the stellar structure and for the nucleosynthetic yields. In this first study, we focus on the pre-explosive evolution and we evolved the models to the final phases of carbon burning. This set of models is the first to implement envelope mixing based on internal gravity waves for the entire evolution of the star. Methods. We used the MESA stellar evolution code to simulate stellar models with an initial mass of 20 M⊙ from zero-age main sequence up to a central core temperature of 109 K, which corresponds to the final phases of carbon burning. We varied the convective boundary mixing, implemented as ‘step-overshoot’, with the overshoot parameter (αov) in the range 0.05−0.4. We varied the amount of envelope mixing (log(Denv/cm2 s−1)) in the range 0−6 with a mixing profile based on internal gravity waves. To study the nucleosynthesis taking place in these stars in great detail, we used a large nuclear network of 212 isotopes from 1H to 66Zn. Results. Enhanced mixing according to the asteroseismology of main-sequence stars, both at the convective core boundary and in the envelope, has significant effects on the nucleosynthetic wind yields. This is especially the case for 36Cl and 41Ca, whose wind yields increase by ten orders of magnitude compared to those of the models without enhance envelope mixing. Our evolutionary models beyond the main sequence diverge in yields from models based on rotational mixing, having longer helium-burning lifetimes and lighter helium-depleted cores. Conclusions. We find that the asteroseismic ranges of internal mixing calibrated from core hydrogen-burning stars lead to similar wind yields as those resulting from the theory of rotational mixing. Adopting the seismic mixing levels beyond the main sequence, we find earlier transitions to radiative carbon burning compared to models based on rotational mixing because they have lower envelope mixing in that phase. This influences the compactness and the occurrence of shell mergers, which may affect the supernova properties and explosive nucleosynthesis.

Benchmarking the effective temperature scale of red giant branch stellar models: The case of the metal-poor halo giant HD 122563

Context. There is plenty of evidence in the literature of significant discrepancies between the observations and models of metal-poor red giant branch stars, in particular regarding the effective temperature Teff scale. Aims. We revisit the benchmark star HD 122563 using the most recent observations from Gaia Data Release 3, to investigate if these new constraints may help in resolving this discrepancy. Methods. We review the most recent spectroscopic determinations of the metallicity [Fe/H] of HD 122563, and provide a new assessment of its fundamental parameters, specifically, bolometric luminosity, Teff, surface gravity, plus a photometric determination of its metal content. Using these constraints, we compare the position of the star in the Hertzsprung-Russell (H–R) diagram with various recent sets of stellar evolution tracks. Results. The H-R diagram analysis reveals a significant disagreement between observed and theoretical Teff values, when adopting the most recent spectroscopic estimate of [Fe/H]. On the other hand, by using the photometric determination of [Fe/H], some of the selected sets of stellar tracks appear in fair agreement with observations. The sets with discrepant Teff can be made to agree with observations either by modifying the prescription adopted to calculate the models’ outer boundary conditions, and/or by reducing the adopted value of the mixing length parameter with respect to the solar-calibration. Conclusions. A definitive assessment of whether the Teff scale of metal-poor stellar red giant branch models is consistent with observations requires a more robust determination of the fundamental parameters of HD 122563 and also a larger sample of calibrators. From the theoretical side, it is crucial to minimise the current uncertainties in the treatment (boundary conditions, temperature gradient) of the outer layers of stellar models with convective envelopes.

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.

Physical properties of strong 1 &lt; z &lt; 3 Balmer and Paschen line emitters observed with JWST

Context. The ultraviolet continuum traces young stars while the near-infrared unveils older stellar populations and dust-obscured regions. Balmer emission lines provide insights into gas properties and young stellar objects but are highly affected by dust attenuation. The near-infrared Paschen lines suffer less dust attenuation and can be used to measure star formation rates (SFRs) in star-forming regions obscured by dust clouds. Aims. We present a new way of combining spectro-photometric data in order to test the robustness of the SFRs and stellar mass estimates of star-forming sources observed with JWST. We also aim to quantify the amount of differential attenuation between the interstellar medium and the birth clouds with the use of Paschen emission lines. Methods. We select 13 sources between redshifts 1 and 3 observed with HST, JWST/NIRCam and NIRSpec based on the availability of at least one Balmer and one Paschen line with S/N ≥ 5. With a newly developed version of CIGALE, we fit their hydrogen line equivalent widths (EWs) and photometric data. We assess the impacts of the removal of spectroscopic data by comparing the quality of the fits of the spectro-photometric data to those with photometric data only. We compare the single (BC03) vs binary (BPASS) stellar population models in the fitting process of spectro-photometric data. We derive the differential attenuation and explore different attenuation recipes by fitting spectro-photometric data with BC03. For each stellar model and for each input dataset (with and without EWs), we quantify the deviation on the SFRs and stellar masses from the “standard” choice. Results. The combination of spectro-photometric data provides robust constraints on the physical properties of galaxies, with a significant reduction in the uncertainties compared to using only photometric data. On average, the SFRs are overestimated and the stellar masses are underestimated when EWs are not included as input data. We find a major contribution of the Hα emission line to the broadband photometric measurements of our sources, and a trend of increasing contribution with specific SFR. Using the BPASS models has a significant impact on the derived SFRs and stellar masses, with SFRs being higher by an average of 0.13 dex and stellar masses being lower by an average of 0.18 dex compared to BC03. We show that a flexible attenuation recipe provides more accurate estimates of the dust attenuation parameters, especially the differential attenuation. Finally, we reconstruct the total effective attenuation curves of the most dust-obscured galaxies in our sample.

Pop-cosmos: A Comprehensive Picture of the Galaxy Population from COSMOS Data

We present pop-cosmos: a comprehensive model characterizing the galaxy population, calibrated to 140,938 (r < 25 selected) galaxies from the Cosmic Evolution Survey (COSMOS) with photometry in 26 bands from the ultraviolet to the infrared. We construct a detailed forward model for the COSMOS data, comprising: a population model describing the joint distribution of galaxy characteristics and its evolution (parameterized by a flexible score-based diffusion model); a state-of-the-art stellar population synthesis model connecting galaxies’ intrinsic properties to their photometry; and a data model for the observation, calibration, and selection processes. By minimizing the optimal transport distance between synthetic and real data, we are able to jointly fit the population and data models, leading to robustly calibrated population-level inferences that account for parameter degeneracies, photometric noise and calibration, and selection. We present a number of key predictions from our model of interest for cosmology and galaxy evolution, including the mass function and redshift distribution; the mass–metallicity-redshift and fundamental metallicity relations; the star-forming sequence; the relation between dust attenuation and stellar mass, star formation rate, and attenuation-law index; and the relation between gas-ionization and star formation. Our model encodes a comprehensive picture of galaxy evolution that faithfully predicts galaxy colors across a broad redshift (z < 4) and wavelength range.

The true nature of HE 0057-5959, the most metal-poor, Li-rich star

The Li-rich stars are a class of rare objects with a surface lithium abundance, A(Li), that exceeds that of other stars in the same evolutionary stage. The origin of these stars is still debated, and valuable routes to look at include the Cameron-Fowler mechanism, a mass-transfer process in a binary system, or the engulfment of rocky planets or brown dwarfs. Metal-poor ([Fe/H]&lt;−1 dex) stars are only a small fraction of the entire population of Li-rich stars. We observed the metal-poor ([Fe/H]=−3.95±0.11 dex) giant star HE 0057–5959 with MIKE at the Magellan Telescope, deriving A(Li)NLTE=+2.09±0.07 dex. Such an Li abundance is significantly higher, by about 1 dex, than that of other stars at the same evolutionary stage. A previous analysis of the same target suggested that its high A(Li) reflects an ongoing first-dredge-up process. We revised the nature of HE 0057-5959 by comparing its stellar parameters and A(Li) with appropriate stellar evolution models describing Li depletion due to the deepening of the convective envelope. This comparison rules out that HE 0057-5959 is caught during its first dredge-up, the latter having already ended according to the parameters of this star. Its A(Li), remarkably higher than the typical lithium plateau drawn by similar giant stars, demonstrates that HE 0057-5959 joins the class of the rare metal-poor, Li-rich stars. HE 0057-5959 is the most metal-poor, Li-rich star discovered so far. We considered different scenarios to explain this star also comparing it with the other metal-poor, Li-rich stars. No internal mixing able to activate the Cameron-Fowler mechanism is known for metal-poor stars at this evolutionary stage. The engulfment of planets is also disfavoured because such metal-poor stars should not host planets. Finally, HE 0057-5959 is one of the most Na-rich among the Li-rich stars, and we found that a strong excess of Na abundance is common to all three Li-rich stars with [Fe/H]&lt;–3 dex. This finding could support the scenario of mass transfer from a massive companion star (able to simultaneously produce large amounts of both elements) in a binary system, even if we found no evidence of radial velocity variations.

Primordial Lithium from Globular Cluster Turnoff Stars: M13 and M71

During Big Bang nucleosynthesis (BBN) in the first 15 minutes of the Universe, some 7Li was created along with isotopes of H and He. The determination of that primordial value of Li can help constrain the conditions at that time. The oldest stars with known ages can be found in globular clusters which have well-determined ages through stellar evolution models. High-resolution spectra of Li have been obtained with the Keck I Telescope and HIRES in several unevolved stars in the clusters M13 and M71 with V magnitudes of 17.6–17.9. Abundances of Li have been determined with spectrum synthesis techniques and show a range of a factor of 4. We attribute that spread to differences in initial angular momentum resulting in different amounts of spin-down, related mixing, and destruction of Li. Our results are compared with similar results for main-sequence and turnoff stars in other globular clusters. The range in age of these clusters is 11.2–14.2 Gyr for an age span of 3 Gyr. These clusters range in [Fe/H] from −0.75 to −2.24 corresponding to a factor of 30 in metallicity. The maximum in the Li abundance for these unevolved stars in all eight clusters is the same corresponding to Li/H = 3.16 × 10−10, while the predicted Li abundance, based on the deuterium abundance and the BBN predictions, is 5.24 × 10−10.

ISHELL K-band Survey of Class I and Flat Spectrum Sources: Magnetic Field Measurements in the Protostellar Phase

We perform the first magnetic field strength survey of Class I and Flat Spectrum (FS) sources using K-band observations with iSHELL. We obtained new observations of 42 Class I and FS sources and additionally included 10 sources from the archive. We detect photospheric lines in 44 of the sources, in addition to Brγ, H2, and CO emission in several objects. We model the photospheric absorption lines of 32 Class I and FS sources and measure their magnetic field strengths, K-band temperatures, gravities, projected rotational velocities, and infrared veiling values. We put the physical properties of Class I and FS sources in context by comparing them to the values derived for a sample of Class II sources. We find that (a) the average magnetic field strength of Class I and FS sources 〈B〉 = 2.0 ± 0.15 kG is consistent with the average magnetic field strength of Class II sources 〈B〉 = 1.8 ± 0.15 kG, and (b) the average gravity of Class I and FS objects logg=3.43±0.07 is lower than the average gravity of Class II sources logg=3.75±0.04 , although there is significant overlap between both gravity distributions. Furthermore, using stellar evolutionary models, we deduce that Class I and FS sources have a median age of ∼0.6 Myr, and are, as a group, younger than the Class II stars with a median age of ∼3 Myr. Our results confirm that Class I and FS sources host strong magnetic fields on their photospheres. Thus, it is likely that these sources accrete disk material through a magnetosphere similar to the more evolved T Tauri stars.

Modeling Compact Object Mergers GW190814 and GW200210 and Other Self-bound Compact Stars with Dark Matter Induced by Gravitational Decoupling and Its Significance to the Mass Gap

We present a rigorous study of compact objects within f(Q) gravity where electrical fields and dark matter are studied and provide novel mass–radius relations with models falling in the mass gap of the events GW190814 and GW200210. After formulating the basic equations and finding their relevant solutions, we impose the boundary conditions on the system under treatment. The decoupled solution for the strange stellar model with the dark matter density profile is obtained. The distribution patterns of the effective energy density and the radial as well as the tangential pressure and anisotropy in the system are intensively examined. The stability properties of the stellar configuration and the influence of dark matter are studied. The recent observations of supermassive compact star candidates such as PSR J1614–2230 and PSR J0952–0607 with observed masses greater than or equal to 2 M ⊙ have been employed in our study. Interestingly, the present study predicts the constraints on mass–radius measurements of the observed stars satisfying the equation of state based on the MIT bag model and employing the condition of mimicking, i.e., ρ θ = χ1ρPI in f(Q) gravity. Our graphical results exhibit that for a particular M–R curve with fixed values of the parameters, neutron stars having mass less than M max exist with larger radii within the context of the f(Q) formalism.

Element abundances of galactic RGB stars in the APO-K2 catalogue

Aims. We conducted an investigation on the chemical abundances of 4316 stars in the red giant branch (RGB) phase from the recently released APO-K2 catalogue. Our aim was to characterize the abundance trends of the single elements with [α/Fe], mainly focusing on C, N, and O, which are the most relevant for the estimation of stellar ages. Methods. The chemical analysis of the RGB sample involved cross-matching data from the APO-K2 catalogue with individual element abundances from APOGEE DR17. Results. The analysis detected a statistically significant difference in the [(C+N+O)/Fe]–[α/Fe] trend with respect to the simple α-enhancement scenario. This difference remained robust across different choices for the reference solar mixture and potential zero-point calibrations of C and N abundances. The primary discrepancy was a steeper increase in [O/Fe] with [α/Fe], reaching a 0.1 dex difference at [α/Fe] = 0.3. Notably, the impact on the evolutionary timescale of such oxygen over-abundance with respect to the commonly adopted uniform α-enhancement is rather limited. We verified that stellar models computed using an ad hoc O-rich mixture sped up the evolution by only 1% at [α/Fe] = 0.3, due to the counterbalancing effects of O enrichment on both the evolutionary timescale and the Z-to-[Fe/H] relationship.

Imaging of I Zw 18 by JWST

Context. The blue compact dwarf (BCD) galaxy I Zw 18 is one of the most metal-poor (Z ∼ 3% Z⊙) star-forming galaxies known in the local Universe. Since its discovery, the evolutionary status of this system has been at the center of numerous debates within the astronomical community. Aims. We aim to probe and resolve the stellar populations of I Zw 18 in the near-IR using JWST/NIRCam’s unprecedented imaging spatial resolution and sensitivity. Additionally, our goal is to derive the spatially resolved star formation history (SFH) of the galaxy within the last 1 Gyr, and to provide constraints for older epochs. Methods. We used the point spread function fitting photometry package DOLPHOT to measure positions and fluxes of point sources in the F115W and F200W filters’ images of I Zw 18, acquired as part of the JWST GTO ID 1233 (PI: Meixner). Furthermore, to derive I Zw 18’s SFH, we applied a state-of-the-art color-magnitude diagram (CMD) fitting technique (SFERA 2.0), using two independent sets of stellar models: PARSEC-COLIBRI and MIST. Results. Our analysis of I Zw 18’s CMD reveal three main stellar populations: one younger than ∼30 Myr, mainly associated with the northwest star-forming (SF) region; an intermediate-age population (∼100 − 800 Myr), associated with the southeast SF region; and a red and faint population, linked to the underlying halo of the galaxy, older than 1 Gyr and possibly as old as 13.8 Gyr. The main body of the galaxy shows a very low star formation rate (SFR) of ~10-4M⊙ yr−1 between 1 and 13.8 Gyr ago. In the last billion years, I Zw 18 shows an increasing trend, culminating in two strong bursts of SF around ∼10 and ∼100 Myr ago. Notably, I Zw 18 Component C mimics the evolution of the main body, but with lower SFRs on average. Conclusions. Our results confirm that I Zw 18 is populated by stars of all ages, without any major gaps. Thus, I Zw 18 is not a truly young galaxy, but rather a system characterized by an old underlying stellar halo, in agreement with what has been found in other BCDs by similar studies. The low SF activity exhibited at epochs older than 1 Gyr is in agreement with the “slow cooking" dwarf scenario proposed in the literature, and could have contributed to its low metal content. The galaxy is now experiencing its strongest episode of star formation (∼ 0.6 M⊙ yr−1) mainly located in the northwest region. A recent gravitational interaction between the main body and Component C is the most likely explanation for this starburst.

The cosmic rate of pair-instability supernovae

Abstract Pair-instability supernovae (PISNe) have crucial implications for many astrophysical topics, including the search for very massive stars, the black hole mass spectrum, and galaxy chemical enrichment. To this end, we need to understand where PISNe are across cosmic time, and what are their favourable galactic environments. We present a new determination of the PISN rate as a function of redshift, obtained by combining up-to-date stellar evolution tracks from the PARSEC and FRANEC codes, with an up-to-date semi-empirical determination of the star formation rate and metallicity evolution of star-forming galaxies throughout cosmic history. We find the PISN rate to exhibit a huge dependence on the model assumptions, including the criterion to identify stars unstable to pair production, and the upper limit of the stellar initial mass function. Remarkably, the interplay between the maximum metallicity at which stars explode as PISNe, and the dispersion of the galaxy metallicity distribution, dominates the uncertainties, causing a ∼ seven-orders-of-magnitude PISN rate range. Furthermore, we show a comparison with the core-collapse supernova rate, and study the properties of the favourable PISN host galaxies. According to our results, the main contribution to the PISN rate comes from metallicities between ∼10−3 and 10−2, against the common assumption that views very-low-metallicity, Population III stars as exclusive or dominant PISN progenitors. The strong dependencies we find offer the opportunity to constrain stellar and galaxy evolution models based on possible future (or the lack of) PISN observations.

A Study of the Abundance of Aluminium in Massive Stars

Abstract: This study presents the results of the calculation of the mass fractions of isotopes from hydrogen to nickel during the entire core burning phases of a 25 M⊙ star. A simple stellar model was used, considering four main parameters: temperature, density, initial composition, and time, without accounting for mass loss or mixing. The mass fractions of the isotopes were calculated using the open-source package NucNet Tools and the updated reaction rates from the JINA Reaclib database. A comparison of our results with existing similar data is performed and acceptable agreement is conspicuous. Hydrodynamic conditions in massive stars favor the production of 26Al, therefore particular attention was given to the isotopes 26Al and 27Al and the ratio of their mass fractions R_Al for comparison with the literature. The averaged value of the controversial R_Al during the star's lifetime is found to be 1.68×10-4. Keywords: Elemental abundances in stars, Evolution, Stellar, Nuclear astrophysics, Hydrostatic stellar nucleosynthesis, Nucleosynthesis: stellar, Stars formation. PACS: 97.10.Tk, 87.23.-n, 97.10.Cv, 26.20.-f, 7.10.Cv, 97.10.Bt.

Two Young Eclipsing Binaries in Orion with Temperatures and Radii Affected by Spots and Third Bodies

In this work we present a model of two young eclipsing binaries in the Orion complex. Both heavily spotted, they present radii and temperatures that are in disagreement with the predictions of standard stellar models. 2M05−06 consists of two stars with different masses (∼0.52 and ∼0.42 M ⊙) but with very similar radii (∼0.9 R ⊙), and with the less massive star having a highly spotted surface that causes it to have a hotter (unspotted) photosphere than the higher-mass star. The other system, 2M05−00, consists of two stars of very similar masses (∼0.34 M ⊙), but very different radii (∼0.7 and ∼1.0 R ⊙), which creates an appearance of the two eclipsing stars being noncoeval. 2M05−00 appears to have a tertiary companion that could offer an explanation for the unusual properties of the eclipsing stars, as has been seen in some other young triple systems. Comparing the empirically measured properties of these eclipsing binaries to the predictions of stellar models, both standard and magnetic, we find that only the magnetic models correctly predict the observed relationship between mass and effective temperature. However, standard (nonmagnetic) models better predict the temperatures of the unspotted photospheres. These observations represent an important step in improving our understanding of pre-main-sequence stellar evolution and the roles of spots and tertiaries on fundamental stellar properties.