Search

The earliest stages of star and cluster formation are hidden within dense cocoons of gas and dust, limiting their detection at optical wavelengths. With the unprecedented infrared capabilities of JWST, we can now observe dust-enshrouded star formation with $\sim$10 pc resolution out to $\sim$20 Mpc. Early findings from PHANGS-JWST suggest that 3.3 $\mu$m polycyclic aromatic hydrocarbon (PAH) emission can identify star clusters in their dust-embedded phases. Here, we extend this analysis to 19 galaxies from the PHANGS-JWST Cycle 1 Treasury Survey, providing the first characterization of compact sources exhibiting 3.3$\mu$m PAH emission across a diverse sample of nearby star-forming galaxies. We establish selection criteria, a median color threshold of F300M-F335M=0.67 at F335M=20, and identify of 1816 sources. These sources are predominantly located in dust lanes, spiral arms, rings, and galaxy centers, with $\sim$87% showing concentration indices similar to optically detected star clusters. Comparison with the PHANGS-HST catalogs suggests that PAH emission fades within $\sim$3 Myr. The H$\alpha$ equivalent width of PAH emitters is 1-2.8 times higher than that of young PHANGS-HST clusters, providing evidence that PAH emitters are on average younger. Analysis of the bright portions of luminosity functions (which should not suffer from incompleteness) shows that young dusty clusters may increase the number of optically visible $\leq$ 3 Myr-old clusters in PHANGS-HST by a factor between $\sim$1.8x-8.5x., Comment: Submitted to ApJ

The CO-to-H$_2$ conversion factor ($\alpha_\mathrm{CO}$) is expected to vary with dust abundance and grain size distribution through the efficiency of shielding gas from CO-dissociation radiation. We present a comprehensive analysis of $\alpha_\mathrm{CO}$ and grain size distribution for nearby galaxies, using the PAH fraction ($q_\mathrm{PAH}$) as an observable proxy of grain size distribution. We adopt the resolved observations at 2-kpc resolution in 42 nearby galaxies, where $\alpha_\mathrm{CO}$ is derived from measured metallicity and surface densities of dust and HI assuming a fixed dust-to-metals ratio. We use an analytical model for the evolution of H$_2$ and CO, in which the evolution of grain size distribution is controlled by the dense gas fraction ($\eta$). We find that the observed level of $q_\mathrm{PAH}$ is consistent with the diffuse-gas-dominated model ($\eta=0.2$) where dust shattering is more efficient. Meanwhile, the slight decreasing trend of observed $q_\mathrm{PAH}$ with metallicity is more consistent with high-$\eta$ predictions, likely due to the more efficient loss of PAHs by coagulation. We discuss how grain size distribution (indicated by $q_\mathrm{PAH}$) and metallicity impact $\alpha_\mathrm{CO}$; we however did not obtain conclusive evidence that the grain size distribution affects $\alpha_\mathrm{CO}$. Observations and model predictions show similar anti-correlation between $\alpha_\mathrm{CO}$ and 12+log(O/H). Meanwhile, there is a considerable difference in how resolved $\alpha_\mathrm{CO}$ behaves with $q_\mathrm{PAH}$. The observed $\alpha_\mathrm{CO}$ has a positive correlation with $q_\mathrm{PAH}$, while the model-predicted $\alpha_\mathrm{CO}$ does not have a definite correlation with $q_\mathrm{PAH}$. This difference is likely due to the limitation of one-zone treatment in the model., Comment: 12 pages, 7 figures, accepted for publication in MNRAS

The James Webb Space Telescope (JWST) has revealed unexpectedly rapid galaxy assembly in the early universe, in tension with models of star and galaxy formation. In the gas conditions typical of early galaxies, particularly their low abundances of heavy elements (metals) and dust, the star-formation process is poorly understood. Some models predict that stars form in atomic gas at low metallicity, in contrast to forming in molecular gas as observed in higher-metallicity galaxies. To understand the very high star-formation rates at early epochs, it is necessary to determine whether molecular gas formation represents a bottleneck to star formation, or if it is plentiful even at extremely low metallicity. Despite repeated searches, star-forming molecular gas has not yet been observed in any galaxy below 7% of the Solar metallicity, leaving the question of how stars form at lower metallicities unresolved. Here, we report the detection of rotationally excited emission from molecular hydrogen in the star-forming region of the nearby, 3% Solar metallicity galaxy Leo P with the MIRI-MRS instrument onboard JWST. These observations place a lower limit on the molecular gas content of Leo P and, combined with our upper limit on carbon monoxide emission from a deep search of this galaxy, demonstrate that MIRI-MRS is sensitive to much smaller molecular gas masses at extremely low metallicity compared to the traditional observational tracer. This discovery pushes the maximum metallicity at which purely atomic gas may fuel star formation a factor of two lower, providing crucial empirical guidance for models of star formation in the early universe., Comment: Under review

By analyzing the spectral energy distributions (SEDs) of resolved stars in nearby galaxies, we can constrain their stellar properties and line-of-sight dust extinction. From the Scylla survey, we obtain ultraviolet to near-infrared photometry from Wide Field Camera 3 onboard the {\it Hubble Space Telescope} for more than 1.5 million stars in the SMC and LMC. We use the Bayesian Extinction and Stellar Tool (BEAST) to analyze the multi-band SEDs of these sources and characterize their initial masses, ages, metallicities, distances, and line-of-sight extinction properties (e.g.~$A_V$, $R_V$). We apply quality cuts and perform validation simulations to construct a catalog of over 550,000 stars with high-reliability SED fits, which we use to analyze the stellar content and extinction properties of the SMC and LMC. We detect stars with masses as low as 0.6 $M_{\odot}$. BEAST stellar age distributions show a jump in observed stars around 6 Gyrs ago, which agrees with star-formation histories. Extinctions ($A_V$) in both galaxies follow a log-normal distribution. We compare $A_V$ with ancillary gas and dust tracers like $HI$, $H_\alpha$, and far infrared (FIR) dust emission and find positive correlations on a field-by-field basis. We convert observed $A_V$ to predicted dust surface densities using the Draine et. al. (2014) model and find $A_V$-based dust surface densities are a factor of $\sim$2.5 lower than observed FIR-based dust surface densities, a correction factor similar to other studies., Comment: Accepted to ApJ, updated to 38 pages, 20 figures

The proximity of the Large and Small Magellanic Clouds (LMC and SMC) provides the opportunity to study the impact of dwarf-dwarf interactions on their mass assembly with a unique level of detail. To this end, we analyze two-filter broadband imaging of 83 Hubble Space Telescope (HST) pointings covering 0.203 deg$^2$ towards the SMC, extending out to $\sim$3.5 kpc in projection from its optical center. Lifetime star formation histories (SFHs) fit to each pointing independently reveal an outside-in age gradient such that fields in the SMC outskirts are older on average. We measure radial gradients of the lookback time to form 90%, 75% and 50% of the cumulative stellar mass for the first time, finding $\delta$($\tau_{90}$, $\tau_{75}$, $\tau_{50}$)/$\delta$R = (0.61$^{+0.08}_{-0.07}$, 0.65$^{+0.09}_{-0.08}$, 0.82$^{+0.12}_{-0.16}$) Gyr/kpc assuming PARSEC evolutionary models and a commonly used elliptical geometry of the SMC, although our results are robust to these assumptions. The wing of the SMC deviates from this trend, forming 25\% of its cumulative mass over the most recent 3 Gyr due to a best-fit star formation rate that remains approximately constant. Our results are consistent with chemodynamical evidence of a tidally stripped SMC component in the foreground, and imply contributions to the observed SFH from multiple previous LMC-SMC interactions. We also compare our SMC SFH with results from a companion study of the LMC, finding that while the two galaxies present different internal, spatially resolved SFH trends, both the LMC and SMC have similar near-constant lifetime SFHs when viewed globally., Comment: ApJ in press. 40 pages, 18 figures, 5 tables including Appendices

Scylla is a deep Hubble Space Telescope survey of the stellar populations, interstellar medium and star formation in the LMC and SMC. As a pure-parallel complement to the Ultraviolet Legacy Library of Young Stars as Essential Standards (ULLYSES) survey, Scylla obtained 342 orbits of ultraviolet (UV) through near-infrared (IR) imaging of the LMC and SMC with Wide Field Camera 3. In this paper, we describe the science objectives, observing strategy, data reduction procedure, and initial results from our photometric analysis of 96 observed fields. Although our observations were constrained by ULYSSES primary exposures, we imaged all fields in at least two filters (F475W and F814W), and 64% of fields in at least three and as many as seven WFC3 filters spanning the UV to IR. Overall, we reach average 50% completeness of $m_{\rm F225W}=26.0$, $m_{\rm F275W}=26.2$, $m_{\rm F336W}=26.9$, $m_{\rm F475W}=27.8$, $m_{\rm F814W}=25.5$, $m_{\rm F110W}=24.7$, and $m_{\rm F160W}=24.0$ Vega magnitudes in our photometric catalogs, which is faintward of the ancient main sequence turnoff in all filters. The primary science goals of Scylla include characterizing the structure and properties of dust in the MCs, as well as their spatially-resolved star formation and chemical enrichment histories. Our images and photometric catalogs, which represent the widest-area coverage of MCs with HST photometry to date, are available as a high-level science product at the Barbara A. Mikulski Archive for Space Telescopes., Comment: 25 pages, 16 figures, 8 tables. Accepted for publication in ApJS

We present resolved $3.6-250~\mu$m dust spectral energy distribution (SED) fitting for $\sim 800$ nearby galaxies. We measure the distribution of radiation field intensities heating the dust, the dust mass surface density ($\Sigma_{\rm d}$), and the fraction of dust in the form of polycyclic aromatic hydrocarbons (PAHs; $q_{\rm PAH}$). We find that the average interstellar radiation field ($\overline{U}$) is correlated both with stellar mass surface density ($\Sigma_{\star}$) and star formation rate surface density ($\Sigma_{\rm SFR}$), while more intense radiation fields are only correlated with $\Sigma_{\rm SFR}$. We show that $q_{\rm PAH}$ is a steeply decreasing function of $\Sigma_{\rm SFR}$, likely reflecting PAH destruction in H II regions. Galaxy integrated $q_{\rm PAH}$ is strongly, negatively correlated with specific star formation rate (sSFR) and offset from the star-forming ``main sequence'' ($\Delta$MS), suggesting that both metallicity and star formation intensity play a role in setting the global $q_{\rm PAH}$. We also find a nearly constant M$_{\rm d}$/M$_\star$ ratio for galaxies on the main sequence, with a lower ratio for more quiescent galaxies, likely due to their lower gas fractions. From these results, we construct prescriptions to estimate the radiation field distribution in both integrated and resolved galaxies. We test these prescriptions by comparing our predicted $\overline{U}$ to results of SED fitting for stacked "main sequence" galaxies at $0

Over the last decade, the Atacama Large Millimeter Array has revealed massive, dusty star-forming galaxies at $z\gtrsim5$, and the James Webb Space Telescope is primed to uncover even more information about them. These observations need dust evolution theory to provide context and are excellent benchmarks to test this theory. Here, we investigate the evolution of galactic dust budget at cosmic dawn using a suite of cosmological zoom-in simulations of moderately massive, high-redshift ($M_{\rm star}\gtrsim10^9 M_{\odot}$; $z\gtrsim5$) galaxies from the FIRE project, the highest resolution ($m_{\rm b} \approx 7100\, M_{\odot}$) of such simulations to date. Our simulations incorporate a dust evolution model that accounts for the dominant sources of dust production, growth, and destruction and follows the evolution of specific dust species, allowing it to replicate a wide range of present-day observations. We find, similar to other theoretical works, that dust growth via gas-dust accretion is the dominant producer of dust mass for these massive, $z\gtrsim 5$ galaxies. However, our fiducial model produces $M_{\rm dust}$ that fall ${\gtrsim}1$ dex below observations at any given $M_{\rm star}$ (typical uncertainties are ${\sim}1$ dex), which we attribute to reduced accretion efficiencies caused by a combination of low galactic metallicities and extremely bursty star formation. Modest enhancements (i.e., within observational/theoretical uncertainties) to accretion and SNe II dust creation raise $M_{\rm dust}$ by ${\lesssim}1$ dex, but this still falls below observations which assume $T_{\rm dust}\sim25$ K. One possibility is that inferred dust masses for $z\gtrsim4$ galaxies are overestimated, and recent observational/analytical works that find $T_{\rm dust}\sim50$ K along with metallicity constraints tentatively support this., Comment: 25 pages, 13 figures, accepted for publication in MNRAS. New mock-JWST visualizations and numerous minor revisions and additional citations

We present new JWST observations of the nearby, prototypical edge-on, spiral galaxy NGC 891. The northern half of the disk was observed with NIRCam in its F150W and F277W filters. Absorption is clearly visible in the mid-plane of the F150W image, along with vertical dusty plumes that closely resemble the ones seen in the optical. A $\sim 10 \times 3~{\rm kpc}^2$ area of the lower circumgalactic medium (CGM) was mapped with MIRI F770W at 12 pc scales. Thanks to the sensitivity and resolution of JWST, we detect dust emission out to $\sim 4$ kpc from the disk, in the form of filaments, arcs, and super-bubbles. Some of these filaments can be traced back to regions with recent star formation activity, suggesting that feedback-driven galactic winds play an important role in regulating baryonic cycling. The presence of dust at these altitudes raises questions about the transport mechanisms at play and suggests that small dust grains are able to survive for several tens of million years after having been ejected by galactic winds in the disk-halo interface. We lay out several scenarios that could explain this emission: dust grains may be shielded in the outer layers of cool dense clouds expelled from the galaxy disk, and/or the emission comes from the mixing layers around these cool clumps where material from the hot gas is able to cool down and mix with these cool cloudlets. This first set of data and upcoming spectroscopy will be very helpful to understand the survival of dust grains in energetic environments, and their contribution to recycling baryonic material in the mid-plane of galaxies., Comment: Accepted for publication in Astronomy & Astrophysics; 16 pages, 8 figures

We measure the high-mass stellar initial mass function (IMF) from resolved stars in M33 young stellar clusters. Leveraging \textit{Hubble Space Telescope's} high resolving power, we fully model the IMF probabilistically. We first model the optical CMD of each cluster to constrain its power-law slope $\Gamma$, marginalized over other cluster parameters in the fit (e.g., cluster age, mass, and radius). We then probabilistically model the distribution of MF slopes for a highly strict cluster sample of 9 clusters more massive than log(Mass/M$_{\odot}$)=3.6; above this mass, all clusters have well-populated main sequences of massive stars and should have accurate recovery of their MF slopes, based on extensive tests with artificial clusters. We find the ensemble IMF is best described by a mean high-mass slope of $\overline{\Gamma} = 1.49\pm0.18$, with an intrinsic scatter of $\sigma^{2}_{\Gamma} = 0.02^{+0.16}_{0.00}$, consistent with a universal IMF. We find no dependence of the IMF on environmental impacts such as the local star formation rate or galactocentric radius within M33, which serves as a proxy for metallicity. This $\overline{\Gamma}$ measurement is consistent with similar measurements in M31, despite M33 having a much higher star formation rate intensity. While this measurement is formally consistent with the canonical Kroupa ($\Gamma = 1.30$) IMF, as well as the Salpeter ($\Gamma = 1.35)$) value, it is the second Local Group cluster sample to show evidence for a somewhat steeper high-mass IMF slope. We explore the impacts a steeper IMF slope has on a number of astronomical sub-fields., Comment: Accepted for publication in ApJ. 9 Figures, 1 Table

We present Cloudy modeling of infrared emission lines in the Wolf-Rayet (WR) nebula N76 caused by one of the most luminous and hottest WR stars in the low metallicity Small Magellanic Cloud. We use spatially resolved mid-infrared Spitzer/IRS and far-infrared Herschel/PACS spectroscopy to establish the physical conditions of the ionized gas. The spatially resolved distribution of the emission allows us to constrain properties much more accurately than using spatially integrated quantities. We construct models with a range of constant hydrogen densities between n$_H$ = 4 - 10 cm$^{-3}$ and a stellar wind-blown cavity of 10 pc which reproduces the intensity and shape of most ionized gas emission lines, including the high ionization lines [OIV] and [NeV], as well as [SIII], [SIV], [OIII], and [NeIII]. Our models suggest that the majority of [SiII] emission (91%) is produced at the edge of the HII region around the transition between ionized and atomic gas while very little of the [CII] emission (<5%) is associated with the ionized gas. The physical conditions of N76 are characterized by a hot HII region with a maximum electron temperature of T$_e$ ~ 24,000 K, electron densities that range from n$_e$ ~ 4 to 12 cm$^{-3}$, and high ionization parameters of log(U) ~ -1.15 to -1.77. By analyzing a low metallicty WR nebula with a single ionization source, this work gives valuable insights on the impact WR stars have to the galaxy-integrated ionized gas properties in nearby dwarf galaxies., Comment: 27 pages, 10 figures, 2 tables. Accepted in the Astrophysical Journal

The PHANGS survey uses ALMA, HST, VLT, and JWST to obtain an unprecedented high-resolution view of nearby galaxies, covering millions of spatially independent regions. The high dimensionality of such a diverse multi-wavelength dataset makes it challenging to identify new trends, particularly when they connect observables from different wavelengths. Here we use unsupervised machine learning algorithms to mine this information-rich dataset to identify novel patterns. We focus on three of the PHANGS-JWST galaxies, for which we extract properties pertaining to their stellar populations; warm ionized and cold molecular gas; and Polycyclic Aromatic Hydrocarbons (PAHs), as measured over 150 pc-scale regions. We show that we can divide the regions into groups with distinct multiphase gas and PAH properties. In the process, we identify previously-unknown galaxy-wide correlations between PAH band and optical line ratios and use our identified groups to interpret them. The correlations we measure can be naturally explained in a scenario where the PAHs and the ionized gas are exposed to different parts of the same radiation field that varies spatially across the galaxies. This scenario has several implications for nearby galaxies: (i) The uniform PAH ionized fraction on 150 pc scales suggests significant self-regulation in the ISM, (ii) the PAH 11.3/7.7 \mic~ band ratio may be used to constrain the shape of the non-ionizing far-ultraviolet to optical part of the radiation field, and (iii) the varying radiation field affects line ratios that are commonly used as PAH size diagnostics. Neglecting this effect leads to incorrect or biased PAH sizes., Comment: Main results in figures 6 and 12. Submitted to ApJ, and comments are welcome!

We present NIRCam and NIRISS modules for DOLPHOT, a widely-used crowded field stellar photometry package. We describe details of the modules including pixel masking, astrometric alignment, star finding, photometry, catalog creation, and artificial star tests (ASTs). We tested these modules using NIRCam and NIRISS images of M92 (a Milky Way globular cluster), Draco II (an ultra-faint dwarf galaxy), and WLM (a star-forming dwarf galaxy). DOLPHOT's photometry is highly precise and the color-magnitude diagrams are deeper and have better definition than anticipated during original program design in 2017. The primary systematic uncertainties in DOLPHOT's photometry arise from mismatches in the model and observed point spread functions (PSFs) and aperture corrections, each contributing $\lesssim0.01$ mag to the photometric error budget. Version 1.2 of WebbPSF models, which include charge diffusion and interpixel capacitance effects, significantly reduced PSF-related uncertainties. We also observed minor ($\lesssim0.05$ mag) chip-to-chip variations in NIRCam's zero points, which will be addressed by the JWST flux calibration program. Globular cluster observations are crucial for photometric calibration. Temporal variations in the photometry are generally $\lesssim0.01$ mag, although rare large misalignment events can introduce errors up to 0.08 mag. We provide recommended DOLPHOT parameters, guidelines for photometric reduction, and advice for improved observing strategies. Our ERS DOLPHOT data products are available on MAST, complemented by comprehensive online documentation and tutorials for using DOLPHOT with JWST imaging data., Comment: 30 pages, 17 figures. Accepted to ApJS. Data products to be hosted on MAST. For DOLPHOT/JWST tutorials, see https://dolphot-jwst.readthedocs.io/en/latest/ . For more program and DOLPHOT info, see https://ers-stars.github.io

We present new observations of the central 1 kpc of the M 82 starburst obtained with the James Webb Space Telescope (JWST) near-infrared camera (NIRCam) instrument at a resolution ~0.05"-0.1" (~1-2 pc). The data comprises images in three mostly continuum filters (F140M, F250M, and F360M), and filters that contain [FeII] (F164N), H2 v=1-0 (F212N), and the 3.3 um PAH feature (F335M). We find prominent plumes of PAH emission extending outward from the central starburst region, together with a network of complex filamentary substructure and edge-brightened bubble-like features. The structure of the PAH emission closely resembles that of the ionized gas, as revealed in Paschen alpha and free-free radio emission. We discuss the origin of the structure, and suggest the PAHs are embedded in a combination of neutral, molecular, and photoionized gas., Comment: Submitted to The Astrophysical Journal

We present an investigation of evolved stars in the nearby star-forming galaxy WLM, using NIRCam imaging from the JWST resolved stellar populations early-release science (ERS) program. We find that various combinations of the F090W, F150W, F250M, and F430M filters can effectively isolate red supergiants (RSGs) and thermally-pulsing asymptotic giant branch (TP-AGB) stars from one another, while also providing a reasonable separation of the primary TP-AGB subtypes: carbon-rich C-type stars and oxygen-rich M-type stars. The classification scheme we present here agrees very well with the well-established Hubble Space Telescope (HST) medium-band filter technique. The ratio of C to M-type stars (C/M) is 0.8$\pm$0.1 for both the new JWST and the HST classifications, which is within one sigma of empirical predictions from optical narrow-band CN and TiO filters. The evolved star colors show good agreement with the predictions from the PARSEC$+$COLIBRI stellar evolutionary models, and the models indicate a strong metallicity dependence that makes stellar identification even more effective at higher metallicity. However, the models also indicate that evolved star identification with NIRCam may be more difficult at lower metallicies. We test every combination of NIRCam filters using the models and present additional filters that are also useful for evolved star studies. We also find that $\approx$90\% of the dusty evolved stars are carbon-rich, suggesting that carbonaceous dust dominates the present-day dust production in WLM, similar to the findings in the Magellanic Clouds. These results demonstrate the usefulness of NIRCam in identifying and classifying dust-producing stars without the need for mid-infrared data., Comment: 15 page, 12 figures, submitted to AAS Journals

Observations indicate dust populations vary between galaxies and within them, suggesting a complex life cycle and evolutionary history. Here we investigate the evolution of galactic dust populations across cosmic time using a suite of cosmological zoom-in simulations from the Feedback in Realistic Environments (FIRE) project, spanning $M_{\rm vir}=10^{9-12}M_{\odot};\,M_{*}=10^{6-11}\,M_{\odot}$. Our simulations incorporate a dust evolution model that accounts for the dominant sources of dust production, growth, and destruction and follows the evolution of specific dust species. All galactic dust populations in our suite exhibit similar evolutionary histories, with gas-dust accretion being the dominant producer of dust mass for all but the most metal-poor galaxies. Similar to previous works, we find the onset of efficient gas-dust accretion occurs above a `critical' metallicity threshold ($Z_{\rm crit}$). Due to this threshold, our simulations reproduce observed trends between galactic D/Z and metallicity and element depletion trends in the ISM. However, we find $Z_{\rm crit}$ varies between dust species due to differences in key element abundances, dust physical properties, and life cycle processes resulting in $Z_{\rm crit}\sim0.05Z_{\odot},\,0.2Z_{\odot},\,0.5Z_{\odot}$ for metallic iron, silicates, and carbonaceous dust, respectively. These variations could explain the lack of small carbonaceous grains observed in the Magellanic Clouds. We also find a delay between the onset of gas-dust accretion and when a dust population reaches equilibrium, which we call the equilibrium timescale ($\tau_{\rm eq}$). The relation between $\tau_{\rm eq}$ and the metal enrichment timescale of a galaxy, determined by its recent evolutionary history, can contribute to the scatter in the observed relation between galactic D/Z and metallicity., Comment: 24 pages, 18 figures, accepted for publication in MNRAS. Minor changes

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 James Webb Space Telescope 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 MIRI imaging mosaics at 5.6 and 21 microns that cover a significant portion of one of M33's spiral arms that has existing panchromatic imaging from the Hubble Space Telescope and deep ALMA 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 shows 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., Comment: 12 pages, 10 figures, 1 tables, accepted for publication at MNRAS

We present the Physics at High Angular resolution in Nearby GalaxieS (PHANGS)-AstroSat atlas, which contains ultraviolet imaging of 31 nearby star-forming galaxies captured by the Ultraviolet Imaging Telescope (UVIT) on the AstroSat satellite. The atlas provides a homogeneous data set of far- and near-ultraviolet maps of galaxies within a distance of 22 Mpc and a median angular resolution of 1.4 arcseconds (corresponding to a physical scale between 25 and 160 pc). After subtracting a uniform ultraviolet background and accounting for Milky Way extinction, we compare our estimated flux densities to GALEX observations, finding good agreement. We find candidate extended UV disks around the galaxies NGC 6744 and IC 5332. We present the first statistical measurements of the clumping of the UV emission and compare it to the clumping of molecular gas traced with ALMA. We find that bars and spiral arms exhibit the highest degree of clumping, and the molecular gas is even more clumped than the FUV emission in galaxies. We investigate the variation of the ratio of observed FUV to H$\alpha$ in different galactic environments and kpc-sized apertures. We report that $\sim 65 \%$ varation of the $\log_{10}$(FUV/H$\alpha$) can be described through a combination of dust attenuation with star formation history parameters. The PHANGS-AstroSat atlas enhances the multi-wavelength coverage of our sample, offering a detailed perspective on star formation. When integrated with PHANGS data sets from ALMA, VLT-MUSE, HST and JWST, it develops our comprehensive understanding of attenuation curves and dust attenuation in star-forming galaxies., Comment: 35 pages, 16 figures. The survey webpage is available at https://sites.google.com/view/phangs/home/data/astrosat and the data archive can be accessed at https://www.canfar.net/storage/vault/list/phangs/RELEASES/PHANGS-AstroSat/v1p0

We present the first star formation history (SFH) and age-metallicity relation (AMR) derived from resolved stellar populations imaged with the JWST NIRCam instrument. The target is the Local Group star-forming galaxy WLM at 970 kpc. The depth of the color-magnitude diagram (CMD) reaches below the oldest main sequence turn-off with a SNR=10 at M_F090W=+4.6 mag; this is the deepest CMD for any galaxy that is not a satellite of the Milky Way. We use Hubble Space Telescope (HST) optical imaging that overlaps with the NIRCam observations to directly evaluate the SFHs derived based on data from the two great observatories. The JWST and HST-based SFHs are in excellent agreement. We use the metallicity distribution function measured from stellar spectra to confirm the trends in the AMRs based on the JWST data. Together, these results confirm the efficacy of recovering a SFH and AMR with the NIRCam F090W-F150W filter combination and provide validation of the sensitivity and accuracy of stellar evolution libraries in the near-infrared relative to the optical for SFH recovery work. From the JWST data, WLM shows an early onset to star formation, followed by an extended pause post-reionization before star formation re-ignites, which is qualitatively similar to what has been observed in the isolated galaxies Leo~A and Aquarius. Quantitatively, 15% of the stellar mass formed in the first Gyr, while only 10% formed over the next ~5 Gyr; the stellar mass then rapidly doubled in ~2.5 Gyr, followed by constant star formation over the last ~5 Gyr., Comment: 27 pages, 10 figures, 1 table

We measure the CO-to-H$_2$ conversion factor ($\alpha_\mathrm{CO}$) in 37 galaxies at 2~kpc resolution, using dust surface density inferred from far-infrared emission as a tracer of the gas surface density and assuming a constant dust-to-metals ratio. In total, we have $\sim790$ and $\sim610$ independent measurements of $\alpha_\mathrm{CO}$ for CO (2-1) and (1-0), respectively. The mean values for $\alpha_\mathrm{CO~(2-1)}$ and $\alpha_\mathrm{CO~(1-0)}$ are $9.3^{+4.6}_{-5.4}$ and $4.2^{+1.9}_{-2.0}~M_\odot~pc^{-2}~(K~km~s^{-1})^{-1}$, respectively. The CO-intensity-weighted mean for $\alpha_\mathrm{CO~(2-1)}$ is 5.69, and 3.33 for $\alpha_\mathrm{CO~(1-0)}$. We examine how $\alpha_\mathrm{CO}$ scales with several physical quantities, e.g.\ star-formation rate (SFR), stellar mass, and dust-mass-weighted average interstellar radiation field strength ($\overline{U}$). Among them, $\overline{U}$, $\Sigma_{\rm SFR}$, and integrated CO intensity ($W_\mathrm{CO}$) have the strongest anti-correlation with spatially resolved $\alpha_\mathrm{CO}$. We provide linear regression results to \aco for all quantities tested. At galaxy integrated scales, we observe significant correlations between $\alpha_\mathrm{CO}$ and $W_\mathrm{CO}$, metallicity, $\overline{U}$, and $\Sigma_{\rm SFR}$. We also find that the normalized $\alpha_\mathrm{CO}$ decreases with stellar mass surface density ($\Sigma_\star$) in the high surface density regions ($\Sigma_\star\geq100~{\rm M_\odot~pc^{-2}}$), following the power-law relations $\alpha_\mathrm{CO~(2-1)}\propto\Sigma_\star^{-0.5}$ and $\alpha_\mathrm{CO~(1-0)}\propto\Sigma_\star^{-0.2}$. The power-law index is insensitive to the assumed dust-to-metals ratio. (abridged), Comment: 27 pages, 13 figures, 4 tables, submitted to ApJ

Determining how galactic environment, especially the high gas densities and complex dynamics in bar-fed galaxy centers, alters the star formation efficiency (SFE) of molecular gas is critical to understanding galaxy evolution. However, these same physical or dynamical effects also alter the emissivity properties of CO, leading to variations in the CO-to-H$_2$ conversion factor ($\alpha_\rm{CO}$) that impact the assessment of the gas column densities and thus of the SFE. To address such issues, we investigate the dependence of $\alpha_\rm{CO}$ on local CO velocity dispersion at 150-pc scales using a new set of dust-based $\alpha_\rm{CO}$ measurements, and propose a new $\alpha_\rm{CO}$ prescription that accounts for CO emissivity variations across galaxies. Based on this prescription, we estimate the SFE in a sample of 65 galaxies from the PHANGS-ALMA survey. We find increasing SFE towards high surface density regions like galaxy centers, while using a constant or metallicity-based $\alpha_\rm{CO}$ results in a more homogeneous SFE throughout the centers and disks. Our prescription further reveals a mean molecular gas depletion time of 700 Myr in the centers of barred galaxies, which is overall 3-4 times shorter than in non-barred galaxy centers or the disks. Across the galaxy disks, the depletion time is consistently around 2-3 Gyr regardless of the choice of $\alpha_\rm{CO}$ prescription. All together, our results suggest that the high level of star formation activity in barred centers is not simply due to an increased amount of molecular gas but also an enhanced SFE compared to non-barred centers or disk regions., Comment: 12 pages of main text, 4 figures, accepted for publication in ApJ

HII region electron temperatures are a critical ingredient in metallicity determinations and recent observations reveal systematic variations in the temperatures measured using different ions. We present electron temperatures ($T_e$) measured using the optical auroral lines ([NII]$\lambda5756$, [OII]$\lambda\lambda7320,7330$, [SII]$\lambda\lambda4069,4076$, [OIII]$\lambda4363$, and [SIII]$\lambda6312$) for a sample of HII regions in seven nearby galaxies. We use observations from the Physics at High Angular resolution in Nearby Galaxies survey (PHANGS) obtained with integral field spectrographs on Keck (Keck Cosmic Web Imager; KCWI) and the Very Large Telescope (Multi-Unit Spectroscopic Explorer; MUSE). We compare the different $T_e$ measurements with HII region and interstellar medium environmental properties such as electron density, ionization parameter, molecular gas velocity dispersion, and stellar association/cluster mass and age obtained from PHANGS. We find that the temperatures from [OII] and [SII] are likely over-estimated due to the presence of electron density inhomogeneities in HII regions. We observe that differences between [NII] and [SIII] temperatures are weakly correlated with stellar association mass and molecular gas velocity dispersion. We measure high [OIII] temperatures in a subset of regions with high molecular gas velocity dispersion and low ionization parameter, which may be explained by the presence of low-velocity shocks. In agreement with previous studies, the $T_{\rm{e}}$--$T_{\rm{e}}$ between [NII] and [SIII] temperatures have the lowest observed scatter and generally follow predictions from photoionization modeling, which suggests that these tracers reflect HII region temperatures across the various ionization zones better than [OII], [SII], and [OIII]., Comment: Accepted to ApJ, 66 pages, 36 figures

The CO-to-H$_2$ conversion factor ($\alpha_\rm{CO}$) is central to measuring the amount and properties of molecular gas. It is known to vary with environmental conditions, and previous studies have revealed lower $\alpha_\rm{CO}$ in the centers of some barred galaxies on kpc scales. To unveil the physical drivers of such variations, we obtained ALMA Band 3, 6, and 7 observations toward the inner 2 kpc of NGC 3627 and NGC 4321 tracing $^{12}$CO, $^{13}$CO, and C$^{18}$O lines on 100 pc scales. Our multi-line modeling and Bayesian likelihood analysis of these datasets reveal variations of molecular gas density, temperature, optical depth, and velocity dispersion, which are among the key drivers of $\alpha_\rm{CO}$. The central 300 pc nuclei in both galaxies show strong enhancement of temperature $T_\rm{k}>100$ K and density $n_\rm{H_2}>10^3$ cm$^{-3}$. Assuming a CO-to-H$_2$ abundance of $3\times10^{-4}$, we derive 4-15 times lower $\alpha_\rm{CO}$ than the Galactic value across our maps, which agrees well with previous kpc-scale measurements. Combining the results with our previous work on NGC 3351, we find a strong correlation of $\alpha_\rm{CO}$ with low-J $^{12}$CO optical depths ($\tau_\rm{CO}$), as well as an anti-correlation with $T_\rm{k}$. The $\tau_\rm{CO}$ correlation explains most of the $\alpha_\rm{CO}$ variation in the three galaxy centers, whereas changes in $T_\rm{k}$ influence $\alpha_\rm{CO}$ to second order. Overall, the observed line width and $^{12}$CO/$^{13}$CO 2-1 line ratio correlate with $\tau_\rm{CO}$ variation in these centers, and thus they are useful observational indicators for $\alpha_\rm{CO}$ variation. We also test current simulation-based $\alpha_\rm{CO}$ prescriptions and find a systematic overprediction, which likely originates from the mismatch of gas conditions between our data and the simulations., Comment: Accepted for publication in ApJ; 30 pages of main text + 3 appendices

We investigate how the dust temperature is affected by local environmental quantities, especially dust surface density ($\Sigma_\mathrm{dust}$), dust-to-gas ratio (D/G) and interstellar radiation field. We compile multi-wavelength observations in 46 nearby galaxies, uniformly processed with a common physical resolution of $2~$kpc. A physical dust model is used to fit the infrared dust emission spectral energy distribution (SED) observed with WISE and Herschel. The star formation rate (SFR) is traced with GALEX ultraviolet data corrected by WISE infrared. We find that the dust temperature correlates well with the SFR surface density ($\Sigma_{\rm SFR}$), which traces the radiation from young stars. The dust temperature decreases with increasing D/G at fixed $\Sigma_{\rm SFR}$ as expected from stronger dust shielding at high D/G, when $\Sigma_\mathrm{SFR}$ is higher than $\sim 2\times 10^{-3}~\rm M_\odot~yr^{-1}~kpc^{-2}$. These measurements are in good agreement with the dust temperature predicted by our proposed analytical model. Below this range of $\Sigma_\mathrm{SFR}$, the observed dust temperature is higher than the model prediction and is only weakly dependent on D/G, which is possibly due to the dust heating from old stellar population or the variation of SFR within the past $10^{10}~$yr. Overall, the dust temperature as a function of $\Sigma_\mathrm{SFR}$ and $\Sigma_\mathrm{dust}$ predicted by our analytical model is consistent with observations. We also notice that at fixed gas surface density, $\Sigma_{\rm SFR}$ tends to increase with D/G, i.e. we can empirically modify the Kennicutt-Schmidt law with a dependence on D/G to better match observations., Comment: 16 pages, 8 figures, accepted for publication in MNRAS

We present criteria for separately classifying stars and unresolved background galaxies in photometric catalogs generated with the point spread function (PSF) fitting photometry software DOLPHOT from images taken of Draco II, WLM, and M92 with the Near Infrared Camera (NIRCam) on JWST. Photometric quality metrics from DOLPHOT in one or two filters can recover a pure sample of stars. Conversely, colors formed between short-wavelength (SW) and long-wavelength (LW) filters can be used to effectively identify pure samples of galaxies. Our results highlight that the existing DOLPHOT output parameters can be used to reliably classify stars in our NIRCam data without the need to resort to external tools or more complex heuristics., Comment: 5 pages, 1 figure

We present the JWST Resolved Stellar Populations Early Release Science (ERS) science program. We obtained 27.5 hours of NIRCam and NIRISS imaging of three targets in the Local Group (Milky Way globular cluster M92, ultra-faint dwarf galaxy Draco II, star-forming dwarf galaxy WLM), which span factors of $\sim10^5$ in luminosity, $\sim10^4$ in distance, and $\sim10^5$ in surface brightness. We describe the survey strategy, scientific and technical goals, implementation details, present select NIRCam color-magnitude diagrams (CMDs), and validate the NIRCam exposure time calculator (ETC). Our CMDs are among the deepest in existence for each class of target. They touch the theoretical hydrogen burning limit in M92 ($<0.08$ $M_{\odot}$; SNR $\sim5$ at $m_{F090W}\sim28.2$; $M_{F090W}\sim+13.6$), include the lowest-mass stars observed outside the Milky Way in Draco II (0.09 $M_{\odot}$; SNR $=10$ at $m_{F090W}\sim29$; $M_{F090W}\sim+12.1$), and reach $\sim1.5$ magnitudes below the oldest main sequence turnoff in WLM (SNR $=10$ at $m_{F090W}\sim29.5$; $M_{F090W}\sim+4.6$). The PARSEC stellar models provide a good qualitative match to the NIRCam CMDs, though are $\sim0.05$ mag too blue compared to M92 F090W$-$F150W data. The NIRCam ETC (v2.0) matches the SNRs based on photon noise from DOLPHOT stellar photometry in uncrowded fields, but the ETC may not be accurate in more crowded fields, similar to what is known for HST. We release beta versions of DOLPHOT NIRCam and NIRISS modules to the community. Results from this ERS program will establish JWST as the premier instrument for resolved stellar populations studies for decades to come., Comment: 23 pages, 7 Figures, 3 Tables. Submitted to AAS Journals. Comments welcome

PHANGS-JWST mid-infrared (MIR) imaging of nearby spiral galaxies has revealed ubiquitous filaments of dust emission in intricate detail. We present a pilot study to systematically map the dust filament network (DFN) at multiple scales between 25-400 pc in NGC 628. MIRI images at 7.7, 10, 11.3 and 21$\mu$m of NGC 628 are used to generate maps of the filaments in emission, while PHANGS-HST B-band imaging yields maps of dust attenuation features. We quantify the correspondence between filaments traced by MIR thermal continuum / polycyclic aromatic hydrocarbon (PAH) emission and filaments detected via extinction / scattering of visible light; the fraction of MIR flux contained in the DFN; and the fraction of HII regions, young star clusters and associations within the DFN. We examine the dependence of these quantities with the physical scale at which the DFN is extracted. With our highest resolution DFN maps (25 pc filament width), we find that filaments in emission and attenuation are co-spatial in 40% of sight lines, often exhibiting detailed morphological agreement; that ~30% of the MIR flux is associated with the DFN; and that 75-80% of HII regions and 60% of star clusters younger than 5 Myr are contained within the DFN. However, the DFN at this scale is anti-correlated with looser associations of stars younger than 5 Myr identified using PHANGS-HST near-UV imaging. We discuss the impact of these findings for studies of star formation and the ISM, and the broad range of new investigations enabled with multi-scale maps of the DFN., Comment: 21 pages, 14 figures, accepted for publication as part of PHANGS-JWST ApJL Focus Issue

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 MIRI observations of IC 5332, NGC 628, NGC 1365, and NGC 7496 from the PHANGS-JWST Treasury. Using scaling relationships calibrated in Leroy et al. (2022), scaled F1130W provides 10--40 pc resolution and 3$\sigma$ sensitivity of $\Sigma_{\rm gas} \sim 2$ M$_\odot$ pc$^{-2}$. We characterize the surface densities of structures seen at $< 7$ M$_\odot$ pc$^{-2}$ in our targets, where we expect the gas to be HI-dominated. We highlight the existence of filaments, inter-arm emission, and holes in the diffuse ISM at these low surface densities. Below $\sim 10$ M$_\odot$ pc$^{-2}$ 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 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., Comment: 18 pages, 6 figures, accepted as part of a PHANGS-JWST Focus Issue to appear in ApJ

Large-scale bars can fuel galaxy centers with molecular gas, often leading to the development of dense ring-like structures where intense star formation occurs, forming a very different environment compared to galactic disks. We pair ~0.3" (30pc) resolution new JWST/MIRI imaging with archival ALMA CO(2-1) mapping of the central ~5kpc of the nearby barred spiral galaxy NGC1365, to investigate the physical mechanisms responsible for this extreme star formation. The molecular gas morphology is resolved into two well-known bright bar lanes that surround a smooth dynamically cold gas disk (R_gal ~ 475pc) reminiscent of non-star-forming disks in early type galaxies and likely fed by gas inflow triggered by stellar feedback in the lanes. The lanes host a large number of JWST-identified massive young star clusters. We find some evidence for temporal star formation evolution along the ring. The complex kinematics in the gas lanes reveal strong streaming motions and may be consistent with convergence of gas streamlines expected there. Indeed, the extreme line-widths are found to be the result of inter-`cloud' motion between gas peaks; ScousePy decomposition reveals multiple components with line widths of ~ 19km/s and surface densities of ~ 800M_sun/pc^2, similar to the properties observed throughout the rest of the central molecular gas structure. Tailored hydro-dynamical simulations exhibit many of the observed properties and imply that the observed structures are transient and highly time-variable. From our study of NGC1365, we conclude that it is predominantly the high gas inflow triggered by the bar that is setting the star formation in its CMZ., Comment: 24 pages, 8 figures, accepted for publications as part of PHANGS-JWST ApJL Focus Issue

Polycyclic aromatic hydrocarbons (PAHs) play a critical role in the reprocessing of stellar radiation and in balancing the heating and cooling processes in the interstellar medium (ISM), but appear to be destroyed in HII regions. However, the mechanisms driving their destruction are still not completely understood. Using PHANGS-JWST and PHANGS-MUSE observations, we investigate how the PAH fraction changes in about 1500 HII regions across four nearby star-forming galaxies (NGC 628, NGC 1365, NGC 7496, IC 5332). We find a strong anti-correlation between the PAH fraction and the ionization parameter (the ratio between the ionizing photon flux and the hydrogen density) of HII regions. This relation becomes steeper for more luminous HII regions. The metallicity of HII regions has only a minor impact on these results in our galaxy sample. We find that the PAH fraction decreases with the H$\alpha$ equivalent width - a proxy for the age of the HII regions - although this trend is much weaker than the one identified using the ionization parameter. Our results are consistent with a scenario where hydrogen-ionizing UV radiation is the dominant source of PAH destruction in star-forming regions., Comment: 16 pages, 7 figures, 1 table. Accepted for publication in ApJL as part of a PHANGS-JWST First Results Focus issue

JWST/MIRI imaging of the nearby galaxies IC 5332, NGC 628, NGC 1365 and NGC 7496 from PHANGS reveals a richness of gas structures that in each case form a quasi-regular network of interconnected filaments, shells and voids. We examine whether this multi-scale network of structure is consistent with the fragmentation of the gas disk through gravitational instability. We use FilFinder to detect the web of filamentary features in each galaxy and determine their characteristic radial and azimuthal spacings. These spacings are then compared to estimates of the most Toomre-unstable length (a few kpc), the turbulent Jeans length (a few hundred pc) and the disk scale height (tens of pc) reconstructed using PHANGS-ALMA observations of the molecular gas as a dynamical tracer. Our analysis of the four galaxies targeted in this work indicates that Jeans-scale structure is pervasive. Future work will be essential for determining how the structure observed in gas disks impacts not only the rate and location of star formation but also how stellar feedback interacts positively or negatively with the surrounding multi-phase gas reservoir., Comment: 15 pages, 4 figures. Accepted as part of a PHANGS-JWST ApJL Focus Issue

The PHANGS collaboration has been building a reference dataset for the multi-scale, multi-phase study of star formation and the interstellar medium in nearby galaxies. With the successful launch and commissioning of JWST, we can now obtain high-resolution infrared imaging to probe the youngest stellar populations and dust emission on the scales of star clusters and molecular clouds ($\sim$5-50 pc). In Cycle 1, PHANGS is conducting an 8-band imaging survey from 2-21$\mu$m of 19 nearby spiral galaxies. CO(2-1) mapping, optical integral field spectroscopy, and UV-optical imaging for all 19 galaxies have been obtained through large programs with ALMA, VLT/MUSE, and Hubble. PHANGS-JWST enables a full inventory of star formation, accurate measurement of the mass and age of star clusters, identification of the youngest embedded stellar populations, and characterization of the physical state of small dust grains. When combined with Hubble catalogs of $\sim$10,000 star clusters, MUSE spectroscopic mapping of $\sim$20,000 HII regions, and $\sim$12,000 ALMA-identified molecular clouds, it becomes possible to measure the timescales and efficiencies of the earliest phases of star formation and feedback, build an empirical model of the dependence of small dust grain properties on local ISM conditions, and test our understanding of how dust-reprocessed starlight traces star formation activity, all across a diversity of galactic environments. Here we describe the PHANGS-JWST Treasury survey, present the remarkable imaging obtained in the first few months of science operations, and provide context for the initial results presented in the first series of PHANGS-JWST publications., Comment: Re-submitted after addressing minor comments from referee. To be published as part of PHANGS-JWST ApJL Focus Issue

We use PHANGS-JWST data to identify and classify 1271 compact 21 $\mu$m sources in four nearby galaxies using MIRI F2100W data. We identify sources using a dendrogram-based algorithm, and we measure the background-subtracted flux densities for JWST bands from 2 $\mu$m to 21 $\mu$m. Using the SED in JWST as well as HST bands, plus ALMA and MUSE/VLT observations, we classify the sources by eye. Then we use this classification to define regions in color-color space, and so establish a quantitative framework for classifying sources. We identify 1085 sources as belonging to the ISM of the target galaxies with the remainder being dusty stars or background galaxies. These 21 $\mu$m sources are strongly spatially associated with HII regions ($>92\%$ of sources), while 74$\%$ of sources are coincident with a stellar association defined in the HST data. Using SED fitting, we find that the stellar masses of the 21 $\mu$m sources span a range of 10$^{2}$ to 10$^{4}~M_\odot$ with mass-weighted ages down to 2 Myr. There is a tight correlation between attenuation-corrected H$\alpha$ and 21 $\mu$m luminosity for $L_{\nu,\mathrm{F2100W}}>10^{19}~\mathrm{W~Hz}^{-1}$. Young embedded source candidates selected at 21 $\mu$m are found below this threshold and have $M_\star < 10^{3}~M_\odot$., Comment: 15 pages, 5 figures, Accepted as part of a PHANGS-JWST Focus Issue to appear in ApJ

We present a high-resolution view of bubbles within The Phantom Galaxy (NGC 628); a nearby (~10Mpc), star-forming (~2Msun/yr), face-on (i~9deg) grand-design spiral galaxy. With new data obtained as part of the PHANGS-JWST treasury program, we perform a detailed case-study of two regions of interest, one of which contains the largest and most prominent bubble in the galaxy (The Phantom Void; over 1kpc in diameter), and the other being a smaller region that may be the precursor to such a large bubble (The Precursor Phantom Void). When comparing to matched resolution Halpha observations from the Hubble Space Telescope (HST), we see that the ionized gas is brightest in the shells of both bubbles, and is coincident with the youngest (~1Myr) and most massive (~100,000Msun) stellar associations. We also find an older generation (~20Myr) of stellar associations is present within the bubble of The Phantom Void. From our kinematic analysis of the HI, H2 (CO) and HII gas across The Phantom Void, we infer a high expansion speed of around 15 to 50km/s. The large size and high expansion speed of The Phantom Void suggest that the driving mechanism is sustained stellar feedback due to multiple mechanisms, where early feedback first cleared a bubble (as we observe now in The Precursor Phantom Void), and since then SNe have been exploding within the cavity, and have accelerated the shell. Finally, comparison to simulations shows a striking resemblance to our JWST observations, and suggests that such large-scale stellar feedback-driven bubbles should be common within other galaxies., Comment: 12 pages total, 8 figures, and 1 table. Accepted as part of a PHANGS-JWST Focus Issue to appear in ApJ

The first JWST observations of nearby galaxies have unveiled a rich population of bubbles that trace the stellar feedback mechanisms responsible for their creation. Studying these bubbles therefore allows us to chart the interaction between stellar feedback and the interstellar medium, and the larger galactic flows needed to regulate star formation processes globally. We present the first catalog of bubbles in NGC628, visually identified using MIRI F770W PHANGS-JWST observations, and use them to statistically evaluate bubble characteristics. We classify 1694 structures as bubbles with radii between 6-552 pc. Of these, 31% contain at least one smaller bubble at their edge, indicating that previous generations of star formation have a local impact on where new stars form. On large scales, most bubbles lie near a spiral arm, and their radii increase downstream compared to upstream. Furthermore, bubbles are elongated in a similar direction to the spiral arm ridge-line. These azimuthal trends demonstrate that star formation is intimately connected to the spiral arm passage. Finally, the bubble size distribution follows a power-law of index $p=-2.2\pm0.1$, which is slightly shallower than the theoretical value by 1-3.5$\sigma$ that did not include bubble mergers. The fraction of bubbles identified within the shells of larger bubbles suggests that bubble merging is a common process. Our analysis therefore allows us to quantify the number of star-forming regions that are influenced by an earlier generation, and the role feedback processes have in setting the global star formation rate. With the full PHANGS-JWST sample, we can do this for more galaxies., Comment: 25 pages total, 13 Figures and 1 Table. Accepted for publication in ApJL as part of a PHANGS-JWST First Results Focus issue

We combine JWST observations with ALMA CO and VLT-MUSE H$\alpha$ data to examine off-spiral arm star formation in the face-on, grand-design spiral galaxy NGC 628. We focus on the northern spiral arm, around a galactocentric radius of 3-4 kpc, and study two spurs. These form an interesting contrast, as one is CO-rich and one CO-poor, and they have a maximum azimuthal offset in MIRI 21$\mu$m and MUSE H$\alpha$ of around 40$^\circ$ (CO-rich) and 55$^\circ$ (CO-poor) from the spiral arm. The star formation rate is higher in the regions of the spurs near to spiral arms, but the star formation efficiency appears relatively constant. Given the spiral pattern speed and rotation curve of this galaxy and assuming material exiting the arms undergoes purely circular motion, these offsets would be reached in 100-150 Myr, significantly longer than the 21$\mu$m and H$\alpha$ star formation timescales (both <10 Myr). The invariance of the star formation efficiency in the spurs versus the spiral arms indicates massive star formation is not only triggered in spiral arms, and cannot simply occur in the arms and then drift away from the wave pattern. These early JWST results show that in-situ star formation likely occurs in the spurs, and that the observed young stars are not simply the `leftovers' of stellar birth in the spiral arms. The excellent physical resolution and sensitivity that JWST can attain in nearby galaxies will well resolve individual star-forming regions and help us to better understand the earliest phases of star formation., Comment: 11 pages, 5 Figures, accepted as part of a PHANGS-JWST Focus Issue to appear in ApJ. Updated to include missing author and published paper references

We combine archival HST and new JWST imaging data, covering the ultraviolet to mid-infrared regime, to morphologically analyze the nuclear star cluster (NSC) of NGC 628, a grand-design spiral galaxy. The cluster is located in a 200 pc x 400 pc cavity, lacking both dust and gas. We find roughly constant values for the effective radius (r_eff ~ 5 pc) and ellipticity ({\epsilon} ~ 0.05), while the S\'ersic index (n) and position angle (PA) drop from n ~ 3 to ~ 2 and PA ~ 130{\deg} to 90{\deg}, respectively. In the mid-infrared, r_eff ~ 12pc, {\epsilon} ~ 0.4, and n ~ 1-1.5, with the same PA ~ 90{\deg}. The NSC has a stellar mass of log10 (M_nsc / M_Sun) = 7.06 +- 0.31, as derived through B-V, confirmed when using multi-wavelength data, and in agreement with the literature value. Fitting the spectral energy distribution, excluding the mid-infrared data, yields a main stellar population's age of (8 +- 3) Gyr with a metallicity of Z = 0.012 +- 0.006. There is no indication of any significant star formation over the last few Gyr. Whether gas and dust were dynamically kept out or evacuated from the central cavity remains unclear. The best-fit suggests an excess of flux in the mid-infrared bands, with further indications that the center of the mid-infrared structure is displaced with respect to the optical center of the NSC. We discuss five potential scenarios, none of them fully explaining both the observed photometry and structure., Comment: 26 pages, 10 figures, 6 tables. Accepted for publication by ApJL

We use globular cluster data from the Resolved Stellar Populations Early Release Science (ERS) program to validate the flux calibration for the Near Infrared Camera (NIRCam) on the James Webb Space Telescope (JWST). We find a significant flux offset between the eight short wavelength detectors, ranging from 1-23% (about 0.01-0.2 mag) that affects all NIRCam imaging observations. We deliver improved zeropoints for the ERS filters and show that alternate zeropoints derived by the community also improve the calibration significantly. We also find that the detector offsets appear to be time variable by up to at least 0.1 mag., Comment: Accepted for publication in RNAAS

Recent strides have been made developing dust evolution models for galaxy formation simulations but these approaches vary in their assumptions and degree of complexity. Here we introduce and compare two separate dust evolution models (labelled 'Elemental' and 'Species'), based on recent approaches, incorporated into the GIZMO code and coupled with FIRE-2 stellar feedback and ISM physics. Both models account for turbulent dust diffusion, stellar production of dust, dust growth via gas-dust accretion, and dust destruction from time-resolved supernovae, thermal sputtering in hot gas, and astration. The "Elemental" model tracks the evolution of generalized dust species and utilizes a simple, 'tunable' dust growth routine, while the "Species" model tracks the evolution of specific dust species with set chemical compositions and incorporates a physically motivated, two-phase dust growth routine. We test and compare these models in an idealized Milky Way-mass galaxy and find that while both produce reasonable galaxy-integrated dust-to-metals (D/Z) ratios and predict gas-dust accretion as the main dust growth mechanism, a chemically motivated model is needed to reproduce the observed scaling relation between individual element depletions and D/Z with column density and local gas density. We also find the inclusion of theoretical metallic iron and O-bearing dust species are needed in the case of specific dust species in order to match observations of O and Fe depletions, and the integration of a sub-resolution dense molecular gas/CO scheme is needed to both match observed C depletions and ensure carbonaceous dust is not overproduced in dense environments., Comment: 28 pages, 18 figures, accepted for publication in MNRAS. Added appendix testing variations in ICs

The CO-to-H$_2$ conversion factor ($\alpha_\rm{CO}$) is critical to studying molecular gas and star formation in galaxies. The value of $\alpha_\rm{CO}$ has been found to vary within and between galaxies, but the specific environmental conditions that cause these variations are not fully understood. Previous observations on $\sim$kpc scales revealed low values of $\alpha_\rm{CO}$ in the centers of some barred spiral galaxies, including NGC 3351. We present new ALMA Band 3, 6, and 7 observations of $^{12}$CO, $^{13}$CO, and C$^{18}$O lines on 100 pc scales in the inner $\sim$2 kpc of NGC 3351. Using multi-line radiative transfer modeling and a Bayesian likelihood analysis, we infer the H$_2$ density, kinetic temperature, CO column density per line width, and CO isotopologue abundances on a pixel-by-pixel basis. Our modeling implies the existence of a dominant gas component with a density of $2{-}3\times10^3$ $\rm{cm^{-3}}$ in the central ${\sim}$1 kpc and a high temperature of 30$-$60 K near the nucleus and near the contact points that connect to the bar-driven inflows. Assuming a CO/H$_2$ abundance of $3\times10^{-4}$, our analysis yields $\alpha_\rm{CO}{\sim}0.5{-}2.0$ $\rm{M_\odot\,(K~km~s^{-1}~pc^2)^{-1}}$ with a decreasing trend with galactocentric radius in the central $\sim$1 kpc. The inflows show a substantially lower $\alpha_\rm{CO} < 0.1$ $\rm{M_\odot\,(K~km~s^{-1}~pc^2)^{-1}}$, likely due to lower optical depths caused by turbulence or shear in the inflows. Over the whole region, this gives an intensity-weighted $\alpha_\rm{CO}$ of ${\sim}1.5$ $\rm{M_\odot\,(K~km~s^{-1}~pc^2)^{-1}}$, which is similar to previous dust modeling based results at kpc scales. This suggests that low $\alpha_\rm{CO}$ on kpc scales in the centers of some barred galaxies may be due to the contribution of low optical depth CO emission in bar-driven inflows., Comment: 31 pages, 17 figures; accepted for publication in ApJ; font size of text and figures adjusted

We present the PHANGS-MUSE survey, a programme using the MUSE IFS at the ESO VLT to map 19 massive $(9.4 < \log(M_{*}/M_\odot) < 11.0)$ nearby (D < 20 Mpc) star-forming disc galaxies. The survey consists of 168 MUSE pointings (1'x1' each), a total of nearly 15 Million spectra, covering ~1.5 Million independent spectra. PHANGS-MUSE provides the first IFS view of star formation across different local environments (including galaxy centres, bars, spiral arms) in external galaxies at a median resolution of 50~pc, better than the mean inter-cloud distance in the ionised interstellar medium. This `cloud-scale' resolution allows detailed demographics and characterisations of HII regions and other ionised nebulae. PHANGS-MUSE further delivers a unique view on the associated gas and stellar kinematics, and provides constraints on the star formation history. The PHANGS-MUSE survey is complemented by dedicated ALMA CO(2-1) and multi-band HST observations, therefore allowing us to probe the key stages of the star formation process from molecular clouds to HII regions and star clusters. This paper describes the scientific motivation, sample selection, observational strategy, data reduction and analysis process of the PHANGS-MUSE survey. We present our bespoke automated data-reduction framework, which is built on the reduction recipes provided by ESO, but additionally allows for mosaicking and homogenisation of the point spread function. We further present a detailed quality assessment and a brief illustration of the potential scientific applications of the large set of PHANGS-MUSE data products generated by our data analysis framework. The data cubes and analysis data products described in this paper represent the basis for the first PHANGS-MUSE public data release and are available in the ESO archive and via the Canadian Astronomy Data Centre., Comment: 48 pages, 28 figures, 4 tables, accepted for pub in A&A, Dec 8, 2021. Data products released publicly at this address https://archive.eso.org/scienceportal/home?data_collection=PHANGS

We characterize the local (2-kpc sized) environments of Type Ia, II, and Ib/c supernovae (SNe) that have recently occurred in nearby ($d\lesssim50$ Mpc) galaxies. Using ultraviolet (UV, from GALEX) and infrared (IR, from WISE) maps of 359 galaxies and a sample of 472 SNe, we measure the star formation rate surface density ($\Sigma_{\rm SFR}$) and stellar mass surface density ($\Sigma_\star$) in a 2-kpc beam centered on each SN site. We show that core-collapse SNe are preferentially located along the resolved galactic star-forming main sequence, whereas Type Ia SNe are extended to lower values of $\Sigma_{\rm SFR}$ at fixed $\Sigma_\star$, indicative of locations inside quiescent galaxies or quiescent regions of galaxies. We also test how well the radial distribution of each SN type matches the radial distributions of UV and IR light in each host galaxy. We find that, to first order, the distributions of all types of SNe mirror that of both near-IR light (3.4 and 4.5 microns, tracing the stellar mass distribution) and mid-IR light (12 and 22 microns, tracing emission from hot, small grains), and also resemble our best-estimate $\Sigma_{\rm SFR}$. All types of SNe appear more radially concentrated than the UV emission of their host galaxies. In more detail, the distributions of Type II SNe show small statistical differences from that of near-IR light. We attribute this overall structural uniformity to the fact that within any individual galaxy, $\Sigma_{\rm SFR}$ and $\Sigma_\star$ track one another well, with variations in $\Sigma_{\rm SFR}/\Sigma_\star$ most visible when comparing between galaxies.

We present PHANGS-ALMA, the first survey to map CO J=2-1 line emission at ~1" ~ 100pc spatial resolution from a representative sample of 90 nearby (d<~20 Mpc) galaxies that lie on or near the z=0 "main sequence" of star-forming galaxies. CO line emission traces the bulk distribution of molecular gas, which is the cold, star-forming phase of the interstellar medium. At the resolution achieved by PHANGS-ALMA, each beam reaches the size of a typical individual giant molecular cloud (GMC), so that these data can be used to measure the demographics, life-cycle, and physical state of molecular clouds across the population of galaxies where the majority of stars form at z=0. This paper describes the scientific motivation and background for the survey, sample selection, global properties of the targets, ALMA observations, and characteristics of the delivered ALMA data and derived data products. As the ALMA sample serves as the parent sample for parallel surveys with VLT/MUSE, HST, AstroSat, VLA, and other facilities, we include a detailed discussion of the sample selection. We detail the estimation of galaxy mass, size, star formation rate, CO luminosity, and other properties, compare estimates using different systems and provide best-estimate integrated measurements for each target. We also report the design and execution of the ALMA observations, which combine a Cycle~5 Large Program, a series of smaller programs, and archival observations. Finally, we present the first 1" resolution atlas of CO emission from nearby galaxies and describe the properties and contents of the first PHANGS-ALMA public data release., Comment: 76 pages, 33 figures. Accepted for publication in the Astrophysical Journal Supplement series. Full resolution version and the image atlas to appear as a figure set in the published version can be found https://sites.google.com/view/phangs/publications . Data release coming soon to the ALMA archive and CADC temporarily available at http://phangs.org/data

We describe the processing of the PHANGS-ALMA survey and present the PHANGS-ALMA pipeline, a public software package that processes calibrated interferometric and total power data into science-ready data products. PHANGS-ALMA is a large, high-resolution survey of CO J=2-1 emission from nearby galaxies. The observations combine ALMA's main 12-m array, the 7-m array, and total power observations and use mosaics of dozens to hundreds of individual pointings. We describe the processing of the u-v data, imaging and deconvolution, linear mosaicking, combining interferometer and total power data, noise estimation, masking, data product creation, and quality assurance. Our pipeline has a general design and can also be applied to VLA and ALMA observations of other spectral lines and continuum emission. We highlight our recipe for deconvolution of complex spectral line observations, which combines multiscale clean, single scale clean, and automatic mask generation in a way that appears robust and effective. We also emphasize our two-track approach to masking and data product creation. We construct one set of "broadly masked" data products, which have high completeness but significant contamination by noise, and another set of "strictly masked" data products, which have high confidence but exclude faint, low signal-to-noise emission. Our quality assurance tests, supported by simulations, demonstrate that 12-m+7-m deconvolved data recover a total flux that is significantly closer to the total power flux than the 7-m deconvolved data alone. In the appendices, we measure the stability of the ALMA total power calibration in PHANGS--ALMA and test the performance of popular short-spacing correction algorithms., Comment: Accepted for publication in the Astrophysical Journal Supplement series. 65 pages, 33 figures. Software available at https://github.com/akleroy/phangs_imaging_scripts . For a full resolution version see https://sites.google.com/view/phangs/publications

With a metallicity of 12 + Log(O/H) $\approx$ 7.1-7.2, I Zw 18 is a canonical low-metallicity blue compact dwarf (BCD) galaxy. A growing number of BCDs, including I Zw 18, have been found to host strong, narrow-lined, nebular He II ($\lambda$4686) emission with enhanced intensities compared to H$\beta$ (e.g., He II($\lambda$4686)/H$\beta$ > 1%). We present new observations of I Zw 18 using the Keck Cosmic Web Imager. These observations reveal two nebular He II emission regions (or He III regions) northwest and southeast of the He III region in the galaxy's main body investigated in previous studies. All regions exhibit He II($\lambda4686$)/Hbeta greater than 2%. The two newly resolved He III regions lie along an axis that intercepts the position of I Zw 18's ultraluminous X-ray (ULX) source. We explore whether the ULX could power the two He III regions via shock activity and/or beamed X-ray emission. We find no evidence of shocks from the gas kinematics. If the ULX powers the two regions, the X-ray emission would need to be beamed. Another potential explanation is that a class of early-type nitrogen-rich Wolf-Rayet stars with low winds could power the two He III regions, in which case the alignment with the ULX would be coincidental., Comment: 9 pages, 4 figures, accepted for publication in ApJL

Previous work has argued that atomic gas mass estimates of galaxies from 21 cm HI emission are systematically low due to a cold opaque atomic gas component. If true, this opaque component necessitates a ~35% correction factor relative to the mass from assuming optically-thin HI emission. These mass corrections are based on fitting HI spectra with a single opaque component model that produces a distinct "top-hat" shaped line profile. Here, we investigate this issue using deep, high spectral resolution HI VLA observations of M31 and M33 to test if these top-hat profiles are instead superpositions of multiple HI components along the line-of-sight. We fit both models and find that >80% of the spectra strongly prefer a multi-component Gaussian model while <2% prefer the single opacity-corrected component model. This strong preference for multiple components argues against previous findings of lines-of-sight dominated by only cold HI. Our findings are enabled by the improved spectral resolution (0.42 km/s), whereas coarser spectral resolution blends multiple components together. We also show that the inferred opaque atomic ISM mass strongly depends on the goodness-of-fit definition and is highly uncertain when the inferred spin temperature has a large uncertainty. Finally, we find that the relation of the HI surface density with the dust surface density and extinction has significantly more scatter when the inferred HI opacity correction is applied. These variations are difficult to explain without additionally requiring large variations in the dust properties. Based on these findings, we suggest that the opaque HI mass is best constrained by HI absorption studies., Comment: Accepted in MNRAS. 25 pg, 13 figures

We investigate the relationship between the dust-to-metals ratio (D/M) and the local interstellar medium environment at ~2 kpc resolution in five nearby galaxies: IC342, M31, M33, M101, and NGC628. A modified blackbody model with a broken power-law emissivity is used to model the dust emission from 100 to 500 um observed by Herschel. We utilize the metallicity gradient derived from auroral line measurements in HII regions whenever possible. Both archival and new CO rotational line and HI 21 cm maps are adopted to calculate gas surface density, including new wide field CO and HI maps for IC342 from IRAM and the VLA, respectively. We experiment with several prescriptions of CO-to-H$_2$ conversion factor, and compare the resulting D/M-metallicity and D/M-density correlations, both of which are expected to be non-negative from depletion studies. The D/M is sensitive to the choice of the conversion factor. The conversion factor prescriptions based on metallicity only yield too much molecular gas in the center of IC342 to obtain the expected correlations. Among the prescriptions tested, the one that yields the expected correlations depends on both metallicity and surface density. The 1-$\sigma$ range of the derived D/M spans 0.40-0.58. Compared to chemical evolution models, our measurements suggest that the dust growth time scale is much shorter than the dust destruction time scale. The measured D/M is consistent with D/M in galaxy-integrated studies derived from infrared dust emission. Meanwhile, the measured D/M is systematically higher than the D/M derived from absorption, which likely indicates a systematic offset between the two methods., Comment: 22 pages, 8 figures, 4 tables; accepted by ApJ

We examine the three-dimensional structure and dust extinction properties in a ~ 200 pc $\times$ 100 pc region in the southwest bar of the Small Magellanic Cloud (SMC). We model a deep Hubble Space Telescope optical color-magnitude diagram (CMD) of red clump and red giant branch stars to infer the dust extinction and galactic structure. We model the distance distribution of the stellar component with a Gaussian and find a centroid distance of 65.2 kpc (distance modulus $\mu$ = 19.07 mag) with a FWHM $\approx$ 11.3 kpc. This large extent along the line of sight reproduces results from previous studies using variable stars and red clump stars. Additionally, we find an offset between the stellar and dust distributions, with the dust on the near side relative to the stars by 3.22 $^{+1.69}_{-1.44}$ kpc, resulting in a 73% reddened fraction of stars. Modeling the dust layer with a log-normal $A_V$ distribution indicates a mean extinction $\langle A_V \rangle$ = 0.41 $\pm$ 0.09 mag. We also calculate $A_V/N_H$ = 3.2 - 4.2 $\times10^{-23}$ mag cm$^2$ H$^{-1}$ which is significantly lower than the Milky Way value but is comparable to previous SMC dust-to-gas ratio measurements. Our results yield the first joint dust extinction and 3D geometry properties in a key region in the SMC. This study demonstrates that CMD modeling can be a powerful tool to simultaneously constrain dust extinction and geometry properties in nearby galaxies., Comment: 16 pages, 10 figures. Accepted for publication in ApJ

We demonstrate a new method for measuring the escape fraction of ionizing photons using Hubble Space Telescope imaging of resolved stars in NGC 4214, a local analog of high-redshift starburst galaxies that are thought to be responsible for cosmic reionization. Specifically, we forward model the UV through near-IR spectral energy distributions of $\sim$83,000 resolved stars to infer their individual ionizing flux outputs. We constrain the local escape fraction by comparing the number of ionizing photons produced by stars to the number that are either absorbed by dust or consumed by ionizing the surrounding neutral hydrogen in individual star-forming regions. We find substantial spatial variation in the escape fraction (0-40%). Integrating over the entire galaxy yields a global escape fraction of 25% (+16%/-15%). This value is much higher than previous escape fractions of zero reported for this galaxy. We discuss sources of this apparent tension, and demonstrate that the viewing angle and the 3D ISM geometric effects are the cause. If we assume the NGC 4214 has no internal dust, like many high-redshift galaxies, we find an escape fraction of 59% (an upper-limit for NGC 4214). This is the first non-zero escape fraction measurement for UV-faint (M$_{\rm FUV}$ = -15.9) galaxies at any redshift, and supports the idea that starburst UV-faint dwarf galaxies can provide a sufficient amount of ionizing photons to the intergalactic medium., Comment: 29 pages, 18 figures, 3 tables, accepted for publication in ApJ