Young Clusters Research Articles

The role of radiation pressure in lithium deficit for stars in the stage of late accretion: failure of lithium test for young stellar clusters

ABSTRACT We study the behaviour of lithium atoms in the radiation field of main-sequence (MS) stars. The radiation pressure on atoms in the resonant lines of the ${\rm Li\, {\small I}}$ 6708 Å doublet exceeds the gravity at Teff > 3000 K. The lifetime of ${\rm Li}$ atoms before their ionization is sufficient to achieve the escape velocity. Free ${\rm Li}$ atoms cannot approach the pre-main sequence (PMS) and MS stars if the surrounding gas is sufficiently transparent. Accretion of small bodies (SBs) can partially restore the lithium abundance when SB directly falls on to the star in elongated orbits. In the case of accretion from discs consisting of SB fragments, when the parent bodies move in quasi-circular orbits, the boundary of the SB sublimation region plays a decisive role. Near hot stars, the free ${\rm Li}$ atom is rapidly ionized and participates in accretion on to the star in the form of ions, contributing to an increase in the ${\rm Li}$ abundance. Low-mass stars with low temperatures cannot hold lithium atoms. Near stars with Planck energy distribution in the radiation field, lithium deficiency occurs if 3000 K < Teff < 3400 K. New stars with [Fe/H] = 0.0 may indicate a lithium deficiency at Teff < 5000 K. Failure of the lithium test to estimate the age of a stellar cluster may indicate a significant contribution of accretion in the form of debris discs consisting of SBs.

A statistical and multiwavelength photometric analysis of a young embedded open star cluster: IC 1590

ABSTRACT We present a statistical and multiwavelength photometric studies of young open cluster IC 1590. We identified 91 cluster members using Gaia DR3 astrometry data using ensemble-based unsupervised machine learning techniques. From Gaia EDR3 data, we estimate the best-fitting parameters for IC 1590 using the Automated Stellar Cluster Analysis package (asteca) yielding the distance d ∼ 2.87 ± 0.02 kpc, age ∼ 3.54 ± 0.05 Myr, metallicity z ∼ 0.0212 ± 0.003, binarity value of ∼ 0.558, and extinction Av ∼ 1.252 ± 0.4 mag for an Rv value of ∼ 3.322 ± 0.23. We estimate the initial mass function slope of the cluster to be α = 1.081 ± 0.112 for single stars and α = 1.490 ± 0.051 for a binary fraction of ∼ 0.558 in the mass range 1 M⊙ ≤ m (M⊙) ≤ 100 M⊙. The G-band luminosity function slope is estimated to be ∼ 0.33 ± 0.09. We use (J − H) versus (H − Ks) colour–colour diagram to identify young stellar objects (YSOs). We found that all the identified YSOs have ages ≤ 2 Myr and masses ∼ 0.35 – 5.5 M⊙. We also fit the radial surface density profile. Using the galpy, we performed orbit analysis of the cluster. The extinction map for the cluster region has been generated using the PNICER technique, and it is almost similar to the dust structure obtained from the 500 μm dust continuum emissions map of Herschel SPIRE. We finally at the end discussed the star formation scenario in the cluster region.

PHANGS Hubble Space Telescope Treasury Survey: Globular Cluster Systems in 17 Nearby Spiral Galaxies

We present new catalogs of likely globular clusters (GCs) in 17 nearby spiral galaxies studied as part of the PHANGS Hubble Space Telescope (HST) Treasury Survey. The galaxies were imaged in five broadband filters from the near-ultraviolet through the I band. PHANGS-HST has produced catalogs of stellar clusters of all ages by selecting extended sources (from multiple concentration index measurements) followed by morphological classification (centrally concentrated and symmetric or asymmetric, multiple peaks, and contaminants) by visually examining the V-band images and separately by a machine-learning algorithm which classified larger samples to reach fainter limits. From both cluster catalogs, we select an initial list of candidate GCs which have B − V ≥ 0.5 and V − I ≥ 0.73 mag, then remove likely contaminants (including reddened young clusters, background galaxies misclassified by the neural network, and chance superpositions/blends of stars) after a careful visual inspection. We find that ≈86% of the color-selected candidates classified as spherically symmetric and 68% of those classified as centrally concentrated but asymmetric are likely to be GCs. The luminosity functions of the GC candidates in two of our 17 galaxies, NGC 628 and NGC 3627, are atypical, and continue to rise at least 1 mag fainter than the expected turnover near M V ∼ −7.4. These faint candidate GCs have more extended spatial distributions than their bright counterparts, and may reside in the disk rather than the bulge/halo, similar to faint GCs previously discovered in M101. These faint clusters may be somewhat younger since the age–metallicity degeneracy makes it difficult to determine precise cluster ages from integrated colors once they reach ≈1 Gyr.

Search for brown dwarfs in IC 1396 with Subaru HSC: interpreting the impact of environmental factors on substellar population

ABSTRACT Young stellar clusters are predominantly the hub of star formation and hence, ideal to perform comprehensive studies over the least explored substellar regime. Various unanswered questions like the mass distribution in brown dwarf regime and the effect of diverse cluster environment on brown dwarf formation efficiency still plague the scientific community. The nearby young cluster, IC 1396 with its feedback-driven environment, is ideal to conduct such study. In this paper, we adopt a multiwavelength approach, using deep Subaru HSC along with other data sets and machine learning techniques to identify the cluster members complete down to ∼ 0.03 M⊙ in the central 22 arcmin area of IC 1396. We identify 458 cluster members including 62 brown dwarfs which are used to determine mass distribution in the region. We obtain a star-to-brown dwarf ratio of ∼ 6 for a stellar mass range 0.03–1 M⊙ in the studied cluster. The brown dwarf fraction is observed to increase across the cluster as radial distance from the central OB-stars increases. This study also compiles 15 young stellar clusters to check the variation of star-to-brown dwarf ratio relative to stellar density and ultraviolet (UV) flux ranging within 4–2500 stars pc−2 and 0.7–7.3 G0, respectively. The brown dwarf fraction is observed to increase with stellar density but the results about the influence of incident UV flux are inconclusive within this range. This is the deepest study of IC 1396 as of yet and it will pave the way to understand various aspects of brown dwarfs using spectroscopic observations in future.

A survey for variable young stars with small telescopes – VIII. Properties of 1687 Gaia selected members in 21 nearby clusters

ABSTRACT The Hunting Outbursting Young Stars (HOYS) project performs long-term, optical, multifilter, high cadence monitoring of 25 nearby young clusters and star-forming regions. Utilizing Gaia DR3 data, we have identified about 17 000 potential young stellar members in 45 coherent astrometric groups in these fields. Twenty one of them are clear young groups or clusters of stars within 1 kpc and they contain 9143 Gaia selected potential members. The cluster distances, proper motions, and membership numbers are determined. We analyse long-term (≈ 7 yr) V-, R-, and I-band light curves from HOYS for 1687 of the potential cluster members. One quarter of the stars are variable in all three optical filters, and two-thirds of these have light curves that are symmetric around the mean. Light curves affected by obscuration from circumstellar materials are more common than those affected by accretion bursts, by a factor of 2–4. The variability fraction in the clusters ranges from 10 per cent to almost 100 per cent, and correlates positively with the fraction of stars with detectable inner discs, indicating that a lot of variability is driven by the disc. About one in six variables shows detectable periodicity, mostly caused by magnetic spots. Two-thirds of the periodic variables with disc excess emission are slow rotators, and amongst the stars without disc excess two-thirds are fast rotators – in agreement with rotation being slowed down by the presence of a disc.

Tracing the Origins of Mass Segregation in M35: Evidence for Primordially Segregated Binaries

M35 is a young open cluster and home to an extensive binary population. Using Gaia Data Release 3, Pan-STARRS, and Two Micron All Sky Survey photometry with the Bayesian statistical software, BASE-9, we derive precise cluster parameters, identify single and binary cluster members, and extract their masses. We identify 571 binaries down to Gaia G = 20.3 and a lower limit on the binary frequency of f b = 0.41 ± 0.02. We extend the binary demographics by many magnitudes faint-ward of previous (radial-velocity) studies of this cluster and further away from the cluster center (1.°78, roughly 10 core radii). We find the binary stars to be more centrally concentrated than the single stars in the cluster. Furthermore, we find strong evidence for mass segregation within the binary population itself, with progressively more-massive binary samples becoming more and more centrally concentrated. For the single stars, we find weaker evidence for mass segregation; only the most massive single stars (>2.5M ⊙) appear more centrally concentrated. Given the cluster age of ∼200 Myr, and our derived half-mass relaxation time for the cluster of 230 ± 84 Myr, we estimate ∼47% of the binary stars and ∼12% of single stars in the cluster have had time to become dynamically mass segregated. Importantly, when we investigate only stars with mass segregation timescales greater than the cluster age, we still find the binaries to be more centrally concentrated than the singles, suggesting the binaries may have formed with a primordially different spatial distribution from that of the single stars.

EWOCS-I: The catalog of X-ray sources in Westerlund 1 from the Extended Westerlund 1 and 2 Open Clusters Survey

Context. With a mass exceeding several 104 M⊙ and a rich and dense population of massive stars, supermassive young star clusters represent the most massive star-forming environment that is dominated by the feedback from massive stars and gravitational interactions among stars. Aims. In this paper we present the Extended Westerlund 1 and 2 Open Clusters Survey (EWOCS) project, which aims to investigate the influence of the starburst environment on the formation of stars and planets, and on the evolution of both low and high mass stars. The primary targets of this project are Westerlund 1 and 2, the closest supermassive star clusters to the Sun. Methods. The project is based primarily on recent observations conducted with the Chandra and JWST observatories. Specifically, the Chandra survey of Westerlund 1 consists of 36 new ACIS-I observations, nearly co-pointed, for a total exposure time of 1 Msec. Additionally, we included 8 archival Chandra/ACIS-S observations. This paper presents the resulting catalog of X-ray sources within and around Westerlund 1. Sources were detected by combining various existing methods, and photon extraction and source validation were carried out using the ACIS-Extract software. Results. The EWOCS X-ray catalog comprises 5963 validated sources out of the 9420 initially provided to ACIS-Extract, reaching a photon flux threshold of approximately 2 × 10−8 photons cm−2 s−1. The X-ray sources exhibit a highly concentrated spatial distribution, with 1075 sources located within the central 1 arcmin. We have successfully detected X-ray emissions from 126 out of the 166 known massive stars of the cluster, and we have collected over 71 000 photons from the magnetar CXO J164710.20-455217.

Constraining Stellar Rotation at the Zero-age Main Sequence with TESS

The zero-age main sequence (ZAMS) is a critical phase for stellar angular momentum evolution, as stars transition from contraction-dominated spin-up to magnetic wind-dominated spin-down. We present the first robust observational constraints on rotation for FGK stars at ≈40 Myr. We have analyzed TESS light curves for 1410 members of five young open clusters with ages between 25 and 55 Myr: IC 2391, IC 2602, NGC 2451A, NGC 2547, and Collinder 135. In total, we measure 868 rotation periods, including 96 new, high-quality periods for stars around 1 M ⊙. This is an increase of ten times the existing literature sample at the ZAMS. We then use the τ 2 method to compare our data to models for stellar angular momentum evolution. Although the ages derived from these rotation models do not match isochronal ages, we show that these observations can clearly discriminate between different models for stellar wind torques. Finally, τ 2 fits indicate that magnetic braking and/or internal angular momentum transport significantly impact rotational evolution even on the pre-main sequence.

Magnetic Fields in the Central Molecular Zone Influenced by Feedback and Weakly Correlated with Star Formation

Magnetic fields of molecular clouds in the central molecular zone (CMZ) have been relatively under-observed at sub-parsec resolution. Here, we report JCMT/POL2 observations of polarized dust emission in the CMZ, which reveal magnetic field structures in dense gas at ∼0.5 pc resolution. The 11 molecular clouds in our sample include two in the western part of the CMZ (Sgr C and a farside cloud candidate), four around the Galactic longitude 0 (the 50 km s−1 cloud, CO 0.02−0.02, the Stone, and the Sticks and Straw among the Three Little Pigs), and five along the Dust Ridge (G0.253+0.016, clouds b, c, d, and e/f), for each of which we estimate the magnetic field strength using the angular dispersion function method. The morphologies of magnetic fields in the clouds suggest potential imprints of feedback from expanding H ii regions and young massive star clusters. A moderate correlation between the total viral parameter versus the star formation rate (SFR) and the dense gas fraction of the clouds is found. A weak correlation between the mass-to-flux ratio and the SFR, and a weak anticorrelation between the magnetic field and the dense gas fraction are also found. Comparisons between magnetic fields and other dynamic components in clouds suggest a more dominant role of self-gravity and turbulence in determining the dynamical states of the clouds and affecting star formation at the studied scales.

Between the Cosmic-Ray “Knee” and the “Ankle”: Contribution from Star Clusters

We show that massive, young star clusters may be possible candidates that can accelerate Galactic cosmic rays (CRs) in the range of 107–109 GeV (between the “knee” and “ankle”). Various plausible scenarios, such as acceleration at the wind termination shock and supernova shocks inside these young star clusters, have been proposed, since it is difficult to accelerate particles up to the 107–109 GeV range in the standard paradigm of CR acceleration in supernova remnants. We consider a model for the production of different nuclei in CRs from massive stellar winds using the observed distribution of young star clusters in the Galactic plane. We present a detailed calculation of CR transport in the Galaxy, taking into account the effect of diffusion, interaction losses during propagation, and particle reacceleration by old supernova remnants to determine the all-particle CR spectrum. Using the maximum energy estimate from the Hillas criterion, we argue that a young, massive star cluster can accelerate protons up to a few tens of PeV. Upon comparison with the observed data, our model requires a CR source spectrum with an exponential cutoff of 5 × 107 Z GeV (50 Z PeV) from these clusters, together with a CR injection fraction of ∼5% of the wind kinetic energy. We discuss the possibility of achieving these requirements in star clusters, as well as the associated uncertainties, in the context of considering star clusters as the natural accelerator of the “second component” of Galactic CRs.

Testing He ii Emission from Wolf–Rayet Stars as a Dust Attenuation Measure in Eight Nearby Star-forming Galaxies

The ability to determine galaxy properties such as masses, ages, and star formation rates robustly is critically limited by the ability to measure dust attenuation accurately. Dust reddening is often characterized by comparing observations to models of either nebular recombination lines or the UV continuum. Here, we use a new technique to measure dust reddening by exploiting the He ii λ1640 and λ4686 emission lines originating from the stellar winds of Wolf–Rayet stars. The intrinsic line ratio is determined by atomic physics, enabling an estimate of the stellar reddening similar to how the Balmer lines probe gas-emission reddening. The He ii line ratio is measured from UV and optical spectroscopy using the Space Telescope Imaging Spectrograph on board the Hubble Space Telescope for eight nearby galaxies hosting young massive star clusters. We compare our results to dust reddening values estimated from UV spectral slopes and from Balmer line ratios and find tentative evidence for systematic differences. The reddening derived from the He ii lines tends to be higher, whereas that from the UV continuum tends to be lower. A larger sample size is needed to confirm this trend. If confirmed, this may indicate an age sequence probing different stages of dust clearing. Broad He ii lines have also been detected in galaxies more distant than in our sample, providing the opportunity to estimate the dust reddening of the youngest stellar populations out to distances of ∼100 Mpc.

On the incidence of episodic accretion in Class I YSOs from VVV

ABSTRACT Episodic accretion is one of the competing models to explain the observed luminosity spread in young stellar clusters. These short-lived high accretion events could also have a strong impact on planet formation. Observations of high-amplitude variability in young stellar objects (YSOs) due to large changes in the accretion rate provide direct observational evidence for episodic accretion. However, there are still uncertainties in the frequency of these events and if episodic accretion is universal among YSOs. To determine the frequency of outbursts in Class I YSOs, we built a large and robust sample of objects at this evolutionary stage, and searched for high-amplitude near-infrared (ΔKS > 2 mag) variability in the VIRAC2 database of the Vista Variables in the Via Lactea survey. By complementing with near-IR (2MASS and DENIS) and mid-IR (WISE/Neo-WISE) data, we find that from ∼7000 Class I YSOs, 97 objects can be classified as eruptive variable YSOs. The duration of the outbursts vary from a few months to longer than 9 yr, and cover a similar range of amplitudes. Values of ΔKS > 5 mag, however, are only observed in outbursts with duration longer than 9 yr. When considering different effects of completeness and contamination, we estimate that the incidence of episodic accretion in Class I YSOs is between 2 and 3 per cent. Finally, we determine a recurrence time-scale of long-term outbursts (a.k.a FUors) of $\tau =1.75^{+1.12}_{-0.87}$ kyr. The latter value agrees with previous estimates and is in line with the expectations of higher frequency of FUor outbursts during younger stages of evolution.

PeV proton acceleration in gamma-ray binaries

Current generation of ground based gamma-ray telescopes observed dozens of sources of photons above 100 TeV. Supernova remnants, pulsar wind nebulae, young stellar clusters and superbubbles are considered as possible sites of PeV-regime particles producing the radiation. Another possible source of PeV particles could be gamma-ray binary systems. In these systems, a strong relativistic outflow from a compact object (neutron star or black hole) collides with the dense wind from a massive companion early-type star. Gamma-ray binaries are observed from radio to high energy gamma-rays as luminous non-thermal sources. Apart from acceleration of very high energy leptons producing most of the non-thermal radiation, these systems may also efficiently accelerate protons. We present here the results of numerical simulation of the PeV-regime proton acceleration in gamma-ray binaries. The simulation is based on relativistic MHD modeling of local flows of magnetized plasma in the region of interaction of two colliding winds. We then inject 0.1 PeV protons into the system and directly follow their trajectories to demonstrate that they are accelerated to energies above PeV. High magnetization of the wind of the young massive star providing a Gauss range field in the winds interaction region is of paramount importance for the acceleration of protons above PeV. The maximum energies of protons accelerated by colliding winds in gamma-ray binaries can significantly exceed the energy of the pulsar potential’s drop, which limits from above the energy of particles accelerated by an isolated pulsar.

The early evolution of young massive clusters

Context. Characterising the outcome of the star formation process is key to understand and predict the evolution of stellar populations. Especially the fraction of massive stars in young stellar clusters is of importance as they are the dominant sources of both mechanical and radiative feedback, strongly influencing the thermal and dynamical state of their birth environments, and beyond. Their supernovae may trigger the formation of new generations of stars in neighbouring regions. It turns out that a significant fraction of massive stars escape from their parent cluster via dynamical interactions of single stars and/or multiple stellar systems. Aims. M 17 is the nearest giant H II region hosting a very young and massive cluster: NGC 6618. Our aim is to identify stars brighter than G ≲ 21 mag that belong to NGC 6618, including the (massive) stars that may have escaped since its formation, and to determine the cluster distance and age. Methods. The Gaia DR3 database was used to identify members of NGC 6618 based on parallax and proper motion within 9′ from the cluster centre. We searched for nearby stars in a field of 5° around the cluster centre that may have originated from the cluster, and we determined their transverse velocity, kinematic age, and impact parameter. Results. We identified 42 members of NGC 6618 of which eight have a spectral type of O, with a mean distance of 1675−18+19 pc and a (transversal) velocity dispersion of about 3 km s−1, and a radial velocity dispersion of ∼6 km s−1. Another ten O stars are associated with NGC 6618, but they cannot be classified as members due to poor astrometry and/or high extinction. We have also identified six O star runaways. The relative transverse velocity of these runaways ranges from 10 to 70 km s−1 and their kinematic age ranges from about 100 to 750 kyr. Given the already established young age of NGC 6618 (≲1 Myr), this implies that massive stars are being ejected from the cluster already directly after (or during) the cluster formation process. Conclusions. When constructing the initial mass function, one has to take into account the massive stars that have already escaped from the cluster, that is, about 30% of the O stars of the original population of NGC 6618. The trajectories of the O runaways can be traced back to the central 0.2–0.3 pc region of NGC 6618. The good agreement between the evolutionary and kinematic age of the runaways implies that the latter provides an independent way to estimate (a lower limit to) the age of the cluster.

Keck/KCWI spectroscopy of globular clusters in local volume dwarf galaxies

Abstract A number of nearby dwarf galaxies have globular cluster (GC) candidates that require spectroscopic confirmation. Here, we present Keck telescope spectra for 15 known GCs and GC candidates that may be associated with a host dwarf galaxy and an additional 3 GCs in the halo of M31 that are candidates for accretion from a now-disrupted dwarf galaxy. We confirm six star clusters (of intermediate-to-old age) to be associated with NGC 247. The vast bulk of its GC system remains to be studied spectroscopically. We also confirm the GC candidates in F8D1 and DDO190, finding both to be young star clusters. The three M31 halo GCs all have radial velocities consistent with M31 and are old and very metal-poor. Their ages and metallicities are consistent with accretion from a low-mass satellite galaxy. Finally, three objects are found to be background galaxies – two are projected near NGC 247 and one (candidate GCC7) is near the IKN dwarf. The IKN dwarf thus has only five confirmed GCs but still a remarkable specific frequency of 124.

Age-dating the young open cluster UBC 1 with g-mode asteroseismology, gyrochronology, and isochrone fitting

Aims. UBC 1 is an open cluster discovered in Gaia data and located near the edge of the Transiting Exoplanet Survey Satellite’s (TESS) continuous viewing zone. We aim to provide age constraints for this poorly studied open cluster from the combination of gravity-mode (g-mode) asteroseismology, gyrochronology, and isochrone fitting. Methods. We established the members of UBC 1 from a spatial-kinematic filtering and estimate the cluster age and its parameters. Firstly, we fitted rotating isochrones to the single star cluster sequence. Secondly, using TESS time-series photometry, we explored the variability of the upper main sequence members and identified potential g-mode pulsators. For one star, we found a clear period spacing pattern that we used to deduce the buoyancy travel time, the near-core rotation rate, and an asteroseismic age. For a third independent age estimate, we employed the rotation periods of low-mass members of UBC 1. Results. Based on isochrone fitting, we find log t = 8.1 ± 0.4, where the large uncertainty occurs because UBC 1 does not host evolved stars. From asteroseismology of one g-mode pulsator, we find a constrained age of log t = 8.24−0.14+0.43. From gyrochronology based on 17 cool star cluster members, we estimate log t = 8.35−0.25+0.16. Combined, all three methods lead to a consistent age in the range of 150 − 300 Myr. Conclusions. Our results show that even a single cluster member with identified g modes can improve age-dating of young open clusters. Combining the gyrochronology of low-mass members with asteroseismology of intermediate-mass members is a powerful tool for young open cluster modelling, including high-precision age-dating.

The Star Formation History of the Young Open Cluster NGC 2244: Study of the Binary Fraction

A star formation history (SFH) is the stellar record of the various processes that shaped a galaxy and star clusters since its birth to its current physical state. Studying the young open cluster NGC 2244, and in particular its massive stellar population, can improve our understanding of its formation mechanism, stellar evolution, feedback, and provide a basis for a better understanding of the formation of extragalactic star clusters. In this contribution, we present the first steps towards the derivation of the SFH of NGC 2244. We discuss the first results concerning the membership determination and the observed binary fractions of O- and B-type stars among NGC 2244. Finally, an overview of the forthcoming steps is also presented.

Searching for stellar population in the molecular cloud GRSMC 045.49+00.04

Abstract Understanding the characteristics of young stellar populations is essential for deriving insights into star formation processes within parent molecular clouds and the influence of massive stars on these processes. This study primarily aims to investigate the young stellar objects (YSOs) within the molecular cloud G 045.49+00.04, including three ultra-compact HII (UC HII) regions: G 45.48+0.13 (IRAS 19117+1107), G 45.45+0.06 (IRAS 19120+1103), and G 45.47+0.05. We used near-, mid-, and far-infrared photometric data along with radiation transfer models and the modified blackbody fitting to identify and study the YSOs and the interstellar medium (ISM). In total, we identified 1482 YSOs in a 12 arcmin radius covering GRSMC 045.49+00.04, with a mass range from 1.5 M ${}_{\odot}$ to 22 M ${}_{\odot}$ . Of these, 315 objects form relatively dense clusters in the UC HII regions, close to the IRAS 19120+1103 and 19117+1107 sources. In each UC HII region, several high-mass stars have been identified, which in all likelihood are responsible for the ionization. The YSOs with 21.8 M ${}_{\odot}$ and 13.7 ± 0.4 M ${}_{\odot}$ are associated with IRAS 19120+1103 and 19117+1107, respectively. The non-cluster YSOs (1168) are uniformly distributed on the field. The distribution of YSOs from both samples on the colour-magnitude diagram and by the evolutionary ages is different. About 75% of objects in the IRAS clusters are concentrated around the Zero Age Main Sequence and have a well-defined peak at an age of Log(Age[years]) $\approx$ 6.75, with a narrow spread. The non-cluster objects have two concentrations located to the right and left of the 0.1 Myr isochrone and two well-defined peaks at Log(Age) $\approx$ 6.25 and 5.25. The fraction of the near-infrared excess stars, as well as the mass function confirm that the evolutionary age of the cluster is about 1 Myr. The K luminosity functions’ α slopes for the IRAS clusters and non-cluster objects are 0.32 ± 0.04 and 0.72 ± 0.13, respectively. The steeper α slope is suggesting that the non-cluster objects are less evolved, which is well consistent with the evolutionary age. Similar results – including evolutionary age, narrow age spread, and the less evolved nature of non-cluster objects – were also observed for the YSOs in the neighbouring G 45.14+00.14. The both regions (G 045.49+00.04 and G 45.14+00.14) are located and distinguished by their brightness and density at the edge of the bubble around the highly variable X-ray binary GRS 1915+105, which includes a black hole and a K-giant companion. Based on the above, we can assume that the process of star formation in the young IRAS clusters was triggered by the GRS 1915+105-initiated shock front inside the ISM massive condensation, through the process of ‘collecting and collapse’. Most non-cluster objects probably belong to a later generation. Their formation could be triggered by the recurrent activity of GRS 1915+105 and/or through the edge collapse scenario and mass accumulation through the gas flows along the ISM filaments.

ALMA detection of CO rotational line emission in red supergiant stars of the massive young star cluster RSGC1

Context. The fate of stars largely depends on the amount of mass lost during the end stages of evolution. For single stars with an initial mass between ∼8–30 M⊙, most mass is lost during the red supergiant (RSG) phase, when stellar winds deplete the H-rich envelope. However, the RSG mass-loss rate (Ṁ) is poorly understood theoretically, and so stellar evolution models rely on empirically derived mass-loss rate prescriptions. However, it has been shown that these empirical relations differ largely, with differences up to 2 orders of magnitude. Aims. We aim to derive a new mass-loss rate prescription for RSGs that is not afflicted with some uncertainties inherent in preceding studies. Methods. We have observed CO rotational line emission towards a sample of RSGs in the open cluster RSGC1 that all are of a similar initial mass. The ALMA CO(2–1) line detections allowed us to retrieve the gas mass-loss rates (ṀCO). In contrast to mass-loss rates derived from the analysis of dust spectral features (ṀSED), the data allowed us a direct determination of the wind velocity and no uncertain dust-to-gas correction factor was needed. Results. Five RSGs in RSGC1 have been detected in CO(2–1). The retrieved ṀCO values are systematically lower than ṀSED. Although only five RSGs in RSGC1 have been detected, the data allow us to propose a new mass-loss rate relation for M-type red supergiants with effective temperatures between ∼3200 and 3800 K that is dependent on the luminosity and initial mass, and that is valid during the phase where nuclear burning determines the evolution along the RSG branch. The new mass-loss rate relation is based on the new ṀCO values for the RSGs in RSGC1 and on prior ṀSED values for RSGs in four clusters, including RSGC1. The new Ṁ-prescription yields a good prediction for the mass-loss rate of some well-known Galactic RSGs that are observed in multiple CO rotational lines, including α Ori, μ Cep and VX Sgr. Moreover, there are indications that a stronger, potentially eruptive, mass-loss process is occurring during some fraction of the RSG lifetime, suggesting that RSGs might experience a phase change in mass loss leading to the wind mass-loss rate dominating the RSG evolution at that stage. Conclusions. Implementing a lower mass-loss rate in evolution codes for massive stars has important consequences as to the nature of their end-state. A reduction of the RSG mass-loss rate implies that quiescent RSG mass loss is not enough to strip a single star’s hydrogen-rich envelope. Upon core collapse such single stars would explode as RSGs. Mass-loss rates of order ∼6 times higher would be needed to strip the H-rich envelope and produce a Wolf-Rayet star while evolving back to the blue side of the Hertzsprung–Russell diagram. Future observations of a larger sample of RSGs in open clusters should allow a more stringent determination of the ṀCO–luminosity relation and a sharper diagnostic as to when the phase change in mass loss is occurring.

Multiwavelength detection of an ongoing FUOr-type outburst on a low-mass YSO

ABSTRACT During the pre-main-sequence evolution, Young Stellar Objects (YSOs) assemble most of their mass during the episodic accretion process. The rarely seen FUOr-type events are valuable laboratories to investigate the outbursting nature of YSOs. Here, we present multiwavelength detection of a high-amplitude eruptive source in the young open cluster VdBH 221 with an ongoing outburst, including optical to mid-infrared time series and near-infrared spectra. The initial outburst has an exceptional amplitude of >6.3 mag in Gaia and 4.6 mag in Ks, with a peak luminosity up to 16 L⊙ and a peak mass accretion rate of 1.4 × 10−5 M⊙ yr−1. The optical to infrared spectral energy distribution of this object is consistent with a low-mass star (0.2 M⊙) with a modest extinction (AV < 2 mag). A 100-d delay between optical and infrared rising stages is detected, suggesting an outside-in origin of the instability. The spectroscopic features of this object reveal a self-luminous accretion disc, very similar to FU Orionis, with a low line-of-sight extinction. Most recently, there has been a gradual increase in brightness throughout the wavelength range, possibly suggesting an enhancement of the mass accretion rate.