Showing 1–50 of 517 results for author: Klessen, R S

On the relation between magnetic field strength and gas density in the interstellar medium: A multiscale analysis

Authors: David J. Whitworth, Sundar Srinivasan, Ralph E. Pudritz, Mordecai M. Mac Low, Rowan J. Smith, Aina Palau, Kate Pattle, Gwendoline Eadie, Hector Robinson, Rachel Pillsworth, James Wadsley, Noe Brucy, Ugo Lebreuilly, Patrick Hennebelle, Philipp Girichidis, Fred A. Gent, Jessy Marin, Lylon Sánchez Valido, Vianey Camacho, Ralf S. Klessen, Enrique Vázquez-Semadeni

Abstract: 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… ▽ More 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. △ Less

Submitted 25 July, 2024; originally announced July 2024.

Comments: 23 figures, 28 pages, submitted to MNRAS, Comments welcome

The radiative torque spin-up efficiency of ballistic dust grain aggregates

Authors: Jonathan A. Jäger, Stefan Reissl, Ralf S. Klessen

Abstract: 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 t… ▽ More 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 $ω_{\mathrm{RAT}}$ caused by RATs of grain aggregates, and to resolve this problem we propose a new parametrization which predicts $ω_{\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 $ω_{\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. △ Less

Submitted 13 July, 2024; originally announced July 2024.

Comments: 15 pages, 12 figures

Gravitational waves from mergers of Population III binary black holes: roles played by two evolution channels

Authors: Boyuan Liu, Tilman Hartwig, Nina S. Sartorio, Irina Dvorkin, Guglielmo Costa, Filippo Santoliquido, Anastasia Fialkov, Ralf S. Klessen, Volker Bromm

Abstract: 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),… ▽ More 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 $Ω^{\rm peak}_{\rm GW}\sim 10^{-11}-8\times 10^{-11}$ around observer-frame frequencies $ν\sim 10-100\ \rm Hz$. The Pop III contribution can be a non-negligible ($\sim 2-32\%$) component in the total SGWB at $ν\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. △ Less

Submitted 25 June, 2024; originally announced June 2024.

Comments: 28 + 5 pages, 25 figures, submitted to MNRAS

NEATH III: a molecular line survey of a simulated star-forming cloud

Authors: F. D. Priestley, P. C. Clark, S. C. O. Glover, S. E. Ragan, O. Fehér, L. R. Prole, R. S. Klessen

Abstract: We present synthetic line observations of a simulated molecular cloud, utilising a self-consistent treatment of the dynamics and time-dependent chemical evolution. We investigate line emission from the three most common CO isotopologues ($^{12}$CO, $^{13}$CO, C$^{18}$O) and six supposed tracers of dense gas (NH$_3$, HCN, N$_2$H$^+$, HCO$^+$, CS, HNC). Our simulation produces a range of line intens… ▽ More We present synthetic line observations of a simulated molecular cloud, utilising a self-consistent treatment of the dynamics and time-dependent chemical evolution. We investigate line emission from the three most common CO isotopologues ($^{12}$CO, $^{13}$CO, C$^{18}$O) and six supposed tracers of dense gas (NH$_3$, HCN, N$_2$H$^+$, HCO$^+$, CS, HNC). Our simulation produces a range of line intensities consistent with that observed in real molecular clouds. The HCN-to-CO intensity ratio is relatively invariant with column density, making HCN (and chemically-similar species such as CS) a poor tracer of high-density material in the cloud. The ratio of N$_2$H$^+$ to HCN or CO, on the other hand, is highly selective of regions with densities above $10^{22} \, {\rm cm^{-2}}$, and the N$_2$H$^+$ line is a very good tracer of the dynamics of high volume density ($>10^4 \, {\rm cm^{-3}}$) material. Focusing on cores formed within the simulated cloud, we find good agreement with the line intensities of an observational sample of prestellar cores, including reproducing observed CS line intensities with an undepleted elemental abundance of sulphur. However, agreement between cores formed in the simulation, and models of isolated cores which have otherwise-comparable properties, is poor. The formation from and interaction with the large-scale environment has a significant impact on the line emission properties of the cores, making isolated models unsuitable for interpreting observational data. △ Less

Submitted 10 June, 2024; originally announced June 2024.

Comments: 14 pages, 14 figures. MNRAS accepted

Magnetized compressible turbulence with a fluctuation dynamo and Reynolds numbers over a million

Authors: James R. Beattie, Christoph Federrath, Ralf S. Klessen, Salvatore Cielo, Amitava Bhattacharjee

Abstract: 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 do… ▽ More 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. △ Less

Submitted 29 May, 2024; v1 submitted 26 May, 2024; originally announced May 2024.

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

The Fraction of Dust Mass in the Form of PAHs on 10-50 pc Scales in Nearby Galaxies

Authors: Jessica Sutter, Karin Sandstrom, Jérémy Chastenet, Adam K. Leroy, Eric W. Koch, Thomas G. Williams, Ryan Chown, Francesco Belfiore, Frank Bigiel, Médéric Boquien, Yixian Cao, Mélanie Chevance, Daniel A. Dale, Oleg V. Egorov, Simon C. O. Glover, Brent Groves, Ralf S. Klessen, Kathryn Kreckel, Kirsten L. Larson, Elias K. Oakes, Debosmita Pathak, Lise Ramambason, Erik Rosolowsky, Elizabeth J. Watkins

Abstract: 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… ▽ More 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%. △ Less

Submitted 23 May, 2024; originally announced May 2024.

Comments: Accepted at ApJ, 39 pages, 25 figures

Massive Star Cluster Formation II. Runaway Stars as Fossils of Sub-Cluster Mergers

Authors: Brooke Polak, Mordecai-Mark Mac Low, Ralf S. Klessen, Simon Portegies Zwart, Eric P. Andersson, Sabrina M. Appel, Claude Cournoyer-Cloutier, Simon C. O. Glover, Stephen L. W. McMillan

Abstract: 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 star… ▽ More 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. △ Less

Submitted 20 May, 2024; originally announced May 2024.

Comments: 11 pages, 10 figures. Submitted to A&A. Comments welcome!

Self-consistent modelling of the Milky Way structure using live potentials

Authors: Eva Durán-Camacho, Ana Duarte-Cabral, Alex R. Pettitt, Robin G. Treß, Paul C. Clark, Ralf S. Klessen, Kamran R. J. Bogue, Rowan J. Smith, Mattia C. Sormani

Abstract: 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 observation… ▽ More 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. △ Less

Submitted 11 June, 2024; v1 submitted 15 May, 2024; originally announced May 2024.

Comments: Accepted for publication in MNRAS. 24 pages, 23 figures, 3 tables

Discovery of $\sim$2200 new supernova remnants in 19 nearby star-forming galaxies with MUSE spectroscopy

Authors: Jing Li, K. Kreckel, S. Sarbadhicary, Oleg V. Egorov, B. Groves, K. S. Long, Enrico Congiu, Francesco Belfiore, Simon C. O. Glover, Ashley . T Barnes, Frank Bigiel, Guillermo A. Blanc, Kathryn Grasha, Ralf S. Klessen, Adam Leroy, Laura A. Lopez, J. Eduardo Méndez-Delgado, Justus Neumann, Eva Schinnerer, Thomas G. Williams, PHANGS collaborators

Abstract: 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 ~… ▽ More 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$α$; (2) line ratio maps of [OI]/H$α$; (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$α$ or [OI]/H$α$ in excess of what is expected at fixed H$α$ 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$α$ line ratio is responsible for selecting the most objects (1,368; 61%), (abridged). △ Less

Submitted 14 May, 2024; originally announced May 2024.

Comments: 35 pages, 24 figures,6 tables, submitted to A&A

Do spiral arms enhance star formation efficiency?

Authors: Miguel Querejeta, Adam K. Leroy, Sharon E. Meidt, Eva Schinnerer, Francesco Belfiore, Eric Emsellem, Ralf S. Klessen, Jiayi Sun, Mattia Sormani, Ivana Bešlic, Yixian Cao, Mélanie Chevance, Dario Colombo, Daniel A. Dale, Santiago García-Burillo, Simon C. O. Glover, Kathryn Grasha, Brent Groves, Eric. W. Koch, Lukas Neumann, Hsi-An Pan, Ismael Pessa, Jérôme Pety, Francesca Pinna, Lise Ramambason , et al. (10 additional authors not shown)

Abstract: 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 fr… ▽ More 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) △ Less

Submitted 8 May, 2024; originally announced May 2024.

Comments: 26 pages, 16 figures. Accepted for publication in A&A

JWST Observations of Starbursts: Cold Clouds and Plumes Launching in the M82 Outflow

Authors: Deanne B. Fisher, Alberto D. Bolatto, John Chisholm, Drummond Fielding, Rebecca C. Levy, Elizabeth Tarantino, Martha L. Boyer, Serena A. Cronin, Laura A. Lopez, J. D. Smith, Danielle A. Berg, Sebastian Lopez, Sylvain Veilleux, Paul P. van der Werf, Torsten Böker, Leindert A. Boogaard, Laura Lenkić, Simon C. O. Glover, Vicente Villanueva, Divakara Mayya, Thomas S. -Y. Lai, Daniel A. Dale, Kimberly L. Emig, Fabian Walter, Monica Relaño , et al. (6 additional authors not shown)

Abstract: In this paper we study the filamentary substructure of 3.3 $μ$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 (… ▽ More In this paper we study the filamentary substructure of 3.3 $μ$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$α$ emission (from HST/NICMOS data), at the scale of both plumes and clouds, though the ratio of PAH-to-Pa$α$ 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. △ Less

Submitted 6 May, 2024; originally announced May 2024.

Comments: Submitted to MNRAS

The SDSS-V Local Volume Mapper (LVM): Scientific Motivation and Project Overview

Authors: Niv Drory, Guillermo A. Blanc, Kathryn Kreckel, Sebastian F. Sanchez, Alfredo Mejia-Narvaez, Evelyn J. Johnston, Amy M. Jones, Eric W. Pellegrini, Nicholas P. Konidaris, Tom Herbst, Jose Sanchez-Gallego, Juna A. Kollmeier, Florence de Almeida, Jorge K. Barrera-Ballesteros, Dmitry Bizyaev, Joel R. Brownstein, Mar Canal i Saguer, Brian Cherinka, Maria-Rosa L. Cioni, Enrico Congiu, Maren Cosens, Bruno Dias, John Donor, Oleg Egorov, Evgeniia Egorova , et al. (26 additional authors not shown)

Abstract: 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 resolution… ▽ More 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. △ Less

Submitted 2 May, 2024; originally announced May 2024.

Comments: 29 pages, 12 figures, accepted for publication in The Astronomical Journal

Inefficient star formation in high Mach number environments II. Numerical simulations and comparison with analytical models

Authors: Noé Brucy, Patrick Hennebelle, Tine Colman, Ralf S. Klessen, Corentin Le Yhuelic

Abstract: 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… ▽ More 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. △ Less

Submitted 25 July, 2024; v1 submitted 26 April, 2024; originally announced April 2024.

Comments: 13 pages, 14 figures, accepted for publication in A&A

H-alpha emission and HII regions at the locations of recent supernovae in nearby galaxies

Authors: Ness Mayker Chen, Adam K. Leroy, Sumit K. Sarbadhicary, Laura A. Lopez, Todd A. Thompson, Ashley T. Barnes, Eric Emsellem, Brent Groves, Rupali Chandar, Mélanie Chevance, Ryan Chown, Daniel A. Dale, Oleg V. Egorov, Simon C. O. Glover, Kathryn Grasha, Ralf S. Klessen, Kathryn Kreckel, Jing Li, J. Eduardo Méndez-Delgado, Eric J. Murphy, Debosmita Pathak, Eva Schinnerer, David A. Thilker, Leonardo Úbeda, Thomas G. Williams

Abstract: 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… ▽ More 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. △ Less

Submitted 16 April, 2024; originally announced April 2024.

Comments: Accepted for publication in ApJ; 33 pages, 13 figures, 3 tables in two-column AASTEX63 format

Simulating nearby disc galaxies on the main star formation sequence I. Bar formation and the building of the central gas reservoir

Authors: Pierrick Verwilghen, Eric Emsellem, Florent Renaud, Milena Valentini, Jiayi Sun, Sarah Jeffreson, Ralf S. Klessen, Mattia C. Sormani, Ashley. T. Barnes, Klaus Dolag, Kathryn Grasha, Fu-Heng Liang, Sharon Meidt, Justus Neumann, Miguel Querejeta, Eva Schinnerer, Thomas G. Williams

Abstract: 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 bui… ▽ More 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. △ Less

Submitted 15 April, 2024; originally announced April 2024.

Comments: 22 pages, 17 figures

Journal ref: A&A 687, A53 (2024)

A broad linewidth, compact, millimeter-bright molecular emission line source near the Galactic Center

Authors: Adam Ginsburg, John Bally, Ashley T. Barnes, Cara Battersby, Nazar Budaiev, Natalie O. Butterfield, Paola Caselli, Laura Colzi, Katarzyna M. Dutkowska, Pablo García, Savannah Gramze, Jonathan D. Henshaw, Yue Hu, Desmond Jeff, Izaskun Jiménez-Serra, Jens Kauffmann, Ralf S. Klessen, Emily M. Levesque, Steven N. Longmore, Xing Lu, Elisabeth A. C. Mills, Mark R. Morris, Francisco Nogueras-Lara, Tomoharu Oka, Jaime E. Pineda , et al. (15 additional authors not shown)

Abstract: 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 $α\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 w… ▽ More 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 $α\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. △ Less

Submitted 1 May, 2024; v1 submitted 11 April, 2024; originally announced April 2024.

Comments: Accepted to ApJL

Testing kinematic distances under a realistic Galactic potential

Authors: Glen H. Hunter, Mattia C. Sormani, Jan P. Beckmann, Eugene Vasiliev, Simon C. O. Glover, Ralf S. Klessen, Juan D. Soler, Noé Brucy, Philipp Girichidis, Junia Göller, Loke Ohlin, Robin Tress, Sergio Molinari, Ortwin Gerhard, Milena Benedettini, Rowan Smith, Patrick Hennebelle, Leonardo Testi

Abstract: 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 t… ▽ More 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. △ Less

Submitted 26 March, 2024; originally announced March 2024.

Comments: 23 pages, 22 figures, 6 tables, submitted to A&A

PHANGS-HST: Globular Cluster Systems in 17 Nearby Spiral Galaxies

Authors: Matthew Floyd, Rupali Chandar, Bradley C. Whitmore, David A. Thilker, Janice C. Lee, Rachel E. Pauline, Zion L. Thomas, William J. Berschback, Kiana F. Henny, Daniel A. Dale, Ralf S. Klessen, Eva Schinnerer, Kathryn Grasha, Mederic Boquien, Kirsten L. Larson, Sinan Deger, Ashley T. Barnes, Adam K. Leroy, Erik Rosolowsky, Thomas G. Williams, Leonardo Ubeda

Abstract: 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 morpho… ▽ More 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. △ Less

Submitted 20 March, 2024; originally announced March 2024.

Comments: 21 pages, 12 figures

Journal ref: The Astronomical Journal, Volume 167, Issue 3, id.95, 16 pp. March 2024

The properties and kinematics of HCN emission across the closest starburst galaxy NGC 253 observed with ALMA

Authors: Ivana Beslic, Ashley T. Barnes, Frank Bigiel, Maria Jesus Jimenez-Donaire, Antonio Usero, Jonathan D. Henshaw, Christopher Faesi, Adam K. Leroy, Erik Rosolowsky, Jakob S. den Brok, Melanie Chevance, Cosima Eibensteiner, Kathryn Grasha, Ralf S. Klessen, J. M. Diedrerik Kruijssen, Daizhong Liu, Sharon Meidt, Justus Neumann, Lukas Neumann, Hsi-An Pan, Johannes Puschnig, Miguel Querejeta, Eva Schinnerer, Thomas G. Williams

Abstract: 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 - an… ▽ More 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. △ Less

Submitted 20 March, 2024; originally announced March 2024.

Comments: Accepted for publication to Astronomy and Astrophysics

Magnetic field morphology and evolution in the Central Molecular Zone and its effect on gas dynamics

Authors: R. G. Tress, M. C. Sormani, P. Girichidis, S. C. O. Glover, R. S. Klessen, R. J. Smith, E. Sobacchi, L. Armillotta, A. T. Barnes, C. Battersby, K. R. J. Bogue, N. Brucy, L. Colzi, C. Federrath, P. García, A. Ginsburg, J. Göller, H P. Hatchfield, C. Henkel, P. Hennebelle, J. D. Henshaw, M. Hirschmann, Y. Hu, J. Kauffmann, J. M. D. Kruijssen , et al. (12 additional authors not shown)

Abstract: The interstellar medium in the Milky Way's Central Molecular Zone (CMZ) is known to be strongly magnetised, but its large-scale morphology and impact on the gas dynamics are not well understood. We explore the impact and properties of magnetic fields in the CMZ using three-dimensional non-self gravitating magnetohydrodynamical simulations of gas flow in an external Milky Way barred potential. We f… ▽ More The interstellar medium in the Milky Way's Central Molecular Zone (CMZ) is known to be strongly magnetised, but its large-scale morphology and impact on the gas dynamics are not well understood. We explore the impact and properties of magnetic fields in the CMZ using three-dimensional non-self gravitating magnetohydrodynamical simulations of gas flow in an external Milky Way barred potential. We find that: (1) The magnetic field is conveniently decomposed into a regular time-averaged component and an irregular turbulent component. The regular component aligns well with the velocity vectors of the gas everywhere, including within the bar lanes. (2) The field geometry transitions from parallel to the Galactic plane near $z=0$ to poloidal away from the plane. (3) The magneto-rotational instability (MRI) causes an in-plane inflow of matter from the CMZ gas ring towards the central few parsecs of $0.01-0.1$ M$_\odot$ yr$^{-1}$ that is absent in the unmagnetised simulations. However, the magnetic fields have no significant effect on the larger-scale bar-driven inflow that brings the gas from the Galactic disc into the CMZ. (4) A combination of bar inflow and MRI-driven turbulence can sustain a turbulent vertical velocity dispersion of $σ_z \simeq 5$ km s$^{-1}$ on scales of $20$ pc in the CMZ ring. The MRI alone sustains a velocity dispersion of $σ_z \simeq 3$ km s$^{-1}$. Both these numbers are lower than the observed velocity dispersion of gas in the CMZ, suggesting that other processes such as stellar feedback are necessary to explain the observations. (5) Dynamo action driven by differential rotation and the MRI amplifies the magnetic fields in the CMZ ring until they saturate at a value that scales with the average local density as $B \simeq 102 (n/10^3 cm^{-3})^{0.33}$ $μ$G. Finally, we discuss the implications of our results within the observational context in the CMZ. △ Less

Submitted 19 March, 2024; originally announced March 2024.

Comments: submitted to A&A, comments welcome

PHANGS-HST catalogs for $\sim$100,000 star clusters and compact associations in 38 galaxies: I. Observed properties

Authors: Daniel Maschmann, Janice C. Lee, David A. Thilker, Bradley C. Whitmore, Sinan Deger, Mederic Boquien, Rupali Chandar, Daniel A. Dale, Aida Wofford, Stephen Hannon, Kirsten L. Larson, Adam K. Leroy, Eva Schinnerer, Erik W. Rosolowsky, Leonardo Ubeda, Ashley Barnes, Eric Emsellem, Kathryn Grasha, Brent Groves, Hwihyun Kim, Ralf S. Klessen, Kathryn Kreckel, Rebecca C. Levy, Francesca Pinna, Jimena Rodriguez , et al. (2 additional authors not shown)

Abstract: 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 ha… ▽ More 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. △ Less

Submitted 7 March, 2024; originally announced March 2024.

Comments: 48 pages, 26 figures, article in review at APJ

Cloud properties across spatial scales in simulations of the interstellar medium

Authors: Tine Colman, Noé Brucy, Philipp Girichidis, Simon C. O Glover, Milena Benedettini, Juan D. Soler, Robin G. Tress, Alessio Traficante, Patrick Hennebelle, Ralf S. Klessen, Sergio Molinari, Marc-Antoine Miville-Deschênes

Abstract: 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… ▽ More 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 $σ\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 $α_\mathrm{vir}\propto M^{-0.4}$, which is slightly flatter compared to observed scaling, but in agreement given the large scatter. △ Less

Submitted 1 March, 2024; originally announced March 2024.

Comments: 22 pages, 16 figures, proposed for acceptance in A&A

Journal ref: A&A 686, A155 (2024)

PHANGS-ML: dissecting multiphase gas and dust in nearby galaxies using machine learning

Authors: Dalya Baron, Karin M. Sandstrom, Erik Rosolowsky, Oleg V. Egorov, Ralf S. Klessen, Adam K. Leroy, Médéric Boquien, Eva Schinnerer, Francesco Belfiore, Brent Groves, Jérémy Chastenet, Daniel A. Dale, Guillermo A. Blanc, José E. Méndez-Delgado, Eric W. Koch, Kathryn Grasha, Mélanie Chevance, David A. Thilker, Dario Colombo, Thomas G. Williams, Debosmita Pathak, Jessica Sutter, Toby Brown, John F. Wu, J. E. G. Peek , et al. (3 additional authors not shown)

Abstract: 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 algorit… ▽ More 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. △ Less

Submitted 6 February, 2024; originally announced February 2024.

Comments: Main results in figures 6 and 12. Submitted to ApJ, and comments are welcome!

JWST Observations of Starbursts: Polycyclic Aromatic Hydrocarbon Emission at the Base of the M 82 Galactic Wind

Authors: Alberto D. Bolatto, Rebecca C. Levy, Elizabeth Tarantino, Martha L. Boyer, Deanne B. Fisher, Adam K. Leroy, Serena A. Cronin, Ralf S. Klessen, J. D. Smith, Dannielle A. Berg, Torsten Boeker, Leindert A. Boogaard, Eve C. Ostriker, Todd A. Thompson, Juergen Ott, Laura Lenkic, Laura A. Lopez, Daniel A. Dale, Sylvain Veilleux, Paul P. van der Werf, Simon C. O. Glover, Karin M. Sandstrom, Evan D. Skillman, John Chisholm, Vicente Villanueva , et al. (15 additional authors not shown)

Abstract: 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 promine… ▽ More 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. △ Less

Submitted 21 April, 2024; v1 submitted 29 January, 2024; originally announced January 2024.

Comments: Submitted to The Astrophysical Journal

PHANGS-JWST: Data Processing Pipeline and First Full Public Data Release

Authors: Thomas G. Williams, Janice C. Lee, Kirsten L. Larson, Adam K. Leroy, Karin Sandstrom, Eva Schinnerer, David A. Thilker, Francesco Belfiore, Oleg V. Egorov, Erik Rosolowsky, Jessica Sutter, Joseph DePasquale, Alyssa Pagan, Travis A. Berger, Gagandeep S. Anand, Ashley T. Barnes, Frank Bigiel, Médéric Boquien, Yixian Cao, Jérémy Chastenet, Mélanie Chevance, Ryan Chown, Daniel A. Dale, Sinan Deger, Cosima Eibensteiner , et al. (33 additional authors not shown)

Abstract: 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… ▽ More 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 $μ$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. △ Less

Submitted 9 May, 2024; v1 submitted 26 January, 2024; originally announced January 2024.

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

Hidden Gems on a Ring: Infant Massive Clusters and Their Formation Timeline Unveiled by ALMA, HST, and JWST in NGC 3351

Authors: Jiayi Sun, Hao He, Kyle Batschkun, Rebecca C. Levy, Kimberly Emig, M. Jimena Rodriguez, Hamid Hassani, Adam K. Leroy, Eva Schinnerer, Eve C. Ostriker, Christine D. Wilson, Alberto D. Bolatto, Elisabeth A. C. Mills, Erik Rosolowsky, Janice C. Lee, Daniel A. Dale, Kirsten L. Larson, David A. Thilker, Leonardo Ubeda, Bradley C. Whitmore, Thomas G. Williams, Ashley. T. Barnes, Frank Bigiel, Melanie Chevance, Simon C. O. Glover , et al. (16 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 10 April, 2024; v1 submitted 25 January, 2024; originally announced January 2024.

Comments: 27 pages, 12 figures; ApJ accepted

Accuracy of ALMA estimates of young disk radii and masses. Predicted observations from numerical simulations

Authors: Ngo-Duy Tung, Leonardo Testi, Ugo Lebreuilly, Patrick Hennebelle, Anaëlle Maury, Ralf S. Klessen, Luca Cacciapuoti, Matthias González, Giovanni Rosotti, Sergio Molinari

Abstract: 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 condit… ▽ More 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. △ Less

Submitted 29 January, 2024; v1 submitted 22 January, 2024; originally announced January 2024.

Comments: Accepted for publication in Astronomy & Astrophysics, 32 pages, 28 figures

Multiple power-law tails in the density and column-density distribution in contracting star-forming clumps

Authors: Todor V. Veltchev, Philipp Girichidis, Lyubov Marinkova, Sava Donkov, Orlin Stanchev, Ralf S. Klessen

Abstract: We present a numerical study of the evolution of power-law tails (PLTs) in the (column-)density distributions ($N$-PDF, $ρ$-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 $ρ$-PDF is a stable featu… ▽ More We present a numerical study of the evolution of power-law tails (PLTs) in the (column-)density distributions ($N$-PDF, $ρ$-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 $ρ$-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 $ρ\propto l^{-2}$ in the classical Larson-Penston collapse model, where $l$ is the radius. The second PLT (PLT 2) in the $ρ$-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 $ρ\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 $ρ$-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. △ Less

Submitted 4 January, 2024; originally announced January 2024.

Comments: 14 pages, 12 figures; accepted for publication in MNRAS (main journal)

Heavy Black Hole Seed Formation in High-z Atomic Cooling Halos

Authors: Lewis R. Prole, John A. Regan, Simon C. O. Glover, Ralf S. Klessen, Felix D. Priestley, Paul C. Clark

Abstract: Halos with masses in excess of the atomic limit are believed to be ideal environments in which to form heavy black hole seeds with masses above 10^3 Msun. In cases where the H_2 fraction is suppressed this is expected to lead to reduced fragmentation of the gas and the generation of a top heavy initial mass function. In extreme cases this can result in the formation of massive black hole seeds. Re… ▽ More Halos with masses in excess of the atomic limit are believed to be ideal environments in which to form heavy black hole seeds with masses above 10^3 Msun. In cases where the H_2 fraction is suppressed this is expected to lead to reduced fragmentation of the gas and the generation of a top heavy initial mass function. In extreme cases this can result in the formation of massive black hole seeds. Resolving the initial fragmentation scale and the resulting protostellar masses has, until now, not been robustly tested. Cosmological simulations were performed with the moving mesh code Arepo using a primordial chemistry network until z = 11. Three haloes with masses in excess of the atomic cooling mass were then selected for detailed examination via zoom-ins. The highest resolution simulations resolve densities up to 10^-6 g cm^-3 (10^18 cm^-3) and capture a further 100 yr of fragmentation behaviour at the center of the halo. Our simulations show intense fragmentation in the central region of the halos, leading to a large number of near-solar mass protostars. Despite the increased fragmentation the halos produce a protostellar mass spectrum that peaks at higher masses relative to standard Population III star forming halos. The most massive protostars have accretion rates of 10^-3-10^-1 Msun yr^-1 after the first 100 years of evolution, while the total mass of the central region grows at 1 Msun yr^-1. Lower resolution zoom-ins show that the total mass of the system continues to accrete at 1 Msun yr^-1 for at least 10^4 yr, although how this mass is distributed amongst the rapidly growing number of protostars is unclear. However, assuming that a fraction of stars can continue to accrete rapidly the formation of a sub-population of stars with masses in excess of 10^3 Msun is likely in these halos. △ Less

Submitted 11 December, 2023; originally announced December 2023.

Comments: Submitted to A&A, comments welcome

Massive Star Cluster Formation I. High Star Formation Efficiency While Resolving Feedback of Individual Stars

Authors: Brooke Polak, Mordecai-Mark Mac Low, Ralf S. Klessen, Jia Wei Teh, Claude Cournoyer-Cloutier, Eric P. Andersson, Sabrina M. Appel, Aaron Tran, Sean C. Lewis, Maite J. C. Wilhelm, Simon Portegies Zwart, Simon C. O. Glover, Long Wang, Stephen L. W. McMillan

Abstract: The mode of star formation that results in the formation of globular clusters and young massive clusters is difficult to constrain through observations. We present models of massive star cluster formation using the Torch framework, which uses AMUSE to couple distinct multi-physics codes that handle star formation, stellar evolution and dynamics, radiative transfer, and magnetohydrodynamics. We upg… ▽ More The mode of star formation that results in the formation of globular clusters and young massive clusters is difficult to constrain through observations. We present models of massive star cluster formation using the Torch framework, which uses AMUSE to couple distinct multi-physics codes that handle star formation, stellar evolution and dynamics, radiative transfer, and magnetohydrodynamics. We upgrade Torch by implementing the N-body code PeTar, thereby enabling Torch to handle massive clusters forming from $10^6\rm\, M_\odot$ clouds with $\ge10^5$ individual stars. We present results from Torch simulations of star clusters forming from $10^4, 10^5$, and $10^6\rm M_\odot$ turbulent, spherical gas clouds (named M4, M5, M6) of radius $R=11.7$ pc. We find that star formation is highly efficient and becomes more so at higher cloud mass and surface density. For M4, M5, and M6 with initial surface densities $2.325\times 10^{1,2,3}\rm\, M_\odot\, pc^{-2}$, after a free-fall time of $t_{ff}=6.7,2.1,0.67$ Myr, we find that $\sim30\%$, 40%, and 60% of the cloud mass has formed into stars, respectively. The final integrated star formation efficiency is $32\%,\, 65\%$, and 85\% for M4, M5, and M6. Observations of nearby clusters similar to M4 have similar integrated star formation efficiencies of $\leq30\%$. The M5 and M6 models represent a different regime of cluster formation that is more appropriate for the conditions in starburst galaxies and gas-rich galaxies at high redshift, and that leads to a significantly higher efficiency of star formation. We argue that young massive clusters build up through short efficient bursts of star formation in regions that are sufficiently dense ($\ge 10^2 \rm\,M_\odot\,pc^{-2}$) and massive ($\ge10^5\rm\, M_\odot$). In such environments, the dynamical time of the cloud becomes short enough that stellar feedback cannot act quickly enough to slow star formation. △ Less

Submitted 11 December, 2023; originally announced December 2023.

Comments: Submitted to A&A

The PHANGS-AstroSat Atlas of Nearby Star Forming Galaxies

Authors: Hamid Hassani, Erik Rosolowsky, Eric W. Koch, Joseph Postma, Joseph Nofech, Harrisen Corbould, David Thilker, Adam K. Leroy, Eva Schinnerer, Francesco Belfiore, Frank Bigiel, Mederic Boquien, Melanie Chevance, Daniel A. Dale, Oleg V. Egorov, Eric Emsellem, Simon C. O. Glover, Kathryn Grasha, Brent Groves, Kiana Henny, Jaeyeon Kim, Ralf S. Klessen, Kathryn Kreckel, J. M. Diederik Kruijssen, Janice C. Lee , et al. (7 additional authors not shown)

Abstract: 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 a… ▽ More 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$α$ in different galactic environments and kpc-sized apertures. We report that $\sim 65 \%$ varation of the $\log_{10}$(FUV/H$α$) 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. △ Less

Submitted 10 December, 2023; originally announced December 2023.

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

A Two-Component Probability Distribution Function Describes the mid-IR Emission from the Disks of Star-Forming Galaxies

Authors: Debosmita Pathak, Adam K. Leroy, Todd A. Thompson, Laura A. Lopez, Francesco Belfiore, Mederic Boquien, Daniel A. Dale, Simon C. O. Glover, Ralf S. Klessen, Eric W. Koch, Erik Rosolowsky, Karin M. Sandstrom, Eva Schinnerer, Rowan Smith, Jiayi Sun, Jessica Sutter, Thomas G. Williams, Frank Bigiel, Yixian Cao, Jeremy Chastenet, Melanie Chevance, Ryan Chown, Eric Emsellem, Christopher M. Faesi, Kirsten L. Larson , et al. (6 additional authors not shown)

Abstract: High-resolution JWST-MIRI images of nearby spiral galaxies reveal emission with complex substructures that trace dust heated both by massive young stars and the diffuse interstellar radiation field. We present high angular (0."85) and physical resolution (20-80 pc) measurements of the probability distribution function (PDF) of mid-infrared (mid-IR) emission (7.7-21 $μ$m) from 19 nearby star-formin… ▽ More High-resolution JWST-MIRI images of nearby spiral galaxies reveal emission with complex substructures that trace dust heated both by massive young stars and the diffuse interstellar radiation field. We present high angular (0."85) and physical resolution (20-80 pc) measurements of the probability distribution function (PDF) of mid-infrared (mid-IR) emission (7.7-21 $μ$m) from 19 nearby star-forming galaxies from the PHANGS-JWST Cycle-1 Treasury. The PDFs of mid-IR emission from the disks of all 19 galaxies consistently show two distinct components: an approximately log-normal distribution at lower intensities and a high-intensity power-law component. These two components only emerge once individual star-forming regions are resolved. Comparing with locations of HII regions identified from VLT/MUSE H$α$-mapping, we infer that the power-law component arises from star-forming regions and thus primarily traces dust heated by young stars. In the continuum-dominated 21 $μ$m band, the power-law is more prominent and contains roughly half of the total flux. At 7.7-11.3 $μ$m, the power-law is suppressed by the destruction of small grains (including PAHs) close to HII regions while the log-normal component tracing the dust column in diffuse regions appears more prominent. The width and shape of the log-normal diffuse emission PDFs in galactic disks remain consistent across our sample, implying a log-normal gas column density $N$(H)$\approx10^{21}$cm$^{-2}$ shaped by supersonic turbulence with typical (isothermal) turbulent Mach numbers $\approx5-15$. Finally, we describe how the PDFs of galactic disks are assembled from dusty HII regions and diffuse gas, and discuss how the measured PDF parameters correlate with global properties such as star-formation rate and gas surface density. △ Less

Submitted 29 November, 2023; originally announced November 2023.

Comments: 30 pages without appendix, 17 figures, (with appendix images of full sample: 56 pages, 39 figures), accepted in AJ

Protostellar chimney flues: are jets and outflows lifting submillimetre dust grains from discs into envelopes?

Authors: L. Cacciapuoti, L. Testi, L. Podio, C. Codella, A. J. Maury, M. De Simone, P. Hennebelle, U. Lebreuilly, R. S. Klessen, S. Molinari

Abstract: Low dust opacity spectral indices ($β< 1$) measured in the inner envelopes of class 0/I young stellar objects (age $\sim 10^{4-5}$ yr) have been interpreted as the presence of (sub-)millimetre dust grains in these environments. The density conditions and the lifetimes of collapsing envelopes have proven unfavorable for the growth of solids up to millimetre sizes. As an alternative, magneto-hydrody… ▽ More Low dust opacity spectral indices ($β< 1$) measured in the inner envelopes of class 0/I young stellar objects (age $\sim 10^{4-5}$ yr) have been interpreted as the presence of (sub-)millimetre dust grains in these environments. The density conditions and the lifetimes of collapsing envelopes have proven unfavorable for the growth of solids up to millimetre sizes. As an alternative, magneto-hydrodynamical simulations suggest that protostellar jets and outflows might lift grains from circumstellar discs and diffuse them in the envelope. We reframe available data for the CALYPSO sample of Class 0/I sources and show tentative evidence for an anti-correlation between the value of $β_{1-3mm}$ measured in the inner envelope and the mass loss rate of their jets and outflows, supporting a connection between the two. We discuss the implications that dust transport from the disc to the inner envelope might have for several aspects of planet formation. Finally, we urge for more accurate measurements of both correlated quantities and extension of this work to larger samples, necessary to further test the transport scenario. △ Less

Submitted 27 November, 2023; originally announced November 2023.

Comments: Accepted for publication on ApJ, 12 pages, 3 figures, 3 tables

A dusty streamer infalling onto the disk of a class I protostar. ALMA dual-band constraints on grain properties and mass infall rate

Authors: L. Cacciapuoti, E. Macias, A. Gupta, L. Testi, A. Miotello, C. Espaillat, M. Kuffmeier, S. van Terwisga, J. Tobin, S. Grant, C. F. Manara, D. Segura-Cox, J. Wendeborn, R. S. Klessen, A. J. Maury, U. Lebreuilly, P. Hennebelle, S. Molinari

Abstract: Observations of interstellar material infalling onto star- and planet-forming systems have become increasingly common thanks to recent advancements in radio interferometry. These structures replenish disks with fresh material, have the potential to significantly alter their dynamics, trigger the formation of substructures, induce shocks, and modify their physical and chemical properties. In this s… ▽ More Observations of interstellar material infalling onto star- and planet-forming systems have become increasingly common thanks to recent advancements in radio interferometry. These structures replenish disks with fresh material, have the potential to significantly alter their dynamics, trigger the formation of substructures, induce shocks, and modify their physical and chemical properties. In this study, we combine new ALMA band 3 and archival band 6 observations to characterize the dust content and mass infall rate of a 4,000 au arc-like structure infalling onto M512, a class I young stellar object located in the Lynds 1641 region of the Orion A molecular cloud. We measure for the first time spectral index maps and derive a dust opacity index profile along a streamer, constraining grain properties and its dust mass. We measure a spectral index $α\sim$ 3.2 across the entire structure, and a dust opacity index $β\sim$ 1.6. Given grain properties consistent with the measured $β$, the structure can host up to 245 M$_{\oplus}$ of dust, being comparable or even exceeding the mass of the inner, unresolved 600 au, which contains the protoplanetary disk of M512. Such a massive streamer can strongly affect the evolution of the star- and planet-forming inner system. Assuming typical ISM dust-to-gas ratio of 1%, free-fall timescales (50 kyr) imply total mass infall rates up to 1.5 $\cdot$ 10$^{-6}$ M$_{\odot}$/yr. M512 has been classified as an outbursting source with multi-epoch photometry, thus representing an interesting case study to explore the possible connection between infalling streamers and accretion outbursts. △ Less

Submitted 22 November, 2023; originally announced November 2023.

Comments: Accepted for publication in A&A, 13 pages, 5 figures, 3 tables

Not All Subhaloes Are Created Equal: Modelling the Diversity of Subhalo Density Profiles in TNG50

Authors: Felix M. Heinze, Giulia Despali, Ralf S. Klessen

Abstract: In this work, we analyse the density profiles of subhaloes with masses $M_\mathrm{sh} \geq 1.4 \times 10^8$ M$_\odot$ in the TNG50 simulation, with the aim of including baryonic effects. We evaluate the performance of frequently used models, such as the standard NFW, the Einasto, and a smoothly truncated version of the NFW profile. We find that these models do not perform well for the majority of… ▽ More In this work, we analyse the density profiles of subhaloes with masses $M_\mathrm{sh} \geq 1.4 \times 10^8$ M$_\odot$ in the TNG50 simulation, with the aim of including baryonic effects. We evaluate the performance of frequently used models, such as the standard NFW, the Einasto, and a smoothly truncated version of the NFW profile. We find that these models do not perform well for the majority of subhaloes, with the NFW profile giving the worst fit in most cases. This is primarily due to mismatches in the inner and outer logarithmic slopes, which are significantly steeper for a large number of subhaloes in the presence of baryons. To address this issue, we propose new three-parameter models and show that they significantly improve the goodness of fit independently of the subhalo's specific properties. Our best-performing model is a modified version of the NFW profile with an inner log-slope of -2 and a variable truncation that is sharper and steeper than the slope transition in the standard NFW profile. Additionally, we investigate how both the parameter values of the best density profile model and the average density profiles vary with subhalo mass, $V_\mathrm{max}$, distance from the host halo centre, baryon content and infall time, and we also present explicit scaling relations for the mean parameters of the individual profiles. The newly proposed fit and the scaling relations are useful to predict the properties of realistic subhaloes in the mass range $10^8$ M$_\odot$ $\leq M_\mathrm{sh} \leq$ $10^{13}$ M$_\odot$ that can be influenced by the presence of baryons. △ Less

Submitted 22 November, 2023; originally announced November 2023.

Comments: 15+5 pages, 17+13 figures, submitted to MNRAS, comments welcome

Journal ref: MNRAS, Volume 527, Issue 4, February 2024, 11996-12015

Density profile of a self-gravitating polytropic turbulent fluid in a rotating disk near to the cloud core

Authors: S. Donkov, I. Zh. Stefanov, T. V. Veltchev, R. S. Klessen

Abstract: We obtain two equations (following from two different approaches) for the density profile in a self-gravitating polytropic cylindrically symmetric and rotating turbulent gas disk. The adopted physical picture is appropriate to describe the conditions near to the cloud core where the equation of state of the gas changes from isothermal (in the outer cloud layers) to one of "hard polytrope", and the… ▽ More We obtain two equations (following from two different approaches) for the density profile in a self-gravitating polytropic cylindrically symmetric and rotating turbulent gas disk. The adopted physical picture is appropriate to describe the conditions near to the cloud core where the equation of state of the gas changes from isothermal (in the outer cloud layers) to one of "hard polytrope", and the symmetry changes from spherical to cylindrical. On the assumption of steady state, as the accreting matter passes through all spatial scales, we show that the total energy per unit mass is an invariant with respect to the fluid flow. The obtained equation describes the balance of the kinetic, thermal and gravitational energy of a fluid element. We also introduce a method for approximating density profile solutions (in a power-law form), leading to the emergence of three different regimes. We apply, as well, dynamical analysis of the motion of a fluid element. Only one of the regimes is in accordance with the two approaches (energy and force balance). It corresponds to a density profile of a slope -2, polytropic exponent 3/2, and sub-Keplerian rotation of the disk, when the gravity is balanced by the thermal pressure. It also matches with some observations and numerical works and, in particular, leads to a second power-law tail (of a slope approx. -1) of the density distribution function in dense, self-gravitating cloud regions. △ Less

Submitted 6 December, 2023; v1 submitted 18 November, 2023; originally announced November 2023.

Comments: 11 pages, 1 figure

Journal ref: MNRAS 527, 2790-2798 (2024)

A sensitive, high-resolution, wide-field IRAM NOEMA CO(1-0) survey of the very nearby spiral galaxy IC 342

Authors: M. Querejeta, J. Pety, A. Schruba, A. K. Leroy, C. N. Herrera, I-D. Chiang, S. E. Meidt, E. Rosolowsky, E. Schinnerer, K. Schuster, J. Sun, K. A. Herrmann, A. T. Barnes, I. Beslic, F. Bigiel, Y. Cao, M. Chevance, C. Eibensteiner, E. Emsellem, C. M. Faesi, A. Hughes, J. Kim, R. S. Klessen, K. Kreckel, J. M. D. Kruijssen , et al. (9 additional authors not shown)

Abstract: We present a new wide-field 10.75 x 10.75 arcmin^2 (~11x11 kpc^2), high-resolution (theta = 3.6" ~ 60 pc) NOEMA CO(1-0) survey of the very nearby (d=3.45 Mpc) spiral galaxy IC 342. The survey spans out to about 1.5 effective radii and covers most of the region where molecular gas dominates the cold interstellar medium. We resolved the CO emission into >600 individual giant molecular clouds and ass… ▽ More We present a new wide-field 10.75 x 10.75 arcmin^2 (~11x11 kpc^2), high-resolution (theta = 3.6" ~ 60 pc) NOEMA CO(1-0) survey of the very nearby (d=3.45 Mpc) spiral galaxy IC 342. The survey spans out to about 1.5 effective radii and covers most of the region where molecular gas dominates the cold interstellar medium. We resolved the CO emission into >600 individual giant molecular clouds and associations. We assessed their properties and found that overall the clouds show approximate virial balance, with typical virial parameters of alpha_vir=1-2. The typical surface density and line width of molecular gas increase from the inter-arm region to the arm and bar region, and they reach their highest values in the inner kiloparsec of the galaxy (median Sigma_mol~80, 140, 160, and 1100 M_sun/pc^2, sigma_CO~6.6, 7.6, 9.7, and 18.4 km/s for inter-arm, arm, bar, and center clouds, respectively). Clouds in the central part of the galaxy show an enhanced line width relative to their surface densities and evidence of additional sources of dynamical broadening. All of these results agree well with studies of clouds in more distant galaxies at a similar physical resolution. Leveraging our measurements to estimate the density and gravitational free-fall time at 90 pc resolution, averaged on 1.5 kpc hexagonal apertures, we estimate a typical star formation efficiency per free-fall time of 0.45% with a 16-84% variation of 0.33-0.71% among such 1.5 kpc regions. We speculate that bar-driven gas inflow could explain the large gas concentration in the central kiloparsec and the buildup of the massive nuclear star cluster. This wide-area CO map of the closest face-on massive spiral galaxy demonstrates the current mapping power of NOEMA and has many potential applications. The data and products are publicly available. △ Less

Submitted 10 October, 2023; originally announced October 2023.

Comments: 16 pages, 6 figures, accepted for publication in A&A

Journal ref: A&A 680, A4 (2023)

NEATH II: N$_2$H$^+$ as a tracer of imminent star formation in quiescent high-density gas

Authors: F. D. Priestley, P. C. Clark, S. C. O. Glover, S. E. Ragan, O. Fehér, L. R. Prole, R. S. Klessen

Abstract: Star formation activity in molecular clouds is often found to be correlated with the amount of material above a column density threshold of $\sim 10^{22} \, {\rm cm^{-2}}$. Attempts to connect this column density threshold to a ${\it volume}$ density above which star formation can occur are limited by the fact that the volume density of gas is difficult to reliably measure from observations. We po… ▽ More Star formation activity in molecular clouds is often found to be correlated with the amount of material above a column density threshold of $\sim 10^{22} \, {\rm cm^{-2}}$. Attempts to connect this column density threshold to a ${\it volume}$ density above which star formation can occur are limited by the fact that the volume density of gas is difficult to reliably measure from observations. We post-process hydrodynamical simulations of molecular clouds with a time-dependent chemical network, and investigate the connection between commonly-observed molecular species and star formation activity. We find that many molecules widely assumed to specifically trace the dense, star-forming component of molecular clouds (e.g. HCN, HCO$^+$, CS) actually also exist in substantial quantities in material only transiently enhanced in density, which will eventually return to a more diffuse state without forming any stars. By contrast, N$_2$H$^+$ only exists in detectable quantities above a volume density of $10^4 \, {\rm cm^{-3}}$, the point at which CO, which reacts destructively with N$_2$H$^+$, begins to deplete out of the gas phase onto grain surfaces. This density threshold for detectable quantities of N$_2$H$^+$ corresponds very closely to the volume density at which gas becomes irreversibly gravitationally bound in the simulations: the material traced by N$_2$H$^+$ never reverts to lower densities, and quiescent regions of molecular clouds with visible N$_2$H$^+$ emission are destined to eventually form stars. The N$_2$H$^+$ line intensity is likely to directly correlate with the star formation rate averaged over timescales of around a Myr. △ Less

Submitted 9 October, 2023; originally announced October 2023.

Comments: 10 pages, 10 figures. MNRAS accepted

Turbulent Structure In Supernova Remnants G46.8-0.3 And G39.2-0.3 From THOR Polarimetry

Authors: Russell Shanahan, Jeroen Stil, Loren Anderson, Henrik Beuther, Paul Goldsmith, Ralf S. Klessen, Michael Rugel, Juan D. Soler

Abstract: We present the continued analysis of polarization and Faraday rotation for the supernova remnants (SNRs) G46.8-0.3 and G39.2-0.3 in L-band (1-2 GHz) radio continuum in The HI/OH/Recombination line (THOR) survey. In this work, we present our investigation of Faraday depth fluctuations from angular scales comparable to the size of the SNRs down to scales less than our 16" beam (<~0.7 pc) from Farada… ▽ More We present the continued analysis of polarization and Faraday rotation for the supernova remnants (SNRs) G46.8-0.3 and G39.2-0.3 in L-band (1-2 GHz) radio continuum in The HI/OH/Recombination line (THOR) survey. In this work, we present our investigation of Faraday depth fluctuations from angular scales comparable to the size of the SNRs down to scales less than our 16" beam (<~0.7 pc) from Faraday dispersion (sigma_phi). From THOR, we find median sigma_phi of 15.9 +/- 3.2 rad m^-2 for G46.8-0.3 and 17.6 +/- 1.6 rad m^-2 for G39.2-0.3. When comparing to polarization at 6cm, we find evidence for sigma_phi >~ 30 rad m^-2 in localized regions where we detect no polarization in THOR. We combine Faraday depth dispersion with the rotation measure (RM) structure function (SF) and find evidence for a break in the SF on scales less than the THOR beam. We estimate the RM SF of the foreground interstellar medium (ISM) using the SF of extra-Galactic radio sources (EGRS) and pulsars to find that the RM fluctuations we measure originate within the SNRs for all but the largest angular scales. △ Less

Submitted 3 October, 2023; originally announced October 2023.

Comments: Has been accepted by ApJ for publication. Figures 3 and 4 did do not render well when using an internet browser to view them, but when the pdf file is downloaded these figures look as they should

Cold atomic gas identified by HI self-absorption. Cold atomic clouds toward giant molecular filaments

Authors: J. Syed, H. Beuther, P. F. Goldsmith, Th. Henning, M. Heyer, R. S. Klessen, J. M. Stil, J. D. Soler, L. D. Anderson, J. S. Urquhart, M. R. Rugel, K. G. Johnston, A. Brunthaler

Abstract: Stars form in the dense interiors of molecular clouds. The dynamics and physical properties of the atomic interstellar medium (ISM) set the conditions under which molecular clouds and eventually stars will form. It is, therefore, critical to investigate the relationship between the atomic and molecular gas phase to understand the global star formation process. Using the high angular resolution dat… ▽ More Stars form in the dense interiors of molecular clouds. The dynamics and physical properties of the atomic interstellar medium (ISM) set the conditions under which molecular clouds and eventually stars will form. It is, therefore, critical to investigate the relationship between the atomic and molecular gas phase to understand the global star formation process. Using the high angular resolution data from The HI/OH/Recombination line survey of the Milky Way (THOR), we aim to constrain the kinematic and physical properties of the cold atomic hydrogen gas phase toward the inner Galactic plane. HI self-absorption (HISA) has proven to be a viable method to detect cold atomic hydrogen clouds in the Galactic plane. With the help of a newly developed self-absorption extraction routine (astroSABER), we build upon previous case studies to identify HI self-absorption toward a sample of Giant Molecular Filaments (GMFs). We find the cold atomic gas to be spatially correlated with the molecular gas on a global scale. The column densities of the cold atomic gas traced by HISA are usually of the order of $10^{20}\rm\,cm^{-2}$ whereas those of molecular hydrogen traced by $\rm^{13}CO$ are at least an order of magnitude higher. The HISA column densities are attributed to a cold gas component that accounts for a fraction of $\sim$5% of the total atomic gas budget within the clouds. The HISA column density distributions show pronounced log-normal shapes that are broader than those traced by HI emission. The cold atomic gas is found to be moderately supersonic with Mach numbers of a $\sim$few. In contrast, highly supersonic dynamics drive the molecular gas within most filaments. △ Less

Submitted 3 October, 2023; originally announced October 2023.

Comments: 41 pages, 28 figures, accepted for publication in A&A

Journal ref: A&A 679, A130 (2023)

Investigating the Drivers of Electron Temperature Variations in HII Regions with Keck-KCWI and VLT-MUSE

Authors: Ryan J. Rickards Vaught, Karin M. Sandstrom, Francesco Belfiore, Kathryn Kreckel, J. Eduardo Méndez-Delgado, Eric Emsellem, Brent Groves, Guillermo A. Blanc, Daniel A. Dale, Oleg V. Egorov, Simon C. O. Glover, Kathryn Grasha, Ralf S. Klessen, Justus Neumann, Thomas G. Williams

Abstract: 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]$λ\lambda7320,7330$, [SII]$λ\lambda4069,4076$, [OIII]$\lambda4363$, and [SIII]$\lambda6312$) for a… ▽ More 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]$λ\lambda7320,7330$, [SII]$λ\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]. △ Less

Submitted 7 March, 2024; v1 submitted 29 September, 2023; originally announced September 2023.

Comments: Accepted to ApJ, 66 pages, 36 figures

A new tool to derive simultaneously exponent and extremes of power-law distributions

Authors: S. Pezzuto, A. Coletta, R. S. Klessen, E. Schisano, M. Benedettini, D. Elia, S. Molinari, J. D. Soler, A. Traficante

Abstract: Many experimental quantities show a power-law distribution $p(x)\propto x^{-α}$. In astrophysics, examples are: size distribution of dust grains or luminosity function of galaxies. Such distributions are characterized by the exponent $α$ and by the extremes $x_\text{min}$ $x_\text{max}$ where the distribution extends. There are no mathematical tools that derive the three unknowns at the same time.… ▽ More Many experimental quantities show a power-law distribution $p(x)\propto x^{-α}$. In astrophysics, examples are: size distribution of dust grains or luminosity function of galaxies. Such distributions are characterized by the exponent $α$ and by the extremes $x_\text{min}$ $x_\text{max}$ where the distribution extends. There are no mathematical tools that derive the three unknowns at the same time. In general, one estimates a set of $α$ corresponding to different guesses of $x_\text{min}$ $x_\text{max}$. Then, the best set of values describing the observed data is selected a posteriori. In this paper, we present a tool that finds contextually the three parameters based on simple assumptions on how the observed values $x_i$ populate the unknown range between $x_\text{min}$ and $x_\text{max}$ for a given $α$. Our tool, freely downloadable, finds the best values through a non-linear least-squares fit. We compare our technique with the maximum likelihood estimators for power-law distributions, both truncated and not. Through simulated data, we show for each method the reliability of the computed parameters as a function of the number $N$ of data in the sample. We then apply our method to observed data to derive: i) the slope of the core mass function in the Perseus star-forming region, finding two power-law distributions: $α=2.576$ between $1.06\,M_{\sun}$ and $3.35\,M_{\sun}$, $α=3.39$ between $3.48\,M_{\sun}$ and $33.4\,M_{\sun}$; ii) the slope of the $γ$-ray spectrum of the blazar J0011.4+0057, extracted from the Fermi-LAT archive. For the latter case, we derive $α=2.89$ between 1,484~MeV and 28.7~GeV; then we derive the time-resolved slopes using subsets 200 photons each. △ Less

Submitted 19 September, 2023; v1 submitted 28 August, 2023; originally announced August 2023.

Comments: Corrected few typos found at proof-reading stage. The only important modification is in Table 4 where "x_M not constrained" is now "x_M=40"

Journal ref: MNRAS, 525, 4744 (2023)

A deep-learning approach to the 3D reconstruction of dust density and temperature in star-forming regions

Authors: Victor F. Ksoll, Stefan Reissl, Ralf S. Klessen, Ian W. Stephens, Rowan J. Smith, Juan D. Soler, Alessio Traficante, Leonardo Testi, Patrick Hennebelle, Sergio Molinari

Abstract: Aims: We introduce a new deep-learning approach for the reconstruction of 3D dust density and temperature distributions from multi-wavelength dust emission observations on the scale of individual star-forming cloud cores (<0.2pc). Methods: We construct a training data set by processing cloud cores from the Cloud Factory simulations with the POLARIS radiative transfer code to produce synthetic du… ▽ More Aims: We introduce a new deep-learning approach for the reconstruction of 3D dust density and temperature distributions from multi-wavelength dust emission observations on the scale of individual star-forming cloud cores (<0.2pc). Methods: We construct a training data set by processing cloud cores from the Cloud Factory simulations with the POLARIS radiative transfer code to produce synthetic dust emission observations at 23 wavelengths between 12 and 1300 $μ$m. We simplify the task by reconstructing the cloud structure along individual lines of sight and train a conditional invertible neural network (cINN) for this purpose. The cINN belongs to the group of normalising flow methods and is able to predict full posterior distributions for the target dust properties. We test different cINN setups, ranging from a scenario that includes all 23 wavelengths down to a more realistically limited case with observations at only seven wavelengths. We evaluate the predictive performance of these models on synthetic test data. Results: We report an excellent reconstruction performance for the 23-wavelengths cINN model, achieving median absolute relative errors of about 1.8% in $\log(n/m^{-3})$ and 1% in $\log(T_{dust}/K)$, respectively. We identify trends towards overestimation at the low end of the density range and towards underestimation at the high end of both the density and temperature values, which may be related to a bias in the training data. Limiting our coverage to a combination of only seven wavelengths, we still find a satisfactory performance with average absolute relative errors of about 3.3% and 2.5% in $\log(n/m^{-3})$ and $\log(T_{dust}/K)$. Conclusions: This proof-of-concept study shows that the cINN-based approach for 3D reconstruction of dust density and temperature is very promising and even compatible with a more realistically constrained wavelength coverage. △ Less

Submitted 9 February, 2024; v1 submitted 18 August, 2023; originally announced August 2023.

Comments: 38 pages, 25 figures, accepted for publication in Astronomy & Astrophysics on 26.01.2024

A constant N$_2$H$^+$(1-0)-to-HCN(1-0) ratio on kiloparsec scales

Authors: M. J. Jiménez-Donaire, A. Usero, I. Bešlić, M. Tafalla, A. Chacón-Tanarro, Q. Salomé, C. Eibensteiner, A. García-Rodríguez, A. Hacar, A. T. Barnes, F. Bigiel, M. Chevance, D. Colombo, D. A. Dale, T. A. Davis, S. C. O. Glover, J. Kauffmann, R. S. Klessen, A. K. Leroy, L. Neumann, H. Pan, J. Pety, M. Querejeta, T. Saito, E. Schinnerer , et al. (2 additional authors not shown)

Abstract: Nitrogen hydrides such as NH$_3$ and N$_2$H$^+$ are widely used by Galactic observers to trace the cold dense regions of the interstellar medium. In external galaxies, because of limited sensitivity, HCN has become the most common tracer of dense gas over large parts of galaxies. We provide the first systematic measurements of N$_2$H$^+$(1-0) across different environments of an external spiral gal… ▽ More Nitrogen hydrides such as NH$_3$ and N$_2$H$^+$ are widely used by Galactic observers to trace the cold dense regions of the interstellar medium. In external galaxies, because of limited sensitivity, HCN has become the most common tracer of dense gas over large parts of galaxies. We provide the first systematic measurements of N$_2$H$^+$(1-0) across different environments of an external spiral galaxy, NGC6946. We find a strong correlation ($r>0.98,p<0.01$) between the HCN(1-0) and N$_2$H$^+$(1-0) intensities across the inner $\sim8\mathrm{kpc}$ of the galaxy, at kiloparsec scales. This correlation is equally strong between the ratios N$_2$H$^+$(1-0)/CO(1-0) and HCN(1-0)/CO(1-0), tracers of dense gas fractions ($f_\mathrm{dense}$). We measure an average intensity ratio of N$_2$H$^+$(1-0)/HCN(1-0)$=0.15\pm0.02$ over our set of five IRAM-30m pointings. These trends are further supported by existing measurements for Galactic and extragalactic sources. This narrow distribution in the average ratio suggests that the observed systematic trends found in kiloparsec-scale extragalactic studies of $f_\mathrm{dense}$ and the efficiency of dense gas (SFE$_\mathrm{dense}$) would not change if we employed N$_2$H$^+$(1-0) as a more direct tracer of dense gas. At kiloparsec scales our results indicate that the HCN(1-0) emission can be used to predict the expected N$_2$H$^+$(1-0) over those regions. Our results suggest that, even if HCN(1-0) and N$_2$H$^+$(1-0) trace different density regimes within molecular clouds, subcloud differences average out at kiloparsec scales, yielding the two tracers proportional to each other. △ Less

Submitted 2 August, 2023; originally announced August 2023.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 676, L11 (2023)

A comparison of the Milky Way's recent star formation revealed by dust thermal emission and high-mass stars

Authors: J. D. Soler, E. Zari, D. Elia, S. Molinari, C. Mininni, E. Schisano, A. Traficante, R. S. Klessen, S. C. O. Glover, P. Hennebelle, T. Colman, N. Frankel, T. Wenger

Abstract: We present a comparison of the Milky Way's star formation rate (SFR) surface density ($Σ_{\rm SFR}$) obtained with two independent state-of-the-art observational methods. The first method infers $Σ_{\rm SFR}$ from observations of the dust thermal emission from interstellar dust grains in far-infrared wavelengths registered in the Herschel infrared Galactic Plane Survey (Hi-GAL). The second method… ▽ More We present a comparison of the Milky Way's star formation rate (SFR) surface density ($Σ_{\rm SFR}$) obtained with two independent state-of-the-art observational methods. The first method infers $Σ_{\rm SFR}$ from observations of the dust thermal emission from interstellar dust grains in far-infrared wavelengths registered in the Herschel infrared Galactic Plane Survey (Hi-GAL). The second method determines $Σ_{\rm SFR}$ by modeling the current population of O-, B-, and A-type stars in a 6 kpc $\times$ 6 kpc area around the Sun. We find an agreement between the two methods within a factor of two for the mean SFRs and the SFR surface density profiles. Given the broad differences between the observational techniques and the independent assumptions in the methods for computing the SFRs, this agreement constitutes a significant advance in our understanding of the star formation of our Galaxy and implies that the local SFR has been roughly constant over the past 10\,Myr. △ Less

Submitted 18 September, 2023; v1 submitted 2 August, 2023; originally announced August 2023.

Comments: Accepted for publication at Astronomy & Astrophysics (11SEP2023). 8 pages, 7 figures

Journal ref: A&A 678, A95 (2023)

Non-Equilibrium Abundances Treated Holistically (NEATH): the molecular composition of star-forming clouds

Authors: F. D. Priestley, P. C. Clark, S. C. O. Glover, S. E. Ragan, O. Fehér, L. R. Prole, R. S. Klessen

Abstract: Much of what we know about molecular clouds, and by extension star formation, comes from molecular line observations. Interpreting these correctly requires knowledge of the underlying molecular abundances. Simulations of molecular clouds typically only model species that are important for the gas thermodynamics, which tend to be poor tracers of the denser material where stars form. We construct a… ▽ More Much of what we know about molecular clouds, and by extension star formation, comes from molecular line observations. Interpreting these correctly requires knowledge of the underlying molecular abundances. Simulations of molecular clouds typically only model species that are important for the gas thermodynamics, which tend to be poor tracers of the denser material where stars form. We construct a framework for post-processing these simulations with a full time-dependent chemical network, allowing us to model the behaviour of observationally-important species not present in the reduced network used for the thermodynamics. We use this to investigate the chemical evolution of molecular gas under realistic physical conditions. We find that molecules can be divided into those which reach peak abundances at moderate densities ($10^3 \, {\rm cm^{-3}}$) and decline sharply thereafter (such as CO and HCN), and those which peak at higher densities and then remain roughly constant (e.g. NH$_3$, N$_2$H$^+$). Evolving the chemistry with physical properties held constant at their final values results in a significant overestimation of gas-phase abundances for all molecules, and does not capture the drastic variations in abundance caused by different evolutionary histories. The dynamical evolution of molecular gas cannot be neglected when modelling its chemistry. △ Less

Submitted 24 July, 2023; originally announced July 2023.

Comments: 14 pages, 13 figures. MNRAS accepted

Quantifying the energy balance between the turbulent ionised gas and young stars

Authors: Oleg V. Egorov, Kathryn Kreckel, Simon C. O. Glover, Brent Groves, Francesco Belfiore, Eric Emsellem, Ralf S. Klessen, Adam K. Leroy, Sharon E. Meidt, Sumit K. Sarbadhicary, Eva Schinnerer, Elizabeth J. Watkins, Brad C. Whitmore, Ashley T. Barnes, Enrico Congiu, Daniel A. Dale, Kathryn Grasha, Kirsten L. Larson, Janice C. Lee, J. Eduardo Méndez-Delgado, David A. Thilker, Thomas G. Williams

Abstract: 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 fee… ▽ More 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$α$ line, we select 1484 regions of locally elevated velocity dispersion ($σ$(H$α$)>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. △ Less

Submitted 17 August, 2023; v1 submitted 18 July, 2023; originally announced July 2023.

Comments: Accepted by A&A. The abstract is abridged. 31 pages (including 8 pages in appendix), 20 figures, 2 tables

Journal ref: A&A 678, A153 (2023)

A reproduction of the Milky Way's Faraday rotation measure map in galaxy simulations from global to local scales

Authors: Stefan Reissl, Ralf S. Klessen, Eric W. Pellegrini, Daniel Rahner, Rüdiger Pakmor, Robert Grand, Facundo Gomez, Federico Marinacci, Volker Springel

Abstract: Magnetic fields are of critical importance for our understanding of the origin and long-term evolution of the Milky Way. This is due to their decisive role in the dynamical evolution of the interstellar medium (ISM) and their influence on the star-formation process. Faraday rotation measures (RM) along many different sightlines across the Galaxy are a primary means to infer the magnetic field topo… ▽ More Magnetic fields are of critical importance for our understanding of the origin and long-term evolution of the Milky Way. This is due to their decisive role in the dynamical evolution of the interstellar medium (ISM) and their influence on the star-formation process. Faraday rotation measures (RM) along many different sightlines across the Galaxy are a primary means to infer the magnetic field topology and strength from observations. However, the interpretation of the data has been hampered by the failure of previous attempts to explain the observations in theoretical models and to synthesize a realistic multi-scale all-sky RM map. We here utilize a cosmological magnetohydrodynamic (MHD) simulation of the formation of the Milky Way, augment it with a novel star cluster population synthesis model for a more realistic structure of the local interstellar medium, and perform detailed polarized radiative transfer calculations on the resulting model. This yields a faithful first principles prediction of the Faraday sky as observed on Earth. The results reproduce the observations of the Galaxy not only on global scales, but also on local scales of individual star-forming clouds. They also imply that the Local Bubble containing our Sun dominates the RM signal over large regions of the sky. Modern cosmological MHD simulations of the Milky Way's formation, combined with a simple and plausible model for the fraction of free electrons in the ISM, explain the RM observations remarkably well, thus indicating the emergence of a firm theoretical understanding of the genesis of magnetic fields in our Universe across cosmic time. △ Less

Submitted 11 July, 2023; originally announced July 2023.

Comments: 14 pages, 3 figures

Calibrating mid-infrared emission as a tracer of obscured star formation on HII-region scales in the era of JWST

Authors: Francesco Belfiore, Adam K. Leroy, Thomas G. Williams, Ashley T. Barnes, Frank Bigiel, Médéric Boquien, Yixian Cao, Jérémy Chastenet, Enrico Congiu, Daniel A. Dale, Oleg V. Egorov, Cosima Eibensteiner, Eric Emsellem, Simon C. O. Glover, Brent Groves, Hamid Hassani, Ralf S. Klessen, Kathryn Kreckel, Lukas Neumann, Justus Neumann, Miguel Querejeta, Erik Rosolowsky, Patricia Sanchez-Blazquez, Karin Sandstrom, Eva Schinnerer , et al. (3 additional authors not shown)

Abstract: 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 dus… ▽ More 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 μm$, and the emission in the PAH-tracing bands at 3.3, 7.7, 10, and 11.3$\rm μ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$α$ flux as the sum of the observed H$α$ 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$α$ 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 μ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 μm$ decreases, signalling destruction of the PAH molecules. For fainter HII regions all bands suffer from an increasing contamination from the diffuse infrared background. △ Less

Submitted 1 September, 2023; v1 submitted 20 June, 2023; originally announced June 2023.

Comments: accepted for publication in A&A

Journal ref: A&A 678, A129 (2023)

Constraining the LyC escape fraction from LEGUS star clusters with SIGNALS HII region observations: A pilot study of NGC 628

Authors: J. W. Teh, K. Grasha, M. R. Krumholz, A. Battisti, D. Calzetti, L. Rousseau-Nepton, C. Rhea, A. Adamo, R. C. Kennicutt, E. K. Grebel, D. O. Cook, F. Combes, M. Messa, S. Linden, R. S. Klessen, J. M. Vilchez, M. Fumagalli, A. F. McLeod, L. J. Smith, L. Chemin, J. Wang, E. Sabbi, E. Sacchi, A. Petric, L. Della Bruna , et al. (1 additional authors not shown)

Abstract: The ionising radiation of young and massive stars is a crucial form of stellar feedback. Most ionising (Lyman-continuum; LyC, $λ< 912A$) photons are absorbed close to the stars that produce them, forming compact HII regions, but some escape into the wider galaxy. Quantifying the fraction of LyC photons that escape is an open problem. In this work, we present a semi-novel method to estimate the esc… ▽ More The ionising radiation of young and massive stars is a crucial form of stellar feedback. Most ionising (Lyman-continuum; LyC, $λ< 912A$) photons are absorbed close to the stars that produce them, forming compact HII regions, but some escape into the wider galaxy. Quantifying the fraction of LyC photons that escape is an open problem. In this work, we present a semi-novel method to estimate the escape fraction by combining broadband photometry of star clusters from the Legacy ExtraGalactic UV Survey (LEGUS) with HII regions observed by the Star formation, Ionized gas, and Nebular Abundances Legacy Survey (SIGNALS) in the nearby spiral galaxy NGC 628. We first assess the completeness of the combined catalogue, and find that 49\% of HII regions lack corresponding star clusters as a result of a difference in the sensitivities of the LEGUS and SIGNALS surveys. For HII regions that do have matching clusters, we infer the escape fraction from the difference between the ionising power required to produce the observed HII luminosity and the predicted ionising photon output of their host star clusters; the latter is computed using a combination of LEGUS photometric observations and a stochastic stellar population synthesis code SLUG (Stochastically Lighting Up Galaxies). Overall, we find an escape fraction of $f_{esc} = 0.09^{+0.06}_{-0.06}$ across our sample of 42 HII regions; in particular, we find HII regions with high $f_{esc}$ are predominantly regions with low H$α$-luminosity. We also report possible correlation between $f_{esc}$ and the emission lines [O ii]/[N ii] and [O ii]/H$β$. △ Less

Submitted 8 June, 2023; originally announced June 2023.

Comments: Accepted for publication at MNRAS