Giant Stars Research Articles

Insights into the Galactic Bulge Chemodynamical Properties from Gaia Data Release 3

We explore the chemodynamical properties of the Galaxy in the azimuthal velocity V ϕ and metallicity [Fe/H] space using red giant stars from Gaia Data Release 3. The row-normalized V ϕ –[Fe/H] maps form a coherent sequence from the bulge to the outer disk, clearly revealing the thin/thick disk and the Splash. The metal-rich stars display bar-like kinematics, while the metal-poor stars show dispersion-dominated kinematics. The intermediate-metallicity population (−1 < [Fe/H]< − 0.4) can be separated into two populations, one that is bar-like, i.e., dynamically cold ( σVR∼80 km s−1) and fast-rotating (V ϕ ≳ 100 km s−1), and the Splash, which is dynamically hot ( σVR∼110 km s−1) and slow-rotating (V ϕ ≲ 100 km s−1). We compare the observations in the bulge region with an Auriga simulation where the last major merger event occurred ∼10 Gyr ago: only stars born around the time of the merger reveal a Splash-like feature in the V ϕ –[Fe/H] space, suggesting that the Splash is likely merger-induced, predominantly made up of heated disk stars and the starburst associated with the last major merger. Since the Splash formed from the proto-disk, its lower metallicity limit coincides with that of the thick disk. The bar formed later from the dynamically hot disk with [Fe/H] > − 1 dex, with the Splash not participating in the bar formation and growth. Moreover, with a set of isolated evolving N-body disk simulations, we confirm that a nonrotating classical bulge can be spun up by the bar and develop cylindrical rotation, consistent with the observations for the metal-poor stars.

An empirical view of the extended atmosphere and inner envelope of the asymptotic giant branch star R Doradus

Context. The mass loss experienced on the asymptotic giant branch (AGB) at the end of the lives of low- and intermediate-mass stars is widely accepted to rely on radiation pressure acting on newly formed dust grains. Dust formation happens in the extended atmospheres of these stars, where the density, velocity, and temperature distributions are strongly affected by convection, stellar pulsation, and heating and cooling processes. The interaction between these processes and how that affects dust formation and growth is complex. Hence, characterising the extended atmospheres empirically is paramount to advance our understanding of the dust formation and wind-driving processes. Aims. We aim to determine the density, temperature, and velocity distributions of the gas in the extended atmosphere of the AGB star R Dor. Methods. We acquired observations using ALMA towards R Dor to study the gas through molecular line absorption and emission. We modelled the observed 12CO v = 0, J = 2 − 1, v = 1, J = 2 − 1, and 3 − 2 and 13CO v = 0, J = 3 − 2 lines using the 3D radiative transfer code LIME to determine the density, temperature, and velocity distributions up to a distance of four times the radius of the star at sub-millimetre wavelengths. Results. The high angular resolution of the sub-millimetre maps allows for even the stellar photosphere to be spatially resolved. By analysing the absorption against the star, we infer that the innermost layer in the near-side hemisphere is mostly falling towards the star, while gas in the layer above that seems to be mostly outflowing. Interestingly, the high angular resolution of the ALMA Band 7 observations reveal that the velocity field of the gas seen against the star is not homogenous across the stellar disc. The gas temperature and density distributions have to be very steep close to the star to fit the observed emission and absorption, but they become shallower for radii larger than ∼1.6 times the stellar sub-millimetre radius. While the gas mass in the innermost region is hundreds of times larger than the mass lost on average by R Dor per pulsation cycle, the gas densities just above this region are consistent with those expected based on the mass-loss rate and expansion velocity of the large-scale outflow. Our fits to the line profiles require the velocity distribution on the far side of the envelope to be mirrored, on average, with respect to that on the near side. Using a sharp absorption feature seen in the CO v = 0, J = 2 − 1 line, we constrained the standard deviation of the stochastic velocity distribution in the large-scale outflow to be ≲0.4 km s−1. We characterised two blobs detected in the CO v = 0, J = 2 − 1 line and found densities substantially larger than those of the surrounding gas. The two blobs also display expansion velocities that are high relative to that of the large-scale outflow. Monitoring the evolution of these blobs will lead to a better understanding of the role of these structures in the mass-loss process of R Dor.

Heavy element abundances in galactic globular clusters

Context. Globular clusters are considered key objects for understanding the formation and evolution of the Milky Way. In this sense, the characterisation of their chemical and orbital parameters can provide constraints on chemical evolution models of the Galaxy. Aims. We use the heavy element abundances of globular clusters to trace their overall behaviour in the Galaxy, with the aim to analyse potential relations between the hot H-burning and s-process elements. Methods. We measured the content of Cu I and s- and r-process elements (Y II, Ba II, La II, and Eu II) in a sample of 210 giant stars in 18 galactic globular clusters from high-quality UVES spectra. These clusters span a broad metallicity range and the sample is the largest that has been uniformly analysed to date, with respect to heavy elements in globular clusters. Results. The Cu abundances did not show a considerable spread in the sample, nor any correlation with Na, indicating that the Na nucleosynthesis process does not affect the Cu abundance. Most GCs closely follow the Cu, Y, Ba, La, and Eu field stars’ distribution, revealing a similar chemical evolution. The Y abundances in mid-metallicity regime GCs (−1.10 dex &lt; [Fe/H] &lt; −1.80 dex) display a mildly significant correlation with the Na abundance, which ought to be further investigated. Finally, we do not find any significant difference between the n-capture abundances among GCs with either Galactic and extragalactic origins.

PENELLOPE. VI. Searching the PENELLOPE/UVES sample with spectro-astrometry: Two new microjets of Sz 103 and XXCha

The physical mechanism leading to the formation of the blue loop in the Hertzsprung-Russell (HR) diagram is not satisfactorily explained by the evolutionary track of single stars. Rapid rotation and low metallicity drastically modify the internal structures and surface compositions of stars. Therefore, they provide a very significant pattern to investigate the evolutionary properties of the blue loop. In this paper, we mainly explore how rapid rotation and low metallicity have an important impact on the occurrence and extension of the blue loop. To this end, we implemented the rotating stellar evolution model, including the angular momentum transportation and chemical element mixing. We incorporated several initial rotational velocities and two characteristic metallicities in various models to explore the blue loop extension. The blue loop can occur when the hydrogen burning shell merges with the hydrogen--helium abundance discontinuity. We find that the blue loop extension strongly depends on the amplitude and gradient of the hydrogen--helium discontinuity. The hydrogen--helium discontinuity is created by the intermediate convective region or the convective dredge-up. A steeper hydrogen gradient in association with a greater amplitude of the hydrogen abundance discontinuity may favour a hotter star. Both the low metallicity and rapid rotation tend to restrain the development of the outer convective envelope and thus disfavour the occurrence and extension of the blue loop. There are three main reasons for this occurrence. Firstly, the helium core and its core potential can be enlarged by rotational mixing or low metallicity. Secondly, rapid rotation reduces the convective dredge-up depth in the star with $ Z=0.014$ and the mass extension of the intermediate convective region in the star with $ Z=0.0008$. Both of these phenomena lead to a reduction of the amplitude of the hydrogen abundance gradient. Thirdly, strong rotational mixing in the model (i.e. $ ini =350$ Km/s) with $ Z=0.0008$ reduces the energy generation rate from the hydrogen burning shell. Without bending towards higher effective temperature in the HR diagram, the additional helium brought near the H-burning shell associated with the larger He core can cause the star to expand towards becoming a red giant star directly after the core hydrogen burning. Rapid rotation and low metallicity tend to produce surface enrichment of the ratio of nitrogen to carbon and reduce the $ C$ left in the core; this has an important influence on the stellar compactness of the supernovae progenitor.

Tip of the Red Giant Branch Distances with JWST: An Absolute Calibration in NGC 4258 and First Applications to Type Ia Supernova Hosts

The tip of the red giant branch (TRGB) allows for the measurement of precise and accurate distances to nearby galaxies based on the brightest ascent of low-mass red giant branch stars before they undergo the helium flash. With the advent of JWST, there is great promise to utilize the technique to measure galaxy distances out to at least 50 Mpc, significantly further than the Hubble Space Telescope's (HST's) reach of 20 Mpc. However, with any standard candle, it is first necessary to provide an absolute reference. Here, we use Cycle 1 data to provide an absolute calibration in the F090W filter. F090W is most similar to the F814W filter commonly used for TRGB measurements with HST, which had been adopted by the community due to its minimal dependence on the underlying metallicities and ages of stars. The imaging we use was taken in the outskirts of NGC 4258, which has a direct geometrical distance measurement from the Keplerian motion of its water megamaser. Utilizing several measurement techniques, we find MTRGBF090W = −4.362 ± 0.033 (stat) ± 0.045 (sys) mag (Vega) for the metal-poor TRGB. We also perform measurements of the TRGB in two Type Ia supernova hosts, NGC 1559 and NGC 5584. We find good agreement between our TRGB distances and previous determinations of distances to these galaxies from Cepheids (Δ = 0.01 ± 0.06 mag), with these differences being too small to explain the Hubble tension (∼0.17 mag). In addition, we showcase the serendipitous discovery of a faint dwarf galaxy near NGC 5584.

Expansion of Accreting Main-sequence Stars during Rapid Mass Transfer

Accreting main-sequence stars expand significantly when the mass accretion timescale is much shorter than their thermal timescales. This occurs during mass transfer from an evolved giant star onto a main-sequence companion in a binary system and is an important phase in the formation of compact binaries including X-ray binaries, cataclysmic variables, and gravitational-wave sources. In this study, we compute 1D stellar models of main-sequence accretors with different initial masses and accretion rates. The calculations are used to derive semianalytical approximations to the maximum expansion radius. We assume that mass transfer remains fully conservative as long as the inflated accretor fits within its Roche lobe, leading stars to behave like hamsters, stuffing excess material behind their expanding cheeks. We suggest a physically motivated prescription for the mass growth of such “hamstars,” which can be used to determine mass-transfer efficiency in rapid binary population synthesis models. With this prescription, we estimate that progenitors of high-mass X-ray binaries and gravitational-wave sources may have experienced highly nonconservative mass transfer. In contrast, for low-mass accretors, the accretion timescale can exceed the thermal timescale by a larger factor without causing significant radial expansion.

The Star Formation History in Local Group Galaxies. I. Ten Dwarf Galaxies

The star formation histories (SFHs) of galaxies provide valuable insights into galaxy evolution and stellar physics. Understanding the SFHs enables the study of chemical enrichment of galaxies, star formation triggered by interactions, and the behavior of various stellar populations. This work investigates the SFHs of ten dwarf galaxies in the Local Group (LG), which spans a wide range of types, masses, luminosities, and metallicities. The analysis is based on our new sample of the member stars in the LG after removing foreground dwarf stars by using information from the near-infrared color–color diagram and Gaia astrometry. The sample includes the most complete and pure red supergiants and asymptotic giant branch stars to gain valuable insights into the recent SFHs of the galaxies. A fitting method based on the color–magnitude diagram (CMD) is introduced to measure the SFH. The Padova isochrones are used to generate initial model CMDs, accounting for photometric errors and completeness through star-field simulations to match the completeness and error distributions of the observed CMDs. Subsequently, the SFHs, distance modulus, and metallicity of the ten dwarf galaxies are determined by fitting the CMDs. The results indicate that the star formation rates of dwarf irregulars show a gradual increase, while those of dwarf ellipticals exhibit a gradual decrease from the past to the present. Furthermore, this work shows that the star formation activity in dwarf ellipticals persisted up to 30 Myr ago. A significant increasing feature in the SFH of NGC 6822 reveals star formation activity triggered by an interaction event.

AspGap: Augmented Stellar Parameters and Abundances for 37 Million Red Giant Branch Stars from Gaia XP Low-resolution Spectra

We present AspGap, a new approach to inferring stellar labels from the low-resolution Gaia XP spectra, including precise [α/M] estimates—the first time these are obtained by such an approach. AspGap is a neural-network-based regression model trained on APOGEE spectra. In the training step, AspGap learns to use not only XP spectra to predict stellar labels but also the high-resolution APOGEE spectra that lead to the same stellar labels. The inclusion of this last model component—dubbed the hallucinator—creates a more physically motivated mapping and significantly improves the prediction of stellar labels in the validation, particularly that of [α/M]. For giant stars, we find cross-validated rms accuracies for T eff, log g, [M/H], and [α/M] of ∼1%, 0.12 dex, 0.07 dex, and 0.03 dex, respectively. We also validate our labels through comparison with external data sets and through a range of astrophysical tests that demonstrate that we are indeed determining [α/M] from the XP spectra, rather than just inferring it indirectly from correlations with other labels. We publicly release the AspGap codebase, along with our stellar parameter catalog for all giants observed by Gaia XP. AspGap enables the discovery of new insights into the formation and chemodynamics of our Galaxy by providing precise [α/M] estimates for 37 million giant stars, including 14 million with radial velocities from Gaia.

The APO-K2 Catalog. II. Accurate Stellar Ages for Red Giant Branch Stars across the Milky Way

We present stellar age determinations for 4661 red giant branch stars in the APO-K2 catalog, derived using mass estimates from K2 asteroseismology from the K2 Galactic Archaeology Program and elemental abundances from the Apache Point Galactic Evolution Experiment survey. Our sample includes 17 of the 19 fields observed by K2, making it one of the most comprehensive catalogs of accurate stellar ages across the Galaxy in terms of the wide range of populations spanned by its stars, enabling rigorous tests of Galactic chemical evolution models. Taking into account the selection functions of the K2 sample, the data appear to support the age−chemistry morphology of stellar populations predicted by both inside-out and late-burst scenarios. We also investigate trends in age versus stellar chemistry and Galactic position, which are consistent with previous findings. Comparisons against APOKASC-3 asteroseismic ages show agreement to within ∼3%. We also discuss offsets between our ages and spectroscopic ages. Finally, we note that ignoring the effects of α-enhancement on stellar opacity (either directly or with the Salaris metallicity correction) results in an ∼10% offset in age estimates for the most α-enhanced stars, which is an important consideration for continued tests of Galactic models with this and other asteroseismic age samples.

High mass function ellipsoidal variables in the Gaia Focused Product Release: searching for black hole candidates in the binary zoo

The recent Gaia Focused Product Release contains radial velocity time-series for more than 9,000 Gaia long-period photometric variables. Here we search for binary systems with large radial velocity amplitudes to identify candidates with massive, unseen companions. Eight targets have binary mass function f(M)&gt;1M⊙, three of which are eclipsing binaries. The remaining five show evidence of ellipsoidal modulations. We fit spectroscopic orbit models to the Gaia radial velocities, and fit the spectral energy distributions of three targets. For the two systems most likely to host dark companions, J0946 and J1640, we use PHOEBE to fit the ASAS-SN light curves and Gaia radial velocities. The derived companion masses are &gt;3M⊙, but the high Galactic dust extinctions towards these objects limit our ability to rule out main sequence companions or subgiants hotter than the photometric primaries. These systems are similar to other stellar-mass black hole impostors, notably the Unicorn (V723 Mon) and the Giraffe (2M04123153+6738486). While it is possible that J1640 and J0946 are similar examples of stripped giant star binaries, high-resolution spectra can be used to determine the nature of their companions.

Fine structure in the Sigma Orionis cluster revealed by Gaia DR3

Sigma Orionis is an open cluster in the nearest giant star formation site -- Orion. Its youth (3-5 Myr), low reddening, and relative vicinity make it an important benchmark cluster to study stellar and substellar formation and evolution. Young star-forming sites are complex and hierarchical. Precision astrometry from Gaia \,DR3 enables the exploration of their fine structure. We used the modified convergent point technique to kinematically re-evaluate the members in the Sigma Orionis cluster and its vicinity. We present clear evidence for three kinematically distinct groups in the Sigma Orionis region. The second group, the RV Orionis association, is adjacent to the Sigma Orionis cluster and is composed only of low-mass stars. The third group, the Flame association, whose age is comparable to that of Sigma Orionis, overlaps with the younger NGC\,2024 in the Flame Nebula. In total, we have discovered 105 members of this complex not previously found in the literature (82 in Sigma Orionis, 19 in the Flame association, and 4 in the RV Orionis association).

12C/13C of Kepler giant stars: The missing piece of the mixing puzzle

Context. Despite a rich observational background, few spectroscopic studies have dealt with the measurement of the carbon isotopic ratio in giant stars. However, it is a key element in understanding the mixing mechanisms that occur in the interiors of giant stars. Aims. We present the CNO and 12C/13C abundances derived for 71 giant field stars. Then, using this new catalogue and complementary data from the Kepler and Gaia satellites, we study the efficiency of mixing occurring in the giant branch as a function of the stellar properties of the stars (e.g. mass, age, metallicity). Methods. We determined the abundances of CNO and more specifically the carbon isotopic ratio using the high-resolution FIbre-fed Echelle Spectrograph on the Nordic Optical Telescope, for 71 giant field stars. In addition, asteroseismology from the Kepler satellite is available for all stars, providing the stellar masses, ages, and evolutionary states. Finally, astrometry from the Gaia data is also available for most of the sample. We compare these new determinations with stellar evolution models taking into account the effects of transport processes. To exploit the complete potential of our extensive catalogue, and considering both the Milky Way evolution and the impact of stellar evolution, we built mock catalogues using the Besançon Galaxy model in which stellar evolution models taking into account the effects of thermohaline instability are included. Results. We confirm that the carbon isotopic ratio at the surface of core He-burning stars is lower than that of first-ascent RGB stars. The carbon isotopic ratio measured at the surface of the core He-burning stars increases with [Fe/H] and stellar mass, while it decreases with stellar age. These trends are all nicely explained by the thermohaline mixing that occurs in red giants. We show that our models can explain the behaviour of 12C/13C versus N/O, although the observations seem to show a lower N/O than the models. We also note that more constraints on the thick disc core He-burning stars are needed to understand this difference. Conclusions. Overall, the current model including thermohaline mixing is able to reproduce very well the 12C/13C with the stellar metallicity and with the stellar mass and age.

Evolved Massive Stars at Low Metallicity. VII. The Lower Mass Limit of the Red Supergiant Population in the Large Magellanic Cloud

The precise definition of the lower mass limit of red supergiant stars (RSGs) is an open question in astrophysics and does not attract much attention. Here, we assemble a spectroscopic evolved cool star sample with 6602 targets, including RSGs, asymptotic giant branch stars, and red giant branch stars, in the Large Magellanic Cloud based on Gaia DR3 and Sloan Digital Sky Survey IV/APOGEE-2. The reference spectrum of each stellar population is built according to the quantile range of relative intensity (1% ∼ 99%). Five different methods, e.g., χ 2, cosine similarity, machine learning (ML), equivalent width, and line ratio, are used in order to separate different stellar populations. ML and χ 2 provide the best and relatively consistent prediction of a certain population. The derived lower limit of the RSG population is able to reach the K s -band tip of the red giant branch (K s ≈12.0 mag), indicating a luminosity as low as about 103.5 L ⊙, which corresponds to a stellar radius of only about 100 R ⊙. Given the mass–luminosity relation of L/L⊙=f(M/M⊙)3 with f ≈ 15.5 ± 3 and taking into account the mass loss of faint RSGs up to now, the minimal initial mass of the RSG population would be about 6.1 ± 0.4 M ⊙, which is much lower than the traditional threshold of 8 M ⊙ for the massive stars. This is the first spectroscopic evidence, indicating that the lower mass limit of the RSG population is around 6 M ⊙. However, the destinies of such faint RSGs are still elusive and may have a large impact on stellar evolutionary and supernova models.

On the Performances of Estimating Stellar Atmospheric Parameters from CSST Broad-band Photometry

Deriving atmospheric parameters of a large sample of stars is of vital importance to understand the formation and evolution of the Milky Way. Photometric surveys, especially those with near-ultraviolet filters, can offer accurate measurements of stellar parameters, with the precision comparable to that from low/medium resolution spectroscopy. In this study, we explore the capability of measuring stellar atmospheric parameters from Chinese Space Station Telescope (CSST) broad-band photometry (particularly in the near-ultraviolet bands), based on synthetic colors derived from model spectra. We find that colors from the optical and near-ultraviolet filter systems adopted by CSST show significant sensitivities to the stellar atmospheric parameters, especially the metallicity. According to our mock data tests, the precision of the photometric metallicity is quite high, with typical values of 0.17 and 0.20 dex for dwarf and giant stars, respectively. The precision of the effective temperature estimated from broad-band colors are within 50 K.

CO emission survey of asymptotic giant branch stars with ultraviolet excesses

Context. The transition from the spherically symmetric envelopes around asymptotic giant branch (AGB) stars to the asymmetric morphologies observed in planetary nebulae is still not well understood, and the shaping mechanisms are a subject of debate. Even though binarity is widely accepted as a promising option, it is limited by the complication of identifying binary AGB stars observationally. Recently, the presence of ultraviolet excesses in AGB stars has been suggested as a potential indicator of binarity. Aims. Our main goals are to characterise the properties of the circumstellar envelopes (CSEs) around candidate AGB binary stars, specifically those selected based on their UV excess emission, and to compare these properties with those derived from previous CO-based studies of AGB stars. Methods. We observed the 12CO (J=1–0) and 12CO (J=2–1) millimetre-wavelength emission in a sample of 29 AGB binary candidates with the IRAM-30 m antenna. We measured the systemic velocities and the terminal expansion velocities from their line profiles. Population diagrams were used to interpret the results, enabling the estimation of excitation temperatures (Tex), mass-loss rates (Ṁ), and the characteristic sizes of the envelope layers where the CO millimetre emission originates (Rs). We explored different trends between the envelope parameters deduced, multiwavelength flux measurements, and other properties of our sample, and compared them with those previously derived from larger samples of AGB stars found in the literature. Results. We detected 12CO emission in 15 sources, of which 5 are first detections. We found relatively low expansion velocities (3 km s−1 ≲ Vexp ≲ 20 km s−1) in our sample. We derived the average excitation temperature and column density of the CO-emitting layers, which we used to estimate self-consistently the average mass-loss rate (10−8 M⊙ yr−1 ≲ Ṁ ≲ 10−5 M⊙ yr−1) and the CO pho-todissociation radius (5 × 1015 cm ≲ RCO ≲ 2 × 1017 cm) of our targets. We find a correlation between CO intensity and IRAS 60 µm fluxes, revealing a CO-to-IRAS 60 µm ratio lower than for AGB stars and closer to that found for pre-planetary nebulae (pPNe). An anti-correlation is observed between 12CO (and IRAS 60 µm) and the near-ultraviolet (NUV), but no such correlation is observed with the far-ultraviolet (FUV). It is also worth noting that there is no correlation between bolometric luminosity and NUV or FUV. Conclusions. For the first time we have studied the mass-loss properties of UV-excess AGB binary candidates and estimated their main CSE parameters. Our sample of uvAGB stars shows similarities with the broader category of AGB stars, except for a distinct CO-to-IRAS 60 µm trend suggesting enhanced CO photodissociation. Our findings, based on single-dish low-J CO line emission observations, support the dust-driven wind scenario and indicate that alternative mass-loss mechanisms are not necessary (in principle) to explain the ~200–2000 yr old mass-loss ejecta in uvAGBs. The different relationships between 12CO and IRAS 60 µm, with NUV and FUV are consistent with an intrinsic origin of NUV emission, but potential dominance of an extrinsic process (e.g. presence of a binary companion) in FUV emission.

Fossil Signatures of Main-sequence Convective Core Overshoot Estimated through Asteroseismic Analyses

Some physical processes that occur during a star's main-sequence evolution also affect its post-main-sequence evolution. It is well known that stars with masses above approximately 1.1 M ⊙ have well-mixed convective cores on the main sequence; however, the structure of the star in the neighborhood of the convective core regions is currently underconstrained. We use asteroseismology to study the properties of the stellar core, in particular convective boundary mixing through convective overshoot, in such intermediate-mass stars. These core regions are poorly constrained by the acoustic (p) mode oscillations observed for cool main-sequence stars. Consequently, we seek fossil signatures of main-sequence core properties during the subgiant and early first-ascent red giant phases of evolution. During these stages of stellar evolution, modes of mixed character that sample the deep interior can be observed. These modes sample the parts of the stars that are affected by the main-sequence structure of these regions. We model the global and near-core properties of 62 subgiant and early first-ascent red giant branch stars observed by the Kepler, K2, and TESS space missions. We find that the effective overshoot parameter, α ov,eff, increases from M = 1.0 M ⊙ to M = 1.2 M ⊙ before flattening out, although we note that the relationship between α ov,eff and mass will depend on the incorporated modeling choices of internal physics and nuclear reaction network. We also situate these results within existing studies of main-sequence convective core boundaries.

Lithium Abundances from the LAMOST Medium-resolution Survey Data Release 9

Lithium is a fragile but crucial chemical element in the Universe, and exhibits interesting and complex behaviors. Thanks to the mass of spectroscopic data from the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) medium-resolution survey (MRS), we can investigate the lithium abundances in a large and diverse sample of stars, which could bring vital help in studying the origin and evolution of lithium. In this work, we use the Li i 6707.8 Å line to derive the lithium abundance through a template-matching method. A catalog of precise lithium abundance is presented for 795,384 spectra corresponding to 455,752 stars from the LAMOST MRS Data Release 9. Comparing our results with those of external high-resolution references, we find good consistency with a typical deviation of σ A(Li) ∼ 0.2 dex. We also analyze the internal errors using stars that have multiple LAMOST MRS observations, which will reach as low as 0.1 dex when the signal-to-noise ratio of the spectra is >20. Besides, our result indicates that a small fraction of giant stars still exhibit a surprisingly high lithium content, and 967 stars are identified as Li-rich giants with A(Li) > 1.5 dex, accounting for ∼2.6% of our samples. If one takes into account the fact that nearly all stars deplete lithium during the main sequence, then the fraction of Li-rich stars may far exceed 2.6%. This new catalog covers a wide range of stellar evolutionary stages from pre-main sequence to giants, and will provide help to the further study of the chemical evolution of lithium.

A comprehensive search for hot subdwarf stars using Gaia and TESS

Hot subdwarf B (sdB) stars are evolved, subluminous, helium-burning stars that most likely form when red giant stars loose their hydrogen envelope via interactions with close companions. They play an important role in our understanding of binary evolution, stellar atmospheres, and interiors. Only a small fraction of the sdB population is known to exhibit pulsations. Pulsating sdBs have typically been discovered serendipitously in various photometric surveys because specific selection criteria for the sample are lacking. Consequently, while individual properties of these stars are well known, a comprehensive understanding of the entire population remains elusive, and many related questions remain unanswered. The Gaia mission has presented an exceptional chance to create an unbiased sample by employing precise criteria and ensuring a high degree of completeness. The progression of high-precision and high-duty cycle photometric monitoring facilitated by space missions such as Kepler/K2 and the Transiting Exoplanet Survey Satellite (TESS) has yielded an unparalleled wealth of data for pulsating sdBs. We created a dataset of confirmed pulsating sdB stars by combining information from various ground- and space-based photometric surveys. With this dataset, we present a thorough approach to search for pulsating sdB stars based on the current Gaia DR3 sample. Based on TESS photometry, we discovered 61 new pulsating sdB stars and 20 variable sdBs whose source of variability remains to be determined through future spectroscopic follow-up observations.

Charting circumstellar chemistry of carbon-rich asymptotic giant branch stars

Context. Asymptotic giant branch (AGB) stars are major contributors to the chemical enrichment of the interstellar medium through nucleosynthesis and extensive mass loss. Direct measures of both processes can be obtained by studying their circumstellar envelopes in molecular line emission. Most of our current knowledge of circumstellar chemistry, in particular in a C-rich environment, is based on observations of the carbon star IRC +10216. Aims. We aim to obtain a more generalised understanding of the chemistry in C-rich AGB circumstellar envelopes by studying a sample of three carbon stars, IRAS 15194–5115, IRAS 15082–4808, and IRAS 07454–7112, and to observationally test the archetypal status often attributed to IRC +10216. Methods. We performed spatially resolved, unbiased spectral surveys in ALMA Band 3 (85–116 GHz). We estimated the sizes of the molecular emitting regions using azimuthally averaged radial profiles of the line brightness distributions. We derived abundance estimates, using a population diagram analysis for molecules with multiple detected lines, and using single-line analytical calculations for the others. Results. We identify a total of 132 rotational transitions from 49 molecular species. There are two main morphologies of the brightness distributions: centrally peaked (CS, SiO, SiS, HCN) and shell-like (CN, HNC, C2H, C3H, C4H, C3N, HC5N, c-C3H2). The brightness distributions of HC3N and SiC2 have both a central and a shell component. The qualitative behaviour of the brightness distributions of all detected molecules, in particular their relative locations with respect to the central star, is the same for all three stars, and consistent with those observed towards IRC +10216. Of the shell distributions, the cyanopolyynes peak at slightly smaller radii than the hydrocarbons, and CN and HNC show the most extended emission. The emitting regions for each species are the smallest for IRAS 07454–7112, consistent with this object having the lowest circumstellar density within our sample. We find that, within the uncertainties of the analysis, the three stars present similar abundances for most species, and also compared to IRC +10216. We find, tentatively, that SiO is more abundant in our three stars compared to IRC+10216, and that the hydrocarbons are under-abundant in IRAS 07454–7112 compared to the other stars and IRC +10216. Our estimated 12C/13C ratios match well the literature values for the three sources and our estimated silicon and sulphur isotopic ratios are very similar across the three stars and IRC +10216. Conclusions. The observed circumstellar chemistry appears very similar across our sample and compared to that of IRC +10216, both in terms of the relative location of the emitting regions and molecular abundances. This implies that, to a first approximation, the chemical models tailored to IRC +10216 are, at least, able to reproduce the observed chemistry in C-rich envelopes across roughly an order of magnitude in wind density.

Three-dimensional time-dependent convection model for asteroseismology

Context. Due to an ill-depicting model of the convective layers below the photosphere in 1D stellar models (structural contribution) and/or a misrepresentation of the coupling between convection and oscillations (modal contribution), a well-known deviation appears between observed and theoretical frequencies, which grows towards high frequencies; the so-called surface effects. While satisfying solutions have been found regarding the structural contribution, the accurate modeling of the modal effect still represents a challenge. Alongside the frequency, the interaction between convection and oscillations also impacts the damping rate of the modes and forms an important part of the driving mechanism behind the stellar oscillations of low-mass stars. With increasing observational capabilities at our disposal with Kepler and TESS, shortcomings in modeling constitute the main limitation to accurate seismic probing of solar-like and red giant stars. Aims. We present the formalism of an approach that changes the current paradigm by addressing three-dimensional space. This new formalism consists in an original nonadiabatic 3D time-dependent convection model for asteroseismology. Methods. We aim to keep the entire 3D structure of the astrophysical flow in these superficial layers in order to fully account for the nature of turbulence in our model via the use of advanced hydrodynamic simulation. We use the perturbative approach and introduce a spectral decomposition approach that results in an entirely new formalism describing standing waves in 3D. This formalism is set to solve the quasi-radial global nonadiabatic oscillation equations in a full 3D framework. Results. Based on physical assumptions, we establish an eigenvalue problem describing the 3D quasi-radial global nonadiabatic stellar oscillation. We also provide a prescription for its numerical resolution alongside a proposed iteration method for our formalism. Finally, we derive the peculiar 3D work integral and establish the expression of the damping rate. We show how our formalism offers the possibility to probe the complex structure of stars and is able to precisely locate regions of the driving and damping of the modes as well as their physical origin.