Search for: All records

We analyze image and spectral data from ≈365 ks of observations from the Chandra X-ray Observatory of the nearby, edge-on starburst galaxy NGC 253 to constrain properties of the hot phase of the outflow. We focus our analysis on the −1.1 to +0.63 kpc region of the outflow and define several regions for spectral extraction where we determine best-fit temperatures and metal abundances. We find that the temperatures and electron densities peak in the central ∼250 pc region of the outflow and decrease with distance. These temperature and density profiles are in disagreement with an adiabatic spherically expanding starburst wind model and suggest the presence of additional physics such as mass loading and nonspherical outflow geometry. Our derived temperatures and densities yield cooling times in the nuclear region of a few million years, which may imply that the hot gas can undergo bulk radiative cooling as it escapes along the minor axis. Our metal abundances of O, Ne, Mg, Si, S, and Fe all peak in the central region and decrease with distance along the outflow, with the exception of Ne, which maintains a flat distribution. The metal abundances indicate significant dilution outside of the starburst region. We also find estimates of the mass outflow rates, which are 2.8M_⊙yr⁻¹in the northern outflow and 3.2M_⊙yr⁻¹in the southern outflow. Additionally, we detect emission from charge exchange and find it makes a significant contribution (20%–42%) to the total broadband (0.5–7 keV) X-ray emission in the central and southern regions of the outflow.

 

The analytic galactic wind model derived by Chevalier and Clegg in 1985 (CC85) assumes uniform energy and mass-injection within the starburst galaxy nucleus. However, the structure of nuclear star clusters, bulges, and star-forming knots are non-uniform. We generalize to cases with spherically-symmetric energy/mass injection that scale as r−Δ within the starburst volume R, providing solutions for Δ = 0, 1/2, 1, 3/2, and 2. In marked contrast with the CC85 model (Δ = 0), which predicts zero velocity at the centre, for a singular isothermal sphere profile (Δ = 2), we find that the flow maintains a constant Mach number of $\mathcal {M}=\sqrt{3/5} \simeq 0.77$ throughout the volume. The fast interior flow can be written as $v_{r \lt R} = (\dot{E}_T/3\dot{M}_T)^{1/2} \simeq 0.41 \, v_\infty$, where v∞ is the asymptotic velocity, and $\dot{E}_T$ and $\dot{M}_T$ are the total energy and mass injection rates. For $v_\infty \simeq 2000 \, \mathrm{km \, s^{-1}}$, $v_{r\lt R} \simeq 820 \, \mathrm{km\, s^{-1}}$ throughout the wind-driving region. The temperature and density profiles of the non-uniform models may be important for interpreting spatially-resolved maps of starburst nuclei. We compute velocity resolved spectra to contrast the Δ = 0 (CC85) and Δ = 2 models. Next generation X-ray space telescopes such as XRISM may assess these kinematic predictions.

 

Abstract We present the ALMA detection of molecular outflowing gas in the central regions of NGC 4945, one of the nearest starbursts and also one of the nearest hosts of an active galactic nucleus (AGN). We detect four outflow plumes in CO J = 3 − 2 at ∼0.″3 resolution that appear to correspond to molecular gas located near the edges of the known ionized outflow cone and its (unobserved) counterpart behind the disk. The fastest and brightest of these plumes has emission reaching observed line-of-sight projected velocities of over 450 km s −1 beyond systemic, equivalent to an estimated physical outflow velocity v ≳ 600 km s −1 for the fastest emission. Most of these plumes have corresponding emission in HCN or HCO + J = 4 − 3. We discuss a kinematic model for the outflow emission where the molecular gas has the geometry of the ionized gas cone and shares the rotation velocity of the galaxy when ejected. We use this model to explain the velocities we observe, constrain the physical speed of the ejected material, and account for the fraction of outflowing gas that is not detected due to confusion with the galaxy disk. We estimate a total molecular mass outflow rate M ̇ mol ∼ 20 M ⊙ yr −1 flowing through a surface within 100 pc of the disk midplane, likely driven by a combination of the central starburst and AGN. 

We measure the molecular gas environment near recent (<100 yr old) supernovae (SNe) using ∼1″ or ≤150 pc resolution CO (2–1) maps from the PHANGS–Atacama Large Millimeter/submillimeter Array (ALMA) survey of nearby star-forming galaxies. This is arguably the first such study to approach the scales of individual massive molecular clouds (M_mol≳ 10^5.3M_⊙). Using the Open Supernova Catalog, we identify 63 SNe within the PHANGS–ALMA footprint. We detect CO (2–1) emission near ∼60% of the sample at 150 pc resolution, compared to ∼35% of map pixels with CO (2–1) emission, and up to ∼95% of the SNe at 1 kpc resolution, compared to ∼80% of map pixels with CO (2–1) emission. We expect the ∼60% of SNe within the same 150 pc beam, as a giant molecular cloud will likely interact with these clouds in the future, consistent with the observation of widespread SN–molecular gas interaction in the Milky Way, while the other ∼40% of SNe without strong CO (2–1) detections will deposit their energy in the diffuse interstellar medium, perhaps helping drive large-scale turbulence or galactic outflows. Broken down by type, we detect CO (2–1) emission at the sites of ∼85% of our 9 stripped-envelope SNe (SESNe), ∼40% of our 34 Type II SNe, and ∼35% of our 13 Type Ia SNe, indicating that SESNe are most closely associated with the brightest CO (2–1) emitting regions in our sample. Our results confirm that SN explosions are not restricted to only the densest gas, and instead exert feedback across a wide range of molecular gas densities.

 

Abstract We compare mid-infrared (mid-IR), extinction-corrected H α , and CO (2–1) emission at 70–160 pc resolution in the first four PHANGS–JWST targets. We report correlation strengths, intensity ratios, and power-law fits relating emission in JWST’s F770W, F1000W, F1130W, and F2100W bands to CO and H α . At these scales, CO and H α each correlate strongly with mid-IR emission, and these correlations are each stronger than the one relating CO to H α emission. This reflects that mid-IR emission simultaneously acts as a dust column density tracer, leading to a good match with the molecular-gas-tracing CO, and as a heating tracer, leading to a good match with the H α . By combining mid-IR, CO, and H α at scales where the overall correlation between cold gas and star formation begins to break down, we are able to separate these two effects. We model the mid-IR above I ν = 0.5 MJy sr −1 at F770W, a cut designed to select regions where the molecular gas dominates the interstellar medium (ISM) mass. This bright emission can be described to first order by a model that combines a CO-tracing component and an H α -tracing component. The best-fitting models imply that ∼50% of the mid-IR flux arises from molecular gas heated by the diffuse interstellar radiation field, with the remaining ∼50% associated with bright, dusty star-forming regions. We discuss differences between the F770W, F1000W, and F1130W bands and the continuum-dominated F2100W band and suggest next steps for using the mid-IR as an ISM tracer. 

Connecting the gas in H ii regions to the underlying source of the ionizing radiation can help us constrain the physical processes of stellar feedback and how H ii regions evolve over time. With PHANGS–MUSE, we detect nearly 24 000 H ii regions across 19 galaxies and measure the physical properties of the ionized gas (e.g. metallicity, ionization parameter, and density). We use catalogues of multiscale stellar associations from PHANGS–HST to obtain constraints on the age of the ionizing sources. We construct a matched catalogue of 4177 H ii regions that are clearly linked to a single ionizing association. A weak anticorrelation is observed between the association ages and the $\mathrm{H}\, \alpha$ equivalent width $\mathrm{EW}(\mathrm{H}\, \alpha)$, the $\mathrm{H}\, \alpha/\mathrm{FUV}$ flux ratio, and the ionization parameter, log q. As all three are expected to decrease as the stellar population ages, this could indicate that we observe an evolutionary sequence. This interpretation is further supported by correlations between all three properties. Interpreting these as evolutionary tracers, we find younger nebulae to be more attenuated by dust and closer to giant molecular clouds, in line with recent models of feedback-regulated star formation. We also observe strong correlations with the local metallicity variations and all three proposed age tracers, suggestive of star formation preferentially occurring in locations of locally enhanced metallicity. Overall, $\mathrm{EW}(\mathrm{H}\, \alpha)$ and log q show the most consistent trends and appear to be most reliable tracers for the age of an H ii region.

 

We present new 0.3–21μm photometry of SN 2021aefx in the spiral galaxy NGC 1566 at +357 days afterB-band maximum, including the first detection of any Type Ia supernova (SN Ia) at >15μm. These observations follow earlier JWST observations of SN 2021aefx at +255 days after the time of maximum brightness, allowing us to probe the temporal evolution of the emission properties. We measure the fraction of flux emerging at different wavelengths and its temporal evolution. Additionally, the integrated 0.3–14μm decay rate of Δm_0.3–14= 1.35 ± 0.05 mag/100 days is higher than the decline rate from the radioactive decay of⁵⁶Co of ∼1.2 mag/100 days. The most plausible explanation for this discrepancy is that flux is shifting to >14μm, and future JWST observations of SNe Ia will be able to directly test this hypothesis. However, models predicting nonradiative energy loss cannot be excluded with the present data.

 

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-JCO (1–0, 2–1 and 4–3) and [Ci] (1–0) mapping, which we use to trace CO excitation viaR₄₂=I_CO(4−3)/I_CO(2−1)andR₂₁=I_CO(2−1)/I_CO(1−0)and dissociation viaR_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 (lowerR₄₂) and increased signatures of dissociation (higherR_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 compareR₄₂andR_CICOwith 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_CICOexhibits ∼0.1 dex tighter correlations thanR₄₂, suggestingCito 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.

 

JWST observations of polycyclic aromatic hydrocarbon (PAH) emission provide some of the deepest and highest resolution views of the cold interstellar medium (ISM) in nearby galaxies. If PAHs are well mixed with the atomic and molecular gas and illuminated by the average diffuse interstellar radiation field, PAH emission may provide an approximately linear, high-resolution, high-sensitivity tracer of diffuse gas surface density. We present a pilot study that explores using PAH emission in this way based on Mid-Infrared Instrument observations of IC 5332, NGC 628, NGC 1365, and NGC 7496 from the Physics at High Angular resolution in Nearby GalaxieS-JWST Treasury. Using scaling relationships calibrated in Leroy et al., scaled F1130W provides 10–40 pc resolution and 3σsensitivity of Σ_gas∼ 2M_⊙pc⁻². We characterize the surface densities of structures seen at <7M_⊙pc⁻²in our targets, where we expect the gas to be Hi-dominated. We highlight the existence of filaments, interarm emission, and holes in the diffuse ISM at these low surface densities. Below ∼10M_⊙pc⁻²for NGC 628, NGC 1365, and NGC 7496 the gas distribution shows a “Swiss cheese”-like topology due to holes and bubbles pervading the relatively smooth distribution of the diffuse ISM. Comparing to recent galaxy simulations, we observe similar topology for the low-surface-density gas, though with notable variations between simulations with different setups and resolution. Such a comparison of high-resolution, low-surface-density gas with simulations is not possible with existing atomic and molecular gas maps, highlighting the unique power of JWST maps of PAH emission.