Showing 1–50 of 264 results for author: Krumholz, M R

A novel numerical method for mixed-frame multigroup radiation-hydrodynamics with GPU acceleration implemented in the QUOKKA code

Authors: Chong-Chong He, Benjamin D. Wibking, Mark R. Krumholz

Abstract: Mixed-frame formulations of radiation-hydrodynamics (RHD), where the radiation quantities are computed in an inertial frame but matter quantities are in a comoving frame, are advantageous because they admit algorithms that conserve energy and momentum to machine precision and combine more naturally with adaptive mesh techniques, since unlike pure comoving-frame methods they do not face the problem… ▽ More Mixed-frame formulations of radiation-hydrodynamics (RHD), where the radiation quantities are computed in an inertial frame but matter quantities are in a comoving frame, are advantageous because they admit algorithms that conserve energy and momentum to machine precision and combine more naturally with adaptive mesh techniques, since unlike pure comoving-frame methods they do not face the problem that radiation quantities must change frame every time a cell is refined or coarsened. However, implementing multigroup RHD in a mixed-frame formulation presents challenges due to the complexity of handling frequency-dependent interactions and the Doppler shift of radiation boundaries. In this paper, we introduce a novel method for multigroup RHD that integrates a mixed-frame formulation with a piecewise powerlaw approximation for frequency dependence within groups. This approach ensures the exact conservation of total energy and momentum while effectively managing the Lorentz transformation of group boundaries and evaluation of group-averaged opacities. Our method takes advantage of the locality of matter-radiation coupling, allowing the source term for $N_g$ frequency groups to be handled with simple equations with a sparse Jacobian matrix of size $N_g + 1$, which can be inverted with $O(N_g)$ complexity. This results in a computational complexity that scales linearly with $N_g$ and requires no more communication than a pure hydrodynamics update, making it highly efficient for massively parallel and GPU-based systems. We implement our method in the GPU-accelerated RHD code QUOKKA and demonstrate that it passes a wide range of numerical tests. We demonstrate that the piecewise powerlaw method shows significant advantages over traditional opacity averaging methods for handling rapidly variable opacities with modest frequency resolution. △ Less

Submitted 25 July, 2024; originally announced July 2024.

Comments: 17 pages, 8 figures, submitted to MNRAS. Comments are welcome. The code described in this article is available at: https://github.com/quokka-astro/quokka

Suppressed Cosmic Ray Energy Densities in Molecular Clouds From Streaming Instability-Regulated Transport

Authors: Margot Fitz Axen, Stella Offner, Phillip F. Hopkins, Mark R. Krumholz, Michael Y. Grudic

Abstract: Cosmic rays (CRs) are the primary driver of ionization in star forming molecular clouds (MCs). Despite their potential impacts on gas dynamics and chemistry, no simulations of star cluster formation following the creation of individual stars have included explicit cosmic ray transport (CRT) to date. We conduct the first numerical simulations following the collapse of a $2000 M_{\odot}$ MC and the… ▽ More Cosmic rays (CRs) are the primary driver of ionization in star forming molecular clouds (MCs). Despite their potential impacts on gas dynamics and chemistry, no simulations of star cluster formation following the creation of individual stars have included explicit cosmic ray transport (CRT) to date. We conduct the first numerical simulations following the collapse of a $2000 M_{\odot}$ MC and the subsequent star formation including CRT using the STARFORGE framework implemented in the GIZMO code. We show that when CR-transport is streaming-dominated, the CR energy in the cloud is strongly attenuated due to energy losses from the streaming instability. Consequently, in a Milky Way like environment the median CR ionization rate (CRIR) in the cloud is low ($ ζ\lesssim 2 \times 10^{-19} \rm s^{-1}$) during the main star forming epoch of the calculation and the impact of CRs on the star formation in the cloud is limited. However, in high-CR environments, the CR distribution in the cloud is elevated ($ζ\lesssim 6 \times 10^{-18}$), and the relatively higher CR pressure outside the cloud causes slightly earlier cloud collapse and increases the star formation efficiency (SFE) by $50 \%$ to $\sim 13 \%$. The initial mass function (IMF) is similar in all cases except with possible variations in a high-CR environment. Further studies are needed to explain the range of ionization rates observed in MCs and explore star formation in extreme CR environments. △ Less

Submitted 24 July, 2024; originally announced July 2024.

Comments: Accepted to The Astrophysical Journal

Comparing metallicity correlations in nearby non-AGN and AGN-host galaxies

Authors: Song-lin Li, Zefeng Li, Emily Wisnioski, Mark R. Krumholz, Sebastián F. Sánchez

Abstract: The gas-phase metallicity distribution within galaxies records critical information about galactic evolution. In this work we investigate how active galactic nuclei (AGN) influence this distribution by measuring the two-point correlation functions of gas-phase metallicity in 95 non-AGN and 37 AGN-host galaxies from the Calar Alto Legacy Integral Field spectroscopy Area integral field spectrographi… ▽ More The gas-phase metallicity distribution within galaxies records critical information about galactic evolution. In this work we investigate how active galactic nuclei (AGN) influence this distribution by measuring the two-point correlation functions of gas-phase metallicity in 95 non-AGN and 37 AGN-host galaxies from the Calar Alto Legacy Integral Field spectroscopy Area integral field spectrographic survey. We measure metallicity using a novel Bayesian method that properly includes both stellar and AGN contributions to emission line fluxes and allows us to measure metallicities in both AGN-host and non-AGN galaxies in a single, consistent framework. We find that the two-point correlation functions of both AGN-host and non-AGN galaxies are well-fit by a simple injection-diffusion model, and that the correlation lengths $l_\mathrm{corr}$ we derive for the non-AGN galaxies are reasonably consistent with those obtained in earlier work. The AGN-host galaxies generally have smaller $l_\mathrm{corr}$ than non-AGN galaxies at fixed stellar mass, but similar $l_\mathrm{corr}$ at fixed star formation rate (SFR), suggesting that the primary effect of hosting an AGN in this sample is a reduction in SFR at fixed stellar mass, and that this in turn suppresses the correlation length. Our findings further indicate that, while both SFR and stellar mass are positively correlated with metallicity correlation length $l_\mathrm{corr}$, the former is more fundamental, implying that fluctuations in the metallicity distribution within galaxies are driven more by short-term responses to physical processes such as star formation that can change much faster than a Hubble time. △ Less

Submitted 5 July, 2024; originally announced July 2024.

Comments: 13 pages, 10 figures, submitted to MNRAS

TeV $γ$-ray emission near globular cluster Terzan 5 as a probe of cosmic ray transport

Authors: Mark R. Krumholz, Roland M. Crocker, Arash Bahramian, Pol Bordas

Abstract: Cosmic rays travelling through interstellar space have their propagation directions repeatedly scattered by fluctuating interstellar magnetic fields. The nature of this scattering is a major unsolved problem in astrophysics, one that has resisted solution largely due to a lack of direct observational constraints on the scattering rate. Here we show that very high-energy $γ$-ray emission from the g… ▽ More Cosmic rays travelling through interstellar space have their propagation directions repeatedly scattered by fluctuating interstellar magnetic fields. The nature of this scattering is a major unsolved problem in astrophysics, one that has resisted solution largely due to a lack of direct observational constraints on the scattering rate. Here we show that very high-energy $γ$-ray emission from the globular cluster Terzan 5, which has unexpectedly been found to be displaced from the cluster, presents a direct probe of this process. We show that this displacement is naturally explained by cosmic rays accelerated in the bow shock around the cluster propagating a finite distance before scattering processes re-orient enough of them towards Earth to produce a detectable $γ$-ray signal. The angular distance between the cluster and the signal places tight constraints on the scattering rate, which we show are consistent with a model whereby scattering is primarily due to excitation of magnetic waves by the cosmic rays themselves. The analysis method we develop here will make it possible to use sources with similarly displaced non-thermal X-ray and TeV $γ$-ray signals as direct probes of cosmic ray scattering across a range of Galactic environments. △ Less

Submitted 26 June, 2024; originally announced June 2024.

Comments: Accepted to Nature Astronomy; 16 pages, 7 figures; this submission contains the main article, supplementary information, and extended data in a single document

On the limitations of H alpha luminosity as a star formation tracer in spatially resolved observations

Authors: Zipeng Hu, Benjamin D. Wibking, Mark R. Krumholz, Christoph Federrath

Abstract: This study examines the limitations of H$α$ luminosity as a tracer of star formation rates (SFR) in spatially resolved observations. We carry out high-resolution simulations of a Milky Way-like galaxy including both supernova and photoionization feedback, and from these we generate synthetic H$α$ emission maps that we compare to maps of the true distribution of young stellar objects (YSOs) on scal… ▽ More This study examines the limitations of H$α$ luminosity as a tracer of star formation rates (SFR) in spatially resolved observations. We carry out high-resolution simulations of a Milky Way-like galaxy including both supernova and photoionization feedback, and from these we generate synthetic H$α$ emission maps that we compare to maps of the true distribution of young stellar objects (YSOs) on scales from whole-galaxy to individual molecular clouds ($\lesssim 100$ pc). Our results reveal significant spatial mismatches between H$α$ and true YSO maps on sub-100 pc scales, primarily due to ionizing photon leakage, with a secondary contribution from young stars drifting away from their parent molecular clouds. On small scales these effect contribute significantly to the observed anti-correlation between gas and star formation, such that there is noticeably less anti-correlation if we replace an H$α$-based star formation map with a YSO-based one; this in turn implies that previous studies have underestimated the time it takes for young stars to disperse their parent molecular clouds. However, these effects are limited in dense regions with hydrogen columns $N_\mathrm{H} > 3 \times 10^{21}$ cm$^{-2}$, where the H$α$- and YSO-based SFR maps show better agreement. Based on this finding we propose a calibration model that can precisely measure the SFR of large molecular clouds (mean radius > 100 pc) with a combination of H$α$ and CO observations, which provides a foundation for future study of star formation processes in extragalactic molecular clouds. △ Less

Submitted 16 June, 2024; originally announced June 2024.

Formation of Giant Clumps in High-$z$ Disc Galaxies by Compressive Turbulence

Authors: Nir Mandelker, Omry Ginzburg, Avishai Dekel, Frederic Bournaud, Mark R. Krumholz, Daniel Ceverino, Joel Primack

Abstract: We address the formation of giant clumps in violently unstable gas-rich disc galaxies at cosmic noon. While these are commonly thought to originate from gravitational Toomre instability, cosmological simulations have indicated that clumps form even in regions where the Toomre $Q$ parameter is well above unity, which should be stable according to linear Toomre theory (Inoue et al., 2016). Examining… ▽ More We address the formation of giant clumps in violently unstable gas-rich disc galaxies at cosmic noon. While these are commonly thought to originate from gravitational Toomre instability, cosmological simulations have indicated that clumps form even in regions where the Toomre $Q$ parameter is well above unity, which should be stable according to linear Toomre theory (Inoue et al., 2016). Examining one of these cosmological simulations, we find that it exhibits an excess in compressive modes of turbulence with converging motions. The energy in converging motions within proto-clump regions is $\sim 70\%$ of the total turbulent energy, compared to $\sim 17\%$ expected in equipartition. When averaged over the whole disc, $\sim 32\%$ of the turbulent energy is in converging motions, with a further $\sim 8\%$ in diverging motions. Thus, a total of $\sim 40\%$ of the turbulent energy is in compressive modes, with the rest in solenoidal modes, compared to the $(1/3)-(2/3)$ division expected in equipartition. By contrast, we find that in an isolated-disc simulation with similar properties, resembling high-$z$ star-forming galaxies, the energy in the different turbulence modes are in equipartition, both in proto-clump regions and over the whole disc. We conclude that the origin of the excessive converging motions in proto-clump regions is external to the disc, and propose several mechanisms that can induce them. This is an additional mechanism for clump formation, complementary to and possibly preceding gravitational instability. △ Less

Submitted 11 June, 2024; originally announced June 2024.

Comments: 6 pages, 4 figures, Submitted to MNRAS Letters Comments welcome!

Feedback in Emerging Extragalactic Star Clusters (JWST--FEAST): Calibration of Star Formation Rates in the Mid-Infrared with NGC 628

Authors: Daniela Calzetti, Angela Adamo, Sean T. Linden, Benjamin Gregg, Mark R. Krumholz, Varun Bajaj, Arjan Bik, Michele Cignoni, Matteo Correnti, Bruce Elmegreen, Helena Faustino Vieira, John S. Gallagher, Kathryn Grasha, Robert A. Gutermuth, Kelsey E. Johnson, Matteo Messa, Jens Melinder, Goran Ostlin, Alex Pedrini, Elena Sabbi, Linda J. Smith, Monica Tosi

Abstract: New JWST near-infrared imaging of the nearby galaxy NGC 628 from the Cycle 1 program JWST-FEAST is combined with archival JWST mid-infrared imaging to calibrate the 21 $μ$m emission as a star formation rate indicator (SFR) at $\sim$120 pc scales. The Pa$α$ ($λ$1.8756 $μ$m) hydrogen recombination emission line targeted by FEAST provides a reference SFR indicator that is relatively insensitive to du… ▽ More New JWST near-infrared imaging of the nearby galaxy NGC 628 from the Cycle 1 program JWST-FEAST is combined with archival JWST mid-infrared imaging to calibrate the 21 $μ$m emission as a star formation rate indicator (SFR) at $\sim$120 pc scales. The Pa$α$ ($λ$1.8756 $μ$m) hydrogen recombination emission line targeted by FEAST provides a reference SFR indicator that is relatively insensitive to dust attenuation, as demonstrated by combining this tracer with the HST H$α$ imaging. Our analysis is restricted to regions that appear compact in nebular line emission and are sufficiently bright to mitigate effects of both age and stochastic sampling of the stellar initial mass function. We find that the 21 $μ$m emission closely correlates with the nebular line emission, with a power-law with exponent=1.07$\pm$0.01, in agreement with past results. We calibrate a hybrid SFR indicator using a combination of H$α$ and 24 $μ$m (extrapolated from 21 $μ$m) tracers and derive the proportionality constant between the two tracers $b=0.095\pm0.007$, which is $\sim$ 3-5 times larger than previous derivations using large regions/entire galaxies. We model these discrepancies as an increasing contribution to the dust heating by progressively older stellar populations for increasing spatial scales, in agreement with earlier findings that star formation is hierarchically distributed in galaxies. Thus, use of hybrid SFR indicators requires prior knowledge of the mean age of the stellar populations dominating the dust heating, which makes their application uncertain. Conversely, non-linear calibrations of SFRs from L(24) alone are more robust, with a factor $\lesssim$2.5 variation across the entire range of L(24) luminosities from HII regions to galaxies. △ Less

Submitted 3 June, 2024; originally announced June 2024.

Comments: 28 pages, 9 figures. Accepted for publication on the Astrophysical Journal

Feedback in emerging extragalactic star clusters, FEAST: The relation between 3.3 $μ$m PAH emission and Star Formation Rate traced by ionized gas in NGC 628

Authors: Benjamin Gregg, Daniela Calzetti, Angela Adamo, Varun Bajaj, Jenna E. Ryon, Sean T. Linden, Matteo Correnti, Michele Cignoni, Matteo Messa, Elena Sabbi, John S. Gallagher, Kathryn Grasha, Alex Pedrini, Robert A. Gutermuth, Jens Melinder, Ralf Kotulla, Gustavo Pérez, Mark R. Krumholz, Arjan Bik, Göran Östlin, Kelsey E. Johnson, Giacomo Bortolini, Linda J. Smith, Monica Tosi, Subhransu Maji , et al. (1 additional authors not shown)

Abstract: We present maps of ionized gas (traced by Pa$α$ and Br$α$) and 3.3 $μ$m Polycyclic Aromatic Hydrocarbon (PAH) emission in the nearby spiral galaxy NGC 628, derived from new JWST/NIRCam data from the FEAST survey. With this data, we investigate and calibrate the relation between 3.3 $μ$m PAH emission and star formation rate (SFR) in and around emerging young star clusters (eYSCs) on a scale of… ▽ More We present maps of ionized gas (traced by Pa$α$ and Br$α$) and 3.3 $μ$m Polycyclic Aromatic Hydrocarbon (PAH) emission in the nearby spiral galaxy NGC 628, derived from new JWST/NIRCam data from the FEAST survey. With this data, we investigate and calibrate the relation between 3.3 $μ$m PAH emission and star formation rate (SFR) in and around emerging young star clusters (eYSCs) on a scale of ${\sim}40$ pc. We find a tight (correlation coefficient $ρ$${\sim}$0.9) sub-linear (power-law exponent $α$${\sim}$0.75) relation between the 3.3 $μ$m PAH luminosity surface density and SFR traced by Br$α$ for compact, cospatial (within 0.16$''$ or ${\sim}$7 pc) peaks in Pa$α$, Br$α$, and 3.3 $μ$m (eYSC-I). The scatter in the relationship does not correlate well with variations in local interstellar medium (ISM) metallicity due to a radial metallicity gradient, but rather is likely due to stochastic sampling of the stellar initial mass function (IMF) and variations in the PAH heating and age of our sources. The deviation from a linear relation may be explained by PAH destruction in more intense ionizing environments, variations in age, and IMF stochasticity at intermediate to low luminosities. We test our results with various continuum subtraction techniques using combinations of NIRCam bands and find that they remain robust with only minor differences in the derived slope and intercept. An unexpected discrepancy is identified between the relations of hydrogen recombination lines (Pa$α$ versus Br$α$; H$α$ versus Br$α$). △ Less

Submitted 15 May, 2024; originally announced May 2024.

Comments: 28 pages, 10 figures, 3 tables

The Interplay between the IMF and Star Formation Efficiency through Radiative Feedback at High Stellar Surface Densities

Authors: Shyam H. Menon, Lachlan Lancaster, Blakesley Burkhart, Rachel S. Somerville, Avishai Dekel, Mark R. Krumholz

Abstract: The observed rest-UV luminosity function at cosmic dawn ($z \sim 8-14$) measured by JWST revealed an excess of UV-luminous galaxies relative to many pre-launch theoretical predictions. A high star-formation efficiency (SFE) and a top-heavy initial mass function (IMF) are among the mechanisms proposed for explaining this excess. Although a top-heavy IMF has been proposed for its ability to increase… ▽ More The observed rest-UV luminosity function at cosmic dawn ($z \sim 8-14$) measured by JWST revealed an excess of UV-luminous galaxies relative to many pre-launch theoretical predictions. A high star-formation efficiency (SFE) and a top-heavy initial mass function (IMF) are among the mechanisms proposed for explaining this excess. Although a top-heavy IMF has been proposed for its ability to increase the light-to-mass ratio (\(Ψ_{\mathrm{UV}}\)), the resulting enhanced radiative pressure from young stars could decrease the star formation efficiency (SFE), potentially driving galaxy luminosities back down. In this Letter, we use idealized radiation hydrodynamic simulations of star cluster formation to explore the effects of a top-heavy IMF on the SFE of clouds typical of the high pressure conditions found at these redshifts. We find that the SFE in star clusters with solar neighbourhood-like dust abundance decreases with increasingly top-heavy IMF's -- by $\sim 20 \%$ for an increase of factor 4 in $Ψ_{\mathrm{UV}}$, and by $50 \%$ for a factor $ \sim 10$ in $Ψ_{\mathrm{UV}}$. However, we find that an expected decrease in the dust-to-gas ratio ($\sim 0.01 \times \mathrm{Solar}$) at these redshifts can completely compensate for the enhanced light output. This leads to a (cloud-scale; $\sim 10 \, \mathrm{pc}$) SFE that is $\gtrsim 70\%$ even for a factor 10 increase in $Ψ_{\mathrm{UV}}$, implying that highly efficient star formation is unavoidable for high surface density and low metallicity conditions. Our results suggest that a top-heavy IMF, if present, likely coexists with efficient star formation in these galaxies. △ Less

Submitted 1 May, 2024; originally announced May 2024.

Comments: 12 pages, 4 figures. Accepted to the Astrophysical Journal Letters

When are galactic winds molecular?

Authors: Aditi Vijayan, Mark R. Krumholz

Abstract: The molecular phase of supernova-driven outflows originates from the cold, molecular gas in the disc of a star-forming galaxy, and may carry a substantial fraction of the wind mass flux in some galaxies, but it remains poorly understood. Observations of this phase come mostly from very nearby galaxies due its low surface brightness and covering fraction, and simulations often lack the spatial reso… ▽ More The molecular phase of supernova-driven outflows originates from the cold, molecular gas in the disc of a star-forming galaxy, and may carry a substantial fraction of the wind mass flux in some galaxies, but it remains poorly understood. Observations of this phase come mostly from very nearby galaxies due its low surface brightness and covering fraction, and simulations often lack the spatial resolution necessary to resolve it. Here we analytically estimate the survivability of this phase in order to understand under what conditions an galactic wind can contain a significant molecular phase. We show that the molecular content of outflows is primarily determined by two dimensionless numbers: a generalised Eddington ratio describing the strength of the outflow and an ionisation parameter-like quantity describing the strength of the radiation field per baryon. We apply this model to a sample of galaxies and show that, while any molecules entrained in the winds of normal star-forming galaxies should be destroyed close to the galactic disc, the outflows of strong starburst should become increasingly dominated by molecules. △ Less

Submitted 22 April, 2024; originally announced April 2024.

Comments: 8 pages, 5 figures, submitted to MNRAS. Comments welcome

An Asymptotically-Correct Implicit-Explicit Time Integration Scheme for Finite Volume Radiation-Hydrodynamics

Authors: Chong-Chong He, Benjamin D. Wibking, Mark R. Krumholz

Abstract: Numerical radiation-hydrodynamics (RHD) for non-relativistic flows is a challenging problem because it encompasses processes acting over a very broad range of timescales, and where the relative importance of these processes often varies by orders of magnitude across the computational domain. Here we present a new implicit-explicit (IMEX) method for numerical RHD that has a number of desirable prop… ▽ More Numerical radiation-hydrodynamics (RHD) for non-relativistic flows is a challenging problem because it encompasses processes acting over a very broad range of timescales, and where the relative importance of these processes often varies by orders of magnitude across the computational domain. Here we present a new implicit-explicit (IMEX) method for numerical RHD that has a number of desirable properties that have not previously been combined in a single method. Our scheme is based on moments and allows machine-precision conservation of energy and momentum, making it highly suitable for adaptive mesh refinement applications; it requires no more communication than hydrodynamics and includes no non-local iterative steps, making it highly suitable for massively parallel and GPU-based systems where communication is a bottleneck; and we show that it is asymptotically-accurate in the streaming, static diffusion, and dynamic diffusion limits, including in the so-called asymptotic diffusion regime where the computational grid does not resolve the photon mean free path. We implement our method in the GPU-accelerated RHD code QUOKKA and show that it passes a wide range of numerical tests. △ Less

Submitted 12 April, 2024; originally announced April 2024.

Comments: 14 pages, 6 figures, submitted to MNRAS; Code available for download at https://github.com/quokka-astro/quokka

The mass-metallicity and fundamental metallicity relations in non-AGN and AGN-host galaxies

Authors: Song-Lin Li, Kathryn Grasha, Mark R. Krumholz, Emily Wisnioski, Ralph S. Sutherland, Lisa J. Kewley, Yan-Mei Chen, Zefeng Li

Abstract: Galaxies' stellar masses, gas-phase oxygen abundances (metallicity), and star formation rates (SFRs) obey a series of empirical correlations, most notably the mass-metallicity relation (MZR) and fundamental metallicity relation (FZR), which relates oxygen abundance to a combination of stellar mass and SFR. However, due to the difficulty of measuring oxygen abundances and SFRs in galaxies that host… ▽ More Galaxies' stellar masses, gas-phase oxygen abundances (metallicity), and star formation rates (SFRs) obey a series of empirical correlations, most notably the mass-metallicity relation (MZR) and fundamental metallicity relation (FZR), which relates oxygen abundance to a combination of stellar mass and SFR. However, due to the difficulty of measuring oxygen abundances and SFRs in galaxies that host powerful active galactic nuclei (AGN), to date it is unknown to what extent AGN-host galaxies also follow these correlations. In this work, we apply Bayesian methods to the MaNGA integral field spectrographic (IFS) survey that allow us to measure oxygen abundances and SFRs in AGN hosts, and use these measurements to explore how the MZR and FZR differ between galaxies that do and do not host AGN. We find similar MZRs at stellar masses above $10^{10.5} \mathrm{M}_\odot$, but that at lower stellar masses AGN hosts show up to $\sim 0.2$ dex higher oxygen abundances. The offset in the FZR is significantly smaller, suggesting that the larger deviation in the MZR is a result of AGN-host galaxies having systematically lower SFRs at fixed stellar mass. However, within the AGN-host sample there is little correlation between SFR and oxygen abundance. These findings support a scenario in which an AGN can halt efficient gas accretion, which drives non-AGN host galaxies to both higher SFR and lower oxygen abundance, resulting in the galaxy evolving off the star-forming main sequence (SFMS). As a consequence, as the SFR declines for an individual system its metallicity remains mostly unchanged. △ Less

Submitted 24 March, 2024; originally announced March 2024.

Comments: 18 pages, 14 figures, accepted for publication in MNRAS

Cosmological evolution of metallicity correlation functions from the Auriga simulations

Authors: Zefeng Li, Robert J. J. Grand, Emily Wisnioski, J. Trevor Mendel, Mark R. Krumholz, Yuan-Sen Ting, Ruediger Pakmor, Facundo A. Gómez, Federico Marinacci, Ioana Ciucă

Abstract: We study the cosmological evolution of the two-point correlation functions of galactic gas-phase metal distributions using the 28 simulated galaxies from the Auriga Project. Using mock observations of the $z = 0$ snapshots to mimic our past work, we show that the correlation functions of the simulated mock observations are well matched to the correlation functions measured from local galaxy survey… ▽ More We study the cosmological evolution of the two-point correlation functions of galactic gas-phase metal distributions using the 28 simulated galaxies from the Auriga Project. Using mock observations of the $z = 0$ snapshots to mimic our past work, we show that the correlation functions of the simulated mock observations are well matched to the correlation functions measured from local galaxy surveys. This comparison suggests that the simulations capture the processes important for determining metal correlation lengths, the key parameter in metallicity correlation functions. We investigate the evolution of metallicity correlations over cosmic time using the true simulation data, showing that individual galaxies undergo no significant systematic evolution in their metal correlation functions from $z\sim 3$ to today. In addition, the fluctuations in metal correlation length are correlated with but lag ahead fluctuations in star formation rate. This suggests that re-arrangement of metals within galaxies occurs at a higher cadence than star formation activity, and is more sensitive to the changes of environment, such as galaxy mergers, gas inflows / outflows, and fly-bys. △ Less

Submitted 13 February, 2024; originally announced February 2024.

Comments: 12 pages, 10 figures, 1 table, accepted for publication in MNRAS

The Molecular Cloud Lifecycle II: Formation and Destruction of Molecular Clouds Diagnosed via H$_2$ Fluorescent Emission Emission

Authors: Blakesley Burkhart, Shmuel Bialy, Daniel Seifried, Stefanie Walch, Erika Hamden, Thomas J. Haworth, Keri Hoadley, Shuo Kong, Madisen Johnson, Sarah Jeffreson, Mark R. Krumholz, Min-Young Lee, Amiel Sternberg, Neal J. Turner

Abstract: Molecular hydrogen (H$_2$) formation and dissociation are key processes that drive the gas lifecycle in galaxies. Using the SImulating the LifeCycle of Molecular Clouds (SILCC) zoom-in simulation suite, we explore the utility of future observations of H$_2$ dissociation and formation for tracking the lifecycle of molecular clouds. The simulations used in this work include non-equilibrium H$_2$ for… ▽ More Molecular hydrogen (H$_2$) formation and dissociation are key processes that drive the gas lifecycle in galaxies. Using the SImulating the LifeCycle of Molecular Clouds (SILCC) zoom-in simulation suite, we explore the utility of future observations of H$_2$ dissociation and formation for tracking the lifecycle of molecular clouds. The simulations used in this work include non-equilibrium H$_2$ formation, stellar radiation, sink particles, and turbulence. We find that, at early times in the cloud evolution, H$_2$ formation rapidly outpaces dissociation and molecular clouds build their mass from the atomic reservoir in their environment. Rapid H$_2$ formation is also associated with a higher early star formation rate. For the clouds studied here, H$_2$ is strongly out of chemical equilibrium during the early stages of cloud formation but settles into a bursty chemical steady-state about 2 Myrs after the first stars form. At the latest stage of cloud evolution, dissociation outweighs formation and the clouds enter a dispersal phase. We discuss how theories for the molecular cloud lifecycle and the star formation efficiency may be distinguished with observational measurements of H$_2$ fluorescence with a space-based high-resolution FUV spectrometer, such as the proposed Hyperion and Eos NASA Explorer missions. Such missions would enable measurements of the H$_2$ dissociation and formation rates, which we demonstrate can be connected to different phases in a molecular cloud's star-forming life, including cloud building, rapidly star-forming, H$_2$ chemical equilibrium, and cloud destruction. △ Less

Submitted 5 February, 2024; v1 submitted 2 February, 2024; originally announced February 2024.

Comments: Submitted to ApJ, comments welcome

ALMA-LEGUS II: The Influence of Sub-Galactic Environment on Molecular Cloud Properties

Authors: Molly K. Finn, Kelsey E. Johnson, Remy Indebetouw, Allison H. Costa, Angela Adamo, Alessandra Aloisi, Lauren Bittle, Daniela Calzetti, Daniel A. Dale, Clare L. Dobbs, Jennifer Donovan Meyer, Bruce G. Elmegreen, Debra M. Elmegreen, Michele Fumagalli, J. S. Gallagher, Kathryn Grasha, Eva K. Grebel, Robert C. Kennicutt, Mark R. Krumholz, Janice C. Lee, Matteo Messa, Preethi Nair, Elena Sabbi, Linda J. Smith, David A. Thilker , et al. (2 additional authors not shown)

Abstract: We compare the molecular cloud properties in sub-galactic regions of two galaxies, barred spiral NGC 1313, which is forming many massive clusters, and flocculent spiral NGC 7793, which is forming significantly fewer massive clusters despite having a similar star formation rate to NGC 1313. We find that there are larger variations in cloud properties between different regions within each galaxy tha… ▽ More We compare the molecular cloud properties in sub-galactic regions of two galaxies, barred spiral NGC 1313, which is forming many massive clusters, and flocculent spiral NGC 7793, which is forming significantly fewer massive clusters despite having a similar star formation rate to NGC 1313. We find that there are larger variations in cloud properties between different regions within each galaxy than there are between the galaxies on a global scale, especially for NGC 1313. There are higher masses, linewidths, pressures, and virial parameters in the arms of NGC 1313 and center of NGC 7793 than in the interarm and outer regions of the galaxies. The massive cluster formation of NGC 1313 may be driven by its greater variation in environments, allowing more clouds with the necessary conditions to arise, although no one parameter seems primarily responsible for the difference in star formation. Meanwhile NGC 7793 has clouds that are as massive and have as much kinetic energy as clouds in the arms of NGC 1313, but have densities and pressures more similar to the interarm regions and so are less inclined to collapse and form stars. The cloud properties in NGC 1313 and NGC 7793 suggest that spiral arms, bars, interarm regions, and flocculent spirals each represent distinct environments with regard to molecular cloud populations. We see surprisingly little difference in surface densities between the regions, suggesting that the differences in surface densities frequently seen between arm and interarm regions of lower-resolution studies are indicative of the sparsity of molecular clouds, rather than differences in their true surface density. △ Less

Submitted 2 January, 2024; originally announced January 2024.

Comments: 24 pages, 12 figures, accepted for publication in ApJ

ALMA-LEGUS I: The Influence of Galaxy Morphology on Molecular Cloud Properties

Authors: Molly K. Finn, Kelsey E. Johnson, Remy Indebetouw, Allison H. Costa, Angela Adamo, Alessandra Aloisi, Lauren Bittle, Daniela Calzetti, Daniel A. Dale, Clare L. Dobbs, Jennifer Donovan Meyer, Bruce G. Elmegreen, Debra M. Elmegreen, Michele Fumagalli, J. S. Gallagher, Kathryn Grasha, Eva K. Grebel, Robert C. Kennicutt, Mark R. Krumholz, Janice C. Lee, Matteo Messa, Preethi Nair, Elena Sabbi, Linda J. Smith, David A. Thilker , et al. (2 additional authors not shown)

Abstract: We present a comparative study of the molecular gas in two galaxies from the LEGUS sample: barred spiral NGC 1313 and flocculent spiral NGC 7793. These two galaxies have similar masses, metallicities, and star formation rates, but NGC 1313 is forming significantly more massive star clusters than NGC 7793, especially young massive clusters (<10 Myr, >10^4 Msol). Using ALMA CO(2-1) observations of t… ▽ More We present a comparative study of the molecular gas in two galaxies from the LEGUS sample: barred spiral NGC 1313 and flocculent spiral NGC 7793. These two galaxies have similar masses, metallicities, and star formation rates, but NGC 1313 is forming significantly more massive star clusters than NGC 7793, especially young massive clusters (<10 Myr, >10^4 Msol). Using ALMA CO(2-1) observations of the two galaxies with the same sensitivities and resolutions of 13 pc, we directly compare the molecular gas in these two similar galaxies to determine the physical conditions responsible for their large disparity in cluster formation. By fitting size-linewidth relations for the clouds in each galaxy, we find that NGC 1313 has a higher intercept than NGC 7793, implying that its clouds have higher kinetic energies at a given size scale. NGC 1313 also has more clouds near virial equilibrium than NGC 7793, which may be connected to its higher rate of massive cluster formation. However, these virially bound clouds do not show a stronger correlation with young clusters than that of the general cloud population. We find surprisingly small differences between the distributions of molecular cloud populations in the two galaxies, though the largest of those differences are that NGC 1313 has higher surface densities and lower free-fall times. △ Less

Submitted 2 January, 2024; originally announced January 2024.

Comments: 25 pages, 10 figures, accepted for publication in ApJ

The ALMaQUEST Survey XI: A strong but non-linear relationship between star formation and dynamical equilibrium pressure

Authors: Sara L. Ellison, Hsi-An Pan, Asa F. L. Bluck, Mark R. Krumholz, Lihwai Lin, Leslie Hunt, Edvige Corbelli, Mallory D. Thorp, Jorge Barrera-Ballesteros, Sebastian F. Sanchez, Jillian M. Scudder, Salvatore Quai

Abstract: We present the extended ALMA MaNGA QUEnching and STar formation survey, a combination of the original 46 ALMaQUEST galaxies plus new ALMA observations for a further 20 interacting galaxies. Three well-studied scaling relations are fit to the 19,999 star-forming spaxels in the extended sample, namely the resolved Schmidt-Kennicutt (rSK) relation, the resolved star forming main sequence (rSFMS) and… ▽ More We present the extended ALMA MaNGA QUEnching and STar formation survey, a combination of the original 46 ALMaQUEST galaxies plus new ALMA observations for a further 20 interacting galaxies. Three well-studied scaling relations are fit to the 19,999 star-forming spaxels in the extended sample, namely the resolved Schmidt-Kennicutt (rSK) relation, the resolved star forming main sequence (rSFMS) and the resolved molecular gas main sequence (rMGMS). We additionally investigate the relationship between the dynamical equilibrium pressure (PDE) and star formation rate surface density (Sigma_SFR), which we refer to as the resolved PDE (rPDE) relation. Contrary to previous studies that have focussed on normal star-forming galaxies and found an approximately linear rPDE relation, the presence of more vigourously star-forming galaxies in the extended ALMaQUEST sample reveals a marked turnover in the relation at high pressures. Although the scatter around the linear fit to the rPDE relation is similar to the other three relations, a random forest analysis, which can extract non-linear dependences, finds that PDE is unambiguously more important than either Sigma_H2 or Sigma_star for predicting Sigma_SFR. We compare the observed rPDE relation to the prediction of the pressure-regulated feedback-modulated (PRFM) model of star formation, finding that galaxies residing on the global SFMS do indeed closely follow the rPDE relation predicted by the PRFM theory. However, galaxies above and below the global SFMS show significant deviations from the model. Galaxies with high SFR are instead consistent with models that include other contributions to turbulence in addition to the local star formation feedback. △ Less

Submitted 5 December, 2023; originally announced December 2023.

Comments: Accepted for publication in MNRAS

Fundamental MHD scales -- II: the kinematic phase of the supersonic small-scale dynamo

Authors: Neco Kriel, James R. Beattie, Christoph Federrath, Mark R. Krumholz, Justin Kin Jun Hew

Abstract: The small-scale dynamo (SSD) amplifies weak magnetic fields exponentially fast via kinetic motions. While there exist well-established theories for SSDs in incompressible flows, many astrophysical SSDs operate in supersonic turbulence. To understand the impact of compressibility on amplified magnetic fields, we perform an extensive set of visco-resistive SSD simulations, covering a wide range of s… ▽ More The small-scale dynamo (SSD) amplifies weak magnetic fields exponentially fast via kinetic motions. While there exist well-established theories for SSDs in incompressible flows, many astrophysical SSDs operate in supersonic turbulence. To understand the impact of compressibility on amplified magnetic fields, we perform an extensive set of visco-resistive SSD simulations, covering a wide range of sonic Mach number $\mathcal{M}$, hydrodynamic Reynolds number Re, and magnetic Prandtl number Pm. We develop robust methods for measuring kinetic and magnetic energy dissipation scales $\ell_ν$ and $\ell_η$, as well as the scale at which magnetic fields are strongest $\ell_p$ during the kinematic phase of these simulations. We show that $\ell_ν/\ell_η\sim$ Pm$^{1/2}$ is a universal feature in the kinematic phase of Pm $\geq 1$ SSDs, regardless of $\mathcal{M}$ or Re, and we confirm earlier predictions that SSDs operating in incompressible plasmas (either $\mathcal{M} \leq 1$ or Re $<$ Re$_{\rm crit} \approx 100$) concentrate magnetic energy at the smallest scales allowed by magnetic dissipation, $\ell_p \sim \ell_η$, and produce fields organised with field strength and field-line curvature inversely correlated. However, we show that these predictions fail for compressible SSDs ($\mathcal{M} > 1$ and Re $>$ Re$_{\rm crit}$), where shocks concentrate magnetic energy in large-scale, over-dense, coherent structures, with size $\ell_p \sim (\ell_{\rm turb} / \ell_{\rm shock})^{1/3} \ell_η\gg \ell_η$, where $\ell_{\rm shock} \sim \mathcal{M}^2 / [$Re $ (\mathcal{M} - 1)^2]$ is shock width, and $\ell_{\rm turb}$ is the turbulent outer scale; magnetic field-line curvature becomes almost independent of the field strength. We discuss the implications for galaxy mergers and for cosmic-ray transport models in the interstellar medium that are sensitive to field-line curvature statistics. △ Less

Submitted 25 October, 2023; originally announced October 2023.

Comments: 25 pages, 15 figures, submitted to MNRAS, json-file w/ data

CONGRuENTS (COsmic-ray, Neutrino, Gamma-ray and Radio Non-Thermal Spectra). II. Population-level correlations between galactic infrared, radio, and γ-ray emission

Authors: Matt A. Roth, Mark R. Krumholz, Roland M. Crocker, Todd A. Thompson

Abstract: Galaxies obey a number of empirical correlations between their radio, γ-ray, and infrared emission, but the physical origins of these correlations remain uncertain. Here we use the CONGRuENTS model for broadband non-thermal emission from star-forming galaxies, which self-consistently calculates energy-dependent transport and non-thermal emission from cosmic ray hadrons and leptons, to predict radi… ▽ More Galaxies obey a number of empirical correlations between their radio, γ-ray, and infrared emission, but the physical origins of these correlations remain uncertain. Here we use the CONGRuENTS model for broadband non-thermal emission from star-forming galaxies, which self-consistently calculates energy-dependent transport and non-thermal emission from cosmic ray hadrons and leptons, to predict radio and γ-ray emission for a synthetic galaxy population with properties drawn from a large deep-field survey. We show that our synthetic galaxies reproduce observed relations such as the FIR-radio correlation, the FIR-γ correlation, and the distribution of radio spectral indices, and we use the model to explain the physical origins of these relations. Our results show that the FIR-radio correlation arises because the amount of cosmic ray electron power ultimately radiated as synchrotron emission varies only weakly with galaxy star formation rate as a result of the constraints imposed on gas properties by hydrostatic balance and turbulent dynamo action; the same physics dictates the extent of proton calorimetry in different galaxies, and thus sets the FIR-γ-ray correlation. We further show that galactic radio spectral indices result primarily from competition between thermal free-free emission and energy-dependent loss of cosmic ray electrons to bremsstrahlung and escape into galactic halos, with shaping of the spectrum by inverse Compton, synchrotron, and ionisation processes typically playing a sub-dominant role. In addition to explaining existing observations, we use our analysis to predict a heretofore unseen correlation between the curvature of galaxies' radio spectra and their pion-driven γ-ray emission, a prediction that will be testable with upcoming facilities. △ Less

Submitted 15 June, 2024; v1 submitted 9 October, 2023; originally announced October 2023.

Comments: 17 pages, 14 figures

Journal ref: MNRAS, Volume 530, Issue 2, May 2024, Pages 1849-1865

Quokka-based Understanding of Outflows (QED) I. Metal loading, phase structure, and convergence testing for Solar neighbourhood conditions

Authors: Aditi Vijayan, Mark R. Krumholz, Benjamin D. Wibking

Abstract: Multiphase galactic outflows, generated by supernova feedback, are likely to be more metal-rich than the interstellar media from which they are driven due to incomplete mixing between supernova ejecta and the ambient ISM. This enrichment is important for shaping galactic metallicities and metallicity gradients, but measuring it quantitatively from simulations requires resolution high enough to res… ▽ More Multiphase galactic outflows, generated by supernova feedback, are likely to be more metal-rich than the interstellar media from which they are driven due to incomplete mixing between supernova ejecta and the ambient ISM. This enrichment is important for shaping galactic metallicities and metallicity gradients, but measuring it quantitatively from simulations requires resolution high enough to resolve mass, momentum and energy exchanges between the different phases of the outflows. In this context, we present simulations of outflows, driven by SN feedback, conducted using \textsc{Quokka}, a new GPU-optimised AMR radiation-hydrodynamics code. This code allows us to reach combinations of resolution, simulation volume, and simulation duration larger than those that have previously been possible, and to resolve all gas phases from cold neutral medium, $T \sim 100$ K, to hot ionised gas, $T \gtrsim 10^7$ K. In this, a first of a series of papers exploring generation and evolution of multiphase outflows from a wide range of galactic environments and star formation rates, we quantify the extent of selective metal loading in Solar neighbourhood-like environments. We explain the selective metal loading we find as a result of the transport of metals within and between phases, a phenomenon we can study owing to the parsec-scale resolution that our simulations achieve. We also quantify the sensitivity of metal loading studies to numerical resolution, and present convergence criteria for future studies. △ Less

Submitted 14 September, 2023; originally announced September 2023.

Comments: Submitted to MNRAS. Full title="Quokka-based Understanding of Outflows Derived from Extensive, Repeated, Accurate, Thorough, Demanding, Expensive, Memory-consuming, Ongoing Numerical Simulations of Transport, Removal, Accretion, Nucleosynthesis, Deposition, and Uplifting of Metals (QUOD ERAT DEMONSTRANDUM, i.e, QED)"

Universal Upper End of the Stellar Initial Mass Function in the Young and Compact LEGUS clusters

Authors: Dooseok Escher Jung, Daniela Calzetti, Matteo Messa, Mark Heyer, Mattia Sirressi, Sean T. Linden, Angela Adamo, Rupali Chandar, Michele Cignoni, David O. Cook, Clare L. Dobbs, Bruce G. Elmegreen, Aaron S. Evans, Michele Fumagalli, John S. Gallagher III, Deidre A. Hunter, Kelsey E. Johnson, Robert C. Kennicutt Jr., Mark R. Krumholz, Daniel Schaerer, Elena Sabbi, Linda J. Smith, Monica Tosi, Aida Wofford

Abstract: We investigate the variation in the upper end of stellar initial mass function (uIMF) in 375 young and compact star clusters in five nearby galaxies within $\sim 5$ Mpc. All the young stellar clusters (YSCs) in the sample have ages $\lesssim 4$ Myr and masses above 500 $M_{\odot}$, according to standard stellar models. The YSC catalogs were produced from Hubble Space Telescope images obtained as p… ▽ More We investigate the variation in the upper end of stellar initial mass function (uIMF) in 375 young and compact star clusters in five nearby galaxies within $\sim 5$ Mpc. All the young stellar clusters (YSCs) in the sample have ages $\lesssim 4$ Myr and masses above 500 $M_{\odot}$, according to standard stellar models. The YSC catalogs were produced from Hubble Space Telescope images obtained as part of the Legacy ExtraGalactic UV Survey (LEGUS) Hubble treasury program. They are used here to test whether the uIMF is universal or changes as a function of the cluster's stellar mass. We perform this test by measuring the H$α$ luminosity of the star clusters as a proxy for their ionizing photon rate, and charting its trend as a function of cluster mass. Large cluster numbers allow us to mitigate the stochastic sampling of the uIMF. The advantage of our approach relative to previous similar attempts is the use of cluster catalogs that have been selected independently of the presence of H$α$ emission, thus removing a potential sample bias. We find that the uIMF, as traced by the H$α$ emission, shows no dependence on cluster mass, suggesting that the maximum stellar mass that can be produced in star clusters is universal, in agreement with previous findings. △ Less

Submitted 28 July, 2023; originally announced July 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

The metallicity dependence of the stellar initial mass function

Authors: Tabassum S. Tanvir, Mark R. Krumholz

Abstract: Dust is important for star formation because it is the crucial component that couples gas to stellar radiation fields, allowing radiation feedback to influence gas fragmentation and thus the stellar initial mass function (IMF). Variations in dust abundance therefore provide a potential avenue by which variation in galaxy metallicity might affect the IMF. In this paper we present a series of radiat… ▽ More Dust is important for star formation because it is the crucial component that couples gas to stellar radiation fields, allowing radiation feedback to influence gas fragmentation and thus the stellar initial mass function (IMF). Variations in dust abundance therefore provide a potential avenue by which variation in galaxy metallicity might affect the IMF. In this paper we present a series of radiation-magnetohydrodynamic simulations in which we vary the metallicity and thus the dust abundance from 1% of Solar to 3$\times$ Solar, spanning the range from the lowest metallicity dwarfs to the most metal-rich early-type galaxies found in the local Universe. We design the simulations to keep all dimensionless parameters constant so that the interaction between feedback and star-forming environments of varying surface density and metallicity is the only factor capable of breaking the symmetry between the simulations and modifying the IMF, allowing us to cleanly isolate and understand the effects of each environmental parameter. We find that at a fixed surface density more metal-rich clouds tend to form a slightly more bottom-heavy IMF than metal-poor ones, primarily because in metal-poor gas radiation feedback is able to propagate further, heating somewhat larger volumes of gas. However, shifts in IMF with metallicity at a fixed surface density are much smaller than shifts with surface density at fixed metallicity; metallicity-induced IMF variations are too small to explain the variations in mass-to-light ratio reported in galaxies of different mass and metallicity. We, therefore, conclude that metallicity variations are much less important than variations in surface density in driving changes in the IMF and that the latter rather than the former are most likely responsible for the IMF variations found in early-type galaxies. △ Less

Submitted 31 May, 2023; originally announced May 2023.

Comments: 10 pages, 12 figures, submitted in MNRAS

Direct observations of the atomic-molecular phase transition in the Milky Way's nuclear wind

Authors: Karlie A. Noon, Mark R. Krumholz, Enrico M. Di Teodoro, Naomi M. McClure-Griffiths, Felix J. Lockman, Lucia Armillotta

Abstract: Hundreds of high-velocity atomic gas clouds exist above and below the Galactic Centre, with some containing a molecular component. However, the origin of these clouds in the Milky Way's wind is unclear. This paper presents new high-resolution MeerKAT observations of three atomic gas clouds and studies the relationship between the atomic and molecular phases at $\sim 1$ pc scales. The clouds' atomi… ▽ More Hundreds of high-velocity atomic gas clouds exist above and below the Galactic Centre, with some containing a molecular component. However, the origin of these clouds in the Milky Way's wind is unclear. This paper presents new high-resolution MeerKAT observations of three atomic gas clouds and studies the relationship between the atomic and molecular phases at $\sim 1$ pc scales. The clouds' atomic hydrogen column densities, $N_{\mathrm{HI}}$, are less than a $\mbox{few}\times 10^{20}$ cm$^{-2}$, but the two clouds closest to the Galactic Centre nonetheless have detectable CO emission. This implies the presence of H$_{2}$ at levels of $N_{\mathrm{HI}}$ at least a factor of ten lower than in the typical Galactic interstellar medium. For the cloud closest to the Galactic Centre, there is little correlation between the $N_{\mathrm{HI}}$ and the probability that it will harbour detectable CO emissions. In contrast, for the intermediate cloud, detectable CO is heavily biased toward the highest values of $N_{\mathrm{HI}}$. The cloud most distant from the Galactic Centre has no detectable CO at similar $N_{\mathrm{HI}}$ values. Moreover, we find that the two clouds with detectable CO are too molecule-rich to be in chemical equilibrium, given the depths of their atomic shielding layers, which suggests a scenario whereby these clouds consist of pre-existing molecular gas from the disc that the Galactic wind has swept up, and that is dissociating into atomic hydrogen as it flows away from the Galaxy. We estimate that entrained molecular material of this type has a $\sim \mathrm{few}-10$ Myr lifetime before photodissociating. △ Less

Submitted 13 July, 2023; v1 submitted 13 April, 2023; originally announced April 2023.

Comments: 11 pages, 6 figures, 2 tables. Accepted for publication in MNRAS

Journal ref: Volume 524, Issue 1, 2023, Pages 1258 - 1268

The interplay between feedback, accretion, transport and winds in setting gas-phase metal distribution in galaxies

Authors: Piyush Sharda, Omri Ginzburg, Mark R. Krumholz, John C. Forbes, Emily Wisnioski, Matilde Mingozzi, Henry R. M. Zovaro, Avishai Dekel

Abstract: The recent decade has seen an exponential growth in spatially-resolved metallicity measurements in the interstellar medium (ISM) of galaxies. To first order, these measurements are characterised by the slope of the radial metallicity profile, known as the metallicity gradient. In this work, we model the relative role of star formation feedback, gas transport, cosmic gas accretion, and galactic win… ▽ More The recent decade has seen an exponential growth in spatially-resolved metallicity measurements in the interstellar medium (ISM) of galaxies. To first order, these measurements are characterised by the slope of the radial metallicity profile, known as the metallicity gradient. In this work, we model the relative role of star formation feedback, gas transport, cosmic gas accretion, and galactic winds in driving radial metallicity profiles and setting the mass-metallicity gradient relation (MZGR). We include a comprehensive treatment of these processes by including them as sources that supply mass, metals, and energy to marginally unstable galactic discs in pressure and energy balance. We show that both feedback and accretion that can drive turbulence and enhance metal-mixing via diffusion are crucial to reproduce the observed MZGR in local galaxies. Metal transport also contributes to setting metallicity profiles, but it is sensitive to the strength of radial gas flows in galaxies. While the mass loading of galactic winds is important to reproduce the mass metallicity relation (MZR), we find that metal mass loading is more important to reproducing the MZGR. Specifically, our model predicts preferential metal enrichment of galactic winds in low-mass galaxies. This conclusion is robust against our adopted scaling of the wind mass-loading factor, uncertainties in measured wind metallicities, and systematics due to metallicity calibrations. Overall, we find that at $z \sim 0$, galactic winds and metal transport are more important in setting metallicity gradients in low-mass galaxies whereas star formation feedback and gas accretion dominate setting metallicity gradients in massive galaxies. △ Less

Submitted 11 January, 2024; v1 submitted 28 March, 2023; originally announced March 2023.

Comments: 21 pages (+ appendix), 11 figures, 4 tables. Accepted by MNRAS

The sub-critical illusion: synthetic Zeeman effect observations from galactic zoom-in simulations

Authors: Zipeng Hu, Benjamin D. Wibking, Mark R. Krumholz

Abstract: Mass-to-flux ratios measured via the Zeeman effect suggest the existence of a transition from a magnetically sub-critical state in HI clouds to a super-critical state in molecular clouds. However, due to projection, chemical, and excitation effects, Zeeman measurements are subject to a number of biases, and may not reflect the true relations between gravitational and magnetic energies. In this pap… ▽ More Mass-to-flux ratios measured via the Zeeman effect suggest the existence of a transition from a magnetically sub-critical state in HI clouds to a super-critical state in molecular clouds. However, due to projection, chemical, and excitation effects, Zeeman measurements are subject to a number of biases, and may not reflect the true relations between gravitational and magnetic energies. In this paper, we carry out simulations of the formation of magnetised molecular clouds, zooming in from an entire galaxy to sub-pc scales, which we post-process to produce synthetic HI and OH Zeeman measurements. The mass-to-flux ratios we recover from the simulated observations show a transition in magnetic criticality that closely matches observations, but we find that the gravitational-magnetic energy ratios on corresponding scales are mostly super-critical, even in the HI regime. We conclude that HI clouds in the process of assembling to form molecular clouds are already super-critical even before H_2 forms, and that the apparent transition from sub- to super-criticality between HI and H_2 is primarily an illusion created by chemical and excitation biases affecting the Zeeman measurements. △ Less

Submitted 22 February, 2023; originally announced February 2023.

Comments: 12 pages, 10 figures, submitted to MNRAS

Dust Buried Compact Sources in the Dwarf Galaxy NGC 4449

Authors: Daniela Calzetti, Sean T. Linden, Timothy McQuaid, Matteo Messa, Zhiyuan Ji, Mark R. Krumholz, Angela Adamo, Bruce Elmegreen, Kathryn Grasha, Kelsey E. Johnson, Elena Sabbi, Linda Smith, Varun Bajaj

Abstract: Multi-wavelength images from the Hubble Space Telescope covering the wavelength range 0.27-1.6 $μ$m show that the central area of the nearby dwarf galaxy NGC4449 contains several tens of compact sources that are emitting in the hydrogen recombination line Pa$β$ (1.2818 $μ$m) but are only marginally detected in H$α$ (0.6563 $μ$m) and undetected at wavelengths $λ\le$0.55 $μ$m. An analysis of the spe… ▽ More Multi-wavelength images from the Hubble Space Telescope covering the wavelength range 0.27-1.6 $μ$m show that the central area of the nearby dwarf galaxy NGC4449 contains several tens of compact sources that are emitting in the hydrogen recombination line Pa$β$ (1.2818 $μ$m) but are only marginally detected in H$α$ (0.6563 $μ$m) and undetected at wavelengths $λ\le$0.55 $μ$m. An analysis of the spectral energy distributions (SEDs) of these sources indicates that they are likely relatively young star clusters heavily attenuated by dust. The selection function used to identify the sources prevents meaningful statistical analyses of their age, mass, and dust extinction distributions. However, these cluster candidates have ages $\sim$5-6 Myr and A$_V>$6 mag, according to their SED fits, and are extremely compact, with typical deconvolved radii of 1 pc. The dusty clusters are located at the periphery of dark clouds within the galaxy and appear to be partially embedded. Density and pressure considerations indicate that the HII regions surrounding these clusters may be stalled, and that pre-supernova feedback has not been able to clear the clusters of their natal cocoons. These findings are in potential tension with existing models that regulate star formation with pre-supernova feedback, since pre-supernova feedback acts on short timescales, $\lesssim$4 Myr, for a standard Stellar Initial Mass function. The existence of a population of dusty star clusters with ages $>$4 Myr, if confirmed by future observations, paints a more complex picture for the role of stellar feedback in controlling star formation. △ Less

Submitted 21 February, 2023; originally announced February 2023.

Comments: 27 pages, 22 figures; accepted for publication in the Astrophysical Journal

Spatially resolved dust properties and quasar-galaxy decomposition of a HyLIRG at z = 4.4

Authors: Takafumi Tsukui, Emily Wisnioski, Mark R. Krumholz, Andrew Battisti

Abstract: We report spatially resolved dust properties of the quasar host galaxy BRI 1335-0417 at redshift $z = 4.4$ constrained by the ALMA observations. The dust temperature map, derived from a greybody fit to rest-frame 90 and 161~$μ$m continuum images, shows a steep increase towards the centre, reaching $57.1 \pm 0.3$ K and a flat median profile at the outer regions of $\sim$38 K. Image decomposition an… ▽ More We report spatially resolved dust properties of the quasar host galaxy BRI 1335-0417 at redshift $z = 4.4$ constrained by the ALMA observations. The dust temperature map, derived from a greybody fit to rest-frame 90 and 161~$μ$m continuum images, shows a steep increase towards the centre, reaching $57.1 \pm 0.3$ K and a flat median profile at the outer regions of $\sim$38 K. Image decomposition analysis reveals the presence of a point source in both dust continuum images spatially coincident with the highest temperature peak and the optical quasar position, which we attribute to warm dust heated by an active galactic nucleus (AGN). We show that a model including this warm component along with cooler dust heated by star formation describes the global SED better than a single component model, with dust temperatures of 87.1$^{+34.1}_{-18.3}$ K (warm component) and 52.6$^{+10.3}_{-11.0}$ K (cold component). The star formation rate (SFR) estimated from the cold dust component is $1700_{-400}^{+500} M_\odot$ yr$^{-1}$, a factor of three smaller than previous estimates due to a large AGN contribution ($53^{+14}_{-15}$\%). The unresolved warm dust component also explains the steep temperature gradient, as the temperature profile derived after the point source subtraction is flat. The point source subtraction also reduces the estimated central SFR surface density $Σ_{\mathrm{SFR}}$ by over a factor of three. With this correction, spatially resolved measurements of $Σ_{\mathrm{SFR}}$ and the surface gas mass density $Σ_{\mathrm{gas}}$ form a roughly linear sequence in the Kennicutt-Schmidt diagram with a constant gas depletion time of 50-200 Myr. The demonstrated AGN-host galaxy decomposition reveals the importance of spatially resolved data for accurate measurements of quasar host galaxy properties, including dust temperature, star-formation rates, and size. △ Less

Submitted 22 June, 2023; v1 submitted 14 February, 2023; originally announced February 2023.

Comments: 26pages, 27figures, 4tables, Published in MNRAS; comments are warmly welcomed

Clouds of Theseus: long-lived molecular clouds are composed of short-lived H2 molecules

Authors: Sarah M. R. Jeffreson, Vadim A. Semenov, Mark R. Krumholz

Abstract: We use passive gas tracer particles in an Arepo simulation of a dwarf spiral galaxy to relate the Lagrangian evolution of star-forming gas parcels and their H2 molecules to the evolution of their host giant molecular clouds. We find that the median chemical lifetime of H2 is just 4 Myr, independent of the lifetime of its host molecular cloud, which may vary from 1 to 90 Myr, with a substantial por… ▽ More We use passive gas tracer particles in an Arepo simulation of a dwarf spiral galaxy to relate the Lagrangian evolution of star-forming gas parcels and their H2 molecules to the evolution of their host giant molecular clouds. We find that the median chemical lifetime of H2 is just 4 Myr, independent of the lifetime of its host molecular cloud, which may vary from 1 to 90 Myr, with a substantial portion of all star formation in the galaxy occurring in relatively long-lived clouds. The rapid ejection of gas from around young massive stars by early stellar feedback is responsible for this short H2 survival time, driving down the density of the surrounding gas, so that its H2 molecules are dissociated by the interstellar radiation field. This ejection of gas from the H2-dominated state is balanced by the constant accretion of new gas from the galactic environment, constituting a "competition model" for molecular cloud evolution. Gas ejection occurs at a rate that is proportional to the molecular cloud mass, so that the cloud lifetime is determined by the accretion rate, which may be as high as 4 x 10^4 Msol/Myr in the longest-lived clouds. Our findings therefore resolve the conflict between observations of rapid gas ejection around young massive stars and observations of long-lived molecular clouds in galaxies, that often survive up to several tens of Myr. We show that the fastest-accreting, longest-lived, highest-mass clouds drive supernova clustering on sub-cloud scales, which in turn is a key driver of galactic outflows. △ Less

Submitted 24 January, 2023; originally announced January 2023.

Comments: 16 pages, 16 figures. Submitted to MNRAS, comments welcome!

Cluster Population Demographics in NGC 628 Derived from Stochastic Population Synthesis Models

Authors: Jianling Tang, Kathryn Grasha, Mark R. Krumholz

Abstract: The physical properties of star cluster populations offer valuable insights into their birth, evolution, and disruption. However, individual stars in clusters beyond the nearest neighbours of the Milky Way are unresolved, forcing analyses of star cluster demographics to rely on integrated light, a process fraught with uncertainty. Here we infer the demographics of the cluster population in the ben… ▽ More The physical properties of star cluster populations offer valuable insights into their birth, evolution, and disruption. However, individual stars in clusters beyond the nearest neighbours of the Milky Way are unresolved, forcing analyses of star cluster demographics to rely on integrated light, a process fraught with uncertainty. Here we infer the demographics of the cluster population in the benchmark galaxy NGC 628 using data from the Legacy Extra-galactic UV Survey (LEGUS) coupled to a novel Bayesian forward-modelling technique. Our method analyzes all 1178 clusters in the LEGUS catalogue, $\sim4$ times more than prior studies severely affected by completeness cuts. Our results indicate that the cluster mass function is either significantly steeper than the commonly-observed slope of $-2$ or is truncated at $\approx 10^{4.5}$ M$_\odot$; the latter possibility is consistent with proposed relations between truncation mass and star formation surface density. We find that cluster disruption is relatively mild for the first $\approx 200$ Myr of cluster evolution; no evidence for mass-dependent disruption is found. We find suggestive but not incontrovertible evidence that inner galaxy clusters may be more prone to disruption and outer galaxy clusters have a more truncated mass function, but confirming or refuting these findings will require larger samples from future observations of outer galaxy clusters. Finally, we find that current stellar track and atmosphere models, along with common forms for cluster mass and age distributions, cannot fully capture all features in the multidimensional photometric distribution of star clusters. While our forward-modelling approach outperforms earlier backward-modeling approaches, some systematic differences persist between observed and modelled photometric distributions. △ Less

Submitted 22 July, 2024; v1 submitted 14 January, 2023; originally announced January 2023.

Comments: 20 pages, 20 figures, 3 tables, accepted for publication in MNRAS on July 22nd 2024

CONGRuENTS (COsmic-ray, Neutrino, Gamma-ray and Radio Non-Thermal Spectra). I. A predictive model for galactic non-thermal emission

Authors: Matt A. Roth, Mark R. Krumholz, Roland M. Crocker, Todd A. Thompson

Abstract: The total luminosity and spectral shape of the non-thermal emission produced by cosmic rays depends on their interstellar environment, a dependence that gives rise to correlations between galaxies' bulk properties -- star formation rate, stellar mass, and others -- and their non-thermal spectra. Understanding the physical mechanisms of cosmic ray transport, loss, and emission is key to understandi… ▽ More The total luminosity and spectral shape of the non-thermal emission produced by cosmic rays depends on their interstellar environment, a dependence that gives rise to correlations between galaxies' bulk properties -- star formation rate, stellar mass, and others -- and their non-thermal spectra. Understanding the physical mechanisms of cosmic ray transport, loss, and emission is key to understanding these correlations. Here, in the first paper of the series, we present a new method to compute the non-thermal spectra of star-forming galaxies, and describe an open-source software package -- COsmic-ray, Neutrino, Gamma-ray and Radio Non-Thermal Spectra (CONGRuENTS) -- that implements it. As a crucial innovation, our method requires as input only a galaxy's effective radius, star formation rate, stellar mass, and redshift, all quantities that are readily available for large samples of galaxies and do not require expensive, spatially resolved gas measurements. From these inputs we derive individual, galaxy-by-galaxy models for the background gas and radiation field through which cosmic rays propagate, from which we compute steady state cosmic ray spectra for hadronic and leptonic particles in both the galactic disc and halo by solving the full kinetic equation. We invoke modern models for cosmic ray transport and include all significant emission and loss mechanisms. In this paper we describe the model and validate it against non-thermal emission measured in nearby star-forming galaxies that span four orders of magnitude in star formation rate. △ Less

Submitted 16 May, 2023; v1 submitted 19 December, 2022; originally announced December 2022.

Comments: 23 pages, 14 figures, 1 table, accepted for publication in MNRAS

Can pre-supernova winds from massive stars enrich the interstellar medium with nitrogen at high redshift?

Authors: Arpita Roy, Mark R. Krumholz, Michael A. Dopita, Ralph S. Sutherland, Lisa J. Kewley, Alexander Heger

Abstract: Understanding the nucleosynthetic origin of nitrogen and the evolution of the N/O ratio in the interstellar medium is crucial for a comprehensive picture of galaxy chemical evolution at high-redshift because most observational metallicity (O/H) estimates are implicitly dependent on the N/O ratio. The observed N/O at high-redshift shows an overall constancy with O/H, albeit with a large scatter. We… ▽ More Understanding the nucleosynthetic origin of nitrogen and the evolution of the N/O ratio in the interstellar medium is crucial for a comprehensive picture of galaxy chemical evolution at high-redshift because most observational metallicity (O/H) estimates are implicitly dependent on the N/O ratio. The observed N/O at high-redshift shows an overall constancy with O/H, albeit with a large scatter. We show that these heretofore unexplained features can be explained by the pre-supernova wind yields from rotating massive stars (M$\gtrsim 10 \, \mathrm{M}_\odot$, $v/v_{\rm{crit}} \gtrsim 0.4$). Our models naturally produce the observed N/O plateau, as well as the scatter at low O/H. We find the scatter to arise from varying star formation efficiency. However, the models that have supernovae dominated yields produce a poor fit to the observed N/O at low O/H. This peculiar abundance pattern at low O/H suggests that dwarf galaxies are most likely to be devoid of SNe yields and are primarily enriched by pre-supernova wind abundances. △ Less

Submitted 9 December, 2022; originally announced December 2022.

Comments: 6 pages, 1 figure, published in Proceedings of the International Astronomical Union No. 366, 2022 titled "The Origin of Outflows in Evolved Stars"

Journal ref: Proceedings of the International Astronomical Union, 16(S366), 33-38, 2022

The impact of carbon and oxygen abundances on the metal-poor initial mass function

Authors: Piyush Sharda, Anish M. Amarsi, Kathryn Grasha, Mark R. Krumholz, David Yong, Gen Chiaki, Arpita Roy, Thomas Nordlander

Abstract: Star formation models predict that the metal-poor initial mass function (IMF) can be substantially different from that observed in the metal-rich Milky Way. This changeover occurs because metal-poor gas clouds cool inefficiently due to their lower abundance of metals and dust. However, predictions for the metal-poor IMF to date rely on assuming Solar-scaled abundances, that is, [X/O] = 0 at all [O… ▽ More Star formation models predict that the metal-poor initial mass function (IMF) can be substantially different from that observed in the metal-rich Milky Way. This changeover occurs because metal-poor gas clouds cool inefficiently due to their lower abundance of metals and dust. However, predictions for the metal-poor IMF to date rely on assuming Solar-scaled abundances, that is, [X/O] = 0 at all [O/H]. There is now growing evidence that elements such as C and O that dominate metal line cooling in the ISM do not follow Solar scaling at low metallicities. In this work, we extend models that predict the variation in the characteristic (or, the peak) IMF mass as a function of metallicity using [C/O] ratios derived from observations of metal-poor Galactic stars and of H II regions in dwarf galaxies. These data show [C/O] < 0 at sub-Solar [O/H], which leads to a substantially different metal-poor IMF in the metallicity range where C I and C II cooling dominate ISM thermodynamics, resulting in an increase in the characteristic mass by a factor as large as 7. An important consequence of this difference is a shift in the location of the transition from a top- to a bottom-heavy IMF upwards by 0.5 $-$ 1 dex in metallicity. Our findings indicate that the IMF is very sensitive to the assumptions around Solar-scaled ISM compositions in metal-poor systems (e.g., dwarf galaxies, the Galactic halo and metal-poor stars) that are a key focus of JWST. △ Less

Submitted 10 November, 2022; originally announced November 2022.

Comments: 14 pages, 8 figures. Accepted by MNRAS

The cosmic ray ionisation and $γ$-ray budgets of star-forming galaxies

Authors: Mark R. Krumholz, Roland M. Crocker, Stella S. R. Offner

Abstract: Cosmic rays in star-forming galaxies are a dominant source of both diffuse $γ$-ray emission and ionisation in gas too deeply shielded for photons to penetrate. Though the cosmic rays responsible for $γ$-rays and ionisation are of different energies, they are produced by the same star formation-driven sources, and thus galaxies' star formation rates, $γ$-ray luminosities, and ionisation rates shoul… ▽ More Cosmic rays in star-forming galaxies are a dominant source of both diffuse $γ$-ray emission and ionisation in gas too deeply shielded for photons to penetrate. Though the cosmic rays responsible for $γ$-rays and ionisation are of different energies, they are produced by the same star formation-driven sources, and thus galaxies' star formation rates, $γ$-ray luminosities, and ionisation rates should all be linked. In this paper we use up-to-date cross-section data to determine this relationship, finding that cosmic rays in a galaxy of star formation rate $\dot{M}_*$ and gas depletion time $t_\mathrm{dep}$ produce a maximum primary ionisation rate $ζ\approx 1\times 10^{-16} (t_\mathrm{dep}/\mbox{Gyr})^{-1}$ s$^{-1}$ and a maximum $γ$-ray luminosity $L_γ\approx 4\times 10^{39} (\dot{M}_*/\mathrm{M}_\odot\mbox{ yr}^{-1})$ erg s$^{-1}$ in the 0.1 - 100 GeV band. These budgets imply either that the ionisation rates measured in Milky Way molecular clouds include a significant contribution from local sources that elevate them above the Galactic mean, or that CR-driven ionisation in the Milky Way is enhanced by sources not linked directly to star formation. Our results also imply that ionisation rates in starburst systems are only moderately enhanced compared to those in the Milky Way. Finally, we point out that measurements of $γ$-ray luminosities can be used to place constraints on galactic ionisation budgets in starburst galaxies that are nearly free of systematic uncertainties on the details of cosmic ray acceleration. △ Less

Submitted 7 February, 2023; v1 submitted 7 November, 2022; originally announced November 2022.

Comments: 19 pages, 8 figures, accepted for publication in MNRAS

Outflows Driven by Direct and Reprocessed Radiation Pressure in Massive Star Clusters

Authors: Shyam H. Menon, Christoph Federrath, Mark R. Krumholz

Abstract: We use three-dimensional radiation hydrodynamic (RHD) simulations to study the formation of massive star clusters under the combined effects of direct ultraviolet (UV) and dust-reprocessed infrared (IR) radiation pressure. We explore a broad range of mass surface density $Σ\sim 10^2$-$10^5 \, \mathrm{M}_{\odot} \, \mathrm{pc}^{-2}$, spanning values typical of weakly star-forming galaxies to extrem… ▽ More We use three-dimensional radiation hydrodynamic (RHD) simulations to study the formation of massive star clusters under the combined effects of direct ultraviolet (UV) and dust-reprocessed infrared (IR) radiation pressure. We explore a broad range of mass surface density $Σ\sim 10^2$-$10^5 \, \mathrm{M}_{\odot} \, \mathrm{pc}^{-2}$, spanning values typical of weakly star-forming galaxies to extreme systems such as clouds forming super-star clusters, where radiation pressure is expected to be the dominant feedback mechanism. We find that star formation can only be regulated by radiation pressure for $Σ\lesssim 10^3 \, \mathrm{M}_{\odot} \, \mathrm{pc}^{-2}$, but that clouds with $Σ\lesssim 10^5 \, \mathrm{M}_{\odot} \, \mathrm{pc}^{-2}$ become super-Eddington once high star formation efficiencies ($\sim 80 \%$) are reached, and therefore launch the remaining gas in a steady outflow. These outflows achieve mass-weighted radial velocities of $\sim 15$ - $30 \,\mathrm{km} \, \mathrm{s}^{-1}$, which is $\sim 0.5$ - $2.0$ times the cloud escape speed. This suggests that radiation pressure is a strong candidate to explain recently observed molecular outflows found in young super-star clusters in nearby starburst galaxies. We quantify the relative importance of UV and IR radiation pressure in different regimes, and deduce that both are equally important for $Σ\sim 10^3 \, \mathrm{M}_{\odot} \, \mathrm{pc}^{-2}$, whereas clouds with higher (lower) density are increasingly dominated by the IR (UV) component. Comparison with control runs without either the UV or IR bands suggests that the outflows are primarily driven by the impulse provided by the UV component, while IR radiation has the effect of rendering a larger fraction of gas super-Eddington, and thereby increasing the outflow mass flux by a factor of $\sim 2$. △ Less

Submitted 20 March, 2023; v1 submitted 6 October, 2022; originally announced October 2022.

Comments: 15 pages, 11 figures. MNRAS accepted. v2: Minor changes in text made to address referee comments

Cosmic Ray Interstellar Propagation Tool using Itô Calculus (criptic): software for simultaneous calculation of cosmic ray transport and observational signatures

Authors: Mark R. Krumholz, Roland M. Crocker, Matt L. Sampson

Abstract: We present criptic, the Cosmic Ray Interstellar Propagation Tool using Itô Calculus, a new open-source software package to simulate the propagation of cosmic rays through the interstellar medium and to calculate the resulting observable non-thermal emission. Criptic solves the Fokker-Planck equation describing transport of cosmic rays on scales larger than that on which their pitch angles become a… ▽ More We present criptic, the Cosmic Ray Interstellar Propagation Tool using Itô Calculus, a new open-source software package to simulate the propagation of cosmic rays through the interstellar medium and to calculate the resulting observable non-thermal emission. Criptic solves the Fokker-Planck equation describing transport of cosmic rays on scales larger than that on which their pitch angles become approximately isotropic, and couples this to a rich and accurate treatment of the microphysical processes by which cosmic rays in the energy range $\sim$MeV to $\sim$PeV lose energy and produce emission. Criptic is deliberately agnostic as to both the cosmic ray transport model and the state of the background plasma through which cosmic rays travel. It can solve problems where cosmic rays stream, diffuse, or perform arbitrary combinations of both, and the coefficients describing these transport processes can be arbitrary functions of the background plasma state, the properties of the cosmic rays themselves, and local integrals of the cosmic ray field itself (e.g., the local cosmic ray pressure or pressure gradient). The code is parallelised using a hybrid OpenMP-MPI paradigm, allowing rapid calculations exploiting multiple cores and nodes on modern supercomputers. Here we describe the numerical methods used in the code, our treatment of the microphysical processes, and the set of code tests and validations we have performed. △ Less

Submitted 19 September, 2022; v1 submitted 27 July, 2022; originally announced July 2022.

Comments: 29 pages, 15 figures, accepted to MNRAS; very minor changes from first version. The software described is available from https://bitbucket.org/krumholz/criptic/src/master/, and a users' guide can be found at https://criptic.readthedocs.io/

Infrared Radiation Feedback Does Not Regulate Star Cluster Formation

Authors: Shyam H. Menon, Christoph Federrath, Mark R. Krumholz

Abstract: We present 3D radiation-hydrodynamical (RHD) simulations of star cluster formation and evolution in massive, self-gravitating clouds, whose dust columns are optically thick to infrared (IR) photons. We use \texttt{VETTAM} -- a recently developed, novel RHD algorithm, which uses the Variable Eddington Tensor (VET) closure -- to model the IR radiation transport through the cloud. We also use realist… ▽ More We present 3D radiation-hydrodynamical (RHD) simulations of star cluster formation and evolution in massive, self-gravitating clouds, whose dust columns are optically thick to infrared (IR) photons. We use \texttt{VETTAM} -- a recently developed, novel RHD algorithm, which uses the Variable Eddington Tensor (VET) closure -- to model the IR radiation transport through the cloud. We also use realistic temperature ($T$) dependent IR opacities ($κ$) in our simulations, improving upon earlier works in this area, which used either constant IR opacities or simplified power laws ($κ\propto T^2$). We investigate the impact of the radiation pressure of these IR photons on the star formation efficiency (SFE) of the cloud, and its potential to drive dusty winds. We find that IR radiation pressure is unable to regulate star formation or prevent accretion onto the star clusters, even for very high gas surface densities ($Σ> 10^5 M_{\odot} \, \mathrm{pc}^{-2}$), contrary to recent semi-analytic predictions and simulation results using simplified treatments of the dust opacity. We find that the commonly adopted simplifications of $κ\propto T^2$ or constant $κ$ for the IR dust opacities leads to this discrepancy, as those approximations overestimate the radiation force. By contrast, with realistic opacities that take into account the micro-physics of the dust, we find that the impact of IR radiation pressure on star formation is very mild, even at significantly high dust-to-gas ratios ($\sim 3$ times solar), suggesting that it is unlikely to be an important feedback mechanism in controlling star formation in the ISM. △ Less

Submitted 22 September, 2022; v1 submitted 28 June, 2022; originally announced June 2022.

Comments: 28 pages, 19 figures. MNRAS accepted. v2: Minor changes made to address referee comments

Spatial metallicity distribution statistics at $\lesssim 100$ pc scales in the AMUSING++ nearby galaxy sample

Authors: Zefeng Li, Emily Wisnioski, J. Trevor Mendel, Mark R. Krumholz, Lisa J. Kewley, Carlos López-Cobá, Sebastian F. Sánchez, Joseph P. Anderson, Lluís Galbany

Abstract: We analyse the spatial statistics of the 2D gas-phase oxygen abundance distributions in a sample of 219 local galaxies. We introduce a new adaptive binning technique to enhance the signal-to-noise ratio of weak lines, which we use to produce well-filled metallicity maps for these galaxies. We show that the two-point correlation functions computed from the metallicity distributions after removing r… ▽ More We analyse the spatial statistics of the 2D gas-phase oxygen abundance distributions in a sample of 219 local galaxies. We introduce a new adaptive binning technique to enhance the signal-to-noise ratio of weak lines, which we use to produce well-filled metallicity maps for these galaxies. We show that the two-point correlation functions computed from the metallicity distributions after removing radial gradients are in most cases well described by a simple injection-diffusion model. Fitting the data to this model yields the correlation length $l_{\rm corr}$, which describes the characteristic interstellar medium mixing length scale. We find typical correlation lengths $l_{\rm corr} \sim 1$ kpc, with a strong correlation between $l_{\rm corr}$ and stellar mass, star formation rate, and effective radius, a weak correlation with Hubble type, and significantly elevated values of $l_{\rm corr}$ in interacting or merging galaxies. We show that the trend with star formation rate can be reproduced by a simple transport+feedback model of interstellar medium turbulence at high star formation rate, and plausibly also at low star formation rate if dwarf galaxy winds have large mass-loading factors. We also report the first measurements of the injection width that describes the initial radii over which supernova remnants deposit metals. Inside this radius the metallicity correlation function is not purely the product of a competition between injection and diffusion. We show that this size scale is generally smaller than 60 pc. △ Less

Submitted 16 June, 2022; originally announced June 2022.

Comments: 18 pages, 18 figures, 1 table, submitted to MNRAS. Comments are welcome

Environmental variation of the low-mass IMF

Authors: Tabassum S Tanvir, Mark R Krumholz, Christoph Federrath

Abstract: We use a series of magnetohydrodynamic simulations including both radiative and protostellar outflow feedback to study the environmental variation of the initial mass function. The simulations represent a carefully-controlled experiment whereby we keep all dimensionless parameters of the flow constant except for those related to feedback. We show that radiation feedback suppresses the formation of… ▽ More We use a series of magnetohydrodynamic simulations including both radiative and protostellar outflow feedback to study the environmental variation of the initial mass function. The simulations represent a carefully-controlled experiment whereby we keep all dimensionless parameters of the flow constant except for those related to feedback. We show that radiation feedback suppresses the formation of lower mass objects more effectively as the surface density increases, but this only partially compensates for the decreasing Jeans mass in denser environments. Similarly, we find that protostellar outflows are more effective at suppressing the formation of massive stars in higher surface density environments. The combined effect of these two trends is towards an IMF with a lower characteristic mass and a narrower overall mass range in high surface density environments. We discuss the implications of these findings for the interpretation of observational evidence of IMF variation in early-type galaxies. △ Less

Submitted 14 September, 2022; v1 submitted 10 June, 2022; originally announced June 2022.

Comments: 15 pages, 13 figures, accepted by MNRAS

Turbulent diffusion of streaming cosmic rays in compressible, partially ionised plasma

Authors: Matt L. Sampson, James R. Beattie, Mark R. Krumholz, Roland M. Crocker, Christoph Federrath, Amit Seta

Abstract: Cosmic rays (CRs) are a dynamically important component of the interstellar medium (ISM) of galaxies. The $\sim$GeV CRs that carry most CR energy and pressure are likely confined by self-generated turbulence, leading them to stream along magnetic field lines at the ion Alfvén speed. However, the consequences of self-confinement for CR propagation on galaxy scales remain highly uncertain. In this p… ▽ More Cosmic rays (CRs) are a dynamically important component of the interstellar medium (ISM) of galaxies. The $\sim$GeV CRs that carry most CR energy and pressure are likely confined by self-generated turbulence, leading them to stream along magnetic field lines at the ion Alfvén speed. However, the consequences of self-confinement for CR propagation on galaxy scales remain highly uncertain. In this paper, we use a large ensemble of magnetohydrodynamical turbulence simulations to quantify how the basic parameters describing ISM turbulence -- the sonic Mach number, $\mathcal{M}$ (plasma compressibility), Alfvén Mach number, $\mathcal{M}_{A0}$ (strength of the large-scale field with respect to the turbulence), and ionisation fraction by mass, $χ$ -- affect the transport of streaming CRs. We show that the large-scale transport of CRs whose small-scale motion consists of streaming along field lines is well described as a combination of streaming along the mean field and superdiffusion both along (parallel to) and across (perpendicular to) it; $\mathcal{M}_{A0}$ drives the level of anisotropy between parallel and perpendicular diffusion and $χ$ modulates the magnitude of the diffusion coefficients, while in our choice of units, $\mathcal{M}$ is unimportant except in the sub-Alfvénic ($\mathcal{M}_{A0} \lesssim 0.5$) regime. Our finding that superdiffusion is ubiquitous potentially explains the apparent discrepancy between CR diffusion coefficients inferred from measurements close to individual sources compared to those measured on larger, Galactic scales. Finally, we present empirical fits for the diffusion coefficients as a function of plasma parameters that may be used as sub-grid recipes for global interstellar medium, galaxy or cosmological simulations. △ Less

Submitted 3 November, 2022; v1 submitted 17 May, 2022; originally announced May 2022.

Comments: Accepted in MNRAS

The Volatile Carbon to Oxygen Ratio as a Tracer for the Formation Locations of Interstellar Comets

Authors: Darryl Z. Seligman, Leslie A. Rogers, Samuel H. C. Cabot, John W. Noonan, Theodore Kareta, Kathleen E. Mandt, Fred Ciesla, Adam McKay, Adina D. Feinstein, W. Garrett Levine, Jacob L. Bean, Thomas Nordlander, Mark R. Krumholz, Megan Mansfield, Devin J. Hoover, Eric Van Clepper

Abstract: Based on the occurrence rates implied by the discoveries of 1I/`Oumuamua and 2I/Borisov, the forthcoming Rubin Observatory Legacy Survey of Space and Time (LSST) should detect $\ge1$ interstellar objects every year (Hoover et al. 2021). We advocate for future measurements of the production rates of H$_2$O, CO$_2$ and CO in these objects to estimate their carbon to oxygen ratios, which traces forma… ▽ More Based on the occurrence rates implied by the discoveries of 1I/`Oumuamua and 2I/Borisov, the forthcoming Rubin Observatory Legacy Survey of Space and Time (LSST) should detect $\ge1$ interstellar objects every year (Hoover et al. 2021). We advocate for future measurements of the production rates of H$_2$O, CO$_2$ and CO in these objects to estimate their carbon to oxygen ratios, which traces formation locations within their original protoplanetary disks. We review similar measurements for Solar System comets, which indicate formation interior to the CO snowline. By quantifying the relative processing in the interstellar medium and Solar System, we estimate that production rates will not be representative of primordial compositions for the majority of interstellar comets. Preferential desorption of CO and CO$_2$ relative to H$_2$O in the interstellar medium implies that measured C/O ratios represent lower limits on the primordial ratios. Specifically, production rate ratios of ${\rm Q}({\rm CO})/{\rm Q}({\rm H_2O})<.2$ and ${\rm Q}({\rm CO})/{\rm Q}({\rm H_2O})>1$ likely indicate formation interior and exterior to the CO snowline, respectively. The high C/O ratio of 2I/Borisov implies that it formed exterior to the CO snowline. We provide an overview of the currently operational facilities capable of obtaining these measurements that will constrain the fraction of ejected comets that formed exterior to the CO snowline. This fraction will provide key insights into the efficiency of and mechanisms for cometary ejection in exoplanetary systems. △ Less

Submitted 1 June, 2022; v1 submitted 27 April, 2022; originally announced April 2022.

Comments: Accepted for publication at PSJ, 33 pages, 14 figures, preprint for reference at Exoplanets IV Program Number 405.03

Gamma-ray emission from the Sagittarius Dwarf Spheroidal galaxy due to millisecond pulsars

Authors: Roland M. Crocker, Oscar Macias, Dougal Mackey, Mark R. Krumholz, Shin'ichiro Ando, Shunsaku Horiuchi, Matthew G. Baring, Chris Gordon, Thomas Venville, Alan R. Duffy, Rui-Zhi Yang, Felix Aharonian, J. A. Hinton, Deheng Song, Ashley J. Ruiter, Miroslav D. Filipović

Abstract: The Fermi Bubbles are giant, gamma-ray emitting lobes emanating from the nucleus of the Milky Way discovered in ~1-100 GeV data collected by the Large Area Telescope on board the Fermi Gamma-Ray Space Telescope. Previous work has revealed substructure within the Fermi Bubbles that has been interpreted as a signature of collimated outflows from the Galaxy's super-massive black hole. Here we show vi… ▽ More The Fermi Bubbles are giant, gamma-ray emitting lobes emanating from the nucleus of the Milky Way discovered in ~1-100 GeV data collected by the Large Area Telescope on board the Fermi Gamma-Ray Space Telescope. Previous work has revealed substructure within the Fermi Bubbles that has been interpreted as a signature of collimated outflows from the Galaxy's super-massive black hole. Here we show via a spatial template analysis that much of the gamma-ray emission associated to the brightest region of substructure -- the so-called cocoon -- is likely due to the Sagittarius dwarf spheroidal (Sgr dSph) galaxy. This large Milky Way satellite is viewed through the Fermi Bubbles from the position of the Solar System. As a tidally and ram-pressure stripped remnant, the Sgr dSph has no on-going star formation, but we nevertheless demonstrate that the dwarf's millisecond pulsar (MSP) population can plausibly supply the gamma-ray signal that our analysis associates to its stellar template. The measured spectrum is naturally explained by inverse Compton scattering of cosmic microwave background photons by high-energy electron-positron pairs injected by MSPs belonging to the Sgr dSph, combined with these objects' magnetospheric emission. This finding plausibly suggests that MSPs produce significant gamma-ray emission amongst old stellar populations, potentially confounding indirect dark matter searches in regions such as the Galactic Centre, the Andromeda galaxy, and other massive Milky Way dwarf spheroidals. △ Less

Submitted 8 September, 2022; v1 submitted 25 April, 2022; originally announced April 2022.

Comments: Updated to match version accepted for publication in Nature Astronomy (2022). 14 pages main text, 3 main figures, 7 extended data figures. For published version of the paper, see https://rdcu.be/cUZ6X

Journal ref: Nature Astronomy (2022)

Gas Accretion Can Drive Turbulence in Galaxies

Authors: John C. Forbes, Razieh Emami, Rachel S. Somerville, Shy Genel, Dylan Nelson, Annalisa Pillepich, Blakesley Burkhart, Greg L. Bryan, Mark R. Krumholz, Lars Hernquist, Stephanie Tonnesen, Paul Torrey, Viraj Pandya, Christopher C. Hayward

Abstract: The driving of turbulence in galaxies is deeply connected with the physics of feedback, star formation, outflows, accretion, and radial transport in disks. The velocity dispersion of gas in galaxies therefore offers a promising observational window into these processes. However, the relative importance of each of these mechanisms remains controversial. In this work we revisit the possibility that… ▽ More The driving of turbulence in galaxies is deeply connected with the physics of feedback, star formation, outflows, accretion, and radial transport in disks. The velocity dispersion of gas in galaxies therefore offers a promising observational window into these processes. However, the relative importance of each of these mechanisms remains controversial. In this work we revisit the possibility that turbulence on galactic scales is driven by the direct impact of accreting gaseous material on the disk. We measure this effect in a disk-like star-forming galaxy in IllustrisTNG, using the high-resolution cosmological magnetohydrodynamical simulation TNG50. We employ Lagrangian tracer particles with a high time cadence of only a few Myr to identify accretion and other events, such as star formation, outflows, and movement within the disk. The energies of particles as they arrive in the disk are measured by stacking the events in bins of time before and after the event. The average effect of each event is measured on the galaxy by fitting explicit models for the kinetic and turbulent energies as a function of time in the disk. These measurements are corroborated by measuring the cross-correlation of the turbulent energy in the different annuli of the disk with other time series, and searching for signals of causality, i.e. asymmetries in the cross-correlation across zero time lag. We find that accretion contributes to the large-scale turbulent kinetic energy even if it is not the dominant driver of turbulence in this $\sim 5 \times 10^{9} M_\odot$ stellar mass galaxy. Extrapolating this finding to a range of galaxy masses, we find that there are regimes where energy from direct accretion may dominate the turbulent energy budget, particularly in disk outskirts, galaxies less massive than the Milky Way, and at redshift $\sim 2$. △ Less

Submitted 11 April, 2022; originally announced April 2022.

Comments: Submitted to AAS Journals, comments welcome

The observable properties of cool winds from galaxies, AGN, and star clusters -- II. 3D models for the multiphase wind of M82

Authors: Yuxuan Yuan, Mark R. Krumholz, Crystal L. Martin

Abstract: Galactic winds are a crucial player in galaxy formation and evolution, but observations of them have proven extraordinarily difficult to interpret, leaving large uncertainties even in basic quantities such as mass outflow rates. Part of this uncertainty arises from the relatively simplistic models to which complex wind observations are often fit, which inevitably discard much of the available info… ▽ More Galactic winds are a crucial player in galaxy formation and evolution, but observations of them have proven extraordinarily difficult to interpret, leaving large uncertainties even in basic quantities such as mass outflow rates. Part of this uncertainty arises from the relatively simplistic models to which complex wind observations are often fit, which inevitably discard much of the available information. Here we present an analysis of the wind of the nearby dwarf starburst galaxy M82 using a semi-analytic model that is able to take advantage of the full three-dimensional information present in position-position-velocity data cubes measured in the Hi 21 cm line, the CO 2-1 line, and the Ha line. Our best-fitting model produces position-dependent spectra in good agreement with the observations, and shows that the total wind mass flux in the atomic and molecular phases is approx 10 M_sun yr-1 (corresponding to a mass loading factor of about 2 - 3), with less than a factor of two uncertainty; the mass flux in the warm ionised phase is more poorly constrained, and may be comparable to or smaller than this. At least over the few kpc off the plane for which we trace the outflow, it appears to be a wind escaping the galaxy, rather than a fountain that falls back. Our fits require that clouds of cool gas entrained into the wind expand only modestly, suggesting they are magnetically confined. Finally, we demonstrate that attempts to model the wind using simplifying assumptions such as instantaneous acceleration and a constant terminal wind speed can yield significantly erroneous results. △ Less

Submitted 11 April, 2022; originally announced April 2022.

Ion Alfvén velocity fluctuations and implications for the diffusion of streaming cosmic rays

Authors: James R. Beattie, Mark R. Krumholz, Christoph Federrath, Matt Sampson, Roland M. Crocker

Abstract: The interstellar medium (ISM) of star-forming galaxies is magnetized and turbulent. Cosmic rays (CRs) propagate through it, and those with energies from $\sim\,\rm{GeV} - \rm{TeV}$ are likely subject to the streaming instability, whereby the wave damping processes balances excitation of resonant ionic Alfvén waves by the CRs, reaching an equilibrium in which the propagation speed of the CRs is ver… ▽ More The interstellar medium (ISM) of star-forming galaxies is magnetized and turbulent. Cosmic rays (CRs) propagate through it, and those with energies from $\sim\,\rm{GeV} - \rm{TeV}$ are likely subject to the streaming instability, whereby the wave damping processes balances excitation of resonant ionic Alfvén waves by the CRs, reaching an equilibrium in which the propagation speed of the CRs is very close to the local ion Alfvén velocity. The transport of streaming CRs is therefore sensitive to ionic Alfvén velocity fluctuations. In this paper we systematically study these fluctuations using a large ensemble of compressible MHD turbulence simulations. We show that for sub-Alfvénic turbulence, as applies for a strongly magnetized ISM, the ionic Alfvén velocity probability density function (PDF) is determined solely by the density fluctuations from shocked gas forming parallel to the magnetic field, and we develop analytical models for the ionic Alfvén velocity PDF up to second moments. For super-Alfvénic turbulence, magnetic and density fluctuations are correlated in complex ways, and these correlations as well as contributions from the magnetic fluctuations sets the ionic Alfvén velocity PDF. We discuss the implications of these findings for underlying "macroscopic" diffusion mechanisms in CRs undergoing the streaming instability, including modeling the macroscopic diffusion coefficient for the parallel transport in sub-Alfvénic plasmas. We also describe how, for highly-magnetized turbulent gas, the gas density PDF, and hence column density PDF, can be used to access information about ionic Alfvén velocity structure from observations of the magnetized ISM. △ Less

Submitted 7 August, 2022; v1 submitted 25 March, 2022; originally announced March 2022.

Comments: Invited submission accepted for publication in Frontiers in Astronomy and Space Sciences: Multi-scale Magnetic Field Measurements in the Multi-Phase Interstellar Medium; 46 pages; 12 Figures, 2 tables

A radiation hydrodynamics scheme on adaptive meshes using the Variable Eddington Tensor (VET) closure

Authors: Shyam H. Menon, Christoph Federrath, Mark R. Krumholz, Rolf Kuiper, Benjamin D. Wibking, Manuel Jung

Abstract: We present a new algorithm to solve the equations of radiation hydrodynamics (RHD) in a frequency-integrated, two-moment formulation. Novel features of the algorithm include i) the adoption of a non-local Variable Eddington Tensor (VET) closure for the radiation moment equations, computed with a ray-tracing method, ii) support for adaptive mesh refinement (AMR), iii) use of a time-implicit Godunov… ▽ More We present a new algorithm to solve the equations of radiation hydrodynamics (RHD) in a frequency-integrated, two-moment formulation. Novel features of the algorithm include i) the adoption of a non-local Variable Eddington Tensor (VET) closure for the radiation moment equations, computed with a ray-tracing method, ii) support for adaptive mesh refinement (AMR), iii) use of a time-implicit Godunov method for the hyperbolic transport of radiation, and iv) a fixed-point Picard iteration scheme to accurately handle the stiff nonlinear gas-radiation energy exchange. Tests demonstrate that our scheme works correctly, yields accurate rates of energy and momentum transfer between gas and radiation, and obtains the correct radiation field distribution even in situations where more commonly used -- but less accurate -- closure relations like the Flux-limited Diffusion and Moment-1 approximations fail. Our scheme presents an important step towards performing RHD simulations with increasing spatial and directional accuracy, effectively improving their predictive capabilities. △ Less

Submitted 22 March, 2022; originally announced March 2022.

Comments: 7 pages, 2 figures, to appear in the proceedings of IAUS362: The predictive power of computational astrophysics as discovery tool (D. Bisikalo, T. Hanawa, C. Boily, J. Stone, eds.). arXiv admin note: substantial text overlap with arXiv:2202.08778

The Life and Times of Giant Molecular Clouds

Authors: Mélanie Chevance, Mark R. Krumholz, Anna F. McLeod, Eve C. Ostriker, Erik W. Rosolowsky, Amiel Sternberg

Abstract: Giant molecular clouds (GMCs) are the sites of star formation and stellar feedback in galaxies. Their properties set the initial conditions for star formation and their lifecycles determine how feedback regulates galaxy evolution. In recent years, the advent of high-resolution telescopes has enabled systematic GMC-scale studies of the molecular interstellar medium in nearby galaxies, now covering… ▽ More Giant molecular clouds (GMCs) are the sites of star formation and stellar feedback in galaxies. Their properties set the initial conditions for star formation and their lifecycles determine how feedback regulates galaxy evolution. In recent years, the advent of high-resolution telescopes has enabled systematic GMC-scale studies of the molecular interstellar medium in nearby galaxies, now covering a wide range of physical conditions and allowing for the first studies of how GMC properties depend on galactic environment. These observational developments have been accompanied by numerical simulations of improving resolution that are increasingly accurately accounting for the effects of the galactic-scale environment on GMCs, while simultaneously improving the treatment of the small-scale processes of star-formation and stellar feedback within them. The combination of these recent developments has greatly improved our understanding of the formation, evolution, and destruction of GMCs. We review the current state of the field, highlight current open questions, and discuss promising avenues for future studies. △ Less

Submitted 17 March, 2022; originally announced March 2022.

Comments: 37 pages, 5 figures; To appear in Protostars and Planets VII; Editors: Shu-ichiro Inutsuka, Yuri Aikawa, Takayuki Muto, Kengo Tomida, and Motohide Tamura

Energy balance and Alfvén Mach numbers in compressible magnetohydrodynamic turbulence with a large-scale magnetic field

Authors: James R. Beattie, Mark R. Krumholz, Raphael Skalidis, Christoph Federrath, Amit Seta, Roland M. Crocker, Philip Mocz, Neco Kriel

Abstract: Energy equipartition is a powerful theoretical tool for understanding astrophysical plasmas. It is invoked, for example, to measure magnetic fields in the interstellar medium (ISM), as evidence for small-scale turbulent dynamo action, and, in general, to estimate the energy budget of star-forming molecular clouds. In this study we motivate and explore the role of the volume-averaged root-mean-squa… ▽ More Energy equipartition is a powerful theoretical tool for understanding astrophysical plasmas. It is invoked, for example, to measure magnetic fields in the interstellar medium (ISM), as evidence for small-scale turbulent dynamo action, and, in general, to estimate the energy budget of star-forming molecular clouds. In this study we motivate and explore the role of the volume-averaged root-mean-squared (rms) magnetic coupling term between the turbulent, $δ\mathbf{B}$ and large-scale, $\mathbf{B}_0$ fields, $\left< (δ\mathbf{B}\cdot\mathbf{B}_0)^{2} \right>^{1/2}_{\mathcal{V}}$. By considering the second moments of the energy balance equations we show that the rms coupling term is in energy equipartition with the volume-averaged turbulent kinetic energy for turbulence with a sub-Alfvénic large-scale field. Under the assumption of exact energy equipartition between these terms, we derive relations for the magnetic and coupling term fluctuations, which provide excellent, parameter-free agreement with time-averaged data from 280 numerical simulations of compressible MHD turbulence. Furthermore, we explore the relation between the turbulent, mean-field and total Alfvén Mach numbers, and demonstrate that sub-Alfvénic turbulence can only be developed through a strong, large-scale magnetic field, which supports an extremely super-Alfvénic turbulent magnetic field. This means that the magnetic field fluctuations are significantly subdominant to the velocity fluctuations in the sub-Alfvénic large-scale field regime. Throughout our study, we broadly discuss the implications for observations of magnetic fields and understanding the dynamics in the magnetised ISM. △ Less

Submitted 24 July, 2022; v1 submitted 25 February, 2022; originally announced February 2022.

Comments: Accepted in MNRAS. 20 pages. 16 figures. 1 table

The evolution of turbulent galactic discs: gravitational instability, feedback and accretion

Authors: Omri Ginzburg, Avishai Dekel, Nir Mandelker, Mark R. Krumholz

Abstract: We study the driving of turbulence in star-froming disc galaxies of different masses at different epochs, using an analytic "bathtub" model. The disc of gas and stars is assumed to be in marginal Toomre instability. Turbulence is assumed to be sustained via an energy balance between its dissipation and three simultaneous energy sources. These are stellar feedback, inward transport due to disc inst… ▽ More We study the driving of turbulence in star-froming disc galaxies of different masses at different epochs, using an analytic "bathtub" model. The disc of gas and stars is assumed to be in marginal Toomre instability. Turbulence is assumed to be sustained via an energy balance between its dissipation and three simultaneous energy sources. These are stellar feedback, inward transport due to disc instability and clumpy accretion via streams. The transport rate is computed with two different formalisms, with similar results. To achieve the energy balance, the disc self-regulates either the mass fraction in clumps or the turbulent viscous torque parameter. In this version of the model, the efficiency by which the stream kinetic energy is converted into turbulence is a free parameter, $ξ_a$. We find that the contributions of the three energy sources are in the same ball park, within a factor of $\sim\!2$ in all discs at all times. In haloes that evolve to a mass $\leq 10^{12}\,\Msun$ by $z=0$ ($\leq 10^{11.5}\,\Msun$ at $z\!\sim\!2$), feedback is the main driver throughout their lifetimes. Above this mass, the main driver is either transport or accretion for very low or very high values of $ξ_a$, respectively. For an assumed $ξ_a(t)$ that declines in time, galaxies in halos with present-day mass $>\!10^{12}$ M$_\odot$ make a transition from accretion to transport dominance at intermediate redshifts, $z\! \sim\!3$, when their mass was $\geq\!10^{11.5}\,\Msun$. The predicted relation between star-formation rate and gas velocity dispersion is consistent with observations. △ Less

Submitted 7 May, 2022; v1 submitted 24 February, 2022; originally announced February 2022.

Comments: 17 pages, 6 figures. Accepted for publication in MNRAS

Journal ref: 2022MNRAS.513.6177G

VETTAM: A scheme for radiation hydrodynamics with adaptive mesh refinement using the variable Eddington tensor method

Authors: Shyam H. Menon, Christoph Federrath, Mark R. Krumholz, Rolf Kuiper, Benjamin D. Wibking, Manuel Jung

Abstract: We present Variable Eddington Tensor-closed Transport on Adaptive Meshes (\texttt{VETTAM}), a new algorithm to solve the equations of radiation hydrodynamics (RHD) with support for adaptive mesh refinement (AMR) in a frequency-integrated, two-moment formulation. The method is based on a non-local Variable Eddington Tensor (VET) closure computed with a hybrid characteristics scheme for ray tracing.… ▽ More We present Variable Eddington Tensor-closed Transport on Adaptive Meshes (\texttt{VETTAM}), a new algorithm to solve the equations of radiation hydrodynamics (RHD) with support for adaptive mesh refinement (AMR) in a frequency-integrated, two-moment formulation. The method is based on a non-local Variable Eddington Tensor (VET) closure computed with a hybrid characteristics scheme for ray tracing. We use a Godunov method for the hyperbolic transport of radiation with an implicit backwards-Euler temporal update to avoid the explicit timestep constraint imposed by the light-crossing time, and a fixed-point Picard iteration scheme to handle the nonlinear gas-radiation exchange term, with the two implicit update stages jointly iterated to convergence. We also develop a modified wave-speed correction method for AMR, which we find to be crucial for obtaining accurate results in the diffusion regime. We demonstrate the robustness of our scheme with a suite of pure radiation and RHD tests, and show that it successfully captures the streaming, static diffusion, and dynamic diffusion regimes and the spatial transitions between them, casts sharp shadows, and yields accurate results for rates of momentum and energy exchange between radiation and gas. A comparison between different closures for the radiation moment equations, with the Eddington approximation (0th-moment closure) and the $M_1$ approximation (1st-moment closure), demonstrates the advantages of the VET method (2nd-moment closure) over the simpler closure schemes. \texttt{VETTAM} has been coupled to the AMR \texttt{FLASH} (magneto-)hydrodynamics code and we summarize by reporting performance features and bottlenecks of our implementation. △ Less

Submitted 17 February, 2022; originally announced February 2022.

Comments: 23 pages, 15 figures. Accepted for publication in MNRAS