Search

We present a comparison of observed polycyclic aromatic hydrocarbon (PAH) feature ratios in 19 nearby galaxies with a grid of theoretical expectations for near- and mid-infrared dust emission. The PAH feature ratios are drawn from Cycle 1 JWST observations and are measured for 7224 stellar clusters and 29176 stellar associations for which we have robust ages and mass estimates from HST five-band photometry. Though there are galaxy-to-galaxy variations, the observed PAH feature ratios largely agree with the theoretical models, particularly those that are skewed toward more ionized and larger PAH size distributions. For each galaxy we also extract PAH feature ratios for 200 pc-wide circular regions in the diffuse interstellar medium, which serve as a non-cluster/association control sample. Compared to what we find for stellar clusters and associations, the 3.3um/7.7um and 3.3um/11.3um ratios from the diffuse interstellar medium are $\sim 0.10-0.15$ dex smaller. When the observed PAH feature ratios are compared to the radiation field hardness as probed by the [OIII]/H$\beta$ ratio, we find anti-correlations for nearly all galaxies in the sample. These results together suggest that the PAH feature ratios are driven by the shape intensity of the radiation field, and that the smallest PAHs -- observed via JWST F335M imaging -- are increasingly 'processed' or destroyed in regions with the most intense and hard radiation fields., Comment: Accepted for publication in AJ

We introduce the Data Analysis Pipeline (DAP) for the Sloan Digital Sky Survey V (SDSS-V) Local Volume Mapper (LVM) project, referred to as the LVM-DAP. We outline our methods for recovering both stellar and emission line components from the optical integral field spectroscopy, highlighting the developments and changes implemented to address specific challenges of the data set. The observations from the LVM project are unique because they cover a wide range of physical resolutions, from approximately 0.05 pc to 100 pc, depending on the distance to the targets. This, along with the varying number of stars sampled in each aperture (ranging from zero, just one of a few, to thousands), presents challenges in using previous spectral synthesis methods and interpreting the spectral fits. We provide a detailed explanation of how we model the stellar content and separate it from the ionized gas emission lines. To assess the accuracy of our results, we compare them with both idealized and more realistic simulations, highlighting the limitations of our methods. We find that the DAP robustly correct for stellar continuum features and recover emission line parameters (e.g. flux, equivalent width, systemtic velocity and velocity dispersion) with a precision and accuracy that fulfill the requirements of the primary goal of the analysis. In addition, the recovered stellar parameters are reliable for single stars, the recovery of integrated populations is less precise. We conclude with a description of the data products we provide, instructions for downloading and using our software, and a showcase illustrating the quality of the data and the analysis on a deep exposure taken on the Huygens region at the center of the Orion Nebula., Comment: 35 pages, 17 figures, 8 tables, accepted for publication in AJ

Combining Atacama Large Millimeter/sub-millimeter Array CO(2-1) mapping and JWST near- and mid-infrared imaging, we characterize the relationship between CO(2-1) and polycyclic aromatic hydrocarbon (PAH) emission at ~100 pc resolution in 66 nearby star-forming galaxies, expanding the sample size from previous ~100 pc resolution studies by more than an order of magnitude. Focusing on regions of galaxies where most of the gas is likely to be molecular, we find strong correlations between CO(2-1) and 3.3 micron, 7.7 micron, and 11.3 micron PAH emission, estimated from JWST's F335M, F770W, and F1130W filters. We derive power law relations between CO(2-1) and PAH emission, which have indices in the range 0.8-1.2, implying relatively weak variations in the observed CO-to-PAH ratios across the regions that we study. We find that CO-to-PAH ratios and scaling relationships near HII regions are similar to those in diffuse sight lines. The main difference between the two types of regions is that sight lines near HII regions show higher intensities in all tracers. Galaxy centers, on the other hand, show higher overall intensities and enhanced CO-to-PAH ratios compared to galaxy disks. Individual galaxies show 0.19 dex scatter in the normalization of CO at fixed I_PAH and this normalization anti-correlates with specific star formation rate (SFR/M*) and correlates with stellar mass. We provide a prescription that accounts for these galaxy-to-galaxy variations and represents our best current empirical predictor to estimate CO(2-1) intensity from PAH emission, which allows one to take advantage of JWST's excellent sensitivity and resolution to trace cold gas., Comment: 21 pages, 4 figures, 3 tables. Submitted to ApJ

Polycyclic aromatic hydrocarbons (PAHs) are a ubiquitous component of the interstellar medium (ISM) in z~0 massive, star-forming galaxies and play key roles in ISM energy balance, chemistry, and shielding. Wide field of view, high resolution mid-infrared (MIR) images from JWST provides the ability to map the fraction of dust in the form of PAHs and the properties of these key dust grains at 10-50 pc resolution in galaxies outside the Local Group. We use MIR JWST photometric observations of a sample of 19 nearby galaxies from the "Physics at High Angular Resolution in Nearby GalaxieS" (PHANGS) survey to investigate the variations of the PAH fraction. By comparison to lower resolution far-IR mapping, we show that a combination of the MIRI filters (R$_{\rm{PAH}}$ = [F770W+F1130W]/F2100W) traces the fraction of dust by mass in the form of PAHs (i.e., the PAH fraction, or q$_{\rm{PAH}}$). Mapping R$_{\rm{PAH}}$ across the 19 PHANGS galaxies, we find that the PAH fraction steeply decreases in HII regions, revealing the destruction of these small grains in regions of ionized gas. Outside HII regions, we find R$_{\rm{PAH}}$ is constant across the PHANGS sample with an average value of 3.43$\pm$0.98, which, for an illuminating radiation field of intensity 2-5 times that of the radiation field in the solar neighborhood, corresponds to q$_{\rm{PAH}}$ values of 3-6%., Comment: Accepted at ApJ, 39 pages, 25 figures

We present the largest extragalactic survey of supernova remnant (SNR) candidates in nearby star-forming galaxies using exquisite spectroscopic maps from MUSE. Supernova remnants exhibit distinctive emission-line ratios and kinematic signatures, which are apparent in optical spectroscopy. Using optical integral field spectra from the PHANGS-MUSE project, we identify SNRs in 19 nearby galaxies at ~ 100~pc scales. We use five different optical diagnostics: (1) line ratio maps of [SII]/H$\alpha$; (2) line ratio maps of [OI]/H$\alpha$; (3) velocity dispersion map of the gas; (4) and (5) two line ratio diagnostic diagrams from BPT diagrams to identify and distinguish SNRs from other nebulae. Given that our SNRs are seen in projection against HII regions and diffuse ionized gas, in our line ratio maps we use a novel technique to search for objects with [SII]/H$\alpha$ or [OI]/H$\alpha$ in excess of what is expected at fixed H$\alpha$ surface brightness within photoionized gas. In total, we identify 2,233 objects using at least one of our diagnostics, and define a subsample of 1,166 high-confidence SNRs that have been detected with at least two diagnostics. The line ratios of these SNRs agree well with the MAPPINGS shock models, and we validate our technique using the well-studied nearby galaxy M83, where all SNRs we found are also identified in literature catalogs and we recover 51% of the known SNRs. The remaining 1,067 objects in our sample are detected with only one diagnostic and we classify them as SNR candidates. We find that ~ 35% of all our objects overlap with the boundaries of HII regions from literature catalogs, highlighting the importance of using indicators beyond line intensity morphology to select SNRs. [OI]/H$\alpha$ line ratio is responsible for selecting the most objects (1,368; 61%), (abridged)., Comment: 35 pages, 24 figures,6 tables, submitted to A&A

We present a statistical analysis of the local, approximately 50-100 pc scale, H-alpha emission at the locations of recent (less than 125 years) supernovae (SNe) in nearby star-forming galaxies. Our sample consists of 32 SNe in 10 galaxies that are targets of the PHANGS-MUSE survey. We find that 41% (13/32) of these SNe occur coincident with a previously identified HII region. For comparison, HII regions cover 32% of the area within 1 kpc of any recent SN. Contrasting this local covering fraction with the fraction of SNe coincident with HII regions, we find a statistical excess of 7.6% +/- 8.7% of all SNe to be associated with HII regions. This increases to an excess of 19.2% +/- 10.4% when considering only core-collapse SNe. These estimates appear to be in good agreement with qualitative results from new, higher resolution HST H-alpha imaging, which also suggest many CCSNe detonate near but not in HII regions. Our results appear consistent with the expectation that only a modest fraction of stars explode during the first 5 Myr of the life of a stellar population, when H-alpha emission is expected to be bright. Of the HII region associated SNe, 8% (11/13) also have associated detected CO(2-1) emission, indicating the presence of molecular gas. The HII region associated SNe have typical Av extinctions approximately equal to 1 mag, consistent with a significant amount of pre-clearing of gas from the region before the SNe explode., Comment: Accepted for publication in ApJ; 33 pages, 13 figures, 3 tables in two-column AASTEX63 format

We present a study of the nearby low-metallicity dwarf galaxy IC 1613, focusing on the search for massive stars and related feedback processes, as well as for faint supernova remnants (SNR) in late stages of evolution. We obtained the deepest images of IC 1613 in the narrow-band H{\alpha}, He II and [S II] emission lines and new long-slit spectroscopy observations using several facilities (6-m BTA, 2.5m SAI MSU, and 150RTT telescopes), in combination with the multi-wavelength archival data from MUSE/VLT, VLA, XMM-Newton, and Swift/XRT. Our deep narrow-band photometry identifies several faint shells in the galaxy, and we further investigate their physical characteristics with the new long-slit spectroscopy observations and the archival multi-wavelength data. Based on energy balance calculations and assumptions about their possible nature, we propose that one of the shells is a possible remnant of a supernova explosion. We study five out of eight Wolf-Rayet (WR) star candidates previously published for this galaxy using the He ii emission line mapping, MUSE/VLT archival spectra, and new long-slit spectra. Our analysis discards the considered WR candidates and finds no new ones. We found P Cyg profiles in H{\alpha} line in two stars, which we classify as Luminous Blue Variable (LBV) star candidates. Overall, the galaxy IC 1613 may have a lower rate of WR star formation than previously suggested.

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 investigated the ionised and atomic gas kinematics and excitation state in the central region of ongoing star formation of the nearby low-metallicity dwarf galaxy Sextans B. The analysis is based on the new observations performed in Ha emission line with high resolution ($R \sim 16000$) scanning Fabry-Perot interferometer at the 6-m BTA SAO RAS telescope, and on the long-slit spectral observations at the 9.2-m SALT and 2.5-m CMO SAI MSU telescopes. Strong non-circular gas motions detected in the studied regions probably resulted from the off-plane gas motions and impact of stellar feedback. We identified six regions of elevated Ha velocity dispersion, five of which exhibit asymmetric or two-component Ha line profiles. Three of these regions are young ($<1.1$ Myr) expanding ($V_{\rm exp} \sim 25-50\ {\rm km\ s^{-1}}$) superbubbles. We argue that at least three regions in the galaxy could be supernova remnants. We conclude that supernovae feedback is the dominant source of energy for superbubbles in Sextans B, which is expected for such a low metallicity, although we cannot rule out a strong impact of pre-supernova feedback for one superbubble., Comment: 16 pages, 6+6 figures, 4 tables. submitted to MNRAS after revision

The exquisite angular resolution and sensitivity of JWST is opening a new window for our understanding of the Universe. In nearby galaxies, JWST observations are revolutionizing our understanding of the first phases of star formation and the dusty interstellar medium. Nineteen local galaxies spanning a range of properties and morphologies across the star-forming main sequence have been observed as part of the PHANGS-JWST Cycle 1 Treasury program at spatial scales of $\sim$5-50pc. Here, we describe pjpipe, an image processing pipeline developed for the PHANGS-JWST program that wraps around and extends the official JWST pipeline. We release this pipeline to the community as it contains a number of tools generally useful for JWST NIRCam and MIRI observations. Particularly for extended sources, pjpipe products provide significant improvements over mosaics from the MAST archive in terms of removing instrumental noise in NIRCam data, background flux matching, and calibration of relative and absolute astrometry. We show that slightly smoothing F2100W MIRI data to 0.9" (degrading the resolution by about 30 percent) reduces the noise by a factor of $\approx$3. We also present the first public release (DR1.1.0) of the pjpipe processed eight-band 2-21 $\mu$m imaging for all nineteen galaxies in the PHANGS-JWST Cycle 1 Treasury program. An additional 55 galaxies will soon follow from a new PHANGS-JWST Cycle 2 Treasury program., Comment: 49 pages (27 in Appendices), 54 Figures (39 in Appendices), 3 Tables. Accepted for publication in ApJS. Updated to match accepted version. Data available at https://archive.stsci.edu/hlsp/phangs/phangs-jwst

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

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

We investigate the ionised gas morphology, excitation properties, and kinematics in 19 nearby star-forming galaxies from the PHANGS-MUSE survey. We directly compare the kinetic energy of expanding superbubbles and the turbulent motions in the interstellar medium with the mechanical energy deposited by massive stars in the form of winds and supernovae, with the aim to answer whether the stellar feedback is responsible for the observed turbulent motions and to quantify the fraction of mechanical energy retained in the superbubbles. Based on the distribution of the flux and velocity dispersion in the H$\alpha$ line, we select 1484 regions of locally elevated velocity dispersion ($\sigma$(H$\alpha$)>45 km/s), including at least 171 expanding superbubbles. We analyse these regions and relate their properties to those of the young stellar associations and star clusters identified in PHANGS-HST data. We find a good correlation between the kinetic energy of the ionised gas and the total mechanical energy input from supernovae and stellar winds from the stellar associations, with a typical coupling efficiency of 10-20%. The contribution of mechanical energy by the supernovae alone is not sufficient to explain the measured kinetic energy of the ionised gas, which implies that pre-supernova feedback in the form of radiation/thermal pressure and winds is necessary. We find that the gas kinetic energy decreases with metallicity for our sample covering Z=0.5-1.0 Zsun, reflecting the lower impact of stellar feedback. For the sample of superbubbles, we find that about 40% of the young stellar associations are preferentially located in their rims. We also find a slightly higher (by ~15%) fraction of the youngest (<3 Myr) stellar associations in the rims of the superbubbles than in the centres, and the opposite for older associations, which implies possible propagation or triggering of star formation., Comment: Accepted by A&A. The abstract is abridged. 31 pages (including 8 pages in appendix), 20 figures, 2 tables

Measurements of the star formation activity on cloud scales are fundamental to uncovering the physics of the molecular cloud, star formation, and stellar feedback cycle in galaxies. Infrared (IR) emission from small dust grains and polycyclic aromatic hydrocarbons (PAHs) are widely used to trace the obscured component of star formation. However, the relation between these emission features and dust attenuation is complicated by the combined effects of dust heating from old stellar populations and an uncertain dust geometry with respect to heating sources. We use images obtained with NIRCam and MIRI as part of the PHANGS--JWST survey to calibrate dust emission at 21$\rm \mu m$, and the emission in the PAH-tracing bands at 3.3, 7.7, 10, and 11.3$\rm \mu m$ as tracers of obscured star formation. We analyse $\sim$ 20000 optically selected HII regions across 19 nearby star-forming galaxies, and benchmark their IR emission against dust attenuation measured from the Balmer decrement. We model the extinction-corrected H$\alpha$ flux as the sum of the observed H$\alpha$ emission and a term proportional to the IR emission, with $a_{IR}$ as the proportionality coefficient. A constant $a_{IR}$ leads to extinction-corrected H$\alpha$ estimates which agree with those obtained with the Balmer decrement with a scatter of $\sim$ 0.1 dex for all bands considered. Among these bands, 21$\rm \mu m$ emission is demonstrated to be the best tracer of dust attenuation. The PAH-tracing bands underestimate the correction for bright HII regions, since in these environments the ratio of PAH-tracing bands to 21$\rm \mu m$ decreases, signalling destruction of the PAH molecules. For fainter HII regions all bands suffer from an increasing contamination from the diffuse infrared background., Comment: accepted for publication in A&A

Connecting the gas in HII regions to the underlying source of the ionizing radiation can help us constrain the physical processes of stellar feedback and how HII regions evolve over time. With PHANGS$\unicode{x2013}$MUSE we detect nearly 24,000 HII regions across 19 galaxies and measure the physical properties of the ionized gas (e.g. metallicity, ionization parameter, density). We use catalogues of multi-scale stellar associations from PHANGS$\unicode{x2013}$HST to obtain constraints on the age of the ionizing sources. We construct a matched catalogue of 4,177 HII regions that are clearly linked to a single ionizing association. A weak anti-correlation is observed between the association ages and the H$\alpha$ equivalent width EW(H$\alpha$), the H$\alpha$/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, EW(H$\alpha$) and log q show the most consistent trends and appear to be most reliable tracers for the age of an HII region., Comment: 15 pages, 12 figures. Accepted for publication in MNRAS

In this work, we present a new catalogue of >40000 ionised nebulae distributed across the 19 galaxies observed by the PHANGS-MUSE survey. The nebulae have been classified using a new model-comparison-based algorithm that exploits the odds ratio principle to assign a probabilistic classification to each nebula in the sample. The resulting catalogue is the largest catalogue containing complete spectral and spatial information for a variety of ionised nebulae available so far in the literature. We developed this new algorithm to address some of the limitations of the traditional classification criteria, such as their binarity, the sharpness of the involved limits, and the limited amount of data they rely on for the classification. The analysis of the catalogue shows that the algorithm performs well when selecting H II regions. We can recover their luminosity function, and its properties are in line with what is available in the literature. We also identify a rather significant population of shock-ionised regions (mostly composed of supernova remnants), an order of magnitude larger than any other homogeneous catalogue of supernova remnants currently available in the literature. The number of supernova remnants we identify per galaxy is in line with results in our Galaxy and other very nearby sources. However, limitations in the source detection algorithm result in an incomplete sample of planetary nebulae, even though their classification seems robust. Finally, we demonstrate how applying a correction for the contribution of the diffuse ionised gas to the nebulae's spectra is essential to obtain a robust classification of the objects and how a correct measurement of the extinction using DIG-corrected line fluxes prompts the use of a higher theoretical Ha/Hb ratio (3.03) than what is commonly used when recovering the E(B-V) via the Balmer decrement technique in massive star-forming galaxies., Comment: 58 pages, 46 figures. Paper accepted for pubblications in A&A. The catalogue will be available via the CDS or at the following link: http://dx.doi.org/10.11570/23.0006

We present new 0.3-21 micron photometry of SN 2021aefx in the spiral galaxy NGC 1566 at +357 days after B-band maximum, including the first detection of any SN Ia at >15 micron. 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 micron decay rate of $\Delta m_{0.3-14} = 1.35 \pm 0.05$ mag/100 days is higher than the decline rate from the radioactive decay of $^{56}$Co of $\sim 1.2$mag/100 days. The most plausible explanation for this discrepancy is that flux is shifting to >14 micron, and future JWST observations of SNe Ia will be able to directly test this hypothesis. However, models predicting non-radiative energy loss cannot be excluded with the present data., Comment: Accepted for publication in ApJL; 11 pages, 4 figures, 2 tables in two-column AASTEX63 format

We present maps tracing the fraction of dust in the form of polycyclic aromatic hydrocarbons (PAHs) in IC 5332, NGC 628, NGC 1365, and NGC 7496 from JWST/MIRI observations. We trace the PAH fraction by combining the F770W ($7.7~\mu$m) and F1130W ($11.3~\mu$m) filters to track ionized and neutral PAH emission, respectively, and comparing the PAH emission to F2100W which traces small, hot dust grains. We find average $R{\rm_{PAH} = (F770W+F1130W)/F2100W}$ values of 3.3, 4.7, 5.1, and 3.6 in IC 5332, NGC 628, NGC 1365, and NGC 7496, respectively. We find that H II regions traced by MUSE H$\alpha$ show a systematically low PAH fraction. The PAH fraction remains relatively constant across other galactic environments, with slight variations. We use CO + H I + H$\alpha$ to trace the interstellar gas phase and find that the PAH fraction decreases above a value of I$_{H\alpha}/\Sigma_{H~I+H_2}$ $\sim~10^{37.5}$ erg s$^{-1}$ kpc$^{-2}$ (M$_\odot$ pc$^{-2}$)$^{-1}$, in all four galaxies. Radial profiles also show a decreasing PAH fraction with increasing radius, correlated with lower metallicity, in line with previous results showing a strong metallicity dependence to the PAH fraction. Our results suggest that the process of PAH destruction in ionized gas operates similarly across the four targets., Comment: 13 pages, 4 figures. Accepted as part of a PHANGS-JWST Focus Issue to appear in ApJ

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

We present maps of the 3.3 micron polycyclic aromatic hydrocarbon (PAH) emission feature in NGC 628, NGC 1365, and NGC 7496 as observed with the Near-Infrared Camera (NIRCam) imager on JWST from the PHANGS-JWST Cycle 1 Treasury project. We create maps that isolate the 3.3 micron PAH feature in the F335M filter (F335M$_{\rm PAH}$) using combinations of the F300M and F360M filters for removal of starlight continuum. This continuum removal is complicated by contamination of the F360M by PAH emission and variations in the stellar spectral energy distribution slopes between 3.0 and 3.6 micron. We modify the empirical prescription from Lai et al. (2020) to remove the starlight continuum in our highly resolved galaxies, which have a range of starlight- and PAH-dominated lines-of-sight. Analyzing radially binned profiles of the F335M$_{\rm PAH}$ emission, we find that between 5-65% of the F335M intensity comes from the 3.3 micron feature within the inner 0.5 $r_{25}$ of our targets. This percentage systematically varies from galaxy to galaxy, and shows radial trends within the galaxies related to each galaxy's distribution of stellar mass, interstellar medium, and star formation. The 3.3 micron emission is well correlated with the 11.3 micron PAH feature traced with the MIRI F1130W filter, as is expected, since both features arise from C-H vibrational modes. The average F335M$_{\rm PAH}$/F1130W ratio agrees with the predictions of recent models by Draine et al. (2021) for PAHs with size and charge distributions shifted towards larger grains with normal or higher ionization., Comment: submitted to AAS journals and revised according to referee comments, part of a PHANGS-JWST Focus Issue to appear in ApJ

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

We compare mid-infrared (mid-IR), extinction-corrected H$\alpha$, 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$\alpha$. At these scales, CO and H$\alpha$ each correlate strongly with mid-IR emission, and these correlations are each stronger than the one relating CO to H$\alpha$ emission. This reflects that mid-IR emission simultaneously acts as a dust column density tracer, leading to the good match with the molecular gas-tracing CO, and as a heating tracer, leading to the good match with the H$\alpha$. By combining mid-IR, CO, and H$\alpha$ 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_\nu = 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$\alpha$-tracing component. The best-fitting models imply that $\sim 50\%$ of the mid-IR flux arises from molecular gas heated by the diffuse interstellar radiation field, with the remaining $\sim 50\%$ associated with bright, dusty star forming regions. We discuss differences between the F770W, F1000W, F1130W bands and the continuum dominated F2100W band and suggest next steps for using the mid-IR as an ISM tracer., Comment: 49 pages, 17 figures, Section 8 provides a detailed summary, accepted for publication in ApJ Letters, part of a PHANGS-JWST Focus Issue to appear in ApJ

Ratios of polycyclic aromatic hydrocarbon (PAH) vibrational bands are a promising tool for measuring the properties of the PAH population and their effect on star formation. The photometric bands of the MIRI and NIRCam instruments on JWST provide the opportunity to measure PAH emission features across entire galaxy disks at unprecedented resolution and sensitivity. Here we present the first results of this analysis in a sample of three nearby galaxies: NGC 628, NGC 1365, and NGC 7496. Based on the variations observed in the 3.3, 7.7, and 11.3 $\mu$m features, we infer changes to the average PAH size and ionization state across the different galaxy environments. High values of F335M$_{\rm PAH}$/F1130W and low values of F1130W/F770W are measured in H II regions in all three galaxies. This suggests that these regions are populated by hotter PAHs, and/or that the PAH ionization fraction is larger. We see additional evidence of heating and/or changes in PAH size in regions with higher molecular gas content as well as increased ionization in regions with higher H$\alpha$ intensity., Comment: 12 pages, 6 figures. Accepted as part of a PHANGS-JWST Focus Issue to appear in ApJ

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 ALMA multi-J CO (1-0, 2-1 and 4-3) and [CI](1-0) mapping, which we use to trace CO excitation via R42 = I_CO(4-3)/I_CO(2-1) and R21 = I_CO(2-1)/I_CO(1-0) and dissociation via RCICO = I_[CI](1-0)/I_CO(2-1) at 330 pc resolution. We find that the gas flowing into the starburst ring from northeast to southwest appears strongly affected by stellar feedback, showing decreased excitation (lower R42) and increased signatures of dissociation (higher RCICO) in the downstream regions. There, radiative transfer modeling suggests that the molecular gas density decreases and temperature and [CI/CO] abundance ratio increase. We compare R42 and RCICO with local conditions across the regions and find that both correlate with near-IR 2 um emission tracing the YMCs and with both PAH (11.3 um) and dust continuum (21 um) emission. In general, RCICO exhibits ~ 0.1 dex tighter correlations than R42, suggesting CI to be a more sensitive tracer of changing physical conditions in the NGC 1365 starburst than CO (4-3). Our results are consistent with a scenario where gas flows into the two arm regions along the bar, becomes condensed/shocked, forms YMCs, and then these YMCs heat and dissociate the gas., Comment: 19 pages, 7 figures and 2 tables in total (12 pages and 6 figures in main text). Accepted as part of a PHANGS-JWST Focus Issue to appear in ApJL

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 comparison of theoretical predictions of dust continuum and polycyclic aromatic hydrocarbon (PAH) emission with new JWST observations in three nearby galaxies: NGC 628, NGC 1365, and NGC 7496. Our analysis focuses on a total of 1063 compact stellar clusters and 2654 stellar associations previously characterized by HST in the three galaxies. We find that the distributions and trends in the observed PAH-focused infrared colors generally agree with theoretical expectations, and that the bulk of the observations is more aligned with models of larger, ionized PAHs. These JWST data usher in a new era of probing interstellar dust and studying how the intense radiation fields near stellar clusters and associations play a role in shaping the physical properties of PAHs., Comment: 8 pages, accepted as part of a PHANGS-JWST Focus Issue to appear in ApJ

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

The earliest stages of star formation, when young stars are still deeply embedded in their natal clouds, represent a critical phase in the matter cycle between gas clouds and young stellar regions. Until now, the high-resolution infrared observations required for characterizing this heavily obscured phase (during which massive stars have formed, but optical emission is not detected) could only be obtained for a handful of the most nearby galaxies. One of the main hurdles has been the limited angular resolution of the Spitzer Space Telescope. With the revolutionary capabilities of the JWST, it is now possible to investigate the matter cycle during the earliest phases of star formation as a function of the galactic environment. In this Letter, we demonstrate this by measuring the duration of the embedded phase of star formation and the implied time over which molecular clouds remain inert in the galaxy NGC628 at a distance of 9.8Mpc, demonstrating that the cosmic volume where this measurement can be made has increased by a factor of $>100$ compared to Spitzer. We show that young massive stars remain embedded for $5.1_{-1.4}^{+2.7}$Myr ($2.3_{-1.4}^{+2.7}$Myr of which being heavily obscured), representing $\sim20\%$ of the total cloud lifetime. These values are in broad agreement with previous measurements in five nearby ($D < 3.5$Mpc) galaxies and constitute a proof of concept for the systematic characterization of the early phase of star formation across the nearby galaxy population with the PHANGS-JWST survey., Comment: 15 pages, 4 figures, 2 tables; accepted for publication in ApJL as part of a PHANGS-JWST First Results Focus issue

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

Mapping star-formation rates (SFR) within galaxies is key to unveiling their assembly and evolution. Calibrations exist for computing SFR from a combination of ultraviolet and infrared bands for galaxies as integrated systems, but their applicability to sub-galactic (kpc) scales remains largely untested. Here we use integral field spectroscopy of 19 nearby ($D <$ 20 Mpc) galaxies obtained by PHANGS-MUSE to derive accurate Balmer decrements (H$\alpha$/H$\beta$) and attenuation-corrected H$\alpha$ maps. We combine this information with mid-infrared maps from WISE at 22 $\rm \mu m$, and ultraviolet maps from GALEX in the far-UV band, to derive SFR surface densities in nearby galaxies on resolved (kpc) scales. Using the H$\alpha$ attenuation-corrected SFR as a reference, we find that hybrid recipes from the literature overestimate the SFR in regions of low SFR surface density, low specific star-formation rate (sSFR), low attenuation and old stellar ages. We attribute these trends to heating of the dust by old stellar populations (IR cirrus). We calibrate this effect by proposing functional forms for the coefficients in front of the IR term which depend on band ratios sensitive to the sSFR. These calibrations prove reliable as a function of physical scale. In particular, they agree within 10% with the attenuation corrections computed from the Balmer decrement on 100 pc scales. Despite small quantitative differences, our calibrations are also applicable to integrated galaxy scales probed by the MaNGA survey, albeit with a larger scatter (up to 0.22 dex). Observations with JWST open up the possibility to calibrate these relations in nearby galaxies with cloud-scale ($\sim$100 pc) resolution mid-IR imaging., Comment: accepted in A&A

Feedback from massive stars shapes the ISM and affects the evolution of galaxies, but its mechanisms acting at the small scales ($\sim 10$ pc) are still not well constrained observationally, especially in the low-metallicity environments. We present the analysis of the ionized gas (focusing on its kinematics, which were never studied before), and its connection to the massive stars in the nearby ($D \sim 1.4$ Mpc) star-forming very metal-poor ($Z\sim 0.07 Z_\odot$) galaxy Sextans A. The analysis is based on the observations with a scanning Fabry-Perot interferometer, long-slit spectroscopy and imaging in emission lines with narrow-band tunable filters. We found 10 expanding superbubbles of ionized gas with ages of 1-3 Myr. We argue that 3 of them are probable supernovae remnants, while the pre-supernovae feedback is an important source of energy for blowing-out the remaining superbubbles. The two brightest sites of star formation exhibit signs of outflowing ionized gas, which is traced by its ionized and atomic gas kinematics and (in one case) by its emission line flux ratios. Overall, the ionized gas kinematics in Sextans A is highly affected by the feedback from several generations of massive stars and inconsistent with the mere solid-body rotation observed in atomic hydrogen., Comment: 19 pages, 12 figures, 5 tables. Accepted for publication in MNRAS

We report the identification of the unusual emission-line stellar-like object in the nearby low-metallicity (Z ~ 0.1Zsun) dwarf galaxy NGC 4068. Our observations performed with long-slit spectrograph and Fabry-Perot interferometer demonstrate high velocity dispersion in Ha line, presence of HeII4686A line and peculiarly low [SII]/[NII] fluxes ratio for this object. From observational data, we derived that the object represents a single star of high bolometric luminosity L~1.5*10^6 Lsun surrounded by an expanding nebula with kinematical age of t ~ 0.5Myr. The nebula exhibits significant nitrogen overabundance (log(N/O) ~ -0.05, that is by ~1.4dex higher than expected for low-metallicity galaxies). We suggested that this is a massive blue supergiant (BSG) or Wolf-Rayet (WR) star surrounded by its ejecta interacting with the interstellar medium. We calculated the models of the nebula using CLOUDY photoionization code, applying CMFGEN-modelled BSG and WR stars as ionisation sources. We found a best agreement between the modelled and observed spectra for the model assuming ionization by low-metallicity WR star of mass M*~80Msun, ionizing the nebula through the strong wind and enriching the interstellar medium with nitrogen., Comment: Minor revision

We study two M 51-type systems Arp 68 and Arp 58, which strongly differ by their stellar masses, gas content and environment. Long-slit spectral observations obtained at the 6-m telescope BTA were used to trace the distributions of a line-of-sight (LOS) velocity and a gas-phase oxygen abundance along the spectral cuts. Two systems are compared by their observed properties. We found a very strong large-scale non-circular motion of gas in both systems and a kpc-size saw-edged velocity profile along the tidal spiral arm of Arp 68, probably caused by the gas outflow due to the stellar feedback. A deep decrease of LOS velocity is also found in the `hinge' region in Arp 58, where the inner spiral arm transforms into the tidal spiral arm, which was predicted earlier for M 51-type galaxies. Local sites of star formation and the satellites are compared with the evolutionary models at the colour-colour diagrams. Unlike the spiral galaxy Arp 58, the main galaxy in Arp 68 system is experiencing an ongoing burst of star formation. Gas-phase metallicity estimates show that Arp 58 has a higher metal abundance and reveals a shallow negative radial gradient of the gas-phase oxygen abundance. The emission gas in Arp 68 has noticeably lower metallicity than it is expected for a given luminosity of this galaxy, which may be connected with its space position in the local void., Comment: 13 pages, 15 figures, accepted for publication in MNRAS

The study of the dynamic properties of bubbles in the interstellar medium is important for understanding the feedback mechanisms from star-formation processes in galaxies. The ongoing integral field spectroscopy of nearby star-forming galaxies reveals many expanding bubbles and superbubbles identified by the local increase in gas velocity dispersion. The limited angular resolution often prevents bona fide measures of the expansion velocities in galaxies outside the Local Group, even despite sufficient spectroscopic resolution. We present a method that makes it possible to measure the expansion velocity of bubbles surrounding massive stars and clusters based on the data about local variations in gas velocity dispersion. We adapt the method for the Fabry-Perot interferometers used with the 6-m telescope of the Special Astrophysical Observatory of the Russian Academy of Sciences, as well as for any integral field spectrograph with a Gaussian line spread function. We apply the method described to analyze the kinematics of ionized superbubbles gas and the only known supernova remnant in the IC 1613 galaxy. The estimate of the kinematic age of the supernova remnant (on the order of 3100 years) agrees well with the previously obtained independent estimate based on X-ray data., Comment: 13 pages, 9 figures, 1 table Accepted for publication in Astrophysical Bulletin

We present the results of a multiwavelength study of the nearby dwarf galaxy DDO 53 - a relatively isolated member of the M 81 group. We analyse the atomic and ionised gas kinematics (based on the observations with Fabry-Perot interferometer in H$\alpha$ line and archival data in HI 21 cm line), distribution, excitation and oxygen abundance of the ionised gas (based on the long-slit and integral-field spectroscopy and on imaging with narrow-band filters), and their relation with the young massive stars (based on archival HST data). We detect a faint 2-kpc sized supershell of ionised gas surrounding the galaxy. Most probably, this structure represents a large-scale gas outflow, however it could be also created by the ionising quanta leaking from star-forming regions to the marginally detected atomic hydrogen surrounding the galactic disc. We analyse the properties of the anomalous HI in the north part of the galaxy and find that its peculiar kinematics is also traced by ionised gas. We argue that this HI feature is related to the accreting gas cloud captured from the intergalactic medium or remaining after the merger event occurred >1 Gyr ago. The infalling gas produces shocks in the interstellar medium and could support the star formation activity in the brightest region in DDO 53., Comment: 20 pages, 13 figures, 4 tables. Accepted by MNRAS

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

The PHANGS program is building the first dataset to enable the multi-phase, multi-scale study of star formation across the nearby spiral galaxy population. This effort is enabled by large survey programs with ALMA, VLT/MUSE, and HST, with which we have obtained CO(2-1) imaging, optical spectroscopic mapping, and high resolution UV-optical imaging, respectively. Here, we present PHANGS-HST, which has obtained five band NUV-U-B-V-I imaging of the disks of 38 spiral galaxies at distances of 4-23 Mpc, and parallel V and I band imaging of their halos, to provide a census of tens of thousands of compact star clusters, and multi-scale stellar associations. The combination of HST, ALMA, and VLT/MUSE observations will yield an unprecedented joint catalog of the observed and physical properties of ~100,000 star clusters, associations, HII regions, and molecular clouds. With these basic units of star formation, PHANGS will systematically chart the evolutionary cycling between gas and stars, across a diversity of galactic environments found in nearby galaxies. We discuss the design of the PHANGS-HST survey, and provide an overview of the HST data processing pipeline and first results, highlighting new methods for selecting star cluster candidates, morphological classification of candidates with convolutional neural networks, and identification of stellar associations over a range of physical scales with a watershed algorithm. We describe the cross-observatory imaging, catalogs, and software products to be released. These high-level science products will seed a broad range of investigations, in particular, the study of embedded stellar populations and dust with JWST, for which a PHANGS Cycle 1 Treasury program to obtain eight band 2-21 $\mu$m imaging has been approved., Comment: ApJS in press. Accepted draft includes moderate revision and updates to the originally posted submitted version. 27 pages, 11 figures, 1 table. Data products released at https://archive.stsci.edu/hlsp/phangs-hst/

We study a peculiar galaxy NGC 90, a pair member of interacting system Arp 65 (NGC 90/93), using the long-slit spectral observations carried out at the Russian 6m telescope BTA and the available SDSS photometric data. This galaxy demonstrates two tidal tails containing young stellar population, being an extension of its `Grand Design' spiral arms. We obtained the distribution of velocity and oxygen abundance of emission gas (O/H) for two slit orientations. In the central part of the galaxy a significant role belongs to non-photoionization mechanism of line emission probably caused by shocks due to LINER-like activity of the nucleus. The O/H has a shallow abundance gradient, typical for interacting galaxies. The most intriguing peculiarity of the galaxy is the presence of the discovered earlier huge HI `cloud' containing about half of total mass of galaxy gas, which is strongly displaced outwards and has a velocity exceeding at about 340 km/s the central velocity of the main galaxy. We found traces of current star formation in the `cloud', even though the cloud is apparently not gravitationally bound with the galaxy. A possible nature of the `cloud' is discussed. We argue that it presents a flow of gas sweeping by ram pressure and elongated along a line of sight., Comment: accepted to MNRAS

Context: GR290 (M33 V0532=Romano's star) is a luminous M33 object undergoing photometric variability typical for luminous blue variable (LBV) stars. It lies inside Wolf-Rayet region in the Hertzsprung-Russell diagram and possesses a WN8 type spectrum at the light minima. Analysis of Gran Telescopio Canarias (GTC) spectra obtained in 2016 led to the conclusion that it is surrounded by an unresolved HII region formed mostly of ejected material from the central star, and disclosed the presence of a second, more extended asymmetrical emission region. Aims: The aim of this paper is to further explore the structure of the nearby environment of GR290. Methods: Long-slit spectra of GR290 were obtained with three slit orientations in the visual and red spectral regions. The emission-line distribution for each slit was analyzed. Results: We confirm the presence of an asymmetric HII region that extends ~50 pc to the south; ~30pc to the north and southeast; ~20 pc to the east and northwest and ~10pc to the west. We also present the first spectrum to be acquired of a star belonging to the neighboring OB88 association, J013501.87+304157.3, which we classify as a B-type supergiant with a possible binary companion., Comment: 11 pages, 10 figures, accepted for publication in A&A

We present results of long-slit and panoramic spectroscopy of extended gaseous disks in 18 nearby S0 galaxies, mostly in groups. The gas in our S0s is found to be often accreted from outside that is implied by its decoupled kinematics: at least 5 galaxies demonstrate strongly inclined large-scale ionized-gas disks smoothly coupled with their outer HI disks, 7 galaxies reveal circumnuclear polar ionized-gas disks, and in NGC 2551 the ionized gas though confined to the main galactic plane however counterrotates the stellar component. The ionized-gas excitation analysis reveals the gas ionization by young stars in 12 of 18 S0 galaxies studied here; the current star formation in these galaxies is confined to the ring-like zones coinciding with the UV-rings. The gas oxygen abundance estimates in the rings are closely concentrated around the value of 0.7 $Z_\odot$ and do not correlate either with the ring radius nor with the metallicity of the underlying stellar population. By applying the tilted-ring analysis to the 2D velocity fields of the ionized gas, we have traced the orientation of the gas rotation-plane lines of nodes along the radius. We have found that current star formation proceeds usually just where the gas lies strictly in the stellar disk planes and rotates there circularly; the sense of the gas rotation does not matter (the counterrotating gas in NGC 2551 form stars currently). In the galaxies without signs of current star formation the extended gaseous disks are either in steady-state quasi-polar orientation (NGC 2655, NGC 2787, NGC 3414, UGC 9519), or are acquired recently through the highly inclined external filaments provoking probably shock-like excitation (NGC 4026, NGC 7280). Our data imply crucial difference of the external-gas accretion regime in S0s with respect to spiral galaxies: the geometry of the gas accretion in S0s is typically off-plane., Comment: 47 pages, re-submitted to ApJS after second revision. Comments welcome!

We study three Arp' systems of peculiar galaxies: Arp 42 (NGC 5829), Arp 82 (NGC 2535/36) and Arp 58 (UGC 4457), using the long-slit spectral observations carried out at the 6m telescope BTA of Special Astrophysical Observatory. Arp 82 and Arp 58 are the M51-type systems. In the third system - Arp 42, there are two extremely luminous kpc-sized clumps observed at the ends of the bifurcated spiral arm, however the source of perturbations remains unknown. From the emission line measurements we analyzed the distribution of gas line-of-sight (LOS) velocity, velocity dispersion and gas-phase metallicity along the slits estimated by different methods. A special attention is paid to the young regions of star formation at the peripheries of the galaxies and in tidal debris, their connection with gas kinematics and abundance. All three systems show the signs of velocity disturbances and the presence of regions of locally enhanced gas velocity dispersion, as well as the regions of diffuse ionized gas. A distribution of oxygen reveals shallow radial abundance gradients, typical for interacting systems. A faint spiral-like branches of tidal bridges which are observed in NGC 2535 and UGC 4457 may represent the remnants of pre-existing "old" mode of spiral waves., Comment: accepted for publication in MNRAS; 15 pages, 14 figures

The ionization state and oxygen abundance distribution in a sample of polar-ring galaxies (PRGs) were studied from the long-slit spectroscopic observations carried out with the SCORPIO-2 focal reducer at the Russian 6-m telescope. The sample consists of 15 PRGs classified as `the best candidates' in the SDSS-based Polar Ring Catalogue. The distributions of line-of-sight velocities of stellar and gaseous components have given kinematic confirmation of polar structures in 13 galaxies in the sample. We show that ionization by young stars dominates in the external parts of polar discs, while shocks have a significant contribution to gas excitation in the inner parts of polar structures. This picture was predicted earlier in a toy model implying the collision between gaseous clouds on polar orbits with the stellar disc gravitational potential well. The exception is a moderately inclined ring to the host galaxy NGC 5014: the accreted gas in the centre has already settled on the main plane and ionized by young stars, while the gas in the internal part of the ring is excited by shocks. The present study three times increases the number of polar structures with an available oxygen abundance estimation. The measured values of the gas metallicity almost do not depend on the galaxy luminosity. The radial [O/H] gradient in the considered polar rings is shallow or absent. No metal-poor gas was detected. We ruled out the scenario of the formation of polar rings due to cold accretion from cosmic filaments for the considered sample of PRGs., Comment: 13 pages, accepted for publication in MNRAS

The formation scenario for giant low surface brightness (gLSB) galaxies with discs as large as 100 kpc still remains unclear. These stellar systems are rare and very hard to observe, therefore a detailed insight on every additional object helps to understand their nature. Here we present a detailed observational study of the gLSB UGC 1922 performed using deep optical imaging and spectroscopic observations combined with archival ultraviolet data. We derived spatially resolved properties of stellar and ionized gas kinematics and characteristics of stellar populations and interstellar medium. We reveal the presence of a kinematically decoupled central component, which counter rotates with respect to the main disc of UGC 1922. The radial metallicity gradient of the ionised gas is in agreement with that found for moderate-size LSB galaxies. At the same time, a slowly rotating and dynamically hot central region of the galaxy hosts a large number of old metal-rich stars, which creates an appearance of a giant elliptical galaxy, that grew an enormous star forming disc. We reproduce most of the observed features of UGC 1922 in N-body/hydrodynamical simulations of an in-plane merger of giant Sa and Sd galaxies. We also discuss alternative formation scenarios of this unusual system., Comment: 15 pages, 10 figures, accepted for publication in MNRAS

We present the results of observations of the galaxy Holmberg I carried out at the Russian 6-m telescope in the narrow-band imaging, long-slit spectroscopy, and scanning Fabry-Perot interferometer modes. A detailed analysis of gas kinematics, ionization conditions, and metallicity of star-forming regions in the galaxy is presented. The aim of the paper is to analyse the propagation of star formation in the galaxy and to understand the role of the ongoing star formation in the evolution of the central 'galaxy-sized' supergiant HI shell (SGS), where all regions of star formation are observed. We show that star formation in the galaxy occurs in large unified complexes rather than in individual giant HII regions. Evidence of the triggered star formation is observed both on scales of individual complexes and of the whole galaxy. We identified two supernova-remnant candidates and one late-type WN star and analysed their spectrum and surrounding-gas kinematics. We provide arguments indicating that the SGS in Holmberg I is destructing by the influence of star formation occurring on its rims., Comment: Accepted by MNRAS; 25 pages, 16 figures, 5 tables

We present the results of optical long-slit spectral observations of interacting system Arp305 carried out at the 6m telescope BTA at the SAO RAS: the radial variation of gas kinematics and oxygen abundance. This study continues the series of spectral observations of tidal debris in the interacting galaxies. Here we pay special attention to the star-forming region between the interacting galaxies - a tidal dwarf galaxy (TDG) candidate. This star-forming system appears to be gravitationally bound or close to this condition. We show that TDG is metal poor in comparison to the parental galaxy. It can indicate that either the origin of star-forming gas of TDG was very far from the centre of parent galaxy or its gas has been diluted by the accretion of metal poor intergalactic gas. Nevertheless the region of the brightest emission clump in TDG where a current star formation takes place reveals a noticeable colour excess evidencing a local gas concentration. We pay attention that TDG is situated in the region of intersection of two gas flows and is in a process of formation probably caused by the collision of these flows. Low velocity difference between TDG, NGC 4017 and the gaseous bridge between them evidences that the lifetime of TDG is restricted by the time of its fall back onto the parent galaxy., Comment: 7 pages, 5 figures, accepted for publication in MNRAS

We studied the nearby edge-on galaxy NGC4656 and its dwarf low surface brightness companion with the enhanced UV brightness, NGC4656UV, belonging to the interacting system NGC4631/56. Regular photometric structure and relatively big size of NGC4656UV allows to consider this dwarf galaxy as a separate group member rather than a tidal dwarf. Spectral long-slit observations were used to obtain the kinematical parameters and gas-phase metallicity of NGC4656UV and NGC4656. Our rough estimate of the total dynamical mass of NGC4656UV allowed us to conclude that this galaxy is the dark-matter dominated LSB dwarf or ultra diffuse galaxy. Young stellar population of NGC4656UV, as well as strong local non-circular gas motions in NGC4656 and the low oxygen gas abundance in the region of this galaxy adjacent to its dwarf companion, give evidence in favour of the accretion of metal-poor gas onto the discs of both galaxies., Comment: 8 pages, 2 tables, 9 figures, accepted for publication in MNRAS