Tracers Of Star Formation Research Articles

Gravitational-wave Dark Siren Cosmology Systematics from Galaxy Weighting

The detection of GW170817 and the measurement of its redshift from the associated electromagnetic counterpart provided the first gravitational-wave (GW) determination of the Hubble constant (H 0), demonstrating the potential power of standard siren cosmology. In contrast to this “bright siren” approach, the “dark siren” approach can be utilized for GW sources in the absence of an electromagnetic counterpart: One considers all galaxies contained within the localization volume as potential hosts. When statistically averaging over the potential host galaxies, weighting them by physically motivated properties (e.g., tracing star formation or stellar mass) could improve convergence. Using mock galaxy catalogs, we explore the impact of these weightings on the measurement of H 0. We find that incorrect weighting schemes can lead to significant biases due to two effects: the assumption of an incorrect galaxy redshift distribution, and preferentially weighting incorrect host galaxies during the inference. The magnitudes of these biases are influenced by the number of galaxies along each line of sight, the measurement uncertainty in the GW luminosity distance, and correlations in the parameter space of galaxies. We show that the bias may be overcome from improved localization constraints in future GW detectors, a strategic choice of priors or weighting prescription, and by restricting the analysis to a subset of high-signal-to-noise ratio events. We propose the use of hierarchical inference as a diagnostic of incorrectly weighted prescriptions. Such approaches can simultaneously infer the correct weighting scheme and the values of the cosmological parameters, thereby mitigating the bias in dark siren cosmology due to incorrect host-galaxy weighting.

Full disc [CII] mapping of nearby star-forming galaxies. SOFIA FIFI/LS observations of NGC 3627, NGC 4321, and NGC 6946

As a major cooling line of interstellar gas, the far-infrared 158 line from singly ionised carbon is an important tracer of various components of the interstellar medium in galaxies across all spatial and morphological scales. Yet, there is still not a strong constraint on the origins of emission. In this work, we derive the resolved star formation rate relation and aim to unravel the complexity of the origin of CII. We used the Field-Imaging Far-Infrared Line Spectrometer on board the Stratospheric Observatory for Infrared Astronomy to map in three nearby star-forming galaxies at sub-kiloparsec scales, namely, NGC,3627, NGC,4321, and NGC,6946, and we compared these observations to the galactic properties derived from complementary data from the literature. We find that the relationship between the fine structure line and star formation rate shows variations between the galaxies as well as between different environments within each galaxy. Our results show that the use of as a tracer for star formation is much more tangled than has previously been suggested within the extragalactic literature, which typically focuses on small regions of galaxies and/or uses large-aperture sampling of many different physical environments. As found within resolved observations of the Milky Way, the picture obtained from observations is complicated by its local interstellar medium conditions. Future studies will require a larger sample and additional observational tracers, obtained on spatial scales within galaxies, in order to accurately disentangle the origin of and calibrate its use as a star formation tracer.

The ALPINE-ALMA [CII] Survey: Modelling ALMA and JWST lines to constrain the interstellar medium of z∼ 5 galaxies

Aims. We have devised a model for estimating the ultraviolet (UV) and optical line emission (i.e. CIII] 1909 Å, Hβ, [OIII] 5007 Å, Hα, and [NII] 6583 Å) that traces HII regions in the interstellar medium (ISM) of a subset of galaxies at z ~ 4-6 from the ALMA large programme ALPINE. The aim is to investigate the combined impact of binary stars in the stellar population and an abrupt quenching in the star formation history (SFH) on the line emission. This is crucial for understanding the ISM’s physical properties in the Universe’s earliest galaxies and identifying new star formation tracers in high-z galaxies. Methods. The model simulates HII plus PhotoDissociation Region (PDR) complexes by performing radiative transfer through 1D slabs characterised by gas density (n), ionisation parameter (U), and metallicity (Z). The model also takes into account (a) the heating from star formation, whose spectrum has been simulated with Starburst99 and Binary Population and Spectral Synthesis (BPASS) to quantify the impact of binary stars; and (b) a constant, exponentially declining, and quenched SFH. For each galaxy, we selected from our CLOUDY models the theoretical ratios between the [CII] line emission that trace PDRs and nebular lines from HII regions. These ratios were then used to derive the expected optical/UV lines from the observed [CII]. Results. We find that binary stars have a strong impact on the line emission after quenching, by keeping the UV photon flux higher for a longer time. This is relevant in maintaining the free electron temperature and ionised column density in HII regions unaltered up to 5 Myr after quenching. Furthermore, we constrained the ISM properties of our subsample, finding a low ionisation parameter of log U≈ − 3.8 ± 0.2 and high densities of log(n/cm−3)≈2.9 ± 0.6. Finally, we derive UV/optical line luminosity-star formation rate relations (log(Lline/erg s−1) = α log(SFR/M⊙ yr−1) + β) for different burstiness parameter (ks) values. We find that in the fiducial BPASS model, the relations have a negligible SFH dependence but depend strongly on the ks value, while in the SB99 case, the dominant dependence is on the SFH. We propose their potential use for characterising the burstiness of galaxies at high z.

Preferential occurrence of fast radio bursts in massive star-forming galaxies.

Fast radio bursts (FRBs) are millisecond-duration events detected from beyond the Milky Way. FRB emission characteristics favour highly magnetized neutron stars, or magnetars, as the sources1, as evidenced by FRB-like bursts from a galactic magnetar2,3, and the star-forming nature of FRB host galaxies4,5. However, the processes that produce FRB sources remain unknown6. Although galactic magnetars are often linked to core-collapse supernovae (CCSNe)7, it is uncertain what determines which supernovae result in magnetars. The galactic environments of FRB sources can be used to investigate their progenitors. Here, we present the stellar population properties of 30 FRB host galaxies discovered by the Deep Synoptic Array (DSA-110). Our analysis shows a marked deficit of low-mass FRB hosts compared with the occurrence of star formation in the Universe, implying that FRBs are a biased tracer of star formation, preferentially selecting massive star-forming galaxies. This bias may be driven by galaxy metallicity, which is positively correlated with stellar mass8. Metal-rich environments may favour the formation of magnetar progenitors through stellar mergers9,10, as higher-metallicity stars are less compact and more likely to fill their Roche lobes, leading to unstable mass transfer. Although massive stars do not have convective interiors to generate strong magnetic fields by dynamo11, merger remnants are thought to have the requisite internal magnetic-field strengths to result in magnetars11,12. The preferential occurrence of FRBs in massive star-forming galaxies suggests that a core-collapse supernova of merger remnants preferentially forms magnetars.

On the limitations of Hα luminosity as a star formation tracer in spatially resolved observations

On the limitations of Hα luminosity as a star formation tracer in spatially resolved observations

A spatially resolved radio spectral study of the galaxy M 51

Context. Radio continuum emission from galaxies at gigahertz frequencies can be used as an extinction-free tracer of star formation. However, at frequencies of a few hundred megahertz, there is evidence for low-frequency spectral flattening. Aims. We wish to understand the origin of this low-frequency flattening better, and to this end, we performed a spatially resolved study of the nearby spiral galaxy M51. We explored the different effects that can cause a flattening of the spectrum towards lower frequencies, such as free–free absorption and cosmic-ray ionisation losses. Methods. We used radio continuum intensity maps between 54 and 8350 MHz at eight different frequencies, with observations at 240 MHz from the Giant Metrewave Radio Telescope presented for the first time. We corrected for the contribution from thermal free–free emission using an H α map that was corrected for extinction with 24 μm data. We fitted free–free absorption models to the radio spectra to determine the emission measure (EM) as well as polynomial functions to measure the non-thermal spectral curvature. We also obtained a new extinction-corrected H α intensity map from the Metal-THINGS survey using integral field unit spectroscopy. Results. The non-thermal low-frequency radio continuum spectrum between 54 and 144 MHz is very flat and even partially inverted, particularly in the spiral arms; in contrast, the spectrum at higher frequencies is typical for a non-thermal radio continuum spectrum. However, we did not find any correlation between the EMs calculated from radio and from H α observations; instead, the non-thermal spectral curvature weakly correlates with the H I gas-mass surface density. This suggests that cosmic-ray ionisation losses play an important role in the low-frequency spectral flattening. Conclusions. The observed spectral flattening towards low frequencies in M51 is caused by a combination of ionisation losses and free–free absorption. The reasons for this flattening need to be understood in order to use sub-gigahertz frequencies as a tracer of star formation.

Type Ia Supernova Progenitor Properties and their Host Galaxies

We present an eigenfunction method to analyze 161 visual light curves (LCs) of Type Ia supernovae (SNe Ia) obtained by the Carnegie Supernova Project to characterize their diversity and host-galaxy correlations. The eigenfunctions are based on the delayed-detonation (DD) scenario using three parameters: the LC stretch s determined by the amount of deflagration burning governing the 56Ni production, the main-sequence mass M MS of the progenitor white dwarf controlling the explosion energy, and its central density ρ c shifting the 56Ni distribution. Our analysis tool (Supernova Parameter Analysis Tool) extracts the parameters from observations and projects them into physical space using their allowed ranges (M MS ≤ 8 M ⊙, ρ c ≤ 7–8 × 109 g cm−3). The residuals between fits and individual LC points are ≈1%–3% for ≈92% of objects. We find two distinct M MS groups corresponding to a fast (≈4–65 Myr) and a slow(≈200–500 Myr) stellar evolution. Most underluminous SNe Ia have hosts with low star formation but high M MS, suggesting slow evolution times of the progenitor system. 91T-like SNe show very similar LCs and high M MS and are correlated to star formation regions, making them potentially important tracers of star formation in the early Universe out to z ≈ 4–11. Some ∼6% outliers with nonphysical parameters using DD scenarios can be attributed to superluminous SNe Ia and subluminous SNe Ia with hosts of active star formation. For deciphering the SNe Ia diversity and high-precision SNe Ia cosmology, the importance is shown for LCs covering out to ≈60 days past maximum. Finally, our method and results are discussed within the framework of multiple explosion scenarios, and in light of upcoming surveys.

Gas conditions of a star-formation selected sample in the first billion years

ABSTRACT We present Atacama Large Millimetre/submillimetre Array (ALMA) observations of the [O iii] 88 $\mu {\rm m}$ emission of a sample of thirteen galaxies at $z = 6$ to 7.6 selected as [C ii]-emitting companion sources of quasars. To disentangle the origins of the luminous Oxygen line in the $z\, \gt\, 6$ Universe, we looked at emission-line galaxies that are selected through an excellent star-formation tracer [C ii] with star-formation rates between 9 and 162 $\rm {\rm M}_{\odot }\,yr^{-1}$. Direct observations reveal [O iii] emission in just a single galaxy ($L_\mathrm{[O\, {\small III}]}/L_\mathrm{[C\, {\small II}]}$$\, = 2.3$), and a stacked image shows no [O iii] detection, providing deep upper limits on the $L_\mathrm{[O\, {\small III}]}/L_\mathrm{[C\, {\small II}]}$ ratios in the $z \gt 6$ Universe ($L_\mathrm{[O\, {\small III}]}/L_\mathrm{[C\, {\small II}]}$$\, \lt 1.2$ at $3 \sigma$). While the fidelity of this sample is high, no obvious optical/near-infrared counterpart is seen in the JWST imaging available for four galaxies. Additionally accounting for low-z CO emitters, line stacking shows that our sample-wide result remains robust: The enhanced $L_\mathrm{[O\, {\small III}]}/L_\mathrm{[C\, {\small II}]}$ reported in the first billion years of the Universe is likely due to the selection towards bright, blue Lyman-break galaxies with high surface star-formation rates or young stellar populations. The deep upper limit on the rest-frame 90 μm continuum emission ($\lt 141 \mu$Jy at $3 \sigma$), implies a low average dust temperature ($T_\mathrm{dust} \lesssim 30\,$ K) and high-dust mass ($M_\mathrm{dust} \sim 10^8\, \mathrm{M}_\odot$). As more normal galaxies are explored in the early Universe, synergy between JWST and ALMA is fundamental to further investigate the ISM properties of the a broad range of samples of high-z galaxies.

JWST/NIRSpec and MIRI observations of an expanding, jet-driven bubble of warm H2 in the radio galaxy 3C 326 N

The physical link between AGN activity and the suppression of star formation in their host galaxies is one of the major open questions of the AGN feedback scenario. The Spitzer space mission revealed a subset of powerful nearby radio galaxies with unusually bright line emission from warm ($T 100$ K) molecular hydrogen, while typical star-formation tracers such as polycyclic aromatic hydrocarbons (PAHs) or a dust continuum have been exceptionally faint or undetected. Here, we present JWST NIRSpec and MIRI MRS IFU observations of one of the best studied galaxies of this class, 3C 326 N at z=0.09. We identified a total of 19 lines of the S, O, and Q series of ro-vibrational H$_2$ emission with NIRSpec at a 0.11 spatial resolution, probing a small quantity odot $) of gas at temperatures of $T 1000$ K. We also mapped the rotational mid-infrared lines of H$_2$ 0--0 S(3), S(5), and S(6) at a spatial resolution of 0.4 with MIRI/MRS, probing most of the $2 odot $ of warm H$_2$ in this galaxy. The CO band heads show a stellar component consistent with a 'slow-rotator' that is typical of a massive ($3 $\,M$_ galaxy, offering a reliable systemic redshift of $z=0.08979 0.0003$. The extended line emission shows a bipolar bubble expanding through the molecular disk at velocities of up to $, delineated by several bright clumps along the northern outer rim, potentially coming from gas fragmentation. Throughout the disk, H$_2$ is very broadly dispersed, with an FWHM of $ 100-1300$ km $ and complex, dual-component Gaussian line profiles. The extended FeII lambda 1.644 and Paalpha follow the same morphology, however NeIII lambda 15.56 is more symmetric about the nucleus. We show that most of the gas (with the exception of NeIII lambda 15.56) is predominantly heated by shocks driven by the radio jets into the gas, both for the ro-vibrational and rotational H$_2$ lines. In addition, the accompanying line broadening is sufficient to suppress star formation in the molecular gas. We also compared the morphology and kinematics of the rotational and ro-vibrational lines, finding the latter to be a good proxy to the global morphology and kinematic properties of the former in strongly turbulent environments. This demonstrates the potential of using the higher frequency ro-vibrational lines in studying turbulent molecular gas. Provided they are bright enough, they would allow us to examine turbulence in galaxies during the early phases of cosmic history, while most rotational lines are red-shifted out of the MIRI bandpass for $z ge1.5$

Molecular isotopologue measurements toward super star clusters and the relation to their ages in NGC 253 with ALCHEMI

Context. Determining the evolution of the CNO isotopes in the interstellar medium (ISM) of starburst galaxies can yield important constraints on the ages of super star clusters (SSCs), or on other aspects and factors contributing to their evolution, such as the initial mass function (IMF). Due to the time-dependent nature of the abundances of isotopes within the ISM – as they are supplied from processes such as nucleosynthesis or chemical fractionation –, this provides the opportunity to test whether or not isotope ratios trace the ages of highly star-forming regions, such as SSCs. Aims. The goal of this study is to investigate whether the isotopic variations in SSC regions within NGC 253 are correlated with their different ages as derived from stellar population modelling. Methods. We measured abundance ratios of CO, HCN, and HCO+ isotopologues in six regions containing SSCs within NGC 253 using high-spatial-resolution (1.6″, ∼28 pc) data from the ALCHEMI (ALma Comprehensive High-resolution Extragalactic Molecular Inventory) ALMA Large program. We then analysed these ratios using RADEX radiative transfer modelling, with the parameter space sampled using the nested sampling Monte Carlo algorithm MLFriends. These abundance ratios were then compared to ages predicted in each region via the fitting of observed star-formation tracers (such as Brγ) to Starburst99 starburst stellar population evolution models. Results. We determined the isotopic column density ratios across multiple regions of SSC activity in NGC 253 using non-LTE radiative transfer modelling. We do not find any significant trend with age for the CO and HCN isotopologue ratios on timescales of the ages of the SSC* regions observed. However, HCO+ may show a correlation with age over these timescales in 12C/13C. Conclusions. The driving factors of these ratios within SSCs could be the IMF or fractionation effects. To further probe these effects in SSCs over time, a larger sample of SSCs must be observed spanning a larger age range.

Constraints on the z ∼ 5 Star-forming Galaxy Luminosity Function From Hubble Space Telescope Imaging of an Unbiased and Complete Sample of Long Gamma-Ray Burst Host Galaxies

We present rest-frame UV Hubble Space Telescope imaging of the largest and most complete sample of 23 long-duration gamma-ray burst (GRB) host galaxies between redshifts 4 and 6. Of these 23, we present new WFC3/F110W imaging for 19 of the hosts, which we combine with archival WFC3/F110W and WFC3/F140W imaging for the remaining four. We use the photometry of the host galaxies from this sample to characterize both the rest-frame UV luminosity function (LF) and the size–luminosity relation of the sample. We find that when assuming the standard Schechter-function parameterization for the UV LF, the GRB host sample is best fit with α=−1.30−0.25+0.30 and M*=−20.33−0.54+0.44 mag, which are consistent with results based on z ∼ 5 Lyman-break galaxies. We find that ∼68% of our size–luminosity measurements fall within or below the same relation for Lyman-break galaxies at z ∼ 4. This study observationally confirms expectations that at z ∼ 5 Lyman-break and GRB host galaxies should trace the same population and demonstrates the utility of GRBs as probes of hidden star formation in the high-redshift Universe. Under the assumption that GRBs unbiasedly trace star formation at this redshift, our nondetection fraction of 7/23 is consistent at the 95% confidence level with 13%–53% of star formation at redshift z ∼ 5 occurring in galaxies fainter than our detection limit of M 1600Å ≈ −18.3 mag.

The Size–Linewidth Relation and Signatures of Feedback from Quiescent to Active Star Forming Regions in the LMC

To investigate the effects of stellar feedback on the gravitational state of giant molecular clouds (GMCs), we study 12CO and 13CO Atacama Large Millimeter/submillimeter Array maps of nine GMCs distributed throughout the Large Magellanic Cloud (LMC), the nearest star-forming galaxy to our own. We perform noise and resolution matching on the sample, working at a common resolution of 3.″5 (0.85 pc at the LMC distance of 50 kpc), and we use the Spectral Clustering for Molecular Emission Segmentation clustering algorithm to identify discrete substructure, or “clumps.” We supplement these data with three tracers of recent star formation: 8 μm surface brightness, continuum-subtracted Hα flux, and interstellar radiation field energy density inferred from dust emission. The 12CO clumps identified cover a range of 3.6 dex in luminosity-based mass and 2.4 dex in average 8 μm surface brightness, representative of the wide range of conditions of the interstellar medium in the LMC. Our observations suggest evidence for increased turbulence in these clouds. While the turbulent linewidths are correlated with clump surface density, in agreement with previous observations, we find even better correlation with the three star formation activity tracers considered, suggesting stellar energy injection plays a significant role in the dynamical state of the clumps. The excess linewidths we measure do not appear to result from opacity broadening. 12CO clumps are found to be typically less gravitationally bound than 13CO clumps, with some evidence of the kinetic-to-gravitational potential energy ratio increasing with star formation tracers. Further multiline analysis may better constrain the assumptions made in these calculations.

Radio-only and Radio-to-far-ultraviolet Spectral Energy Distribution Modeling of 14 ULIRGs: Insights into the Global Properties of Infrared Bright Galaxies

We present detailed spectral energy distribution (SED) modeling of 14 local ultraluminous infrared galaxies (ULIRGs) with outstanding photometric data from the literature covering the ultraviolet–infrared (FIR) and radio bands (∼50 MHz to ∼30 GHz). We employ the CIGALE SED fitting code to model the ultraviolet–FIR–radio SED. For the radio-only SED modeling, we use the UltraNest package, leveraging its nested sampling algorithm. Combining the results from our previous study on 11 luminous infrared galaxies (LIRGs), we discuss the global astrophysical properties of a sample of 25 starburst galaxies (z < 0.5). Their radio spectra are frequently characterized by bends and turnovers, with no indication of ULIRGs exhibiting more complicated SEDs than LIRGs despite showing more signs of interactions. Including radio measurements in the CIGALE modeling constrained the dust luminosity and star formation rate (SFR) estimates by more than 1 order of magnitude better than previously reported for starburst galaxies. We show that total and nonthermal radio luminosity at 1.4 and 4.8 GHz frequencies can be good estimators of recent SFRs for all LIRGs and those ULIRGS with an insignificant influence of active galactic nuclei. A weaker but still significant correlation is observed between radio SFRs at 1.4 GHz and old (averaged over 100 Myr) SFRs based on SED modeling, indicative of multiple episodes of starburst activity during their lifetime. The thermal radio luminosity at 4.8 GHz is a better tracer of recent star formation than the thermal luminosity at 1.4 GHz. Statistically, our modeled nonthermal radio spectral indices do not significantly correlate with redshift, stellar mass, SFR, specific SFR, and dust mass.

Supernova environments in J-PLUS

We investigated the local environmental properties of 418 supernovae (SNe) of all types using data from the Javalambre Photometric Local Universe Survey (J-PLUS), which includes five broad-band and seven narrow-band imaging filters. Our study involves two independent analyses: (1) the normalized cumulative-rank (NCR) method, which utilizes all 12 single bands along with five continuum-subtracted narrow-band emission and absorption bands, and (2) simple stellar population (SSP) synthesis, where we build spectral energy distributions (SED) of the surrounding 1 kpc2 SN environment using the 12 broad- and narrow-band filters. Improvements on previous works include: (i) the extension of the NCR technique to other filters (broad and narrow) and the use a set of homogeneous data (same telescope and instruments); (ii) a correction for extinction to all bands based on the relation between the g − i color and the color excess E(B − V); and (iii) a correction for the contamination of the [N II] λ6583 line that falls within the Hα filter. All NCR distributions in the broad-band filters, tracing the overall light distribution in each galaxy, are similar to each other. The main difference is that type Ia, II, and IIb SNe are preferably located in redder environments than the other SN types. The radial distribution of the SNe shows that type IIb SNe seem to have a preference for occurring in the inner regions of galaxies, whereas other types of SNe occur throughout the galaxies without a distinct preference for a specific location. For the Hα filter we recover the sequence from SNe Ic, which has the highest NCR, to SNe Ia, which has the lowest; this is interpreted as a sequence in progenitor mass and age. All core-collapse SN types are strongly correlated to the [O II] emission, which also traces star formation rate (SFR), following the same sequence as in Hα. The NCR distributions of the Ca II triplet show a clear division between II-IIb-Ia and Ib-Ic-IIn subtypes, which is interpreted as a difference in the environmental metallicity. Regarding the SSP synthesis, we found that including the seven J-PLUS narrow filters in the fitting process has a more significant effect on the core-collapse SN environmental parameters than for SNe Ia, shifting their values toward more extincted, younger, and more star-forming environments, due to the presence of strong emission lines and stellar absorptions in those narrow bands.

The ALCHEMI Atlas: Principal Component Analysis Reveals Starburst Evolution in NGC 253

Molecular lines are powerful diagnostics of the physical and chemical properties of the interstellar medium (ISM). These ISM properties, which affect future star formation, are expected to differ in starburst galaxies from those of more quiescent galaxies. We investigate the ISM properties in the central molecular zone of the nearby starburst galaxy NGC 253 using the ultrawide millimeter spectral scan survey from the Atacama Large Millimeter/submillimeter Array Large Program ALCHEMI. We present an atlas of velocity-integrated images at a 1.″6 resolution of 148 unblended transitions from 44 species, including the first extragalactic detection of HCNH+ and the first interferometric images of C3H+, NO, and HCS+. We conduct a principal component analysis (PCA) on these images to extract correlated chemical species and to identify key groups of diagnostic transitions. To the best of our knowledge, our data set is currently the largest astronomical set of molecular lines to which PCA has been applied. The PCA can categorize transitions coming from different physical components in NGC 253 such as (i) young starburst tracers characterized by high-excitation transitions of HC3N and complex organic molecules versus tracers of on-going star formation (radio recombination lines) and high-excitation transitions of CCH and CN tracing photodissociation regions, (ii) tracers of cloud-collision-induced shocks (low-excitation transitions of CH3OH, HNCO, HOCO+, and OCS) versus shocks from star formation-induced outflows (high-excitation transitions of SiO), as well as (iii) outflows showing emission from HOC+, CCH, H3O+, CO isotopologues, HCN, HCO+, CS, and CN. Our findings show these intensities vary with galactic dynamics, star formation activities, and stellar feedback.

JOYS: MIRI/MRS spectroscopy of gas-phase molecules from the high-mass star-forming region IRAS 23385+6053

Context. Space-based mid-infrared (IR) spectroscopy is a powerful tool for the characterization of important star formation tracers of warm gas which are unobservable from the ground. The previous mid-IR spectra of bright high-mass protostars with the Infrared Space Observatory (ISO) in the hot-core phase typically show strong absorption features from molecules such as CO2, C2H2, and HCN. However, little is known about their fainter counterparts at earlier stages. Aims. We aim to characterize the gas-phase molecular features in James Webb Space Telescope Mid-Infrared Instrument Medium Resolution Spectrometer (MIRI/MRS) spectra of the young and clustered high-mass star-forming region IRAS 23385+6053. Methods. Spectra were extracted from several locations in the MIRI/MRS field of view, targeting two mid-IR sources tracing embedded massive protostars as well as three H2 bright outflow knots at distances of &gt;8000 au from the multiple. Molecular features in the spectra were fit with local thermodynamic equilibrium (LTE) slab models, with their caveats discussed in detail. Results. Rich molecular spectra with emission from CO, H2, HD, H2O, C2H2, HCN, CO2, and OH are detected towards the two mid-IR sources. However, only CO and OH are seen towards the brightest H2 knot positions, suggesting that the majority of the observed species are associated with disks or hot core regions rather than outflows or shocks. The LTE model fits to 12CO2, C2H2, HCN emission suggest warm 120–200 K emission arising from a disk surface around one or both protostars. The abundances of CO2 and C2H2 of ~10−7 are consistent with previous observations of high-mass protostars. Weak ~500 K H2O emission at ~6–7 µm is detected towards one mid-IR source, whereas 250–1050 K H2O absorption is found in the other. The H2O absorption may occur in the disk atmosphere due to strong accretion-heating of the midplane, or in a disk wind viewed at an ideal angle for absorption. CO emission may originate in the hot inner disk or outflow shocks, but NIRSpec data covering the 4.6 µm band head are required to determine the physical conditions of the CO gas, as the high temperatures seen in the MIRI data may be due to optical depth. OH emission is detected towards both mid-IR source positions and one of the shocks, and is likely excited by water photodissociation or chemical formation pumping in a highly non-LTE manner. Conclusions. The observed molecular spectra are consistent with disks having already formed around two protostars in the young IRAS 23385+6054 system. Molecular features mostly appear in emission from a variety of species, in contrast to the more evolved hot core phase protostars which typically show only absorption; however, further observations of young high-mass protostars are needed to disentangle geometry and viewing angle effects from evolution.

CMZoom. IV. Incipient High-mass Star Formation throughout the Central Molecular Zone

In this work, we constrain the star-forming properties of all possible sites of incipient high-mass star formation in the Milky Way’s Galactic Center. We identify dense structures using the CMZoom 1.3 mm dust continuum catalog of objects with typical radii of ∼0.1 pc, and measure their association with tracers of high-mass star formation. We incorporate compact emission at 8, 21, 24, 25, and 70 μm from the Midcourse Space Experiment, Spitzer, Herschel, and SOFIA, cataloged young stellar objects, and water and methanol masers to characterize each source. We find an incipient star formation rate (SFR) for the Central Molecular Zone (CMZ) of ∼0.08 M ⊙ yr−1 over the next few 105 yr. We calculate upper and lower limits on the CMZ’s incipient SFR of ∼0.45 and ∼0.05 M ⊙ yr−1,respectively, spanning roughly equal to and several times greater than other estimates of CMZ’s recent SFR. Despite substantial uncertainties, our results suggest the incipient SFR in the CMZ may be higher than previously estimated. We find that the prevalence of star formation tracers does not correlate with source volume density, but instead ≳75% of high-mass star formation is found in regions above a column density ratio (N SMA/N Herschel) of ∼1.5. Finally, we highlight the detection of atoll sources, a reoccurring morphology of cold dust encircling evolved infrared sources, possibly representing H ii regions in the process of destroying their envelopes.

Polycyclic aromatic hydrocarbon luminous galaxies in JWST CEERS data

ABSTRACT It has been an unanswered question how many dusty galaxies have been undetected from the state-of-the-art observational surveys. JWST enables us to detect faint infrared (IR) galaxies that have prominent polycyclic aromatic hydrocarbon (PAH) features in the mid-IR wavelengths. PAH is a valuable tracer of star formation and dust properties in the mid-IR wavelength. The JWST Cosmic Evolution Early Release Science (CEERS) fields provide us with wavelength coverage from 7.7 to 21 μm using six photometric bands of the mid-IR instrument (MIRI). We have identified galaxies dominated by mid-IR emission from PAHs, termed PAH galaxies. From our multiband photometry catalogue, we selected 10 PAH galaxies displaying high flux ratios of log (S15/S10) &amp;gt; 0.8. The SED fitting analysis indicates that these galaxies are star-forming galaxies with total IR luminosities of 1010 ∼ 1011.5 L⊙ at z ∼1. The morphology of PAH galaxies does not show any clear signatures of major merging or interaction within the MIRI resolution. The majority of them are on the star-formation main sequence at z ∼ 1. Our result demonstrates that JWST can detect PAH emissions from normal star-forming galaxies at z ∼ 1, in addition to ultra-luminous infrared galaxies or luminous IR galaxies (LIRGs).

JWST reveals star formation across a spiral arm in M33

ABSTRACT Young stellar objects (YSOs) are the gold standard for tracing star formation in galaxies but have been unobservable beyond the Milky Way and Magellanic Clouds. But that all changed when the JWST was launched, which we use to identify YSOs in the Local Group galaxy M33, marking the first time that individual YSOs have been identified at these large distances. We present Mid-Infrared Instrument (MIRI) imaging mosaics at 5.6 and 21 $\mu$m that cover a significant portion of one of M33’s spiral arms that has existing panchromatic imaging from the Hubble Space Telescope and deep Atacama Large Millimeter/submillimeter Array CO measurements. Using these MIRI and Hubble Space Telescope images, we identify point sources using the new dolphot MIRI module. We identify 793 candidate YSOs from cuts based on colour, proximity to giant molecular clouds (GMCs), and visual inspection. Similar to Milky Way GMCs, we find that higher mass GMCs contain more YSOs and YSO emission, which further show YSOs identify star formation better than most tracers that cannot capture this relationship at cloud scales. We find evidence of enhanced star formation efficiency in the southern spiral arm by comparing the YSOs to the molecular gas mass.

Massive star formation in the hub–filament system of RCW 117

ABSTRACT We present a multiwavelength investigation of the hub–filament system RCW 117 (IRAS 17059-4132), which shows intricate filamentary features in the far-infrared, mapped using Herschel images. We obtain the column density and dust temperature maps for the region using the Herschel images and identify 88 cores and 12 filaments from the column density map of the region (18 arcmin × 18 arcmin). We employ the ThrUMMS 13CO (J = 1–0) data for probing the kinematics in RCW 117, and find velocity gradients (∼0.3–1 km s−1 pc−1) with hints of matter inflow along the filamentary structures. Ionized gas emission from the associated H ii region is examined using the Giant Metrewave Radio Telescope (GMRT) at 610 and 1280 MHz and is found to be of extent 5 × 3 pc2 with intensity being brightest towards the hub. We estimate the peak electron density towards the hub to be ∼750 cm−3. Thirty-four Class 0/I young stellar objects (YSOs) have been identified in the region using the Spitzer GLIMPSE colour–colour diagram, with many lying along the filamentary structures. Based on the (i) presence of filamentary structures, (ii) distribution of cores across the region, with ∼39 per cent found along the filamentary structures, (iii) massive star formation tracers in the hub, and (iv) the kinematics, we believe that global hierarchical collapse can plausibly explain the observed features in RCW 117.