Main Sequence Research Articles

Binarity at LOw Metallicity (BLOeM)

Surveys in the Milky Way and Large Magellanic Cloud have revealed that the majority of massive stars will interact with companions during their lives. However, knowledge of the binary properties of massive stars at low metallicity, and therefore in conditions approaching those of the Early Universe, remain sparse. We present the Binarity at LOw Metallicity (BLOeM) campaign, an ESO large programme designed to obtain 25 epochs of spectroscopy for 929 massive stars in the Small Magellanic Cloud, allowing us to probe multiplicity in the lowest-metallicity conditions to date (Z = 0.2 Z⊙). BLOeM will provide (i) the binary fraction, (ii) the orbital configurations of systems with periods of P ≲ 3 yr, (iii) dormant black-hole binary candidates (OB+BH), and (iv) a legacy database of physical parameters of massive stars at low metallicity. Main sequence (OB-type) and evolved (OBAF-type) massive stars are observed with the LR02 setup of the GIRAFFE instrument of the Very Large Telescope (3960–4570 Å resolving power R = 6200; typical signal-to-noise ratio(S/N) ≈70–100). This paper utilises the first nine epochs obtained over a three-month time period. We describe the survey and data reduction, perform a spectral classification of the stacked spectra, and construct a Hertzsprung-Russell diagram of the sample via spectral-type and photometric calibrations. Our detailed classification reveals that the sample covers spectral types from O4 to F5, spanning the effective temperature and luminosity ranges 6.5 ≲ Teff/kK ≲ 45 and 3.7 < log L/L⊙ < 6.1 and initial masses of 8 ≲ Mini ≲ 80 M⊙. The sample comprises 159 O-type stars, 331 early B-type (B0–3) dwarfs and giants (luminosity classes V–III), 303 early B-type supergiants (II–I), and 136 late-type BAF supergiants. At least 82 stars are OBe stars: 20 O-type and 62 B-type (13% and 11% of the respective samples). In addition, the sample includes 4 high-mass X-ray binaries, 3 stars resembling luminous blue variables, 2 bloated stripped-star candidates, 2 candidate magnetic stars, and 74 eclipsing binaries.

JWST imaging of the closest globular clusters

We present observations of the two closest globular clusters, NGC 6121 and NGC 6397, taken with the NIRISS detector of JWST. The combination of our new JWST data with archival Hubble Space Telescope (HST) images allows us to compute proper motions, disentangle cluster members from field objects, and probe the main sequence (MS) of the clusters down to <0.1 M⊙ as well as the brighter part of the white-dwarf sequence. We show that theoretical isochrones fall short in modeling the low-mass MS and discuss possible explanations for the observed discrepancies. Our analysis suggests that the lowest-mass members of both clusters are significantly more metal-rich and oxygen-poor than their higher-mass counterparts. It is unclear whether the difference is caused by a genuine mass-dependent chemical heterogeneity, low-temperature atmospheric processes altering the observed abundances, or systematic shortcomings in the models. We computed the present-day local luminosity and mass functions of the two clusters; our data reveal a strong flattening of the mass function indicative of a significant preferential loss of low-mass stars in agreement with previous dynamical models for these two clusters. We have made our NIRISS astro-photometric catalogs and stacked images publicly available to the community.

Observations of Extremely Metal-poor O Stars: Weak Winds and Constraints for Evolution Models

Metal-poor massive stars drive the evolution of low-mass galaxies, both locally and at high redshift. However, quantifying the feedback they impart to their local surroundings remains uncertain because models of stellar evolution, mass loss, and ionizing spectra are unconstrained by observations below 20% solar metallicity (Z ⊙). We present new Keck Cosmic Web Imager optical spectroscopy of three O-type stars in the nearby dwarf galaxies Leo P, Sextans A, and WLM, which have gas-phase oxygen abundances of 3%–14% Z ⊙. To characterize their fundamental stellar properties and radiation-driven winds, we fit PoWR atmosphere models to the optical spectra simultaneously with Hubble Space Telescope far-ultraviolet (FUV) spectra and multiwavelength photometry. We find that all three stars have effective temperatures consistent with their spectral types and surface gravities typical of main-sequence dwarf stars. Yet, the combination of those inferred parameters and luminosity for the two lower-Z stars is not reproduced by stellar evolution models, even those that include rotation or binary interactions. The scenario of multiple-star systems is difficult to reconcile with all available data, suggesting that these observations pose a challenge to current evolution models. We highlight the importance of validating the relationship between stellar mass, temperature, and luminosity at very low Z for accurate estimates of ionizing photon production and spectral hardness. Finally, all three stars’ FUV wind profiles reveal low mass-loss rates and terminal wind velocities in tension with expectations from widely adopted radiation-driven wind models. These results provide empirical benchmarks for future development of mass-loss and evolution models for metal-poor stellar populations.

Detailed Abundances of the Planet-hosting TOI-1173 A/B System: Possible Evidence of Planet Engulfment in a Very Wide Binary

Over the last decade, studies of large samples of binary systems have identified chemical anomalies and shown that they might be attributed to planet formation or planet engulfment. However, both scenarios have primarily been tested in pairs without known exoplanets. In this work, we explore these scenarios in the newly detected planet-hosting wide binary TOI-1173 A/B (projected separation ∼11,400 au), using high-resolution MAROON-X and ARCES spectra. We determined photospheric stellar parameters both by fitting stellar models and via the spectroscopic equilibrium approach. Both analyses agree and suggest that they are cool main-sequence stars located in the thin disk. A line-by-line differential analysis between the components (B−A) displays an abundance pattern in the condensation temperature plane, where the planet-hosting star TOI-1173 A is enhanced in refractory elements such as iron by more than 0.05 dex. This suggests the engulfment of ∼18 M ⊕ of rocky material in star A. Our hypothesis is supported by the dynamics of the system (detailed in our companion paper), which suggest that the super-Neptune TOI-1173 A b might have been delivered to its current short period (∼7 days) through circularization and von Zeipel–Lidov–Kozai mechanisms, thereby triggering the engulfment of inner rocky exoplanets.

Characterizing the Rapid Hydrogen Disappearance in SN 2022crv: Evidence of a Continuum between Type Ib and IIb Supernova Properties

We present optical and near-infrared (NIR) observations of SN 2022crv, a stripped-envelope supernova in NGC 3054, discovered within 12 hr of explosion by the Distance Less Than 40 Mpc Survey. We suggest that SN 2022crv is a transitional object on the continuum between Type Ib supernovae (SNe Ib) and Type IIb supernovae (SNe IIb). A high-velocity hydrogen feature (∼ −20,000 to −16,000 km s−1) was conspicuous in SN 2022crv at early phases, and then quickly disappeared. We find that a hydrogen envelope of ∼10−3 M ⊙ can reproduce the observed behavior of the hydrogen feature. The lack of early envelope cooling emission implies that SN 2022crv had a compact progenitor with an extremely low amount of hydrogen. A nebular spectral analysis shows that SN 2022crv is consistent with the explosion of a He star with a final mass of ∼4.5–5.6 M ⊙ that evolved from a ∼16 to 22 M ⊙ zero-age main-sequence star in a binary system with ∼1.0–1.7 M ⊙ of oxygen finally synthesized in the core. In order to retain such a small amount of hydrogen, the initial orbital separation of the binary system is likely larger than ∼1000 R ⊙. The NIR spectra of SN 2022crv show a unique absorption feature on the blue side of the He i line at ∼1.005 μm. This is the first time such a feature has been observed in SNe Ib/IIb, and it could be due to Sr II. Further detailed modeling of SN 2022crv can shed light on the progenitor and the origin of the mysterious absorption feature in the NIR.

The YMCA (Yes, Magellanic Clouds Again) survey: Probing the outer regions of the Magellanic system with VST

Context. The Magellanic Clouds (MCs) are the Milky Way’s most massive dwarf satellites. As they also represent the closest pair of galaxies in an ongoing tidal interaction while simultaneously infalling into the Milky Way halo, they provide a unique opportunity to study in detail an ongoing three-body encounter. Aims. We present the YMCA (Yes, Magellanic Clouds Again) survey: Probing the outer regions of the Magellanic system with VST, based on deep optical photometry carried out with the VLT Survey Telescope (VST). Methods. The YMCA survey targeted 110 square degrees, in the g and i filters, in the periphery of both MCs, including a long strip in between the Large Magellanic Cloud (LMC) and the Small Magellanic Cloud (SMC). The photometry of YMCA is sufficiently deep (50% complete down to g ≃ 23.5 − 24.0 mag) to allow for a detailed analysis of main-sequence stars in regions of the MCs that have remained relatively unexplored at these faint magnitudes. Results. The resulting colour–magnitude diagrams reveal that the outskirts of the MCs are predominantly characterised by intermediate-age and old stellar populations, with limited or negligible evidence of recent star formation. The analysis of the age distribution of star clusters (SCs) within the surveyed area, both already known and newly discovered candidates, hints at a close fly-by between the LMC and SMC that occurred ≃2.5 − 3.0 Gyr ago, in agreement with previous results. We also report the discovery of candidate SCs with ages within the so-called age-gap, thus questioning its real existence.

Active Galactic Nuclei in the Green Valley at z ∼ 0.7

We present near-infrared (NIR) spectroscopy using the MMT and Magellan infrared spectrograph for a sample of 29 massive galaxies ( logM*/M⊙≳10 ) at z ∼ 0.7 with optical spectroscopy from the Large Early Galaxy Astrophysics Census (LEGA-C) survey. Having both optical and NIR spectroscopy at this redshift allows us to measure the full suite of rest-optical strong emission lines, enabling the study of ionization sources and the rest-optical selection of active galactic nuclei (AGN), as well as the measurement of dust-corrected Hα-based star formation rates (SFRs). We find that 11 out of 29 galaxies host AGN. We infer the nonparametric star formation histories with the spectral energy distribution fitting code Prospector and classify galaxies as star-forming, green valley, or quiescent based on their most recent specific SFRs. We explore the connection between AGN activity and suppressed star formation, and find that 89% ± 15% of galaxies in the green valley or below host AGN, while only 15% ± 8% of galaxies above the green valley host AGN. We construct the star-forming main sequence (SFMS) and find that the AGN host galaxies are 0.37 dex below the SFMS, while galaxies without detectable AGN are consistent with being on the SFMS. However, when compared to a bootstrapped mass-matched sample, the SFRs of our sample of AGN host galaxies are consistent with the full LEGA-C sample. Based on this mass-matched analysis, we cannot rule out that this suppression of star formation is driven by other processes associated with the higher mass of the AGN sample. We therefore cannot link the presence of AGN activity to the quenching of star formation.

TIC 290061484: A Triply Eclipsing Triple System with the Shortest Known Outer Period of 24.5 Days

We have discovered a triply eclipsing triple-star system, TIC 290061484, with the shortest known outer period, P out, of only 24.5 days. This “eclipses” the previous record set by λ Tauri at 33.02 days, which held for 68 yr. The inner binary, with an orbital period of P in = 1.8 days, produces primary and secondary eclipses and exhibits prominent eclipse timing variations with the same periodicity as the outer orbit. The tertiary star eclipses, and is eclipsed by, the inner binary with pronounced asymmetric profiles. The inclinations of both orbits evolve on observable timescales such that the third-body eclipses exhibit dramatic depth variations in TESS data. A photodynamical model provides a complete solution for all orbital and physical parameters of the triple system, showing that the three stars have masses of 6.85, 6.11, and 7.90 M ⊙, radii near those corresponding to the main sequence, and T eff in the range of 21,000–23,700 K. Remarkably, the model shows that the triple is in fact a subsystem of a hierarchical 2+1+1 quadruple with a distant fourth star. The outermost star has a period of ∼3200 days and a mass comparable to the stars in the inner triple. In ∼20 Myr, all three components of the triple subsystem will merge, undergo a Type II supernova explosion, and leave a single remnant neutron star. At the time of writing, TIC 290061484 is the most compact triple system and one of the tighter known compact triples (i.e., P out/P in = 13.7).

Variable Stars in M31 Stellar Clusters from the Panchromatic Hubble Andromeda Treasury

Variable stars in stellar clusters can offer key constraints on stellar evolution and pulsation models, utilizing estimates of host cluster properties to constrain stellar physical parameters. We present a catalog of 86 luminous (F814W < 19) variable stars in M31 clusters identified by mining the archival Panchromatic Hubble Andromeda Treasury (PHAT) survey using a combination of statistical analysis of sparse PHAT light curves and difference imaging. We determine the evolutionary phases and initial masses of these variable stars by matching them with theoretical isochrones generated using host cluster properties from the literature. We calculate the probability of PHAT photometry being blended due to the highly crowded nature of cluster environments for each cluster-variable star, using these probabilities to inform our level of confidence in the derived properties of each star. Our 86 cluster-variable stars have initial masses between 0.8 and 67 M ⊙. Their evolutionary phases span the main sequence, more evolved hydrogen- and helium-burning phases, and the post–asymptotic giant branch. We identify numerous candidate variable star types: RV Tauri variables, red supergiants, and slowly pulsating B-type supergiants, along with Wolf–Rayet stars, α Cygni and Mira variables, a classical Cepheid, and a possible superasymptotic giant. We characterize 12 cluster-variable stars at higher confidence based on their difference image quality and lower blending probability. Ours is the first systematic study of variable stars in extragalactic stellar clusters leveraging the superior resolution of the Hubble Space Telescope and demonstrating the unique power of stellar clusters in constraining the fundamental properties of variable stars.

Lyman Continuum Leakers at z > 3 in the GOODS-S Field: Starburst or Not?

We investigate the star-forming properties of 23 Lyman continuum (LyC) leakers at z > 3 in the Great Observatories Origins Deep Survey-South field based on a systematic review of LyC observations from the literature. Using data from the Hubble Space Telescope and the James Webb Space Telescope, we construct the spectral energy distributions for these LyC leakers, covering the spectrum from rest-frame ultraviolet to near-infrared. Through the application of a unified modeling approach, we measure the ultraviolet slope (β), star formation rate, and stellar mass for these LyC leakers in a consistent manner. These high-redshift LyC leakers demonstrate statistically blue UV-continuum slopes (β), which is consistent with their high escape fraction of LyC photons. We find that these high-redshift LyC leakers span a wide range of specific star formation rate ( log(sSFR /yr) from −8.6 to −6.7). Ten of these LyC leakers are located on the star formation main sequence, instead of all being in the starburst mode. The results indicate that intense bursts of star formation are not necessarily required for the leakage of LyC photons for galaxies at z > 3.

Giant Eruptions in Massive Stars and their Effect on the Stellar Structure

Giant eruptions (GEs) in luminous blue variables are years-to-decades-long episodes of enhanced mass loss from the outer layers of the star during which the star undergoes major changes in its physical and observed properties. We use the Modules for Experiments in Stellar Astrophysics stellar evolution code to model the evolution of a 70 M ⊙ star that undergoes a GE. We let the star evolve to the termination of the main sequence, and when it reaches T ≃ 19,400 K we emulate a GE by removing mass from its outer layers at a rate of 0.15 M ⊙ yr−1 for 20 yr. As mass is being lost, the star contracts and releases a substantial amount of gravitational energy. The star undergoes an initial ≃3 days of expansion followed by years of contraction. During that time the star tries to reach an equilibrium state, and as a result of loss in gravitational energy, its luminosity drops about 1 order of magnitude. As the GE terminates, we let the star continue to evolve without any further mass loss and track its recovery as it regains its equilibrium by adjusting its internal structure. After ≃87 yr it reaches a state very close to the one where the GE was first initiated. We suggest that at this point another GE or a cycle of GEs may occur.

JADES + JEMS: A Detailed Look at the Buildup of Central Stellar Cores and Suppression of Star Formation in Galaxies at Redshifts 3 < z < 4.5

We present a spatially resolved study of stellar populations in six galaxies with stellar masses M * ∼ 1010 M ☉ at z ∼ 3.7 using 14-filter James Webb Space Telescope (JWST)/NIRCam imaging from the JADES and JEMS surveys. The six galaxies are visually selected to have clumpy substructures with distinct colors over rest frame 3600−4100 Å, including a red, dominant stellar core that is close to their stellar-light centroids. With 23-filter photometry from the Hubble Space Telescope to JWST, we measure the stellar-population properties of individual structural components via spectral energy distribution fitting using Prospector. We find that the central stellar cores are ≳2 times more massive than the Toomre mass, indicating they may not form via single in situ fragmentation. The stellar cores have stellar ages of 0.4−0.7 Gyr that are similar to the timescale of clump inward migration due to dynamical friction, suggesting that they likely instead formed through the coalescence of giant stellar clumps. While they have not yet quenched, the six galaxies are below the star-forming main sequence by 0.2−0.7 dex. Within each galaxy, we find that the specific star formation rate is lower in the central stellar core, and the stellar-mass surface density of the core is already similar to quenched galaxies of the same masses and redshifts. Meanwhile, the stellar ages of the cores are either comparable to or younger than the extended, smooth parts of the galaxies. Our findings are consistent with model predictions of the gas-rich compaction scenario for the buildup of galaxies’ central regions at high redshifts. We are likely witnessing the coeval formation of dense central cores, along with the onset of galaxy-wide quenching at z > 3.

Searching for new cataclysmic variables in the Chandra Source Catalog

Aims. Cataclysmic variables (CVs) are compact binary systems in which a white dwarf accretes matter from a Roche-lobe-filling companion star. For this study we searched for new CVs in the Milky Way in the Chandra Source Catalog v2.0, cross-matched with Gaia Data Release 3 (DR3). Methods. We identified new CV candidates by combining X-ray and optical data in a color-color diagram called the X-ray main sequence. We used two different cuts in this diagram to compile pure and optically variable samples of CV candidates. We undertook optical spectroscopic follow-up observations with the Keck and Palomar Observatories to confirm the nature of these sources. Results. We assembled a sample of 25 887 Galactic X-ray sources and found 14 new CV candidates. Seven objects show X-ray and/or optical variability. All sources show X-ray luminosity in the 1029 − 1032 erg s−1 range, and their X-ray spectra can be approximated by a power-law model with photon indices in the Γ ∼ 1 − 3 range or an optically thin thermal emission model in the kT ∼ 1 − 70 keV range. We spectroscopically confirmed four CVs, discovering two new polars, one low accretion rate polar and a WZ Sge-like low accretion rate CV. X-ray and optical properties of the other nine objects suggest that they are also CVs (likely magnetic or dwarf novae), and one other object could be an eclipsing binary, but revealing their true nature requires further observations. Conclusions. These results show that a joint X-ray and optical analysis can be a powerful tool for finding new CVs in large X-ray and optical catalogs. X-ray observations such as those by Chandra are particularly efficient at discovering magnetic and low accretion rate CVs, which could be missed by purely optical surveys.

The California Legacy Survey. V. Chromospheric Activity Cycles in Main-sequence Stars

We present optical spectroscopy of 710 solar neighborhood stars collected over 20 years to catalog chromospheric activity and search for stellar activity cycles. The California Legacy Survey stars are amenable to exoplanet detection using precise radial velocities, and we present their Ca ii H and K time series as a proxy for stellar and chromospheric activity. Using the High Resolution Echelle Spectrometer at Keck Observatory, we measured stellar flux in the cores of the Ca ii H and K lines to determine S-values on the Mount Wilson scale and the log(RHK′) metric, which is comparable across a wide range of spectral types. From the 710 stars, with 52,372 observations, 285 stars were sufficiently sampled to search for stellar activity cycles with periods of 2–25 yr, and 138 stars showed stellar cycles of varying length and amplitude. S-values can be used to mitigate stellar activity in the detection and characterization of exoplanets. We used them to probe stellar dynamos and to place the Sun's magnetic activity into context among solar neighborhood stars. Using precise stellar parameters and time-averaged activity measurements, we found tightly constrained cycle periods as a function of stellar temperature between log(RHK′) of −4.7 and −4.9, a range of activity in which nearly every star has a periodic cycle. These observations present the largest sample of spectroscopically determined stellar activity cycles to date.

An empirical isochrone archive for nearby open clusters

Context. The ages of star clusters and co-moving stellar groups contain essential information about the Milky Way. Their special properties and placement throughout the galactic disk make them excellent tracers of galactic structure and key components to unlocking its star formation history. Yet, even though the importance of stellar population ages has been widely recognized, their determination remains a challenging task often associated with highly model-dependent and uncertain results. Aims. We propose a new approach to this long-standing problem, which relies on empirical isochrones of known clusters extracted from high-quality observational data. These purely observation-based data products open up the possibility of relative age determination, free of stellar evolution model assumptions. Methods. For the derivation of the empirical isochrones, we used a combination of the statistical analysis tool principal component analysis for preprocessing and the supervised machine learning method support vector regression for curve extraction. To improve the statistical reliability of our result, we defined the empirical isochrone of a color-magnitude diagram (CMD) of a cluster as the median calculated from a set of nboot = 1000 curves derived from bootstrapped data. The algorithm requires no physical priors, is computationally fast, and can easily be generalized over a large range of CMD combinations and evolutionary stages of clusters. Results. We provide empirical isochrones in all Gaia DR2 and DR3 color combinations for 83 nearby clusters (d &lt; 500 pc), which cover an estimated age range of 7 Myr to 3 Gyr. In doing so, we pave the way for a relative comparison between individual stellar populations based on an age-scaling ladder of empirical isochrones of known clusters. Furthermore, due to the exceptional precision of the available observational data, we report accurate lower main sequence empirical isochrones for many clusters in our sample, which are of special interest as this region is known to be especially complex to model. We validate our method and results by comparing the extracted empirical isochrones to cluster ages in the literature. We also investigate the added information that empirical isochrones covering the lower main sequence can provide on case studies of the IC 4665 cluster and the Meingast 1 stream. Conclusions. The archive of empirical isochrones offers a novel approach to validating age estimates and can be used as an age-scaling ladder or age brackets for new populations and serve as calibration data for further constraining stellar evolution models.

The δ Scuti stars of the Cep–Her Complex – I. Pulsator fraction, rotation, asteroseismic large spacings, and the νmax relation

ABSTRACT We identify delta Scuti ($\delta$ Sct) pulsators amongst members of the recently discovered Cep–Her Complex using light curves from the Transiting Exoplanet Survey Satellite (TESS). We use Gaia colours and magnitudes to isolate a subsample of provisional Cep–Her members that are located in a narrow band on the colour–magnitude diagram compatible with the zero-age main sequence. The $\delta$ Sct pulsator fraction amongst these stars peaks at 100 per cent and we describe a trend of higher pulsator fractions for younger stellar associations. We use four methods to measure the frequency of maximum amplitude or power, $\nu _{\rm max}$, to minimize methodological bias, and we demonstrate their sound performance. The $\nu _{\rm max}$ measurements display a correlation with effective temperature, but with a scatter that is too large for the relation to be useful. We find two ridges in the $\nu _{\rm max}$–$T_{\rm eff}$ diagram, one of which appears to be the result of rapid rotation causing stars to pulsate in low-order modes. We measure the $\nu _{\rm max}$ values of $\delta$ Sct stars in four other clusters or associations of similar age (Trumpler 10, the Pleiades, NGC 2516, and Praesepe) and find similar behaviour with $T_{\rm eff}$. Using échelle diagrams, we measure the asteroseismic large spacing, $\Delta \nu$, for 70 stars, and find a correlation between $\Delta \nu$, rotation, and luminosity that allows rapid rotators seen at low inclinations to be distinguished from slow rotators. We find that rapid rotators are more likely than slow rotators to pulsate, but they do so with less regular pulsation patterns. We also investigate the reliability of Gaia’s vbroad measurement for A-type stars, finding that it is mostly accurate but underestimates $v\sin i$ for slow rotators ($v\sin i \lt 50$ km s−1) by 10–15 per cent.

Gaia BH1 and BH2 — Evolutionary Models with Overshooting of the Black Hole Progenitors within the Present-Day Binary Separation

Abstract Three black holes (BHs) in wide binaries — Gaia BH1, BH2 and BH3 — were recently discovered. The likely progenitors of the BHs were massive stars that experienced a supergiant phase, reaching radii of ∼1000R⊙, before collapsing to form the BH. Such radii are difficult to accommodate with the present-day orbits of BH1 and BH2 — with semi-major axes of 1.4 and 3.7 au, respectively. In this letter, we show that the maximal radii of the supergiants are not necessarily so large, and realistic stellar evolution models, with some assumed overshooting above the convective core into the radiative stellar envelope, produce substantially smaller maximal radii. The limited expansion of supergiants is consistent with the empirical Humphreys-Davidson limit — the absence of red supergiants above an upper luminosity limit, notably lower than the highest luminosity of main-sequence stars. We propose that the evolution that led to the formation of Gaia BH1 and BH2 simply did not involve an expansion to the cool supergiant phase.

The shape of convection in 2D and 3D global simulations of stellar interiors

Theoretical descriptions of convective overshooting in stellar interiors often rely on a basic one-dimensional parameterization of the flow called the filling factor for convection. Several different definitions of the filling factor have been developed for this purpose, based on: (1) the percentage of the volume, (2) the mass flux, and (3) the convective flux that moves through the boundary. We examine these definitions of the filling factor with the goal of establishing their ability to explain differences between 2D and 3D global simulations of stellar interiors that include fully compressible hydrodynamics and realistic microphysics for stars. We study convection and overshooting in pairs of identical two-dimensional (2D) and three-dimensional (3D) global simulations of stars produced with a fully compressible, time-implicit hydrodynamics code. We examine pairs of simulations for (1) a $3$ odot $ red giant star near the first dredge-up point, (2) a $1$ M$_ odot $ pre-main-sequence star with a large convection zone, (3) the current sun, and (4) a $20$ M$_ odot $ main-sequence star with a large convective core. Our calculations of the filling factor based on the volume percentage and the mass flux indicate asymmetrical convection near the surface for each star with an outer convection zone. However, near the convective boundary, convective flows achieve inward-outward symmetry for each star that we study; for 2D and 3D simulations, these filling factors are indistinguishable. A filling factor based on the convective flux is contaminated by boundary-layer-like flows, making a theoretical interpretation difficult. We present two possible new alternatives to these frequently used definitions of a filling factor, which instead compare flows at two different radial points. The first alternative is the penetration parameter of anders2022stellar . The second alternative is a new statistic that we call the plume interaction parameter. We demonstrate that both of these parameters captures systematic differences between 2D and 3D simulations around the convective boundary.

Finding the unusual red giant remnants of cataclysmic variable mergers

Mergers between helium white dwarfs and main-sequence stars are likely common, producing red giant-like remnants making up roughly a few percent of all low-mass ( ≲2M⊙) red giants. Through detailed modeling, we show that these merger remnants possess distinctive photometric, asteroseismic, and surface abundance signatures through which they may be identified. During hydrogen shell burning, merger remnants reach higher luminosities and possess pulsations which depart from the usual degenerate sequence on the asteroseismic Δν– ΔΠ diagram for red giant branch stars. For sufficiently massive helium white dwarfs, merger remnants undergo especially violent helium flashes which can dredge up a large amount of core material (up to ∼0.1M⊙) into the envelope. Such post-dredge-up remnants are more luminous than normal red clump stars, are surface carbon-, helium-, and possibly lithium-rich, and possess a wider range of asteroseismic g-mode period spacings and mixed-mode couplings. Recent asteroseismically determined low-mass ( ≲0.8M⊙) red clump stars may be core helium-burning remnants of mergers involving lower-mass helium white dwarfs.

Toward a Comprehensive Grid of Cepheid Models with MESA. I. Uncertainties of the Evolutionary Tracks of Intermediate-mass Stars

Helium-burning stars, in particular Cepheids, are especially difficult to model, as the choice of free parameters can greatly impact the shape of the blue loops—the part of the evolutionary track at which the instability strip is crossed. Contemporary one-dimensional stellar evolution codes, like Modules for Experiments in Stellar Astrophysics (MESA), come with a large number of free parameters that allow us to model the physical processes in stellar interiors under many assumptions. The uncertainties that arise from this freedom are rarely discussed in the literature despite their impact on the evolution of the model. We calculate a grid of evolutionary models with MESA, varying several controls, like solar mixture of heavy elements, mixing-length theory prescription, nuclear reaction rates, the scheme to determine convective boundaries, atmosphere model, and temporal and spatial resolution, and quantify their impact on age and location of the evolutionary track on the H-R diagram from the main sequence until the end of core helium burning. Our investigation was conducted for a full range of masses and metallicities expected for classical Cepheids. The uncertainties are significant, especially during core helium burning, reaching or exceeding the observational uncertainties of logTeff and logL for detached eclipsing binary systems. For ≥9 M ⊙ models, thin convective shells develop and evolve erratically, not allowing the models to converge. A careful inspection of Kippenhahn diagrams and convergence study is advised for a given mass and metallicity, to assess how severe this problem is and to what extent it may affect the evolution.