Young Massive Star Clusters Research Articles

Witnessing an extreme, highly efficient galaxy formation mode with resolved Lyman-α and Lyman-continuum emission

J1316+2614 at z = 3.613 is the UV-brightest (MUV = −24.7) and strongest Lyman continuum-emitting (fescLyC ≈ 90%) star-forming galaxy known; it also shows signatures of inflowing gas from its blue-dominated Lyα profile. We present high-resolution imaging with the Hubble Space Telescope (HST) and the Very Large Telescope (VLT) of the LyC, Lyα, rest-UV, and optical emission of J1316+2614. Detailed analysis of the LyC and UV light distributions reveals compact yet resolved profiles, with LyC and UV morphologies showing identical half-light radii of reff ≃ 220 pc. The continuum-subtracted Lyα emission, obtained with the HST ramp-filter FR551N, reveals an extended filamentary structure of ≃6.0 kpc oriented south to north with only residual flux within the stellar core, suggesting a Lyα ‘hole’. Our spectral energy distribution analysis shows that J1316+2614 is characterised by a young (5.7 ± 1.0 Myr), nearly un-obscured stellar population with a high star-formation rate (SFR = 898 ± 181 M⊙ yr−1) and a stellar mass of M⋆young = (4.8 ± 0.3) × 109 M⊙. Additionally, the spectral energy distribution analysis supports the absence of an underlying old stellar population (M⋆old ≤ 2.8 × 109 M⊙, 3σ). J1316+2614 presents remarkably high SFR and stellar mass surface densities of log(Σ SFR[M⊙ yr−1 kpc−2]) = 3.47 ± 0.11 and log(ΣM⋆[M⊙pc−2]) = 4.20 ± 0.06, respectively, which are among the highest observed in star-forming galaxies and are more typically observed in local young massive star clusters and globular clusters. Our findings indicate that J1316+2614 is a powerful, young, and compact starburst that is leaking a significant amount of LyC photons due to a lack of gas and dust within the starburst. We explored the conditions for gas expulsion using a simple energetic balance and find that, given the strong binding force in J1316+2614, a high star-formation efficiency (ϵSF ≥ 0.7) is necessary to explain the removal of gas and its exposed nature. Our results thus suggest a close link between high ϵSF and high fescLyC. This high efficiency can also naturally explain the remarkably high SFR, UV luminosity, and efficient mass growth of J1316+2614, which acquired at least 62% of its mass in the last 6 Myr. J1316+2614 may exemplify an intense, feedback-free starburst with a high ϵSF, similar to those proposed for UV-bright galaxies at high redshifts.

Detection of extended gamma-ray emission in the vicinity of Cl Danks 1 and 2

Abstract We report the detection of high-energy gamma-ray emission towards the G305 star-forming region. Using almost 15 years of observation data from Fermi Large Area Telescope, we detected an extended gamma-ray source in this region with a significance of ∼13σ. The gamma-ray radiation reveals a clear pion-bump feature and can be fitted with the power law parent proton spectrum with an index of −2.5. The total cosmic ray (CR) proton energy in the gamma-ray production region is estimated to be the order of $10^{49}\ \rm erg$. We further derived the CR radial distribution from both the gamma-ray emission and gas distribution and found it roughly obeys the 1/r type profile, though a constant profile is not ruled out. This is consistent with other similar systems and expected from the continuous injection of CRs by the central powerful young massive star cluster Danks 1 or Danks 2 in this region. Together with former detections of similar gamma-ray structures, such as Cygnus cocoon, Westerlund 1, Westerlund 2, NGC 3603, and W40, the detection supports the hypothesis that young massive star clusters are CR accelerators.

GOALS-JWST: Constraining the Emergence Timescale for Massive Star Clusters in NGC 3256

We present the results of a James Webb Space Telescope NIRCam and NIRSpec investigation into the young massive star cluster (YMC) population of NGC 3256, the most cluster-rich luminous infrared galaxy in the Great Observatories All Sky LIRG Survey. We detect 3061 compact YMC candidates with a signal-to-noise ratio ≥3 at F150W, F200W, and F335M. Based on yggdrasil stellar population models, we identify 116/3061 sources with F150W – F200W > 0.47 and F200W – F355M > −1.37 colors, suggesting that they are young (t ≤ 5 Myr), dusty (A V = 5−15), and massive (M ⊙ > 105). This increases the sample of dust-enshrouded YMCs detected in this system by an order of magnitude relative to previous Hubble Space Telescope studies. With NIRSpec integral field unit pointings centered on the northern and southern nucleus, we extract the Paα and 3.3 μm polycyclic aromatic hydrocarbon (PAH) equivalent widths for eight bright and isolated YMCs. Variations in both the F200W – F335M color and 3.3 μm PAH emission with the Paα line strength suggest a rapid dust clearing (<3−4 Myr) for the emerging YMCs in the nuclei of NGC 3256. Finally, with both the age and dust emission accurately measured, we use yggdrasil to derive the color excess (E(B − V)) for all eight YMCs. We demonstrate that YMCs with strong 3.3 μm PAH emission (F200W – F335M > 0) correspond to sources with E(B − V) > 3, which are typically missed in UV-optical studies. This underscores the importance of deep near-infrared imaging for finding and characterizing these very young and dust-embedded sources.

Using 26Al to detect ongoing self-enrichment in young massive star clusters

Abstract Self-enrichment is one of the leading explanations for chemical anomalies in globular clusters. In this scenario, various candidate polluter stars have been proposed to eject gas with altered chemical composition during the self-enrichment process. Most of the proposed polluters will also eject radioactive 26Al into the surroundings. Hence, any detection of 26Al in young massive star clusters (YMCs) would support the self-enrichment scenario if YMCs were indeed the progenitors of globular clusters. Observations of gamma-ray data from COMPTEL and INTEGRAL, as well as detections of 26AlF molecules by the Atacama Large Millimeter-submillimeter Array (ALMA), indicate the maturing of 26Al detection methods. Detection possibilities will be enhanced in the short- to mid-term by the upcoming launch of the Compton Spectrometer and Imager (COSI). The Square Kilometer Array (SKA) could in principle also detect radio recombination lines of the positronium formed from the decay products of 26Al. Here, we show for a sample of YMCs in the nearby Universe, where self-enrichment could plausibly take place. For some nearby galaxies, this could enhance 26Al by an order of one magnitude. Detecting 26AlF with ALMA appears feasible for many candidate self-enrichment clusters, although significant challenges remain with other detection methods. The Large Magellanic Cloud, with its YMC R136, stands out as the most promising candidate. Detecting a 1.8 MeV radioactive decay line of 26Al here would require at least 15 months of targeted observation with COSI, assuming ongoing self-enrichment in R136.

On the Correlation between Young Massive Star Clusters and Gamma-Ray Unassociated Sources

Star clusters (SCs) are potential cosmic-ray accelerators and therefore are expected to emit high-energy radiation. However, a clear detection of gamma-ray emission from this source class has only been possible for a handful of cases. This could in principle result from two different reasons: either detectable SCs are limited to a small fraction of the total number of Galactic SCs, or gamma-ray-emitting SCs are not recognized as such and therefore are listed in the ensemble of unidentified sources. In this Letter we investigate this latter scenario by comparing available catalogs of SCs and H ii regions, obtained from Gaia and Wide-field Infrared Survey Explorer observations, to the gamma-ray GeV and TeV catalogs built from Fermi Large Area Telescope (LAT), H.E.S.S., and LHAASO data. The significance of the correlation between catalogs is evaluated by comparing the results with simulations of synthetic populations. A strong correlation emerges between Fermi-LAT-unidentified sources and H ii regions that trace massive SCs in the earliest (≲1–2 Myr) phase of their life, where no supernova explosions have happened yet, confirming that winds of massive stars can alone accelerate particles and produce gamma-ray emission at least up to GeV energies. The association with TeV energy sources is less evident. Similarly, no significant association is found between Gaia SCs and GeV nor TeV sources. We ascribe this fact to the larger extension of these objects but also to an intrinsic bias in the Gaia selection toward SCs surrounded by a lower target gas density, which would otherwise hinder the detection in the optical wave band.

Very-high-energy γ-Ray Emission from Young Massive Star Clusters in the Large Magellanic Cloud

The Tarantula Nebula in the Large Magellanic Cloud is known for its high star formation activity. At its center lies the young massive star cluster R136, providing a significant amount of the energy that makes the nebula shine so brightly at many wavelengths. Recently, young massive star clusters have been suggested to also efficiently produce very high-energy cosmic rays, potentially beyond PeV energies. Here, we report the detection of very-high-energy γ-ray emission from the direction of R136 with the High Energy Stereoscopic System, achieved through a multicomponent, likelihood-based modeling of the data. This supports the hypothesis that R136 is indeed a very powerful cosmic-ray accelerator. Moreover, from the same analysis, we provide an updated measurement of the γ-ray emission from 30 Dor C, the only superbubble detected at TeV energies presently. The γ-ray luminosity above 0.5 TeV of both sources is (2–3) × 1035 erg s−1. This exceeds by more than a factor of 2 the luminosity of HESS J1646−458, which is associated with the most massive young star cluster in the Milky Way, Westerlund 1. Furthermore, the γ-ray emission from each source is extended with a significance of >3σ and a Gaussian width of about 30 pc. For 30 Dor C, a connection between the γ-ray emission and the nonthermal X-ray emission appears likely. Different interpretations of the γ-ray signal from R136 are discussed.

Cygnus OB2 as a test case for particle acceleration in young massive star clusters

In this paper, we focus on the scientific case of Cygnus OB2, a young massive star cluster (YMSC) located in the northern sky towards the Cygnus X star-forming complex. We consider a model that assumes cosmic-ray acceleration occurring only at the termination shock of the collective wind of the YMSC and address the question of whether or not, and under what hypotheses, hadronic emission by the accelerated particles can account for the observations of Cygnus OB2 obtained by Fermi-LAT and HAWC. To this end, we carefully review the available information on this source, also confronting different estimates of the relevant parameters with ad hoc developed simulations. Once other model parameters are fixed, the spectral and spatial properties of the emission are found to be very sensitive to the unknown properties of the turbulent magnetic field. Comparison with the data shows that our suggested scenario is incompatible with Kolmogorov turbulence. Assuming Kraichnan- or Bohm-type turbulence spectra, the model accounts well for the very high energy (VHE) data, but fails to reproduce the centrally peaked morphology observed by Fermi-LAT, suggesting that additional effects might be important for lower energy γ-ray emission. We discuss how additional progress can be made with more detailed and extended knowledge of the spectral and morphological properties of the emission.

GOALS-JWST: Mid-infrared Molecular Gas Excitation Probes the Local Conditions of Nuclear Star Clusters and the Active Galactic Nucleus in the LIRG VV 114

The enormous increase in mid-IR sensitivity and spatial and spectral resolution provided by the JWST spectrographs enables, for the first time, detailed extragalactic studies of molecular vibrational bands. This opens an entirely new window for the study of the molecular interstellar medium in luminous infrared galaxies (LIRGs). We present a detailed analysis of rovibrational bands of gas-phase CO, H2O, C2H2, and HCN toward the heavily obscured eastern nucleus of the LIRG VV 114, as observed by NIRSpec and the medium resolution spectrograph on the Mid-InfraRed Instrument (MIRI MRS). Spectra extracted from apertures of 130 pc in radius show a clear dichotomy between the obscured active galactic nucleus (AGN) and two intense starburst regions. We detect the 2.3 μm CO bandheads, characteristic of cool stellar atmospheres, in the star-forming regions, but not toward the AGN. Surprisingly, at 4.7 μm, we find highly excited CO (T ex ≈ 700–800 K out to at least rotational level J = 27) toward the star-forming regions, but only cooler gas (T ex ≈ 200 K) toward the AGN. We conclude that only mid-infrared pumping through the rovibrational lines can account for the equilibrium conditions found for CO and H2O in the deeply embedded starbursts. Here, the CO bands probe regions with an intense local radiation field inside dusty young massive star clusters or near the most massive young stars. The lack of high-excitation molecular gas toward the AGN is attributed to geometric dilution of the intense radiation from the bright point source. An overview of the relevant excitation and radiative transfer physics is provided in an appendix.

Is the RSGC4 (Alicante 8) Cluster a Real Star Cluster? Peculiar Radial Velocities of Red Supergiant Stars

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.

GeV γ-ray emission in the field of young massive star cluster RCW 38

ABSTRACT We report the detection of γ-ray emission by the Fermi Large Area Telescope towards the young massive star cluster RCW 38 in the 1–500 GeV photon energy range. We found spatially extended GeV emission towards the direction of RCW 38, which is best modelled by a Gaussian disc of 0.23° radius with a significance of the extension of ∼11.4σ. Furthermore, the spatial correlation with the ionized and molecular gas content favours the hadronic origin of the γ-ray emission. The γ-ray spectrum of RCW 38 has a relatively hard photon index of 2.44 ± 0.03, which is similar to other young massive star clusters. We argue that the diffuse GeV γ-ray emission in this region likely originates from the interaction of accelerated protons in the stellar cluster with the ambient gas.

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.

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.

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.

Diffuse gamma-ray emission around the Rosette Nebula

ABSTRACT The Rosette Nebula is a young stellar cluster and molecular cloud complex, located at the edge of the southern shell of a middle-aged supernova remnant Monoceros Loop (G205.5+0.5). We revisited the GeV gamma-ray emission towards the Rosette Nebula using more than 13 yr of Fermi-LAT data. We tested several spatial models and found that compared to the result using the CO gas template only, the inclusion of the H ii gas template can significantly improve the likelihood fit. We performed spectral analysis using the new spatial template. With both the gamma-ray observation and CO+H ii gas data, we derived the cosmic ray spectrum of different components in the vicinity of the Rosette Nebula. We found the gamma-ray emissions from Rosette Nebula are substantially harder than previously reported, which may imply that Rosette Nebula is another example of a gamma-ray emitting young massive star cluster.

Cosmic Ray Enrichment with 22Ne in Young Massive Star Clusters

—We developed the hypothesis that the 22Ne overabundance in cosmic rays (CRs) is generated in massive star clusters which contain populations of Wolf–Rayet stars. Winds of Wolf–Rayet stars are considered to have high content of 22Ne. We present a model of cosmic ray enrichment with 22Ne in young massive star clusters, adding isotopic yields from supernovae and considering the acceleration efficiency during the lifetime of the stars. The impact of the parameters (the initial mass function in the cluster, rotation velocity, black hole cut-off mass) is discussed. The luminosity of the sources is found to be sufficient to provide the observed CR flux. The local source hypothesis is investigated.

Clues on the presence and segregation of very massive stars in the Sunburst Lyman-continuum cluster at z = 2.37

We report the identification of very massive stars (VMS; mass &gt; 100 M⊙) that may be segregated in the center of the young massive star cluster at z = 2.37 hosted in the lensed galaxy called Sunburst galaxy. This result is based on two pieces of evidence: (1) VLT/MUSE spectra of several multiple images of the same star cluster show key spectral signatures of VMS, such as the He IIλ1640 broad emission, N IVλ1486 emission, and an N IVλ1720 P Cygni profile. In particular, He IIλ1640 is broad (∼1610 ± 300 km s−1), with an equivalent width of 3 Å, and asymmetric profile. These features require an extremely young (∼2.5 Myr) stellar population component in which the masses of the stars exceed 100 M⊙. When a Salpeter initial mass function and BPASS models for normal massive stars are assumed, the observed spectral features require ∼400 VMS. (2) The same star cluster is detected at a signal-to-noise ratio of ∼100 in the Lyman continuum domain (λ &lt; 900 Å). The Lyman continuum emission emerges from a region with a radius that is at least twice smaller than what is observed at 1700 Å (independently of magnification) and is located in the center of the cluster. After delensing, the effective radii in absolute scales are Reff[LyC] ∼ 4.7 ± 1.5 pc and Reff[1700] = 7.8 ± 1.4 pc. The Lyman continuum radiation is mainly produced by hot and massive stars, which implies that their spatial distribution (including that of VMS) is preferentially more confined in the central parts of the cluster. Approximately 400 VMS hosted by a cluster of ∼107 M⊙ produce ∼15% of the escaping Lyman continuum photons, and the remaining photons are produced by other massive early-type stars.

JWST/NIRCam Probes Young Star Clusters in the Reionization Era Sunrise Arc

Star cluster formation in the early universe and its contribution to reionization remains largely unconstrained to date. Here we present JWST/NIRCam imaging of the most highly magnified galaxy known at z ∼ 6, the Sunrise arc. We identify six young massive star clusters (YMCs) with measured radii spanning from ∼20 down to ∼1 pc (corrected for lensing magnification), estimated stellar masses of ∼106–7 M ⊙, and ages of 1–30 Myr based on SED fitting to photometry measured in eight filters extending to rest frame 7000 Å. The resulting stellar mass surface densities are higher than 1000 M ⊙ pc−2 (up to a few 105 M ⊙ pc−2), and their inferred dynamical ages qualify the majority of these systems as gravitationally bound stellar clusters. The star cluster ages map the progression of star formation along the arc, with two evolved systems (≳10 Myr old) followed by very young clusters. The youngest stellar clusters (<5 Myr) show evidence of prominent Hβ+[O iii] emission based on photometry with equivalent widths larger than >1000 Å rest frame and are hosted in a 200 pc sized star-forming complex. Such a region dominates the ionizing photon production with a high efficiency log(ξion[Hzerg−1])∼25.7 . A significant fraction of the recently formed stellar mass of the galaxy (10%–30%) occurred in these YMCs. We speculate that such sources of ionizing radiation boost the ionizing photon production efficiency, which eventually carves ionized channels that might favor the escape of Lyman continuum radiation. The survival of some of the clusters would make them the progenitors of massive and relatively metal-poor globular clusters in the local universe.

Understanding the TeV γ-ray emission surrounding the young massive star cluster Westerlund 1

Context. Young massive star clusters (YMCs) have increasingly become the focus of discussions on the origin of galactic cosmic rays (CRs). The proposition that CRs are accelerated inside superbubbles (SBs) blown by the strong winds of these clusters avoids issues faced by the standard paradigm of acceleration at supernova remnant shocks. Aims. We provide an interpretation of the latest TeV γ-ray observations of the region around the YMC Westerlund 1 taken with the High Energy Stereoscopic System (H.E.S.S.) in terms of diffusive shock acceleration at the cluster wind termination shock, taking the spectrum and morphology of the emission into account. As Westerlund 1 is a prototypical example of a YMC, this study is relevant to the general question about the role of YMCs for the Galactic CR population. Methods. We generated model γ-ray spectra, characterised particle propagation inside the SB based on the advection, diffusion, and cooling timescales, and constrained key parameters of the system. We considered hadronic emission from proton-proton interaction and subsequent pion decay and leptonic emission from inverse Compton scattering on all relevant photon fields, including the cosmic microwave background, diffuse and dust-scattered starlight, and the photon field of Westerlund 1 itself. The effect of the magnetic field on cooling and propagation is discussed. Klein-Nishina effects are found to be important in determining the spectral evolution of the electron population. Results. A leptonic origin of the bulk of the observed γ-rays is preferable. The model is energetically plausible, consistent with the presence of a strong shock, and allows for the observed energy-independent morphology. The hadronic model faces two main issues: confinement of particles to the emission region, and an unrealistic energy requirement.

PHANGS–JWST First Results: Stellar-feedback-driven Excitation and Dissociation of Molecular Gas in the Starburst Ring of NGC 1365?

We compare embedded young massive star clusters (YMCs) to (sub-)millimeter line observations tracing the excitation and dissociation of molecular gas in the starburst ring of NGC 1365. This galaxy hosts one of the strongest nuclear starbursts and richest populations of YMCs within 20 Mpc. Here we combine near-/mid-IR PHANGS–JWST imaging with new Atacama Large Millimeter/submillimeter Array multi-J CO (1–0, 2–1 and 4–3) and [C i] (1–0) mapping, which we use to trace CO excitation via R 42 = I CO(4−3)/I CO(2−1) and R 21 = I CO(2−1)/I CO(1−0) and dissociation via R CICO = I [CI](1−0)/I CO(2−1) at 330 pc resolution. We find that the gas flowing into the starburst ring from northeast to southwest appears strongly affected by stellar feedback, showing decreased excitation (lower R 42) and increased signatures of dissociation (higher R CICO) in the downstream regions. There, radiative-transfer modeling suggests that the molecular gas density decreases and temperature and [CI/CO] abundance ratio increase. We compare R 42 and R CICO with local conditions across the regions and find that both correlate with near-IR 2 μm emission tracing the YMCs and with both polycyclic aromatic hydrocarbon (11.3 μm) and dust continuum (21 μm) emission. In general, R CICO exhibits ∼0.1 dex tighter correlations than R 42, suggesting C i to be a more sensitive tracer of changing physical conditions in the NGC 1365 starburst than CO (4–3). Our results are consistent with a scenario where gas flows into the two arm regions along the bar, becomes condensed/shocked, forms YMCs, and then these YMCs heat and dissociate the gas.

GOALS-JWST: Revealing the Buried Star Clusters in the Luminous Infrared Galaxy VV 114

We present the results of a James Webb Space Telescope NIRCam investigation into the young massive star cluster (YMC) population in the luminous infrared galaxy VV 114. We identify 374 compact YMC candidates with signal-to-noise ratios ≥ 3, 5, and 5 at F150W, F200W, and F356W, respectively. A direct comparison with our HST cluster catalog reveals that ∼20% of these sources are undetected at optical wavelengths. Based on yggdrasil stellar population models, we identify 17 YMC candidates in our JWST imaging alone with F150W – F200W and F200W – F356W colors suggesting they are all very young, dusty (A V = 5–15), and massive (105.8 < M ⊙ < 106.1). The discovery of these “hidden” sources, many of which are found in the “overlap” region between the two nuclei, quadruples the number of t < 3 Myr clusters and nearly doubles the number of t < 6 Myr clusters detected in VV 114. Now extending the cluster age distribution () to the youngest ages, we find a slope of γ = −1.30 ± 0.39 for 106 < τ(yr) < 107, which is consistent with the previously determined value from 107 < τ(yr) < 108.5, and confirms that VV 114 has a steep age distribution slope for all massive star clusters across the entire range of cluster ages observed. Finally, the consistency between our JWST- and HST-derived age distribution slopes indicates that the balance between cluster formation and destruction has not been significantly altered in VV 114 over the last 0.5 Gyr.