Massive Star Clusters Research Articles

Dynamics of star associations in an SMBH-IMBH system: The case of IRS13 in the Galactic centre

The existence of intermediate-mass black holes (IMBHs) still poses challenges to theoretical and observational astronomers. Several candidates have been proposed, including the one in the IRS13 cluster in the Galactic centre, where the evidence is based on the velocity dispersion of its members; however, none have been confirmed to date. We aim to gain insights into the presence of an IMBH in the Galactic centre through a numerical study of the dynamical interplay between an IMBH and star clusters (SCs) in the vicinity of a supermassive black hole (SMBH). We used high-precision $N$-body models of IRS13-like SCs in the Galactic centre, and of more massive SCs that fall into the centre of the Galaxy from larger distances. We find that at IRS13's physical distance of $0.4\ an IRS13-sized SC cannot remain gravitationally bound even if it contains an IMBH of thousands of Msun . Thus, IRS13 appears to be an incidental present-day clumping of stars. Furthermore, we show that the velocity dispersion of tidally disrupted SCs (the likely origin of IRS13) can be fully accounted for by the tidal forces of the central SMBH; the IMBH's influence is not essential.

Magnetic Field Alignment Relative to Multiple Tracers in the High-mass Star-forming Region RCW 36

We use polarization data from SOFIA HAWC+ to investigate the interplay between magnetic fields and stellar feedback in altering gas dynamics within the high-mass star-forming region RCW 36, located in Vela C. This region is of particular interest as it has a bipolar H ii region powered by a massive star cluster, which may be impacting the surrounding magnetic field. To determine if this is the case, we apply the histogram of relative orientations (HRO) method to quantify the relative alignment between the inferred magnetic field and elongated structures observed in several data sets such as dust emission, column density, temperature, and spectral line intensity maps. The HRO results indicate a bimodal alignment trend, where structures observed with dense gas tracers show a statistically significant preference for perpendicular alignment relative to the magnetic field, while structures probed by the photodissociation region (PDR) tracers tend to align preferentially parallel relative to the magnetic field. Moreover, the dense gas and PDR associated structures are found to be kinematically distinct such that a bimodal alignment trend is also observed as a function of line-of-sight velocity. This suggests that the magnetic field may have been dynamically important and set a preferred direction of gas flow at the time that RCW 36 formed, resulting in a dense ridge developing perpendicular to the magnetic field. However, on filament scales near the PDR region, feedback may be energetically dominating the magnetic field, warping its geometry and the associated flux-frozen gas structures, causing the observed preference for parallel relative alignment.

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.

The masses of open star clusters and their tidal tails and the stellar initial mass function

Unresolved binaries have a strong influence on the observed parameters of SC We quantify this influence and compute the resulting mass underestimates and MF N-body simulations of realistic SC were used to investigate the evolution of the binary population in a SC and its tidal tails. Together with an empirically gauged stellar mass-luminosity relation, the results were then used to determine how the presence of binaries changes the photometric mass and MF of the SC and its tails as deduced from observations. T1 which is the tidal tail caused by gas expulsion, contains a larger fraction of binaries than both the SC and T2 which forms after gas expulsion. Additionally T1 has a larger velocity dispersion. Using the luminosity of an unresolved binary, an observer would underestimate its mass. This bias sensitively depends on the companion masses due to the structure of the stellar mass-luminosity relation. Combining the effect of all binaries in the simulation, the total photometric mass of the SC is underestimated by 15<!PCT!>. Dark objects (black holes and neutron stars) increase the difference between the real and observed mass of the SC further. For both the SC and the tails, the observed power-law index of the MF between a stellar mass of 0.3 and 0.7 $M_ is smaller by up to 0.2 than the real one, the real IMF being steeper by this amount. This difference is larger for stars with a larger velocity dispersion or binary fraction. Since the stars formed in SC are the progenitors of the Galactic field stars, this work suggests that the binary fractions of different populations of stars in the Galactic disc will differ as a function of the velocity dispersion. However, the direction of this correlation is currently unclear, and a complete population synthesis will be needed to investigate this effect. Variations in the binary fractions of different clusters can lead to perceived variations of the deduced stellar MF

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.

Binary black hole mergers from Population III star clusters

Binary black holes (BBHs) that are born from the evolution of Population III (Pop. III) stars are one of the main high-redshift targets for next-generation ground-based gravitational-wave (GW) detectors. Their predicted initial mass function and lack of metals make them the ideal progenitors of black holes above the upper edge of the pair-instability mass gap, that is, with a mass higher than ~134 (241) M⊙ for stars that become (or do not become) chemically homogeneous during their evolution. We investigate the effects of cluster dynamics on the mass function of BBHs that are born from Pop. III stars by considering the main uncertainties on the mass function of Pop. III stars, the orbital properties of the binary systems, the star cluster mass, and the disruption time. In our dynamical models, at least ~5% and up to 100% BBH mergers in Pop. III star clusters have a primary mass m1 above the upper edge of the pair-instability mass gap. In contrast, only ≲3% isolated BBH mergers have a primary mass above the gap, unless their progenitors evolved as chemically homogeneous stars. The lack of systems with a primary and/or secondary mass inside the gap defines a zone of avoidance with sharp boundaries in the plane of the primary mass-mass ratio. Finally, we estimate the merger rate density of BBHs. In the most optimistic case, we find a maximum of ℛ ~ 200 Gpc-3 yr-1 at ɀ ~ 15 for BBHs that formed via dynamical capture. For comparison, the merger rate density of isolated Pop. III BBHs is ℛ ≤ 10 Gpc−3 yr−1 for the same model of Pop. III star formation history.

Multiwavelength Galactic center gamma-ray observations explained by a unified cosmic-ray dynamics model

Context. High-energy (HE) and very high-energy (VHE) gamma-ray observations from the Galactic center (GC) detected extended emission correlated with the morphology of the central molecular zone (CMZ). Emission in both bands is expected to be produced by hadronic interaction between cosmic rays (CRs) and ambient gas. Aims. We examine if our three previously proposed scenarios for the CR sources and dynamics, which are consistent with the VHE gamma-ray data (1–100 TeV), also match the HE gamma-ray observations (10–300 GeV). Additionally, we analyze the effect of the isotropic Galactic CR “sea” inside the CMZ. Methods. We generated synthetic gamma-ray maps considering a simplified isotropic diffusion, but more realistic dynamics with two diffusion zones (in and out of the CMZ) and polar advection, for mono-energetic particles of 3 TeV. Additionally, we considered two gas distributions for the CMZ (with and without an inner cavity), and CR populations injected from the clusters of young massive stars (the Arches Cluster, the Quintuplet Cluster, and the nuclear star cluster), plus the supernova Sgr A East. Results. Only the combination of more realistic CR dynamics, the CMZ with an inner cavity, CR injection from all proposed sources, and a CR sea similar to that observed in the Solar System reproduced the current HE and VHE gamma-ray detection from the CMZ and was consistent with the observed gamma-rays from Sagittarius A*. Conclusions. The HE and VHE gamma-rays observations of the GC can be reproduced by a unified model for the CRs.

Massive star cluster formation

The mode of star formation that results in the formation of globular clusters and young massive clusters is difficult to constrain through observations. We present models of massive star cluster formation using the TORCH framework, which uses the Astrophysical MUltipurpose Software Environment (AMUSE) to couple distinct multi-physics codes that handle star formation, stellar evolution and dynamics, radiative transfer, and magnetohydrodynamics. We upgraded TORCH by implementing the N-body code PETAR, thereby enabling TORCH to handle massive clusters forming from 106 M⊙ clouds with ≥105 individual stars. We present results from TORCH simulations of star clusters forming from 104, 105, and 106 M⊙ turbulent spherical gas clouds (named M4, M5, M6) of radius R = 11.7 pc. We find that star formation is highly efficient and becomes more so at a higher cloud mass and surface density. For M4, M5, and M6 with initial surface densities 2.325 × 101,2,3 M⊙ pc−2, after a free-fall time of tff = 6.7,2.1,0.67 Myr, we find that ∼30%, 40%, and 60% of the cloud mass has formed into stars, respectively. The end of simulation-integrated star formation efficiencies for M4, M5, and M6 are ϵ⋆ = M⋆/Mcloud = 36%, 65%, and 85%. Observations of nearby clusters similar in mass and size to M4 have instantaneous star formation efficiencies of ϵinst ≤ 30%, which is slightly lower than the integrated star formation efficiency of M4. The M5 and M6 models represent a different regime of cluster formation that is more appropriate for the conditions in starburst galaxies and gas-rich galaxies at high redshift, and that leads to a significantly higher efficiency of star formation. We argue that young massive clusters build up through short efficient bursts of star formation in regions that are sufficiently dense (Σ ≥ 102 M⊙ pc−2) and massive (Mcloud ≥ 105 M⊙). In such environments, stellar feedback from winds and radiation is not strong enough to counteract the gravity from gas and stars until a majority of the gas has formed into stars.

Massive star cluster formation

Two main mechanisms have classically been proposed for the formation of runaway stars. In the binary supernova scenario (BSS), a massive star in a binary explodes as a supernova, ejecting its companion. In the dynamical ejection scenario, a star is ejected during a strong dynamical encounter between multiple stars. We propose a third mechanism for the formation of runaway stars: the subcluster ejection scenario (SCES), where a subset of stars from an infalling subcluster is ejected out of the cluster via a tidal interaction with the contracting gravitational potential of the assembling cluster. We demonstrate the SCES in a star-by-star simulation of the formation of a young massive cluster from a 106 M⊙ gas cloud using the TORCH framework. This star cluster forms hierarchically through a sequence of subcluster mergers determined by the initial turbulent, spherical conditions of the gas. We find that these mergers drive the formation of runaway stars in our model. Late-forming subclusters fall into the central potential, where they are tidally disrupted, forming tidal tails of runaway stars that are distributed highly anisotropically. Runaways formed in the same SCES have similar ages, velocities, and ejection directions. Surveying observations, we identify several SCES candidate groups with anisotropic ejection directions. The SCES is capable of producing runaway binaries: two wide dynamical binaries in infalling subclusters were tightened through ejection. This allows for another velocity kick via subsequent via a subsequent BSS ejection. An SCES-BSS ejection is a possible avenue for the creation of hypervelocity stars unbound to the Galaxy. The SCES occurs when subcluster formation is resolved. We expect nonspherical initial gas distributions to increase the number of calculated runaway stars, bringing it closer to observed values. The observation of groups of runaway stars formed via the SCES can thus reveal the assembly history of their natal clusters.

JWST Observations of Starbursts: Massive Star Clusters in the Central Starburst of M82

We present a near-infrared (NIR) candidate star cluster catalog for the central kiloparsec of M82 based on new JWST NIRCam images. We identify star cluster candidates using the F250M filter, finding 1357 star cluster candidates with stellar masses >104 M ⊙. Compared to previous optical catalogs, nearly all (87%) of the candidates we identify are new. The star cluster candidates have a median intrinsic cluster radius of ≈1 pc and stellar masses up to 106 M ⊙. By comparing the color–color diagram to dust-free yggdrasil stellar population models, we estimate that the star cluster candidates have A V ∼ 3−24 mag, corresponding to A 2.5μm ∼ 0.3−2.1 mag. There is still appreciable dust extinction toward these clusters into the NIR. We measure the stellar masses of the star cluster candidates, assuming ages of 0 and 8 Myr. The slope of the resulting cluster mass function is β = 1.9 ± 0.2, in excellent agreement with studies of star clusters in other galaxies.

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.

Chemical evolution of a young super star cluster at the Sunburst Arc

ABSTRACT Recent observations of high-redshift galaxies have revealed starburst galaxies with excessive amounts of nitrogen, well above that expected in standard evolutionary models. The Sunburst Arc galaxy, particularly its young and massive star cluster, represents the closest ($z=2.4$) and brightest of these as a strongly lensed object. In this work, we study the chemical history of this star cluster to determine the origin of the elevated gas-phase nitrogen using a chemical evolution model. Our model includes the enrichment of OB stars through stellar winds and core-collapse supernovae assuming that massive stars ($M\gt 25$ $\mathrm{ M}_\odot$) collapse directly into black holes at the end of their lives. We fit the model parameters to the observed chemical abundances of the Sunburst Arc cluster: O/H, C/O, and N/O. We find that the observed chemical abundances can be explained by models featuring intense star formation events, characterized by rapid gas accretion and high star formation efficiencies. Additionally, the stellar population contributing to the gas enrichment must exclude Wolf–Rayet stars. These conditions might be present in other nitrogen-rich objects as their similar chemical abundances suggest a common history. As previous studies have proposed the presence of Wolf–Rayet stars in the new nitrogen-rich objects, further research using chemodynamic modeling is necessary to ascertain the true nature of these objects.

Impacts of stellar wind and supernovae on star cluster formation: Origins of extremely high N/O ratios and multiple stellar populations

Abstract We study metal enrichment originating from stellar wind and supernovae in low-metallicity clouds by performing three-dimensional radiation hydrodynamics simulations. We find that metals ejected from stellar wind are accumulated, leading to subsequent star formation in the nitrogen-enriched gas. During this early phase, the ${\rm N/O}$ ratios are similar to observed nitrogen-enriched galaxies (${\rm [N/O]}\gtrsim 0.5$). Then, once supernovae occur, the ${\rm N/O}$ ratios decrease significantly. If the duration of star formation is comparable to the time-scale of supernovae, the mass fraction of nitrogen-enriched stars reaches half the mass of star clusters. We suggest that the mass of the star cluster needs to exceed $\sim \!10^6$ M$_{\odot }$ to have multiple populations due to stellar wind, considering the condition for massive star cluster formation and the timescales of stellar evolution.

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.

Hubble Space Telescope proper motions of Large Magellanic Cloud star clusters

We present proper motion (PM) measurements for a sample of 23 massive star clusters within the Large Magellanic Cloud using multi-epoch data from the Hubble Space Telescope (HST). We combined archival data from the ACS/WFC and WFC3/UVIS instruments with observations from a dedicated HST programme, resulting in time baselines between 4.7 and 18.2 yr available for PM determinations. For bright well-measured stars, we achieved nominal PM precisions of 55 μas yr−1 down to 11 μas yr−1. To demonstrate the potential and limitations of our PM data set, we analysed the cluster NGC 1850 and showcase a selection of different science applications. The precision of the PM measurements allows us to disentangle the kinematics of the various stellar populations that are present in the HST field. The cluster has a centre-of-mass motion that is different from the surrounding old field stars and also differs from the mean motion of a close-by group of very young stars. We determined the velocity dispersion of field stars to be 0.128 ± 0.003 mas yr−1 (corresponding to 30.3 ± 0.7 km s−1). The velocity dispersion of the cluster inferred from the PM data set most probably overestimates the true value, suggesting that the precision of the measurements at this stage is not sufficient for a reliable analysis of the internal kinematics of extra-galactic star clusters. Finally, we exploit the PM-cleaned catalogue of likely cluster members to determine any radial segregation between fast and slowly-rotating stars, finding that the former are more centrally concentrated. With this paper, we also release the astro-photometric catalogues for each cluster.

Hydrodynamic simulation of Cygnus OB2: the absence of a cluster wind termination shock

ABSTRACT We perform a large-scale hydrodynamic simulation of a massive star cluster whose stellar population mimics that of the Cygnus OB2 association. The main-sequence stars are first simulated during 1.6 Myr, until a quasi-stationary state is reached. At this time, the three Wolf–Rayet stars observed in Cygnus OB2 are added to the simulation, which continues to 2 Myr. Using a high-resolution grid in the centre of the domain, we can resolve the most massive stars individually, which allows us to probe the kinetic structures at small (parsec) scales. We find that, although the cluster excavates a spherical ‘superbubble’ cavity, the stellar population is too loosely distributed to blow a large-scale cluster wind termination shock, and that collective effects from wind–wind interactions are much less efficient than usually assumed. This challenges our understanding of the ultra-high energy emission observed from the region.

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.

Violent Starbursts and Quiescence Induced by Far-ultraviolet Radiation Feedback in Metal-poor Galaxies at High Redshift

JWST observations of galaxies at z ≳ 8 suggest that they are more luminous and clumpier than predicted by most models, prompting several proposals on the physics of star formation and feedback in the first galaxies. In this paper, we focus on the role of ultraviolet (UV) radiation in regulating star formation by performing a set of cosmological radiation hydrodynamics simulations of one galaxy at subparsec resolution with different radiative feedback models. We find that the suppression of cooling by far-UV (FUV) radiation (i.e., H2 dissociating radiation) from Population II stars is the main physical process triggering the formation of compact and massive star clusters and is responsible for the bursty star formation observed in metal-poor galaxies at z ≳ 10. Indeed, artificially suppressing FUV radiation leads to a less intense continuous mode of star formation distributed into numerous but low-mass open star clusters. Due to the intense FUV field, low-metallicity clouds remain warm (∼104 K) until they reach a relatively high density (≳103 cm−3), before becoming self-shielded and transitioning to a colder (∼100 K), partially molecular phase. As a result, star formation is delayed until the clouds accumulate enough mass to become gravitationally unstable. At this point, the clouds undergo rapid star formation, converting gas into stars with high efficiency. We therefore observe exceptionally bright galaxies (10 times brighter than for continuous star formation) and subsequent quenched “dead” galaxies that did not form stars for tens of Myr.

Gas Kinematics and Dynamics of Carina Pillars: A Case Study of G287.76-0.87

We study the kinematics of a pillar, namely G287.76-0.87, using three rotational lines of 12CO(5-4), 12CO(8-7), 12CO(11-10), and a fine structure line of [O i] 63 μm in southern Carina observed by SOFIA/GREAT. This pillar is irradiated by the associated massive star cluster Trumpler 16, which includes η Carina. Our analysis shows that the relative velocity of the pillar with respect to this ionization source is small, ∼1 km s−1, and the gas motion in the tail is more turbulent than in the head. We also performed analytical calculations to estimate the gas column density in local thermal equilibrium (LTE) conditions, which yields N CO as (∼0.2–5) × 1017 cm−2. We further constrain the gas’s physical properties in non-LTE conditions using RADEX. The non-LTE estimations result in nH2≃105cm−3 and N CO ≃ 1016 cm−2. We found that the thermal pressure within the G287.76-0.87 pillar is sufficiently high to make it stable for the surrounding hot gas and radiation feedback if the winds are not active. While they are active, stellar winds from the clustered stars sculpt the surrounding molecular cloud into pillars within the giant bubble around η Carina.