Search

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 present a near infrared (NIR) candidate star cluster catalog for the central kiloparsec of M82 based on new JWST NIRCam images. We identify star cluster candidates using the F250M filter, finding 1357 star cluster candidates with stellar masses $>10^4$ M$_\odot$. Compared to previous optical catalogs, nearly all (87%) of the candidates we identify are new. The star cluster candidates have a median intrinsic cluster radius of $\approx$1 pc and have stellar masses up to $10^6$ M$_\odot$. By comparing the color-color diagram to dust-free yggdrasil stellar population models, we estimate that the star cluster candidates have A$_{\rm V}\sim3-24$ mag, corresponding to A$_{\rm 2.5\mu m}\sim0.3-2.1$ mag. There is still appreciable dust extinction towards these clusters into the NIR. We measure the stellar masses of the star cluster candidates, assuming ages of 0 and 8 Myr. The slope of the resulting cluster mass function is $\beta=1.9\pm0.2$, in excellent agreement with studies of star clusters in other galaxies., Comment: Resubmitted to ApJL

The formation of massive objects via gravitational collapse is relevant both for explaining the origin of the first supermassive black holes and in the context of massive star formation. Here, we analyze simulations of the formation of massive objects pursued by different groups and in various environments, concerning the formation of supermassive black holes, primordial stars, as well as present-day massive stars. We focus particularly on the regime of small virial parameters, i.e., low ratios of the initial kinetic to gravitational energy, low to moderate Mach numbers, and the phase before feedback is very efficient. We compare the outcomes of collapse under different conditions using dimensionless parameters, particularly the star formation efficiency \epsilon_*, the fraction f_* of mass in the most massive object relative to the total stellar mass, and the fraction f_{\rm tot} of mass of the most massive object as a function of the total mass. We find that in all simulations analyzed here, f_{\rm tot} increases as a function of \epsilon_*, although the steepness of the increase depends on the environment. The relation between f_* and \epsilon_* is found to be more complex and also strongly depends on the number of protostars present at the beginning of the simulations. We show that a collision parameter, estimated as the ratio of the system size divided by the typical collision length, allows us to approximately characterize whether collisions will be important. We analyze the statistical correlation between the dimensionless quantities using the Spearman coefficient and confirm via a machine learning analysis that good predictions of f_* can be obtained from \epsilon_* together with a rough estimate of the collision parameter. This suggests that a good estimate of the mass of the most massive object can be obtained once the maximum efficiency for a given environment is known., Comment: 13 pages, 8 figures, 2 tables. Submitted to A&A

The relation between magnetic field strength B and gas density n in the interstellar medium is of fundamental importance to many areas of astrophysics, from protostellar disks to galaxy evolution. We present and compare Bayesian analyses of the B - n relation for a comprehensive observational data set, as well as a large body of numerical MHD simulations. We extend the original Zeeman relation of Crutcher et al. (2010) with a large body of magnetic data that includes 700 observations with the Davis-Chandrasekhar-Fermi method. By using a new multiparameter Bayesian analysis we present a new, more general, time-averaged observational relation: B \propto n^{0.27 \pm 0.017} for n \leq n_0 and B \propto n^{0.54 \pm 0.18} for n \geq n_0, with n_0 = 924^(+145-144) cm^-3. We perform a separate analysis on 19 numerical magnetohydrodynamics simulations that cover a wide range of scales, resolutions, initial conditions, and completed with a variety of codes: arepo, flash, pencil, and ramses. The power law exponents derived from the simulations depend on several physical factors including: dynamo effects, time scales, turbulence, and the initial seed field strength. In particular, early-time simulations where the density, velocity and magnetic fields are unevolved do not match the observational scalings. The simulations that trace the observed density range best, the evolved dwarf galaxy and Milky Way like galaxy simulations, settle into a near consistent exponent of = 0.5 in the dense gas, with variability in the diffuse gas exponent., Comment: 23 figures, 28 pages, submitted to MNRAS, Comments welcome

Aims. Understanding the rotational dynamics of interstellar dust grains is quintessential for the analysis of the observed dust polarization signal. We aim to constrain the set of parameters for an accurate description of the rotational spin-up process of ballistic dust grain aggregates driven by radiative torques (RATs). Methods. We model the dust grains as complex fractal aggregates grown by the ballistic aggregation of uniform spherical particles (monomers) of different sizes. A broad variation of dust materials, shapes, and sizes are studied in the presence of different radiation sources. Results. We find that the canonical parametrization for the torque efficiency overestimates the maximal angular velocity $\omega_{\mathrm{RAT}}$ caused by RATs of grain aggregates, and to resolve this problem we propose a new parametrization which predicts $\omega_{\mathrm{RAT}}$ more accurately. We find that RATs are strongest for larger grains with smaller monomer-density, this manifests as size and monomer-density dependence in the constant part of the parametrization. Following the constant part, the parametrization has two power-laws with different slopes, which retain the universality for all grain sizes. The maximum grain rotation does not scale linearly with radiation strength, due to different drag mechanisms dominating depending on grain material and environment. The angular velocity $\omega_{\mathrm{RAT}}$ of individual single dust grains has a wide distribution and may even differ from the mean by up to two orders of magnitude., Comment: 15 pages, 12 figures

The gravitational wave (GW) signal from binary black hole (BBH) mergers is a promising probe of Population III (Pop III) stars. To fully unleash the power of the GW probe, one important step is to understand the relative importance and features of different BBH evolution channels. We model two channels, isolated binary stellar evolution (IBSE) and nuclear star cluster-dynamical hardening (NSC-DH), in one theoretical framework based on the semi-analytical code A-SLOTH, under various assumptions on Pop III initial mass function (IMF), initial binary statistics and high-$z$ nuclear star clusters (NSCs). The NSC-DH channel contributes $\sim 8-95\%$ of Pop III BBH mergers across cosmic history, with higher contributions achieved by initially wider binary stars, more top-heavy IMFs, and more abundant high-$z$ NSCs. The dimensionless stochastic GW background (SGWB) produced by Pop III BBH mergers has peak values $\Omega^{\rm peak}_{\rm GW}\sim 10^{-11}-8\times 10^{-11}$ around observer-frame frequencies $\nu\sim 10-100\ \rm Hz$. The Pop III contribution can be a non-negligible ($\sim 2-32\%$) component in the total SGWB at $\nu\lesssim 10\ \rm Hz$. The estimated detection rates of Pop III BBH mergers by the Einstein Telescope are $\sim 6-230\ \rm yr^{-1}$ and $\sim 30-1230\ \rm yr^{-1}$ for the NSC-DH and IBSE channels, respectively. Pop III BBH mergers in NSCs are more massive than those from IBSE, so they dominate the Pop III SGWB below $20$ Hz in most cases. Besides, the detection rate of Pop III BBH mergers involving at least one intermediate-mass BH above $100\ \rm M_\odot$ by the Einstein Telescope is $\sim 0.5-200\ \rm yr^{-1}$ in NSCs but remains below $0.1\ \rm yr^{-1}$ for IBSE., Comment: 28 + 5 pages, 25 figures, submitted to MNRAS

Supersonic magnetohydrodynamic (MHD) turbulence is a ubiquitous state for many astrophysical plasmas. However, even the basic statistics for this type of turbulence remains uncertain. We present results from supersonic MHD turbulence simulations at unparalleled resolutions, with plasma Reynolds numbers of over a million. In the kinetic energy spectrum we find a break between the scales that are dominated by kinetic energy, with spectral index $-2$, and those that become strongly magnetized, with spectral index $-3/2$. By analyzing the Helmholtz decomposed kinetic energy spectrum, we find that the compressible modes are not passively mixed through the cascade of the incompressible modes. At high magnetic Reynolds number, above $10^5$, we find a power law in the magnetic energy spectrum with spectral index $-9/5$. On the strongly magnetized, subsonic scales the plasma tends to self-organize into locally relaxed regions, where there is strong alignment between the current density, magnetic field, velocity field and vorticity field, depleting both the nonlinearities and magnetic terms in the MHD equations, which we attribute to plasma relaxation on scales where the magnetic fluctuations evolve on shorter timescales than the velocity fluctuations. This process constrains the cascade to inhomogenous, volume-poor, fractal surfaces between relaxed regions, which has significant repercussions for understanding the nature of magnetized turbulence in astrophysical plasmas and the saturation of the fluctuation dynamo., Comment: 22 pages. 10 figures. Currently the world's largest MHD turbulence simulation on a grid of 10,080^3 (Re and Rm over a million). Comments welcome, email james.beattie@princeton.edu

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

Two main mechanisms have classically been proposed for the formation of runaway stars. The binary supernova scenario (BSS) suggests that a massive star in a binary explodes as a supernova, ejecting its companion. The dynamical ejection scenario suggests that a star is ejected through a strong dynamical encounter between multiple stars. We propose a third mechanism for the formation of runaway stars: the sub-cluster ejection scenario (SCES), where an infalling sub-cluster of stars is ejected out of the cluster by a slingshot interaction with the contracting gravitational potential of the assembling cluster. We demonstrate the SCES in a star-by-star simulation of a young massive cluster forming from a $10^6\rm~M_\odot$ gas cloud using the Torch framework. This star cluster forms hierarchically through a sequence of sub-cluster mergers, determined by the initial turbulent, spherical initial conditions of the gas. We find that these mergers drive the formation of runaway stars in our model. Late-forming sub-clusters fall into the central potential, where they are ejected on a slingshot trajectory, forming groups of runaway stars that are distributed highly anisotropically. Runaways formed by the same SCES share similar ages, velocities, and ejection directions. Surveying observations, we identify several SCES candidate groups with anisotropic ejection directions. The SCES is capable of producing runaway binaries: two wide dynamical binaries in infalling sub-clusters were tightened through ejection. This allows for another velocity kick via subsequent BSS, which is a promising avenue for producing hypervelocity stars unbound to the Galaxy. The SCES occurs when sub-cluster formation is resolved. We expect non-spherical initial gas distributions to increase runaway star numbers. The observation of groups of runaway stars formed via SCES thus reveals the assembly history of their natal clusters., Comment: 11 pages, 10 figures. Submitted to A&A. Comments welcome!

To advance our understanding of the evolution of the interstellar medium (ISM) of our Galaxy, numerical models of Milky Way (MW) type galaxies are widely used. However, most models only vaguely resemble the MW (e.g. in total mass), and often use imposed analytic potentials (which cannot evolve dynamically). This poses a problem in asserting their applicability for the interpretation of observations of our own Galaxy. The goal of this work is to identify a numerical model that is not only a MW-type galaxy, but one that can mimic some of the main observed structures of our Galaxy, using dynamically evolving potentials, so that it can be used as a base model to study the ISM cycle in a galaxy like our own. This paper introduces a suite of 15 MW-type galaxy models developed using the {\sc arepo} numerical code, that are compared to Galactic observations of $^{12}$CO and \ion{H}{I} emission via longitude-velocity plots, from where we extract and compare the skeletons of major galactic features and the terminal gas velocities. We found that our best-fitting model to the overall structure, also reproduces some of the more specific observed features of the MW, including a bar with a pattern speed of $30.0 \pm 0.2$ km\,s$^{-1}$\,kpc$^{-1}$, a bar half-length of $3.2 \pm 0.8$\,kpc. Our model shows large streaming motions around spiral arms, and strong radial motions well beyond the inner bar. This model highlights the complex motions of a dynamic MW-type galaxy and has the potential to offer valuable insight into how our Galaxy regulates the ISM and star formation., Comment: Accepted for publication in MNRAS. 24 pages, 23 figures, 3 tables

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

Spiral arms are some of the most spectacular features in disc galaxies, and also present in our own Milky Way. It has been argued that star formation should proceed more efficiently in spiral arms as a result of gas compression. Yet, observational studies have so far yielded contradictory results. Here we examine arm/interarm surface density contrasts at ~100 pc resolution in 28 spiral galaxies from the PHANGS survey. We find that the arm/interarm contrast in stellar mass surface density (Sigma_*) is very modest, typically a few tens of percent. This is much smaller than the contrasts measured for molecular gas (Sigma_mol) or star formation rate (Sigma_SFR) surface density, which typically reach a factor of ~2-3. Yet, Sigma_mol and Sigma_SFR contrasts show a significant correlation with the enhancement in Sigma_*, suggesting that the small stellar contrast largely dictates the stronger accumulation of gas and star formation. All these contrasts increase for grand-design spirals compared to multi-armed and flocculent systems (and for galaxies with high stellar mass). The median star formation efficiency (SFE) of the molecular gas is 16% higher in spiral arms than in interarm regions, with a large scatter, and the contrast increases significantly (median SFE contrast 2.34) for regions of particularly enhanced stellar contrast (Sigma_* contrast >1.97). The molecular-to-atomic gas ratio (Sigma_mol/Sigma_atom) is higher in spiral arms, pointing to a transformation of atomic to molecular gas. In conclusion, the boost in the star formation efficiency of molecular gas in spiral arms is generally modest or absent, except for locations with exceptionally large stellar contrasts. (abridged), Comment: 26 pages, 16 figures. Accepted for publication in A&A

In this paper we study the filamentary substructure of 3.3 $\mu$m PAH emission from JWST/NIRCam observations in the base of the M82 star-burst driven wind. We identify plume-like substructure within the PAH emission with widths of $\sim$50 pc. Several of the plumes extend to the edge of the field-of-view, and thus are at least 200-300 pc in length. In this region of the outflow, the vast majority ($\sim$70\%) of PAH emission is associated with the plumes. We show that those structures contain smaller scale "clouds" with widths that are $\sim$5-15 pc, and they are morphologically similar to the results of "cloud-crushing" simulations. We estimate the cloud-crushing time-scales of $\sim$0.5-3 Myr, depending on assumptions. We show this time scale is consistent with a picture in which these observed PAH clouds survived break-out from the disk rather than being destroyed by the hot wind. The PAH emission in both the midplane and the outflow is shown to tightly correlate with that of Pa$\alpha$ emission (from HST/NICMOS data), at the scale of both plumes and clouds, though the ratio of PAH-to-Pa$\alpha$ increases at further distances from the midplane. Finally, we show that the outflow PAH emission is suppressed in regions of the M82 wind that are bright in X-ray emission. Overall, our results are broadly consistent with a picture in which cold gas in galactic outflows is launched via hierarchically structured plumes, and those small scale clouds are more likely to survive the wind environment when collected into the larger plume structure., Comment: Submitted to MNRAS

We present the Sloan Digital Sky Survey V (SDSS-V) Local Volume Mapper (LVM). The LVM is an integral-field spectroscopic survey of the Milky Way, Magellanic Clouds, and of a sample of local volume galaxies, connecting resolved pc-scale individual sources of feedback to kpc-scale ionized interstellar medium (ISM) properties. The 4-year survey covers the southern Milky Way disk at spatial resolutions of 0.05 to 1 pc, the Magellanic Clouds at 10 pc resolution, and nearby large galaxies at larger scales totaling $>4300$ square degrees of sky, and more than 55M spectra. It utilizes a new facility of alt-alt mounted siderostats feeding 16 cm refractive telescopes, lenslet-coupled fiber-optics, and spectrographs covering 3600-9800A at R ~ 4000. The ultra-wide field IFU has a diameter of 0.5 degrees with 1801 hexagonally packed fibers of 35.3 arcsec apertures. The siderostats allow for a completely stationary fiber system, avoiding instability of the line spread function seen in traditional fiber feeds. Scientifically, LVM resolves the regions where energy, momentum, and chemical elements are injected into the ISM at the scale of gas clouds, while simultaneously charting where energy is being dissipated (via cooling, shocks, turbulence, bulk flows, etc.) to global scales. This combined local and global view enables us to constrain physical processes regulating how stellar feedback operates and couples to galactic kinematics and disk-scale structures, such as the bar and spiral arms, as well as gas in- and out-flows., Comment: 29 pages, 12 figures, accepted for publication in The Astronomical Journal

Predicting the star formation rate (SFR) in galaxies is crucial to understand their evolution and morphology. To do so requires a fine understanding of how dense structures of gas are created and collapse. In that, turbulence and gravity play a major role. Within the gravo-turbulent framework, we assume that turbulence shapes the ISM, creating density fluctuations that, if gravitationally unstable, will collapse and form stars. The goal of this work is to quantify how different regimes of turbulence, characterized by the strength and compressibility of the driving, shape the density field. We are interested in the outcome in terms of SFR and how it compares with existing analytical models for the SFR. We run a series of hydrodynamical simulations of turbulent gas. The simulations are first conducted without gravity, so that the density and velocity are shaped by the turbulence driving. Gravity is then switched on, and the SFR is measured and compared with analytical models. The physics included in these simulations is very close to the one assumed in the classical gravo-turbulent SFR analytical models, which makes the comparison straightforward. We found that the existing analytical models convincingly agree with simulations at low Mach number, but we measure a much lower SFR in the simulation with a high Mach number. We develop, in a companion paper, an updated physically-motivated SFR model that reproduces well the inefficient high Mach regime of the simulations. Our work demonstrates that accurate estimations of the turbulent-driven replenishment time of dense structures and the dense gas spatial distribution are necessary to correctly predict the SFR in the high Mach regime. The inefficient high-Mach regime is a possible explanation for the low SFR found in dense and turbulent environments such as the centers of our Milky Way and other galaxies., Comment: 13 pages, 14 figures, accepted for publication in 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

Past studies have long emphasised the key role played by galactic stellar bars in the context of disc secular evolution, via the redistribution of gas and stars, the triggering of star formation, and the formation of prominent structures such as rings and central mass concentrations. However, the exact physical processes acting on those structures, as well as the timescales associated with the building and consumption of central gas reservoirs are still not well understood. We are building a suite of hydro-dynamical RAMSES simulations of isolated, low-redshift galaxies that mimic the properties of the PHANGS sample. The initial conditions of the models reproduce the observed stellar mass, disc scale length, or gas fraction, and this paper presents a first subset of these models. Most of our simulated galaxies develop a prominent bar structure, which itself triggers central gas fuelling and the building of an over-density with a typical scale of 100-1000 pc. We confirm that if the host galaxy features an ellipsoidal component, the formation of the bar and gas fuelling are delayed. We show that most of our simulations follow a common time evolution, when accounting for mass scaling and the bar formation time. In our simulations, the stellar mass of $10^{10}$~M$_{\odot}$ seems to mark a change in the phases describing the time evolution of the bar and its impact on the interstellar medium. In massive discs (M$_{\star} \geq 10^{10}$~M$_{\odot}$), we observe the formation of a central gas reservoir with star formation mostly occurring within a restricted starburst region, leading to a gas depletion phase. Lower-mass systems (M$_{\star} < 10^{10}$~M$_{\odot}$) do not exhibit such a depletion phase, and show a more homogeneous spread of star-forming regions along the bar structure, and do not appear to host inner bar-driven discs or rings., Comment: 22 pages, 17 figures

A compact source, G0.02467-0.0727, was detected in ALMA \threemm observations in continuum and very broad line emission. The continuum emission has a spectral index $\alpha\approx3.3$, suggesting that the emission is from dust. The line emission is detected in several transitions of CS, SO, and SO$_2$ and exhibits a line width FWHM $\approx160$ \kms. The line profile appears Gaussian. The emission is weakly spatially resolved, coming from an area on the sky $\lesssim1"$ in diameter ($\lesssim10^4$ AU at the distance of the Galactic Center; GC). The centroid velocity is $v_{LSR}\approx40$-$50$ \kms, which is consistent with a location in the Galactic Center. With multiple SO lines detected, and assuming local thermodynamic equilibrium (LTE) conditions, $T_\mathrm{LTE} = 13$ K, which is colder than seen in typical GC clouds, though we cannot rule out low-density, subthermally excited, warmer gas. Despite the high velocity dispersion, no emission is observed from SiO, suggesting that there are no strong ($\gtrsim10~\mathrm{km~s}^{-1}$) shocks in the molecular gas. There are no detections at other wavelengths, including X-ray, infrared, and radio. We consider several explanations for the Millimeter Ultra-Broad Line Object (MUBLO), including protostellar outflow, explosive outflow, collapsing cloud, evolved star, stellar merger, high-velocity compact cloud, intermediate mass black hole, and background galaxy. Most of these conceptual models are either inconsistent with the data or do not fully explain it. The MUBLO is, at present, an observationally unique object., Comment: Accepted to ApJL

Obtaining reliable distance estimates to gas clouds within the Milky Way is challenging in the absence of certain tracers. The kinematic distance approach has been used as an alternative, derived from the assumption of circular trajectories around the Galactic centre. Consequently, significant errors are expected in regions where gas flow deviates from purely circular motions. We aim to quantify the systematic errors that arise from the kinematic distance method in the presence of a Galactic potential that is non-axisymmetric. We investigate how these errors differ in certain regions of the Galaxy and how they relate to the underlying dynamics. We perform 2D isothermal hydrodynamical simulation of the gas disk with the moving-mesh code Arepo, adding the capability of using an external potential provided by the Agama library for galactic dynamics. We introduce a new analytic potential of the Milky Way, taking elements from existing models and adjusting parameters to match recent observational constraints. We find significant errors in the kinematic distance estimate for gas close to the Sun, along sight lines towards the Galactic centre and anti-centre, and significant deviations associated with the Galactic bar. Kinematic distance errors are low within the spiral arms as gas resides close to local potential minima and the resulting line-of-sight velocity is close to what is expected for an axisymmetric potential. Interarm regions exhibit large deviations at any given Galactic radius. This is caused by the gas being sped up or slowed down as it travels into or out of the spiral arm. We are able to define 'zones of avoidance' in the lv-diagram, where the kinematic distance method is particularly unreliable and should only be used with caution. We report a power law relation between the kinematic distance error and the deviation of the project line-of-sight velocity from circular motion., Comment: 23 pages, 22 figures, 6 tables, submitted to A&A

We present new catalogs of likely globular clusters (GCs) in 17 nearby spiral galaxies studied as part of the PHANGS-HST Treasury Survey. The galaxies were imaged in five broad-band filters from the near-ultraviolet through the $I$ band. PHANGS-HST has produced catalogs of stellar clusters of all ages by selecting extended sources (from multiple concentration index measurements) followed by morphological classification (centrally concentrated and symmetric or asymmetric, multiple peaks, contaminant) by visually examining the V-band image and separately by a machine-learning algorithm which classified larger samples to reach fainter limits. From both cluster catalogs, we select an initial list of candidate GCs to have $B-V \geq 0.5$ and $V-I \geq 0.73$~mag, then remove likely contaminants (including reddened young clusters, background galaxies misclassified by the neural network, and chance superpositions/blends of stars) after a careful visual inspection. We find that $\approx86$ % of the color-selected candidates classified as spherically symmetric, and $\approx68$ of those classified as centrally concentrated but asymmetric are likely to be GCs. The luminosity functions of the GC candidates in 2 of our 17 galaxies, NGC 628 and NGC 3627, are atypical, and continue to rise at least 1~mag fainter than the expected turnover near $M_V \sim -7.4$. These faint candidate GCs have more extended spatial distributions than their bright counterparts, and may reside in the disk rather than the bulge/halo, similar to faint GCs previously discovered in M101. These faint clusters may be somewhat younger since the age-metallicity degeneracy makes it difficult to determine precise cluster ages from integrated colors once they reach $\approx1$~Gyr., Comment: 21 pages, 12 figures

Studying molecular gas in nearby galaxies using hydrogen cyanide (HCN) as a tracer for higher densities than CO emission still poses a significant challenge. Even though several galaxies have HCN maps on a few kpc scales, higher-resolution maps are still required. Our goal is to examine the contrast in intensity between two tracers that probe different density regimes - HCN(1-0)/CO(2-1) ratio - and their kinematics across NGC 253. By utilizing the advanced capabilities of the Atacama Large Millimeter/submillimeter Array (ALMA), we can map these features at high resolution across a large field of view and uncover the nature of such dense gas in extragalactic systems. We present new ALMA Atacama Compact Array and Total Power (ACA+TP) observations of the HCN emission across NGC 253, covering the inner 8.6' of the galaxy disk at 300 pc scales. We analyze the integrated intensity and mean velocity of HCN and CO along each line of sight and use SCOUSE software to perform spectral decomposition, which considers each velocity component separately. Molecular gas traced by HCN piles up in a ring-like structure at a radius of 2 kpc. The HCN emission is enhanced by 2 orders of magnitude in the central 2 kpc regions, beyond which its intensity decreases with increasing galactocentric distance. The number of components in the HCN spectra shows a robust environmental dependence, with multiple velocity features across the center and bar. We have identified an increase in the HCN/CO ratio in these regions, corresponding to a velocity component likely associated with a molecular outflow. We have also discovered that the ratio between the total infrared luminosity and dense gas mass, which indicates the star formation efficiency of dense gas, is anti-correlated with the molecular gas surface density up to approximately 200 Msul/pc^2. In contrast, beyond this point, the ratio starts to increase., Comment: Accepted for publication to Astronomy and Astrophysics

We present the largest catalog to-date of star clusters and compact associations in nearby galaxies. We have performed a V-band-selected census of clusters across the 38 spiral galaxies of the PHANGS-HST Treasury Survey, and measured integrated, aperture-corrected NUV-U-B-V-I photometry. This work has resulted in uniform catalogs that contain $\sim$20,000 clusters and compact associations which have passed human inspection and morphological classification, and a larger sample of $\sim$100,000 classified by neural network models. Here, we report on the observed properties of these samples, and demonstrate that tremendous insight can be gained from just the observed properties of clusters, even in the absence of their transformation into physical quantities. In particular, we show the utility of the UBVI color-color diagram, and the three principal features revealed by the PHANGS-HST cluster sample: the young cluster locus, the middle-age plume, and the old globular cluster clump. We present an atlas of maps of the 2D spatial distribution of clusters and compact associations in the context of the molecular clouds from PHANGS-ALMA. We explore new ways of understanding this large dataset in a multi-scale context by bringing together once-separate techniques for the characterization of clusters (color-color diagrams and spatial distributions) and their parent galaxies (galaxy morphology and location relative to the galaxy main sequence). A companion paper presents the physical properties: ages, masses, and dust reddenings derived using improved spectral energy distribution (SED) fitting techniques., Comment: 48 pages, 26 figures, article in review at APJ

Molecular clouds (MC) are structures of dense gas in the interstellar medium (ISM), that extend from ten to a few hundred parsecs and form the main gas reservoir available for star formation. Hydrodynamical simulations of varying complexity are a promising way to investigate MC evolution and their properties. However, each simulation typically has a limited range in resolution and different cloud extraction algorithms are used, which complicates the comparison between simulations. In this work, we aim to extract clouds from different simulations covering a wide range of spatial scales. We compare their properties, such as size, shape, mass, internal velocity dispersion and virial state. We apply the Hop cloud detection algorithm on (M)HD numerical simulations of stratified ISM boxes and isolated galactic disk simulations that were produced using Flash Ramses and Arepo We find that the extracted clouds are complex in shape ranging from round objects to complex filamentary networks in all setups. Despite the wide range of scales, resolution, and sub-grid physics, we observe surprisingly robust trends in the investigated metrics. The mass spectrum matches in the overlap between simulations without rescaling and with a high-mass slope of $\mathrm{d} N/\mathrm{d}\ln M\propto-1$ in accordance with theoretical predictions. The internal velocity dispersion scales with the size of the cloud as $\sigma\propto R^{0.75}$ for large clouds ($R\gtrsim3\,\mathrm{pc}$). For small clouds we find larger sigma compared to the power-law scaling, as seen in observations, which is due to supernova-driven turbulence. Almost all clouds are gravitationally unbound with the virial parameter scaling as $\alpha_\mathrm{vir}\propto M^{-0.4}$, which is slightly flatter compared to observed scaling, but in agreement given the large scatter., Comment: 22 pages, 16 figures, proposed for acceptance in A&A

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 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

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 study young massive clusters (YMCs) in their embedded "infant" phase with $\sim0.\!^{\prime\prime}1$ ALMA, HST, and JWST observations targeting the central starburst ring in NGC 3351, a nearby Milky Way analog galaxy. Our new ALMA data reveal 18 bright and compact (sub-)millimeter continuum sources, of which 8 have counterparts in JWST images and only 6 have counterparts in HST images. Based on the ALMA continuum and molecular line data, as well as ancillary measurements for the HST and JWST counterparts, we identify 14 sources as infant star clusters with high stellar and/or gas masses (${\sim}10^5\;\mathrm{M_\odot}$), small radii (${\lesssim}\,5\;\mathrm{pc}$), large escape velocities ($6{-}10\;\mathrm{km/s}$), and short free-fall times ($0.5{-}1\;\mathrm{Myr}$). Their multiwavelength properties motivate us to divide them into four categories, likely corresponding to four evolutionary stages from starless clumps to exposed HII region-cluster complexes. Leveraging age estimates for HST-identified clusters in the same region, we infer an evolutionary timeline going from $\sim$1-2 Myr before cluster formation as starless clumps, to $\sim$4-6 Myr after as exposed HII region-cluster complexes. Finally, we show that the YMCs make up a substantial fraction of recent star formation across the ring, exhibit an non-uniform azimuthal distribution without a very coherent evolutionary trend along the ring, and are capable of driving large-scale gas outflows., Comment: 27 pages, 12 figures; ApJ accepted

Protoplanetary disks, which are the natural consequence of the gravitational collapse of the dense molecular cloud cores, host the formation of the known planetary systems in our universe. Substantial efforts have been dedicated to investigating the properties of these disks in the more mature Class II stage, either via numerical simulations of disk evolution from a limited range of initial conditions or observations of their dust continuum and line emission from specific molecular tracers. The results coming from these two standpoints have been used to draw comparisons. However, few studies have investigated the main limitations at work when measuring the embedded Class 0/I disk properties from observations, especially in a statistical fashion. In this study, we provide a first attempt to compare the accuracy of some critical disk parameters in Class 0/I systems, as derived on real ALMA observational data, with the corresponding physical parameters that can be directly defined by theoreticians and modellers in numerical simulations. The approach we follow here is to provide full post-processing of the numerical simulations and apply it to the synthetic observations the same techniques used by observers to derive the physical parameters. We performed 3D Monte Carlo radiative transfer and mock interferometric observations of the disk populations formed in a magnetohydrodynamic (MHD) simulation model of disk formation through the collapse of massive clumps with the tools \textsc{Radmc-3d} and \textsc{Casa}, respectively, to obtain their synthetic observations. With these observations, we re-employed the techniques commonly used in disk modelling from their continuum emissions to infer the properties that would most likely be obtained with real interferometers. We then demonstrated how these properties may vary with respect to the gas kinematics analyses and dust continuum modelling., Comment: Accepted for publication in Astronomy & Astrophysics, 32 pages, 28 figures

We present a numerical study of the evolution of power-law tails (PLTs) in the (column-)density distributions ($N$-PDF, $\rho$-PDF) in contracting star-forming clumps in primordial gas, without and with some initial rotational and/or turbulent support. In all considered runs multiple PLTs emerge shortly after the formation of the first protostar. The first PLT (PLT 1) in the $\rho$-PDF is a stable feature with slope $q_1\simeq -1.3$ which corresponds -- under the condition of preserved spherical symmetry -- to the outer envelope of the protostellar object with density profile $\rho\propto l^{-2}$ in the classical Larson-Penston collapse model, where $l$ is the radius. The second PLT (PLT 2) in the $\rho$-PDF is stable in the pure-infall runs but fluctuates significantly in the runs with initial support against gravity as dozens of protostars form and their mutual tidal forces change the density structure. Its mean slope, $\langle q_2\rangle\simeq -2$, corresponds to a density profile of $\rho\propto l^{-3/2}$ which describes a core in free fall in the classical Larson-Penston collapse model or an attractor solution at scales with dominating protostellar gravity. PLT 1 and PLT 2 in the $N$-PDFs are generally consistent with the observational data of Galactic low-mass star-forming regions from {\it Herschel} data. In the runs with initial support against gravity a third PLT (PLT~3) in the $\rho$-PDFs appears simultaneously with or after the emergence of PLT 2. It is very shallow, with mean slope of $\langle q_3\rangle\simeq -1$, and is associated with the formation of thin protostellar accretion disks., Comment: 14 pages, 12 figures; accepted for publication in MNRAS (main journal)