Search

The granulation of red supergiants (RSGs) in the Magellanic Clouds is systematically investigated by combining the latest RSG samples and light curves from the Optical Gravitational Lensing Experiment and the All-Sky Automated Survey for Supernovae. The present RSG samples are first examined for foreground stars and possible misidentified sources, and the light curves are sequentially checked to remove the outliers by white noise and photometric quality. The Gaussian process (GP) regression is used to model the granulation, and the Markov Chain Monte Carlo is applied to derive the granulation amplitude σ and the period of the undamped oscillator ρ, as well as the damping timescale τ. The dimensionless quality factor Q is then calculated through Q = π τ / ρ. RSGs around Q = 1 / 2 are considered to have significant granulation signals and are used for further analysis. Combining granulation parameters with stellar parameters, robust scaling relations for the timescale ρ are established, while the scaling relations for amplitude σ are represented by a piecewise function, possibly related to the tendency of amplitudes in faint RSGs to converge toward a certain value. Comparing results between the Small Magellanic Cloud and the Large Magellanic Cloud confirms that amplitudes and timescales become larger with metallicity. In examining the scaling relations between the two galaxies, it is found that ρ is nearly independent of metallicity, whereas σ is more significantly affected by metallicity. The GP method is compared with the periodogram fitting of the granulations, and the advantages of either are discussed. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

SN 2023ixf was discovered in M101 within a day of the explosion and rapidly classified as a Type II supernova with flash features. Here we present ultraviolet (UV) spectra obtained with the Hubble Space Telescope 14, 19, 24, and 66 days after the explosion. Interaction between the supernova ejecta and circumstellar material (CSM) is seen in the UV throughout our observations in the flux of the first three epochs and asymmetric Mg ii emission on day 66. We compare our observations to CMFGEN supernova models that include CSM interaction ( $\dot{M}\lt {10}^{-3}$ M _⊙ yr ^−1 ) and find that the power from CSM interaction is decreasing with time, from L _sh ≈ 5 × 10 ^42 erg s ^−1 to L _sh ≈ 1 × 10 ^40 erg s ^−1 between days 14 and 66. We examine the contribution of individual atomic species to the spectra on days 14 and 19, showing that the majority of the features are dominated by iron, nickel, magnesium, and chromium absorption in the ejecta. The UV spectral energy distribution of SN 2023ixf sits between that of supernovae, which show no definitive signs of CSM interaction, and those with persistent signatures assuming the same progenitor radius and metallicity. Finally, we show that the evolution and asymmetric shape of the Mg ii λ λ 2796, 2802 emission are not unique to SN 2023ixf. These observations add to the early measurements of dense, confined CSM interaction, tracing the mass-loss history of SN 2023ixf to ∼33 yr prior to the explosion and the density profile to a radius of ∼5.7 × 10 ^15 cm. They show the relatively short evolution from a quiescent red supergiant wind to high mass loss.

We present high-cadence photometric and spectroscopic observations of supernova (SN) 2024ggi, a Type II SN with flash spectroscopy features, which exploded in the nearby galaxy NGC 3621 at ∼7 Mpc. The light-curve evolution over the first 30 hr can be fit by two power-law indices with a break after 22 hr, rising from M _V ≈ −12.95 mag at +0.66 day to M _V ≈ −17.91 mag after 7 days. In addition, the densely sampled color curve shows a strong blueward evolution over the first few days and then behaves as a normal SN II with a redward evolution as the ejecta cool. Such deviations could be due to interaction with circumstellar material (CSM). Early high- and low-resolution spectra clearly show high-ionization flash features from the first spectrum to +3.42 days after the explosion. From the high-resolution spectra, we calculate the CSM velocity to be 37 ± 4 km s ^−1 . We also see the line strength evolve rapidly from 1.22 to 1.49 days in the earliest high-resolution spectra. Comparison of the low-resolution spectra with CMFGEN models suggests that the pre-explosion mass-loss rate of SN 2024ggi falls in the range of 10 ^−3 –10 ^−2 M _☉ yr ^−1 , which is similar to that derived for SN 2023ixf. However, the rapid temporal evolution of the narrow lines in the spectra of SN 2024ggi ( R _CSM ∼ 2.7 × 10 ^14 cm) could indicate a smaller spatial extent of the CSM than in SN 2023ixf ( R _CSM ∼ 5.4 × 10 ^14 cm), which in turn implies a lower total CSM mass for SN 2024ggi.

We present early-phase good-cadence (hour-to-day) simultaneous multiband ( ugi and vrz bands) imaging of the nearby supernova SN 2024ggi, which exploded in the nearby galaxy NGC 3621. A quick follow-up was conducted within less than a day after the explosion and continued for ∼23 days. The uvg band light curves display a rapid rise (∼1.4 mag day ^−1 ) to maximum in ∼4 days and absolute magnitude M _g ∼ −17.75 mag. The postpeak decay rate in redder bands is ∼0.01 mag day ^−1 . Different colors (e.g., u − g and v − r ) of SN 2024ggi are slightly redder than SN 2023ixf. A significant rise (∼12.5 kK) in blackbody temperature (optical) was noticed within ∼2 days after the explosion, which successively decreased, indicating shock break out inside a dense circumstellar medium surrounding the progenitor. Using semianalytical modeling, the ejecta mass and progenitor radius were estimated as 1.2 M _☉ and ∼550 R _☉ . The archival deep images ( g , r , i , and z bands) from the Dark Energy Camera Legacy Survey were examined, and a possible progenitor was detected in each band (∼22–22.5 mag) and had a mass range of 14–17 M _☉ .

We present a detailed analysis of the progenitor and its local environment for the recently discovered Type II supernova (SN) 2024ggi at a distance of about 6.7 Mpc, by utilizing the pre-explosion images from the Hubble Space Telescope and Spitzer Space Telescope. The progenitor is identified as a red bright variable star, with absolute F814W-band magnitudes being −6.2 mag in 1995 to −7.2 mag in 2003, respectively, consistent with that of a normal red supergiant star. Combining with the historical mid-infrared light curves, a pulsational period of about 379 days can be inferred for the progenitor star. Fitting its spectral energy distribution with stellar spectral models yields the stellar parameters of temperature, radius, and bolometric luminosity as ${T}_{* }={3290}_{-27}^{+19}$ K, ${R}_{* }={887}_{-51}^{+60}$ R _⊙ , and log( L / L _⊙ ) $\,=\,{4.92}_{-0.04}^{+0.05}$ , respectively. The above parameters indicate that the progenitor of SN 2024ggi is consistent with the stellar evolutionary track of a solar-metallicity massive star with an initial mass of ${13}_{-1}^{+1}$ M _⊙ . Moreover, our analysis indicates a relatively low mass-loss rate (i.e.,

Emission lines from Rydberg transitions are detected for the first time from a region close to the surface of Betelgeuse. The H30 α line is observed at 231.905 GHz, with an FWHM ∼42 km s ^−1 and extended wings. A second line at 232.025 GHz (FWHM ∼21 km s ^−1 ), is modeled as a combination of Rydberg transitions of abundant low first ionization potential metals. Both H30 α and the Rydberg combined line X30 α are fitted by Voigt profiles, and collisional broadening with electrons may be partly responsible for the Lorentzian contribution, indicating electron densities of a few 10 ^8 cm ^−3 . X30 α is located in a relatively smooth ring at a projected radius of 0.9× the optical photospheric radius R _⋆ , whereas H30 α is more clumpy, reaching a peak at ∼1.4 R _⋆ . We use a semiempirical thermodynamic atmospheric model of Betelgeuse to compute the 232 GHz (1.29 mm) continuum and line profiles making simple assumptions. Photoionized abundant metals dominate the electron density, and the predicted surface of continuum optical depth unity at 232 GHz occurs at ∼1.3 R _⋆ , in good agreement with observations. Assuming a Saha–Boltzmann distribution for the level populations of Mg, Si, and Fe, the model predicts that the X30 α emission arises in a region of radially increasing temperature and turbulence. Inclusion of ionized C and non-LTE effects could modify the integrated fluxes and location of emission. These simulations confirm the identity of the Rydberg transition lines observed toward Betelgeuse and reveal that such diagnostics can improve future atmospheric models.

The apparent tension between the luminosity functions of red supergiant (RSG) stars and of RSG progenitors of Type II supernovae (SNe) is often referred to as the RSG problem and it motivated some to suggest that many RSGs end their life without an SN explosion. However, the luminosity functions of RSG SN progenitors presented so far were biased to high luminosities, because the sensitivity of the search was not considered. Here, we use limiting magnitudes to calculate a bias-corrected RSG progenitor luminosity function. We find that only (36 ± 11)% of all RSG progenitors are brighter than a bolometric magnitude of −7 mag, a significantly smaller fraction than (56 ± 5)% quoted by Davies & Beasor. The larger uncertainty is due to the relatively small progenitor sample, while uncertainties on measured quantities such as magnitudes, bolometric corrections, extinction, or SN distances, only have a minor impact, as long as they fluctuate randomly for different objects in the sample. The bias-corrected luminosity functions of RSG SN progenitors and Type M supergiants in the Large Magellanic Cloud are consistent with each other, as also found by Davies & Beasor for the uncorrected luminosity function. The RSG progenitor luminosity function, hence, does not imply the existence of failed SNe. The presented statistical method is not limited to progenitor searches, but applies to any situation in which a measurement is done for a sample of detected objects, but the probed quantity or property can only be determined for part of the sample.

The evolved stages of massive stars are poorly understood, but invaluable constraints can be derived from spatially resolved observations of nearby red supergiants, such as Betelgeuse. Atacama Large Millimeter/submillimeter Array (ALMA) observations of Betelgeuse showing a dipolar velocity field have been interpreted as evidence for a projected rotation rate of about 5 km s ^−1 . This is 2 orders of magnitude larger than predicted by single-star evolution, which led to suggestions that Betelgeuse is a binary merger. We propose instead that large-scale convective motions can mimic rotation, especially if they are only partially resolved. We support this claim with 3D CO5BOLD simulations of nonrotating red supergiants that we postprocessed to predict ALMA images and SiO spectra. We show that our synthetic radial velocity maps have a 90% chance of being falsely interpreted as evidence for a projected rotation rate of 2 km s ^−1 or larger for our fiducial simulation. We conclude that we need at least another ALMA observation to firmly establish whether Betelgeuse is indeed rapidly rotating. Such observations would also provide insight into the role of angular momentum and binary interaction in the late evolutionary stages. The data will further probe the structure and complex physical processes in the atmospheres of red supergiants, which are immediate progenitors of supernovae and are believed to be essential in the formation of gravitational-wave sources.

We perform a comprehensive search for optical precursor emission at the position of SN 2023ixf using data from the DLT40, ZTF, and ATLAS surveys. By comparing the current data set with precursor outburst hydrodynamical model light curves, we find that the probability of a significant outburst within 5 yr of explosion is low, and the circumstellar material (CSM) ejected during any possible precursor outburst is likely smaller than ∼0.015 M ⊙. By comparing to a set of toy models, we find that, if there was a precursor outburst, the duration must have been shorter than ∼100 days for a typical brightness of M r ≃ −9 mag or shorter than 200 days for M r ≃ −8 mag; brighter, longer outbursts would have been discovered. Precursor activity like that observed in the normal Type II SN 2020tlf (M r ≃ −11.5) can be excluded in SN 2023ixf. If the dense CSM inferred by early flash spectroscopy and other studies is related to one or more precursor outbursts, then our observations indicate that any such outburst would have to be faint and only last for days to months, or it occurred more than 5 yr prior to the explosion. Alternatively, any dense, confined CSM may not be due to eruptive mass loss from a single red supergiant progenitor. Taken together, the results of SN 2023ixf and SN 2020tlf indicate that there may be more than one physical mechanism behind the dense CSM inferred around some normal Type II supernovae. [ABSTRACT FROM AUTHOR]

We present a detailed investigation of photometric, spectroscopic, and polarimetric observations of the Type II SN 2023ixf. Earlier studies have provided compelling evidence for a delayed shock breakout from a confined dense circumstellar matter (CSM) enveloping the progenitor star. The temporal evolution of polarization in the SN 2023ixf phase revealed three distinct peaks in polarization evolution at 1.4 days, 6.4 days, and 79.2 days, indicating an asymmetric dense CSM, an aspherical shock front and clumpiness in the low-density extended CSM, and an aspherical inner ejecta/He-core. SN 2023ixf displayed two dominant axes, one along the CSM-outer ejecta and the other along the inner ejecta/He-core, showcasing the independent origin of asymmetry in the early and late evolution. The argument for an aspherical shock front is further strengthened by the presence of a high-velocity broad absorption feature in the blue wing of the Balmer features in addition to the P-Cygni absorption post-16 days. Hydrodynamical light-curve modeling indicated a progenitor mass of 10 M _⊙ with a radius of 470 R _⊙ and explosion energy of 2 × 10 ^51 erg, along with 0.06 M _⊙ of ^56 Ni, though these properties are not unique due to modeling degeneracies. The modeling also indicated a two-zone CSM: a confined dense CSM extending up to 5 × 10 ^14 cm with a mass-loss rate of 10 ^−2 M _⊙ yr ^−1 and an extended CSM spanning from 5 × 10 ^14 to at least 10 ^16 cm with a mass-loss rate of 10 ^−4 M _⊙ yr ^−1 , both assuming a wind-velocity of 10 km s ^−1 . The early-nebular phase observations display an axisymmetric line profile of [O i ], redward attenuation of the emission of H α post 125 days, and flattening in the Ks -band, marking the onset of dust formation.

We present detailed multiband photometric and spectroscopic observations and analysis of a rare core-collapse supernova, SN 2021wvw, that includes photometric evolution up to 250 days and spectroscopic coverage up to 100 days postexplosion. A unique event that does not fit well within the general trends observed for Type IIP supernovae, SN 2021wvw shows an intermediate luminosity with a short plateau phase of just about 75 days, followed by a very sharp (∼10 days) transition to the tail phase. Even in the velocity space, it lies at a lower velocity compared to a larger Type II sample. The observed peak absolute magnitude is −16.1 mag in r -band, and the nickel mass is well constrained to 0.020 ± 0.006 M _⊙ . Detailed hydrodynamical modeling using MESA+STELLA suggests a radially compact, low-metallicity, high-mass red supergiant progenitor ( M _ZAMS = 18 M _⊙ ), which exploded with ∼0.2 × 10 ^51 erg s ^−1 leaving an ejecta mass of M _ej ≈ 5 M _⊙ . Significant late-time fallback during the shock propagation phase is also seen in progenitor+explosion models consistent with the light-curve properties. As the faintest short-plateau supernova characterized to date, this event adds to the growing diversity of transitional events between the canonical ∼100 days plateau Type IIP and stripped-envelope events.

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.

We present comprehensive optical observations of SN 2021gmj, a Type II supernova (SN II) discovered within a day of explosion by the Distance Less Than 40 Mpc survey. Follow-up observations show that SN 2021gmj is a low-luminosity SN II (LL SN II), with a peak magnitude M _V = −15.45 and an Fe ii velocity of ∼1800 km s ^−1 at 50 days past explosion. Using the expanding photosphere method, we derive a distance of ${17.8}_{-0.4}^{+0.6}$ Mpc. From the tail of the light curve we obtain a radioactive nickel mass of ${M}_{{}^{56}\mathrm{Ni}}$ = 0.014 ± 0.001 M _⊙ . The presence of circumstellar material (CSM) is suggested by the early-time light curve, early spectra, and high-velocity H α in absorption. Analytical shock-cooling models of the light curve cannot reproduce the fast rise, supporting the idea that the early-time emission is partially powered by the interaction of the SN ejecta and CSM. The inferred low CSM mass of 0.025 M _⊙ in our hydrodynamic-modeling light-curve analysis is also consistent with our spectroscopy. We observe a broad feature near 4600 Å, which may be high-ionization lines of C, N, or/and He ii . This feature is reproduced by radiation-hydrodynamic simulations of red supergiants with extended atmospheres. Several LL SNe II show similar spectral features, implying that high-density material around the progenitor may be common among them.

Simultaneous very-long-baseline interferometry monitoring observations of H _2 O and SiO masers toward VX Sagittarii were conducted from 2014 February to 2019 January. Thirty epochs of observations revealed that the H _2 O and SiO masers had asymmetric and ring-like structures, respectively. However, from 2017 September to 2018 March, the SiO maser transformed from a ring-like structure to a northeast–southwest (NE–SW) extension, and the 43.1 and 86.2 GHz SiO maser components had velocities of 39.48 and 10.65 km s ^−1 in the NE–SW direction, suggesting a possible localized strong shock wave. The H _2 O maser had a double-sided structure oriented in the NE–SW direction with near-stellar velocity components, which aligned with the extended direction of the SiO maser. The nonregular optical brightness and maser intensity variations were speculated to be related to the morphological evolution of the SiO maser. During the stable states attained by regular pulsations, the SiO maser region was presumed to experience radial acceleration, which reverted the SiO maser to a ring-like structure. However, the H _2 O maser region, where the acceleration almost terminates, retained its asymmetric morphology due to the prior influence of external forces. The results suggest that substantial energy transfer can alter the dynamics of the SiO maser and surrounding atmosphere, leading to an asymmetric distribution in the H _2 O maser region.

The granulation of red supergiants (RSGs) in the Magellanic Clouds is systematically investigated by combining the latest RSG samples and light curves from the Optical Gravitational Lensing Experiment and the All-Sky Automated Survey for Supernovae. The present RSG samples are first examined for foreground stars and possible misidentified sources, and the light curves are sequentially checked to remove the outliers by white noise and photometric quality. The Gaussian process (GP) regression is used to model the granulation, and the Markov Chain Monte Carlo is applied to derive the granulation amplitude σ and the period of the undamped oscillator ρ , as well as the damping timescale τ . The dimensionless quality factor Q is then calculated through Q = π τ / ρ . RSGs around $Q=1/\sqrt{2}$ are considered to have significant granulation signals and are used for further analysis. Combining granulation parameters with stellar parameters, robust scaling relations for the timescale ρ are established, while the scaling relations for amplitude σ are represented by a piecewise function, possibly related to the tendency of amplitudes in faint RSGs to converge toward a certain value. Comparing results between the Small Magellanic Cloud and the Large Magellanic Cloud confirms that amplitudes and timescales become larger with metallicity. In examining the scaling relations between the two galaxies, it is found that ρ is nearly independent of metallicity, whereas σ is more significantly affected by metallicity. The GP method is compared with the periodogram fitting of the granulations, and the advantages of either are discussed.

We present a JWST imaging survey of I Zw 18, the archetypal extremely metal-poor, star-forming (SF), blue compact dwarf galaxy. With an oxygen abundance of only ∼3% Z _⊙ , it is among the lowest-metallicity systems known in the local Universe, and is, therefore, an excellent accessible analog for the galactic building blocks which existed at early epochs of ionization and star formation. These JWST data provide a comprehensive infrared (IR) view of I Zw 18 with eight filters utilizing both Near Infrared Camera (F115W, F200W, F356W, and F444W) and Mid-Infrared Instrument (F770W, F1000W, F1500W, and F1800W) photometry, which we have used to identify key stellar populations that are bright in the near- and mid-IR. These data allow for a better understanding of the origins of dust and dust-production mechanisms in metal-poor environments by characterizing the population of massive, evolved stars in the red supergiant (RSG) and asymptotic giant branch (AGB) phases. In addition, it enables the identification of the brightest dust-enshrouded young stellar objects (YSOs), which provide insight into the formation of massive stars at extremely low metallicities typical of the very early Universe. This paper provides an overview of the observational strategy and data processing, and presents first science results, including identifications of dusty AGB, RSG, and bright YSO candidates. These first results assess the scientific quality of JWST data and provide a guide for obtaining and interpreting future observations of the dusty and evolved stars inhabiting compact dwarf SF galaxies in the local Universe.

We observed the Galactic open cluster Berkeley 50 in order to determine its stellar content, distance, and age. We obtained UBV photometry of 1145 stars in a 12.′3 × 12.′3 field, and used Gaia proper motions and parallaxes to identify 64 members, of which we obtained spectra of the 17 brightest members. The majority of the observed population we classified as B dwarfs, with the exception of a newly identified red supergiant star, which our spectroscopy shows has a B-type companion. Our study establishes the distance as 3.8 kpc, with an average color-excess E ( B − V ) = 0.9. Comparison of the physical properties of the cluster with the Geneva evolutionary tracks places the age of the cluster as 50–60 Myr, with its most massive members being ∼7 M _⊙

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.

Young massive star clusters, like the six red supergiant clusters in the Scutum complex, provide valuable insights into star formation and galaxy structures. We investigated the high-resolution near-infrared spectra of 60 RSG candidates in these clusters using the Immersion Grating Infrared Spectrograph. Among the candidates in RSGC4, we found significant scattering in radial velocity (−64 to 115 km s ^−1 ), unlike other clusters with velocities of ∼100 km s ^−1 . Most candidates in RSGC4 have Q _GK _s values larger than 1.7, suggesting that they could be early AGB stars. Four candidates in RSGC4 exhibit infrared excess and distinct absorption features absent in other candidates. Two of these stars exhibit absorption lines resembling those of D-type symbiotic stars, showing radial velocity changes in multiepoch observations. Analysis of relative proper motions revealed no runaway/walkaway stars in RSGC4. The dynamic properties of RSGC4 and RSGC1 differ from the disklike motions of other clusters: RSGC4 has low normalized horizontal action J _hor = J _ϕ / J _tot and vertical action J _ver = ( J _z − J _R )/ J _tot values and high eccentricities, while RSGC1 has vertical motions with high J _ver values and inclinations. We propose that RSGC4 may not be a genuine star cluster but rather a composite of RSGs and AGB stars distributed along the line of sight at similar distances, possibly originating from various environments. Our results suggest a complex and hierarchical secular evolution of star clusters in the Scutum complex, emphasizing the importance of considering factors beyond density crowding when identifying star clusters in the bulge regions.

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 10 ^3.5 L _⊙ , which corresponds to a stellar radius of only about 100 R _⊙ . Given the mass–luminosity relation of $L/{L}_{\odot }=f{(M/{M}_{\odot })}^{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.

We investigate the merger between a 16 M _⊙ star, on its way to becoming a red supergiant (RSG), and a 4 M _⊙ main-sequence companion. Our study employs three-dimensional hydrodynamic simulations using the state-of-the-art adaptive mesh refinement code O cto -T iger . The initially corotating binary undergoes interaction and mass transfer, resulting in the accumulation of mass around the companion and its subsequent loss through the second Lagrangian point (L2). The companion eventually plunges into the envelope of the primary, leading to its spin-up and subsequent merger with the helium core. We examine the internal structural properties of the post-merger star, as well as the merger environment and the outflow driven by the merger. Our findings reveal the ejection of approximately ∼0.6 M _⊙ of material in an asymmetric and somewhat bipolar outflow. We import the post-merger stellar structure into the MESA stellar evolution code to model its long-term nuclear evolution. In certain cases, the post-merger star exhibits persistent rapid equatorial surface rotation as it evolves in the H – R diagram toward the observed location of Betelgeuse. These cases demonstrate surface rotation velocities of a similar magnitude to those observed in Betelgeuse, along with a chemical composition resembling that of Betelgeuse. In other cases, efficient rotationally induced mixing leads to slower surface rotation. This pioneering study aims to model stellar mergers across critical timescales, encompassing dynamical, thermal, and nuclear evolutionary stages.

The red hypergiant VY CMa is famous for its very visible record of high-mass-loss events. Recent CO observations with the Atacama Large Millimeter/submillimeter Array (ALMA) revealed three previously unknown large-scale outflows (Singh et al). In this paper, we use the CO maps to investigate the motions of a cluster of four clumps close to the star, not visible in the optical or infrared images. We present their proper motions measured from two epochs of ALMA images and determine the line-of-sight velocities of the gas in emission at the clumps. We estimate their masses and ages, or time since ejection, and conclude that all four were ejected during VY CMa’s active period in the early 20th century. Together with two additional knots observed with the Hubble Space Telescope, VY CMa experienced at least six massive outflows during a 30 yr period, with a total mass lost ≥0.07 M _⊙ . The position–velocity map of the ^12 CO emission reveals previously unnoticed attributes of the older outer ejecta. In a very narrow range of Doppler velocities, ^12 CO absorption and emission causes some of this outer material to be quite opaque. At those frequencies the inner structure is hidden and we see only emission from an extended outer region. This fact produces a conspicuous but illusory dark spot if one attempts to subtract the continuum in a normal way.

We present high-cadence photometric and spectroscopic observations of SN 2023axu, a classical Type II supernova with an absolute V -band peak magnitude of –17.2 ± 0.1 mag. SN 2023axu was discovered by the Distance Less Than 40 Mpc (DLT40) survey within 1 day of the last nondetection in the nearby galaxy NGC 2283 at 13.7 Mpc. We modeled the early light curve using a recently updated shock cooling model that includes the effects of line blanketing and found the explosion epoch to be MJD 59971.48 ± 0.03 and the probable progenitor to be a red supergiant. The shock cooling model underpredicts the overall UV data, which point to a possible interaction with circumstellar material. This interpretation is further supported by spectral behavior. We see a ledge feature around 4600 Å in the very early spectra (+1.1 and +1.5 days after the explosion), which can be a sign of circumstellar interaction. The signs of circumstellar material are further bolstered by the presence of absorption features blueward of H α and H β at day >40, which is also generally attributed to circumstellar interaction. Our analysis shows the need for high-cadence early photometric and spectroscopic data to decipher the mass-loss history of the progenitor.

Mass loss is a crucial process that affects the observational properties, evolution path, and fate of highly evolved stars. However, the mechanism of mass loss is still unclear, and the mass-loss rate (MLR) of red supergiant stars (RSGs) requires further research and precise evaluation. To address this, we utilized an updated and complete sample of RSGs in the Large Magellanic Cloud (LMC) and employed the 2-DUST radiation transfer model and spectral energy distribution fitting approach to determine the dust-production rates (DPRs) and dust properties of the RSGs. We have fitted 4714 selected RSGs with over 100,000 theoretical templates of evolved stars. Our results show that the DPR range of RSGs in the LMC is 10 ^−11 M _⊙ yr ^−1 –10 ^−7 M _⊙ yr ^−1 , and the total DPR of all RSGs is 1.14 ×10 ^−6 M _⊙ yr ^−1 . We find that 63.3% RSGs are oxygen-rich, and they account for 97.2% of the total DPR. The optically thin RSG, which comprise 30.6% of our sample, contribute only 0.1% of the total DPR, while carbon-rich RSGs (6.1%) produce 2.7% of the total DPR. Overall, 208 RSGs contributed 76.6% of the total DPR. We have established a new relationship between the MLR and luminosity of RSGs in the LMC, which exhibits a positive trend and a clear turning point at $\mathrm{log}L/{L}_{\odot }\approx 4.4$ .

In this work, for the first time, we use seismic data as well as surface abundances to model the supergiant α -Ori, with the goal of setting an upper bound on the axion–photon coupling constant g a γ . We find that, in general, stellar models with g a γ ∈ [0.002; 2.0] × 10−10 GeV−1 agree with the observational data, but beyond that upper limit, we do not find stellar models that are compatible with the observational constraints and the current literature. From g a γ = 3.5 × 10−10 GeV−1 on, the algorithm did not find any fitting models. Even so, all the axionic models considered present distinct internal profiles from the reference case, without axions. Moreover, as the axion energy losses become more significant, the behavior of the stellar models becomes more diversified, even with very similar input parameters. Nonetheless, the consecutive increments of g a γ still show systematic tendencies, resulting from the axion energy losses. Moreover, we establish three important conclusions: (1) the increased luminosity and higher neutrino production are measurable effects, possibly associated with axion energy losses; (2) stellar models with axion energy loss show a quite distinct internal structure; and (3) the importance of future asteroseismic missions in observing low-degree nonradial modes in massive stars is emphasized—as internal gravity waves probe the near-core regions, where axion effects are most intense. Thus, more seismic data will allow us to constrain g a γ better and to prove or dismiss the existence of axion energy loss inside massive stars. [ABSTRACT FROM AUTHOR]

The evolution of granulation is an important mechanism of the light variations of red supergiants (RSGs). Based on pure and complete samples of RSGs in the Magellanic Clouds, the mechanisms and characteristics of the granulation of RSGs are investigated based on time-series data. As predicted by the basic physical process of granulation and previous works, there are tight relations between granulation and stellar parameters of RSGs (i.e., the scaling relations). The scaling relations of RSGs provide a new method to infer stellar parameters by using the characteristic timescale and amplitude of granulations. Some faint sources deviate from the scaling relations, which may be due to the difference in the properties of the granulation of the RSGs before and after the blue loop or contamination by Mira variables. However, both of these possibilities suggest that the scaling relations of granulation is different among different types of stars. [ABSTRACT FROM AUTHOR]

We present the optical spectroscopic evolution of SN 2023ixf seen in subnight cadence spectra from 1.18 to 15 days after explosion. We identify high-ionization emission features, signatures of interaction with material surrounding the progenitor star, that fade over the first 7 days, with rapid evolution between spectra observed within the same night. We compare the emission lines present and their relative strength to those of other supernovae with early interaction, finding a close match to SN 2020pni and SN 2017ahn in the first spectrum and SN 2014G at later epochs. To physically interpret our observations, we compare them to CMFGEN models with confined, dense circumstellar material around a red supergiant (RSG) progenitor from the literature. We find that very few models reproduce the blended N iii ( λλ 4634.0,4640.6)/C iii ( λλ 4647.5,4650.0) emission lines observed in the first few spectra and their rapid disappearance thereafter, making this a unique diagnostic. From the best models, we find a mass-loss rate of 10 ^−3 –10 ^−2 M _⊙ yr ^−1 , which far exceeds the mass-loss rate for any steady wind, especially for an RSG in the initial mass range of the detected progenitor. These mass-loss rates are, however, similar to rates inferred for other supernovae with early circumstellar interaction. Using the phase when the narrow emission features disappear, we calculate an outer dense radius of circumstellar material R _CSM,out ≈ 5 × 10 ^14 cm, and a mean circumstellar material density of ρ = 5.6 × 10 ^−14 g cm ^−3 . This is consistent with the lower limit on the outer radius of the circumstellar material we calculate from the peak H α emission flux, R _CSM,out ≳ 9 × 10 ^13 cm.

We present the discovery of the Type II supernova SN 2023ixf in M101 and follow-up photometric and spectroscopic observations, respectively, in the first month and week of its evolution. Our discovery was made within a day of estimated first light, and the following light curve is characterized by a rapid rise (≈5 days) to a luminous peak ( M _V ≈ − 18.2 mag) and plateau ( M _V ≈ − 17.6 mag) extending to 30 days with a fast decline rate of ≈0.03 mag day ^−1 . During the rising phase, U − V color shows blueward evolution, followed by redward evolution in the plateau phase. Prominent flash features of hydrogen, helium, carbon, and nitrogen dominate the spectra up to ≈5 days after first light, with a transition to a higher ionization state in the first ≈2 days. Both the U − V color and flash ionization states suggest a rise in the temperature, indicative of a delayed shock breakout inside dense circumstellar material (CSM). From the timescales of CSM interaction, we estimate its compact radial extent of ∼(3–7) × 10 ^14 cm. We then construct numerical light-curve models based on both continuous and eruptive mass-loss scenarios shortly before explosion. For the continuous mass-loss scenario, we infer a range of mass-loss history with 0.1–1.0 M _⊙ yr ^−1 in the final 2−1 yr before explosion, with a potentially decreasing mass loss of 0.01–0.1 M _⊙ yr ^−1 in ∼0.7–0.4 yr toward the explosion. For the eruptive mass-loss scenario, we favor eruptions releasing 0.3–1 M _⊙ of the envelope at about a year before explosion, which result in CSM with mass and extent similar to the continuous scenario. We discuss the implications of the available multiwavelength constraints obtained thus far on the progenitor candidate and SN 2023ixf to our variable CSM models.

We present UV and/or optical observations and models of SN 2023ixf, a type II supernova (SN) located in Messier 101 at 6.9 Mpc. Early time ( flash ) spectroscopy of SN 2023ixf, obtained primarily at Lick Observatory, reveals emission lines of H i , He i/ii , C iv , and N iii/iv/v with a narrow core and broad, symmetric wings arising from the photoionization of dense, close-in circumstellar material (CSM) located around the progenitor star prior to shock breakout. These electron-scattering broadened line profiles persist for ∼8 days with respect to first light, at which time Doppler broadened the features from the fastest SN ejecta form, suggesting a reduction in CSM density at r ≳ 10 ^15 cm. The early time light curve of SN 2023ixf shows peak absolute magnitudes (e.g., M _u = −18.6 mag, M _g = −18.4 mag) that are ≳2 mag brighter than typical type II SNe, this photometric boost also being consistent with the shock power supplied from CSM interaction. Comparison of SN 2023ixf to a grid of light-curve and multiepoch spectral models from the non-LTE radiative transfer code CMFGEN and the radiation-hydrodynamics code HERACLES suggests dense, solar-metallicity CSM confined to r = (0.5–1) × 10 ^15 cm, and a progenitor mass-loss rate of $\dot{M}={10}^{-2}\,{M}_{\odot }$ yr ^−1 . For the assumed progenitor wind velocity of v _w = 50 km s ^−1 , this corresponds to enhanced mass loss (i.e., superwind phase) during the last ∼3–6 yr before explosion.

We present the densely sampled early light curve of the Type II supernova (SN) 2023ixf, first observed within hours of explosion in the nearby Pinwheel Galaxy (Messier 101; 6.7 Mpc). Comparing these data to recently updated models of shock-cooling emission, we find that the progenitor likely had a radius of 410 ± 10 R _⊙ . Our estimate is model dependent but consistent with a red supergiant. These models provide a good fit to the data starting about 1 day after the explosion, despite the fact that the classification spectrum shows signatures of circumstellar material around SN 2023ixf during that time. Photometry during the first day after the explosion, provided almost entirely by amateur astronomers, does not agree with the shock-cooling models or a simple power-law rise fit to data after 1 day. We consider the possible causes of this discrepancy, including precursor activity from the progenitor star, circumstellar interaction, and emission from the shock before or after it breaks out of the stellar surface. The very low luminosity (−11 mag > M > −14 mag) and short duration of the initial excess lead us to prefer a scenario related to prolonged emission from the SN shock traveling through the progenitor system.

The direct detection of core-collapse supernova (SN) progenitor stars is a powerful way of probing the last stages of stellar evolution. However, detections in archival Hubble Space Telescope images are limited to about one detection per year. Here, we explore whether we can increase the detection rate by using data from ground-based wide-field surveys. Due to crowding and atmospheric blurring, progenitor stars can typically not be identified in preexplosion images alone. Instead, we combine many pre-SN and late-time images to search for the disappearance of the progenitor star. As a proof of concept, we implement our search of ZTF data. For a few hundred images, we achieve limiting magnitudes of ∼23 mag in the g and r bands. However, no progenitor stars or long-lived outbursts are detected for 29 SNe within z ≤ 0.01, and the ZTF limits are typically several magnitudes less constraining than detected progenitors in the literature. Next, we estimate progenitor detection rates for the Legacy Survey of Space and Time (LSST) with the Vera C. Rubin telescope by simulating a population of nearby SNe. The background from bright host galaxies reduces the nominal LSST sensitivity by, on average, 0.4 mag. Over the 10 yr survey, we expect the detection of ∼50 red supergiant progenitors and several yellow and blue supergiants. The progenitors of Type Ib and Ic SNe will be detectable if they are brighter than −4.7 or −4.0 mag in the LSST i band, respectively. In addition, we expect the detection of hundreds of pre-SN outbursts depending on their brightness and duration.

We perform a comprehensive search for optical precursor emission at the position of SN 2023ixf using data from the DLT40, ZTF, and ATLAS surveys. By comparing the current data set with precursor outburst hydrodynamical model light curves, we find that the probability of a significant outburst within 5 yr of explosion is low, and the circumstellar material (CSM) ejected during any possible precursor outburst is likely smaller than ∼0.015 M _⊙ . By comparing to a set of toy models, we find that, if there was a precursor outburst, the duration must have been shorter than ∼100 days for a typical brightness of M _r ≃ −9 mag or shorter than 200 days for M _r ≃ −8 mag; brighter, longer outbursts would have been discovered. Precursor activity like that observed in the normal Type II SN 2020tlf ( M _r ≃ −11.5) can be excluded in SN 2023ixf. If the dense CSM inferred by early flash spectroscopy and other studies is related to one or more precursor outbursts, then our observations indicate that any such outburst would have to be faint and only last for days to months, or it occurred more than 5 yr prior to the explosion. Alternatively, any dense, confined CSM may not be due to eruptive mass loss from a single red supergiant progenitor. Taken together, the results of SN 2023ixf and SN 2020tlf indicate that there may be more than one physical mechanism behind the dense CSM inferred around some normal Type II supernovae.

Betelgeuse, the nearest red supergiant, dimmed to an unprecedented level in early 2020. The star emerged from this “Great Dimming” episode with its typical, roughly 400 days pulsation cycle halved, and a new dominant period of around 200 days. The dimming event has been attributed to a surface mass ejection, in which rising material drove shocks through the stellar atmosphere and expelled some material, partially obscuring the star as it formed molecules and dust. In this paper, we use hydrodynamic simulations to reveal the connections between Betelgeuse's vigorously convective envelope, the surface mass ejection, and the pulsation mode switching that ensued. An anomalously hot convective plume, generated rarely but naturally in the star's turbulent envelope, can rise and break free from the surface, powering an upwelling that becomes the surface mass ejection. The rising plume also breaks the phase coherence of the star's pulsation, causing the surface to keep expanding even as the deeper layers contract. This drives a switch from the 400 days fundamental mode of pulsation, in which the whole star expands and contracts synchronously, to the 200 days first overtone, where a radial node separates the interior and exterior of the envelope moving in opposite phase. We predict that the star's convective motions will damp the overtone oscillation and Betelgeuse will return to its previous, 400 days fundamental-mode pulsation in the next 5–10 yr. With its resolved surface and unprecedentedly detailed characterization, Betelgeuse opens a window to episodic surface mass ejection in the late-stage evolution of massive stars.

We find that the convective motion in the envelopes of red supergiant (RSG) stars supplies a non-negligible stochastic angular momentum to the mass that a secondary star accretes in a common envelope evolution (CEE), such that jets that the secondary star launches wobble. The orbital motion of the secondary star in a CEE and the density gradient in the envelope impose a nonzero angular momentum to the accreted mass with a constant direction parallel to the orbital angular momentum. From one-dimensional stellar evolution simulations with the numerical code mesa we find that the stochastic convection motion in the envelope of RSG stars adds a stochastic angular momentum component with an amplitude that is about 0.1–1 times that of the constant component due to the orbital motion. We mimic a CEE of the RSG star by removing envelope mass at a high rate and by depositing energy into its envelope. The stochastic angular momentum implies that the accretion disk around the secondary star (which we do not simulate), and therefore the jets that it launches, wobble with angles of up to tens of degrees with respect to the orbital angular momentum axis. This wobbling makes it harder for jets to break out from the envelope and can shape small bubbles in the ejecta that compress filaments that appear as arcs in the ejected nebula, i.e., in planetary nebulae when the giant is an asymptotic giant branch star.

We present an extensive, panchromatic photometric (UV, optical, and near-IR) and low-resolution optical spectroscopic coverage of a Type IIP supernova SN 2018gj that occurred on the outskirts of the host galaxy NGC 6217. From the V- band light curve, we estimate the plateau length to be ∼ 70 ± 2 days, placing it among the very few well-sampled short plateau supernovae (SNe). With V -band peak absolute magnitude M _V ≤ −17.0 ± 0.1 mag, it falls in the middle of the luminosity distribution of the Type II SNe. The color evolution is typical to other Type II SNe except for an early elbow-like feature in the evolution of V − R color owing to its early transition from the plateau to the nebular phase. Using the expanding photospheric method, we present an independent estimate of the distance to SN 2018gj. We report the spectral evolution to be typical of a Type II SNe. However, we see a persistent blueshift in emission lines until the late nebular phase, not ordinarily observed in Type II SNe. The amount of radioactive nickel ( ^56 Ni) yield in the explosion was estimated to be 0.026 ± 0.007 M _⊙ . We infer from semianalytical modeling, nebular spectrum, and 1D hydrodynamical modeling that the probable progenitor was a red supergiant with a zero-age-main-sequence mass ≤13 M _⊙ . In the simulated hydrodynamical model light curves, reproducing the early optical bolometric light curve required an additional radiation source, which could be the interaction with the proximal circumstellar matter.

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

In this work, for the first time, we use seismic data as well as surface abundances to model the supergiant α -Ori, with the goal of setting an upper bound on the axion–photon coupling constant g _a _γ . We find that, in general, stellar models with g _a _γ ∈ [0.002; 2.0] × 10 ^−10 GeV ^−1 agree with the observational data, but beyond that upper limit, we do not find stellar models that are compatible with the observational constraints and the current literature. From g _a _γ = 3.5 × 10 ^−10 GeV ^−1 on, the algorithm did not find any fitting models. Even so, all the axionic models considered present distinct internal profiles from the reference case, without axions. Moreover, as the axion energy losses become more significant, the behavior of the stellar models becomes more diversified, even with very similar input parameters. Nonetheless, the consecutive increments of g _a _γ still show systematic tendencies, resulting from the axion energy losses. Moreover, we establish three important conclusions: (1) the increased luminosity and higher neutrino production are measurable effects, possibly associated with axion energy losses; (2) stellar models with axion energy loss show a quite distinct internal structure; and (3) the importance of future asteroseismic missions in observing low-degree nonradial modes in massive stars is emphasized—as internal gravity waves probe the near-core regions, where axion effects are most intense. Thus, more seismic data will allow us to constrain g _a _γ better and to prove or dismiss the existence of axion energy loss inside massive stars.

Precise interstellar dust extinction laws are important to infer the intrinsic properties of reddened objects and correctly interpret observations. In this work, we attempt to measure the optical–infrared extinction laws of the Large and Small Magellanic Clouds (LMC and SMC) by using red supergiant (RSG) stars and classical Cepheids as extinction tracers. The spectroscopic RSG samples are constructed based on the APOGEE spectral parameters, Gaia astrometric data, and multiband photometry. We establish the effective temperature–intrinsic color relations for RSG stars and determine the color excess ratio (CER) E ( λ − G _RP )/ E ( G _BP − G _RP ) for the LMC and SMC. We use classical Cepheids to derive base relative extinction ${A}_{{G}_{\mathrm{RP}}}/E({G}_{\mathrm{BP}}-{G}_{\mathrm{RP}})$ . The results are 1.589 ± 0.014 and 1.412 ± 0.041 for the LMC and SMC. By combining CERs with ${A}_{{G}_{\mathrm{RP}}}/E({G}_{\mathrm{BP}}-{G}_{\mathrm{RP}})$ , the optical–infrared extinction coefficients ${A}_{\lambda }/{A}_{{G}_{\mathrm{RP}}}$ are determined for 16 bands. We adjust the parameters of R _V -dependent extinction laws and obtain average extinction laws for the LMC and SMC of R _V = 3.40 ± 0.07 and R _V = 2.53 ± 0.10, respectively, which are consistent with Gordon et al. In the optical bands, the adjusted R _V extinction curves agree with the observations with deviations less than 3%.