Search

A fraction of the dense cores that form within a turbulent molecular cloud will eventually collapse, leading to star formation. Identifying the physical criteria for cores to become unstable, and analyzing critical core properties, thus constitutes a necessary step toward the complete theory of star formation. To this end, here we quantify the characteristics of an ensemble of ``critical cores'' that are on the verge of collapse. This critical epoch was identified in a companion paper, which followed the dynamical evolution of prestellar cores in numerical simulations of turbulent, self-gravitating clouds. We find that radial profiles of density and turbulent velocity dispersion constructed for individual critical cores are consistent with our new model for turbulent equilibrium spheres (TESs). While there exists a global linewidth--size relation for a cloud with given size and Mach number, the turbulent scaling relations constructed around each core exhibit significant variations, locally regulating the critical density for a core to become unstable. As a result, there is no single density threshold for collapse, but instead cores collapse at a wide range of densities determined by the local sonic scale, modulated by the local gravitational potential environment, with a distribution expected for TESs with a limited range of turbulent velocity dispersion. The critical cores found in our simulations are mostly transonic; we do not find either purely thermal or highly turbulent cores. We find that the core mass function (CMF) of critical cores peaks around the characteristic mass scale associated with the average properties of a turbulent cloud. We highlight the importance of constructing the CMF at the critical time instead of sink particle mass functions, and derive the resolution requirements to unambiguously identify the peak of the CMF., Comment: 24 pages, 11 figures, submitted to ApJ; updated to include cross-references to the accompanying paper

A fundamental issue in star formation is understanding the precise mechanisms leading to the formation of prestellar cores, and their subsequent gravitationally unstable evolution. To address this question, we carefully construct a suite of turbulent, self-gravitating numerical simulations, and analyze the development and collapse of individual prestellar cores. We show that the numerical requirements for resolving the sonic scale and internal structure of anticipated cores are essentially the same in self-gravitating clouds, calling for the number of cells per dimension to increase quadratically with the cloud's Mach number. In our simulations, we follow evolution of individual cores by tracking the region around each gravitational potential minimum over time. Evolution in nascent cores is towards increasing density and decreasing turbulence, and there is a wide range of critical density for initiating collapse. At given spatial scale the turbulence level also varies widely, and tends to be correlated with density. By directly measuring the radial forces acting within cores, we identify a distinct transition to a state of gravitational runaway. We use our new theory for turbulent equilibrium spheres to predict the onset of each core's collapse. Instability is expected when the critical radius becomes smaller than the tidal radius; we find good agreement with the simulations. Interestingly, the imbalance between gravity and opposing forces is only $\sim 20\%$ during core collapse, meaning that this is a quasi-equilibrium rather than a free-fall process. For most of their evolution, cores exhibit both subsonic contraction and transonic turbulence inherited from core-building flows; supersonic radial velocities accelerated by gravity only appear near the end of the collapse., Comment: 33 pages, 19 figures, submitted to ApJ; updated to include cross-references to the accompanying paper

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

Traditional star formation subgrid models implemented in cosmological galaxy formation simulations, such as that of Springel & Hernquist (2003, hereafter SH03), employ adjustable parameters to satisfy constraints measured in the local Universe. In recent years, however, theory and spatially-resolved simulations of the turbulent, multiphase, star-forming ISM have begun to produce new first-principles models, which when fully developed can replace traditional subgrid prescriptions. This approach has advantages of being physically motivated and predictive rather than empirically tuned, and allowing for varying environmental conditions rather than being tied to local Universe conditions. As a prototype of this new approach, by combining calibrations from the TIGRESS numerical framework with the Pressure-Regulated Feedback-Modulated (PRFM) theory, simple formulae can be obtained for both the gas depletion time and an effective equation of state. Considering galaxies in TNG50, we compare the "native" simulation outputs with post-processed predictions from PRFM. At TNG50 resolution, the total midplane pressure is nearly equal to the total ISM weight, indicating that galaxies in TNG50 are close to satisfying vertical equilibrium. The measured gas scale height is also close to theoretical equilibrium predictions. The slopes of the effective equations of states are similar, but with effective velocity dispersion normalization from SH03 slightly larger than that from current TIGRESS simulations. Because of this and the decrease in PRFM feedback yield at high pressure, the PRFM model predicts shorter gas depletion times than the SH03 model at high densities and redshift. Our results represent a first step towards implementing new, numerically calibrated subgrid algorithms in cosmological galaxy formation simulations., Comment: 33 pages, 12 figures, Accepted for publication in ApJ. This is a Learning the Universe Publication. All codes and data used to produce this work can be found at the following $\href{https://github.com/sultan-hassan/tng50-post-processing-prfm}{GitHub \,Link.}$

We present a suite of six high-resolution chemo-dynamical simulations of isolated galaxies, spanning observed disk-dominated environments on the star-forming main sequence, as well as quenched, bulge-dominated environments. We compare and contrast the physics driving star formation and stellar feedback amongst the galaxies, with a view to modeling these processes in cosmological simulations. We find that the mass-loading of galactic outflows is coupled to the clustering of supernova explosions, which varies strongly with the rate of galactic rotation $\Omega = v_c/R$ via the Toomre length, leading to smoother gas disks in the bulge-dominated galaxies. This sets an equation of state in the star-forming gas that also varies strongly with $\Omega$, so that the bulge-dominated galaxies have higher mid-plane densities, lower velocity dispersions, and higher molecular gas fractions than their main sequence counterparts. The star formation rate in five out of six galaxies is independent of $\Omega$, and is consistent with regulation by the mid-plane gas pressure alone. In the sixth galaxy, which has the most centrally-concentrated bulge and thus the highest $\Omega$, we reproduce dynamical suppression of the star formation efficiency (SFE) in agreement with observations. This produces a transition away from pressure-regulated star formation., Comment: This is a Learning the Universe Publication, accepted for publication in ApJ. 29 pages, 14 figures, and analysis code at: https://github.com/sjeffreson/pressure_regulated_SF_analysis

With numerical simulations that employ adaptive ray-tracing (ART) for radiative transfer at the same time as evolving gas magnetohydrodynamics, thermodynamics, and photochemistry, it is possible to obtain a high resolution view of ultraviolet (UV) fields and their effects in realistic models of the multiphase interstellar medium. Here, we analyze results from TIGRESS-NCR simulations, which follow both far-UV (FUV) wavelengths, important for photoelectric heating and PAH excitation, and the Lyman continuum (LyC), which photoionizes hydrogen. Considering two models, representing solar neighborhood and inner galaxy conditions, we characterize the spatial distribution and time variation of UV radiation fields, and quantify their correlations with gas. We compare four approximate models for the FUV to simulated values to evaluate alternatives when full ART is infeasible. By convolving FUV radiation with density, we produce mock maps of dust emission. We introduce a method to calibrate mid-IR observations, for example from JWST, to obtain high resolution gas surface density maps. We then consider the LyC radiation field, finding most of the gas exposed to this radiation to be in ionization-recombination equilibrium and to have a low neutral fraction. Additionally, we characterize the ionization parameter as a function of environment. Using a simplified model of the LyC radiation field, we produce synthetic maps of emission measure (EM). We show that the simplified model can be used to extract an estimate of the neutral fraction of the photoionized gas and mean free path of ionizing radiation from observed EM maps in galaxies., Comment: 38 pages, 26 figures, accepted for publication in ApJ

Dense cores inherit turbulent motions from the interstellar medium in which they form. As a tool for comparison to both simulations and observations, it is valuable to construct theoretical core models that can relate their internal density and velocity structure while predicting their stability to gravitational collapse. To this end, we solve the angle-averaged equations of hydrodynamics under two assumptions: 1) the system is in a quasi-steady equilibrium; 2) the velocity field consists of radial bulk motion plus isotropic turbulence, with turbulent dispersion increasing as a power-law in the radius. The resulting turbulent equilibrium sphere (TES) solutions form a two-parameter family, characterized by the sonic radius $r_s$ and the power-law index $p$. The TES is equivalent to the Bonnor-Ebert (BE) sphere when $r_s\to \infty$. The density profile in outer regions of the TES is slightly shallower than the BE sphere, but is steeper than the logotropic model. Stability analysis shows that the TESs with size exceeding a certain critical radius are unstable to radial perturbations. The center-to-edge density contrast, mass, and radius of the marginally stable TES all increase with increasing average velocity dispersion. The FWHM of the column density profile is always smaller than the critical radius, by a larger factor at higher velocity dispersion, suggesting that observations need to probe beyond the FWHM to capture the full extent of turbulent cores. When applied to the highly turbulent regime typical of cluster-forming clumps, the critical mass and radius of the TES intriguingly resembles the typical mass and radius of observed star clusters., Comment: 17 pages, 11 figures, accepted for publication in ApJ. A python package that implements the model presented in this paper including a Jupyter notebook for reproducing the figures is available in https://github.com/sanghyukmoon/tesphere

Arkenstone is a new scheme that allows multiphase, stellar feedback-driven winds to be included in coarse resolution cosmological simulations. The evolution of galactic winds and their subsequent impact on the circumgalactic medium are altered by exchanges of mass, energy, momentum, and metals between their component phases. These exchanges are governed by complex, small-scale physical processes that cannot be resolved in cosmological simulations. In this second presentation paper, we describe Arkenstone's novel cloud particle approach for modelling unresolvable cool clouds entrained in hot, fast winds. This general framework allows models of the cloud-wind interaction, derived from state-of-the-art high-resolution simulations, to be applied in a large-scale context. In this work, we adopt a cloud evolution model that captures simultaneous cloud mass loss to and gain from the ambient hot phase via turbulent mixing and radiative cooling, respectively. We demonstrate the scheme using non-cosmological idealized simulations of a galaxy with a realistic circumgalactic medium component, using the Arepo code. We show that the ability of a high-specific energy wind component to perform preventative feedback may be limited by heavy loading of cool clouds coupled into it. We demonstrate that the diverging evolution of clouds of initially differing masses leads to a complex velocity field for the cool phase and a cloud mass function that varies both spatially and temporally in a non-trivial manner. These latter two phenomena can manifest in the simulation because of our choice of a Lagrangian discretisation of the cloud population, in contrast to other proposed schemes. This is a Learning the Universe publication., Comment: Submitted to MNRAS. 27 pages, 11 figures. This is a Learning the Universe publication

The molecular-to-atomic gas ratio is crucial to the evolution of the interstellar medium in galaxies. We investigate the balance between the atomic ($\Sigma_{\rm HI}$) and molecular gas ($\Sigma_{\rm H2}$) surface densities in eight nearby star-forming galaxies using new high-quality observations from MeerKAT and ALMA (for HI and CO, respectively). We define the molecular gas ratio as $R_{\rm mol} = \Sigma_{\rm H2} / \Sigma_{\rm HI}$ and measure how it depends on local conditions in the galaxy disks using multi-wavelength observations. We find that, depending on the galaxy, HI is detected at $>3\sigma$ out to 20-120 kpc in galactocentric radius ($r_{\rm gal}$). The typical radius at which $\Sigma_{\rm HI}$ reaches 1~$\rm M_\odot~pc^{-2}$ is $r_{\rm HI}\approx22$~kpc, which corresponds to 1-3 times the optical radius ($r_{25}$). $R_{\rm mol}$ correlates best with the dynamical equilibrium pressure, P$_{\rm DE}$, among potential drivers studied, with a median correlation coefficient of $<\rho>=0.89$. Correlations between $R_{\rm mol}$ and star formation rate, total gas and stellar surface density, metallicity, and $\Sigma_{\rm SFR}$/P$_{\rm DE}$ are present but somewhat weaker. Our results also show a direct correlation between P$_{\rm DE}$ and $\Sigma_{\rm SFR}$, supporting self-regulation models. Quantitatively, we measure similar scalings as previous works and attribute the modest differences that we find to the effect of varying resolution and sensitivity. At $r_{\rm gal} {\gtrsim}0.4~r_{25}$, atomic gas dominates over molecular gas, and at the balance of these two gas phases, we find that the baryon mass is dominated by stars, with $\Sigma_{*} > 5~\Sigma_{\rm gas}$. Our study constitutes an important step in the statistical investigation of how local galaxy properties impact the conversion from atomic to molecular gas in nearby galaxies., Comment: accepted for publication in A&A; 20 pages, 12 Figures (+4 appendix pages)

We present a new simulation suite for the star-forming interstellar medium (ISM) in galactic disks using the TIGRESS-NCR framework. Distinctive aspects of our simulation suite are: (1) sophisticated and comprehensive numerical treatments of essential physical processes including magnetohydrodynamics, self-gravity, and galactic differential rotation, as well as photochemistry, cooling, and heating coupled with ray-tracing UV radiation transfer and resolved supernova feedback and (2) wide parameter coverage including metallicity over $Z'\equiv Z/Z_\odot\sim0.1-3$, gas surface density $\Sigma_{\rm gas}\sim5-150 M_{\odot}{\rm pc^{-2}}$, and stellar surface density $\Sigma_{\rm star}\sim 1-50 M_{\odot}{\rm pc^{-2}}$. The range of emergent star formation rate surface density is $\Sigma_{\rm SFR}\sim 10^{-4}-0.5 M_{\odot}{\rm kpc^{-2}yr^{-1}}$ and ISM total midplane pressure is $P_{\rm tot}/k_B=10^3-10^6{\rm cm^{-3}K}$, with $P_{\rm tot}$ equal to the ISM weight $W$. For given $\Sigma_{\rm gas}$ and $\Sigma_{\rm star}$, we find $\Sigma_{\rm SFR} \propto Z'^{0.3}$. We provide an interpretation based on the pressure-regulated feedback-modulated (PRFM) star formation theory. We characterize feedback modulation in terms of the yield $\Upsilon$, defined as the ratio of each stress to $\Sigma_{\rm SFR}$. The thermal feedback yield varies sensitively with both weight and metallicity as $\Upsilon_{\rm th}\propto W^{-0.46}Z'^{-0.53}$, while the combined turbulent and magnetic feedback yield shows weaker dependence $\Upsilon_{\rm turb+mag}\propto W^{-0.22}Z'^{-0.18}$. The reduction in $\Sigma_{\rm SFR}$ at low metallicity is due mainly to enhanced thermal feedback yield, resulting from reduced attenuation of UV radiation. With the metallicity-dependent calibrations we provide, PRFM theory can be used for a new subgrid star formation prescription in cosmological simulations where the ISM is unresolved., Comment: Resubmitted to ApJ after minor revision

Bubbles driven by energy and mass injection from small scales are ubiquitous in astrophysical fluid systems and essential to feedback across multiple scales. In particular, O stars in young clusters produce high velocity winds that create hot bubbles in the surrounding gas. We demonstrate that the dynamical evolution of these bubbles is critically dependent upon the geometry of their interfaces with their surroundings and the nature of heat transport across these interfaces. These factors together determine the amount of energy that can be lost from the interior through cooling at the interface, which in turn determines the ability of the bubble to do work on its surroundings. We further demonstrate that the scales relevant to physical dissipation across this interface are extremely difficult to resolve in global numerical simulations of bubbles for parameter values of interest. This means the dissipation driving evolution of these bubbles in numerical simulations is often of a numerical nature. We describe the physical and numerical principles that determine the level of dissipation in these simulations; we use this, along with a fractal model for the geometry of the interfaces, to explain differences in convergence behavior between hydrodynamical and magneto-hydrodynamical simulations presented here. We additionally derive an expression for momentum as a function of bubble radius expected when the relevant dissipative scales are resolved and show that it still results in efficiently-cooled solutions as postulated in previous work., Comment: 41 pages, 20 figures, submitted to ApJ, comments welcome

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

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

We investigate the dynamical interaction between cosmic rays (CRs) and the multiphase interstellar medium (ISM) using numerical magnetohydrodynamic (MHD) simulations with a two-moment CR solver and TIGRESS simulations of star-forming galactic disks. We previously studied transport of CRs within TIGRESS outputs using a "post-processing" approach, and we now assess the effects of the MHD backreaction to CR pressure. We confirm our previous conclusion that there are three quite different regimes of CR transport in multiphase ISM gas, while also finding that simulations with "live MHD" predict a smoother CR pressure distribution. The CR pressure near the midplane is comparable to other pressure components in the gas, but the scale height of CRs is far larger. Next, with a goal of understanding the role of CRs in driving galactic outflows, we conduct a set of controlled simulations of the extraplanar region above $z=500$ pc, with imposed boundary conditions flowing from the midplane into this region. We explore a range of thermal and kinematic properties for the injected thermal gas, encompassing both hot, fast-moving outflows, and cooler, slower-moving outflows. The boundary conditions for CR energy density and flux are scaled from the supernova rate in the underlying TIGRESS model. Our simulations reveal that CRs efficiently accelerate extra-planar material if the latter is mostly warm/warm-hot gas, in which CRs stream at the Alfv\'en speed and the effective sound speed increases as density decreases. In contrast, CRs have very little effect on fast, hot outflows where the Alfv\'en speed is small, even when the injected CR momentum flux exceeds the injected MHD momentum flux., Comment: Resubmitted to ApJ after minor revisions

The internal velocity structure within dense gaseous cores plays a crucial role in providing the initial conditions for star formation in molecular clouds. However, the kinematic properties of dense gas at core scales (~0.01 - 0.1 pc) has not been extensively characterized because of instrument limitations until the unique capabilities of GBT-Argus became available. The ongoing GBT-Argus Large Program, Dynamics in Star-forming Cores (DiSCo) thus aims to investigate the origin and distribution of angular momenta of star-forming cores. DiSCo will survey all starless cores and Class 0 protostellar cores in the Perseus molecular complex down to ~0.01 pc scales with < 0.05 km/s velocity resolution using the dense gas tracer N$_2$H$^+$. Here, we present the first datasets from DiSCo toward the B1 and NGC 1333 regions in Perseus. Our results suggest that a dense core's internal velocity structure has little correlation with other core-scale properties, indicating these gas motions may be originated externally from cloud-scale turbulence. These first datasets also reaffirm the ability of GBT-Argus for studying dense core velocity structure and provided an empirical basis for future studies that address the angular momentum problem with a statistically broad sample., Comment: 17 pages, 12 figures, accepted by MNRAS

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

Chemistry plays a key role in many aspects of astrophysical fluids. Atoms and molecules are agents for heating and cooling, determine the ionization fraction, serve as observational tracers, and build the molecular foundation of life. We present the implementation of a chemistry module in the publicly available magneto-hydrodynamic code Athena++. We implement several chemical networks and heating and cooling processes suitable for simulating the interstellar medium (ISM). A general chemical network framework in the KIDA format is also included, allowing the user to easily implement their own chemistry. Radiation transfer and cosmic-ray ionization are coupled with chemistry and solved with the simple six-ray approximation. The chemical and thermal processes are evolved as a system of coupled ODEs with an implicit solver from the CVODE library. We perform and present a series of tests to ensure the numerical accuracy and convergence of the code. Many tests combine chemistry with gas dynamics, including comparisons with analytic solutions, 1D problems of the photo-dissociation regions and shocks, and realistic 3D simulations of the turbulent ISM. We release the code with the new public version of Athena++, aiming to provide a robust and flexible code for the astrochemical simulation community.

Nuclear rings at the centers of barred galaxies are known to be strongly magnetized. To explore the effects of magnetic fields on star formation in these rings and nuclear gas flows, we run magnetohydrodynamic simulations in which there is a temporally-constant magnetized inflow to the ring, representing a bar-driven inflow. The mass inflow rate is $1\,M_\odot\,\mathrm{yr}^{-1}$, and we explore models with a range of field strength in the inflow. We adopt the TIGRESS framework developed by Kim & Ostriker to handle radiative heating and cooling, star formation, and resulting supernova (SN) feedback. We find that magnetic fields are efficiently amplified in the ring due to rotational shear and SN feedback. Within a few $100\,\mathrm{Myr}$, the turbulent component $B_\mathrm{trb}$ in the ring saturates at $\sim 35\,\mu\mathrm{G}$ (in rough equipartition with the turbulent kinetic energy density), while the regular component $B_\mathrm{reg}$ exceeds $50\,\mu\mathrm{G}$. Expanding superbubbles created by clustered SN explosions vertically drag predominantly-toroidal fields from near the midplane to produce poloidal fields in high-altitude regions. The growth of magnetic fields greatly suppresses star formation at late times. Simultaneously, strong magnetic tension in the ring drives radially inward accretion flows from the ring to form a circumnuclear disk in the central region; this feature is absent in the unmagnetized model., Comment: 29 pages, 17 figures, accepted for publication in ApJ

We measure empirical relationships between the local star formation rate (SFR) and properties of the star-forming molecular gas on 1.5 kpc scales across 80 nearby galaxies. These relationships, commonly referred to as "star formation laws," aim at predicting the local SFR surface density from various combinations of molecular gas surface density, galactic orbital time, molecular cloud free-fall time, and the interstellar medium dynamical equilibrium pressure. Leveraging a multiwavelength database built for the PHANGS survey, we measure these quantities consistently across all galaxies and quantify systematic uncertainties stemming from choices of SFR calibrations and the CO-to-H$_2$ conversion factors. The star formation laws we examine show 0.3-0.4 dex of intrinsic scatter, among which the molecular Kennicutt-Schmidt relation shows a $\sim$10% larger scatter than the other three. The slope of this relation ranges $\beta\approx0.9{-}1.2$, implying that the molecular gas depletion time remains roughly constant across the environments probed in our sample. The other relations have shallower slopes ($\beta\approx0.6{-}1.0$), suggesting that the star formation efficiency (SFE) per orbital time, the SFE per free-fall time, and the pressure-to-SFR surface density ratio (i.e., the feedback yield) may vary systematically with local molecular gas and SFR surface densities. Last but not least, the shapes of the star formation laws depend sensitively on methodological choices. Different choices of SFR calibrations can introduce systematic uncertainties of at least 10-15% in the star formation law slopes and 0.15-0.25 dex in their normalization, while the CO-to-H$_2$ conversion factors can additionally produce uncertainties of 20-25% for the slope and 0.10-0.20 dex for the normalization., Comment: 10 pages main text + 2 appendices. ApJL in press. Data products available at https://www.canfar.net/storage/list/phangs/RELEASES/Sun_etal_2023 . Slides summarizing key results can be found at https://www.dropbox.com/s/5gsegexeo9n0t05/Sun_et_PHANGS_2023.pptx?dl=0

We perform 3D radiation hydrodynamic local shearing box simulations to study the outcome of gravitational instability (GI) in optically thick Active Galactic Nuclei (AGN) accretion disks. GI develops when the Toomre parameter QT \leq 1, and may lead to turbulent heating that balances radiative cooling. However, when radiative cooling is too efficient, the disk may undergo runaway gravitational fragmentation. In the fully gas-pressure-dominated case, we confirm the classical result that such a thermal balance holds when the Shakura-Sunyaev viscosity parameter (alpha) due to the gravitationally-driven turbulence is \sim 0.2, corresponding to dimensionless cooling times Omega tcool \sim 5. As the fraction of support by radiation pressure increases, the disk becomes more prone to fragmentation, with a reduced (increased) critical value of alpha (omega tcool). The effect is already significant when the radiation pressure exceeds 10% of the gas pressure, while fully radiation-pressure-dominated disks fragment at Omega tcool <50 . The latter translates to a maximum turbulence level alpha<0.02, comparable to that generated by Magnetorotational Instability (MRI). Our results suggest that gravitationally unstable (QT \sim 1) outer regions of AGN disks with significant radiation pressure (likely for high/near- Eddington accretion rates) should always fragment into stars, and perhaps black holes., Comment: 26 pages, 19 figures, ApJ in Press

Arkenstone is a new model for multiphase, stellar feedback driven galactic winds designed for inclusion in coarse resolution cosmological simulations. In this first paper of a series, we describe the features that allow Arkenstone to properly treat high specific energy wind components and demonstrate them using idealised non-cosmological simulations of a galaxy with a realistic CGM, using the Arepo code. Hot, fast gas phases with low mass loadings are predicted to dominate the energy content of multiphase outflows. In order to treat the huge dynamic range of spatial scales involved in cosmological galaxy formation at feasible computational expense, cosmological volume simulations typically employ a Lagrangian code or else use adaptive mesh refinement with a quasi-Lagrangian refinement strategy. However, it is difficult to inject a high specific energy wind in a Lagrangian scheme without incurring artificial burstiness. Additionally, the low densities inherent to this type of flow result in poor spatial resolution. Arkenstone addresses these issues with a novel scheme for coupling energy into the ISM/CGM transition region which also provides the necessary level of refinement at the base of the wind. In the absence of our improvements, we show that poor spatial resolution near the sonic point of a hot, fast outflow leads to an underestimation of gas acceleration as the wind propagates. We explore the different mechanisms by which low and high specific energy winds can regulate the SFR of galaxies. In future work, we will demonstrate other aspects of the Arkenstone model., Comment: Published MNRAS, 27 pages, 17 figures; updated for consistency with journal version

We present the first results from a high resolution simulation with a focus on galactic wind driving for an isolated galaxy with a halo mass of $\sim 10^{11}$ M$_{\odot}$ (similar to the Large Magellanic Cloud) and a total gas mass of $\sim 6 \times 10^{8}$ M$_{\odot}$, resulting in $\sim 10^{8}$ gas cells at $\sim 4$ M$_{\odot}$ mass resolution. We adopt a resolved stellar feedback model with non-equilibrium cooling and heating, including photoelectric heating and photo-ionizing radiation, as well as supernovae (SNe), coupled to the second order meshless finite mass (MFM) method for hydrodynamics. These features make this the largest resolved-ISM galaxy model run to date. We find mean star formation rates around $0.05$ M$_{\odot}$ yr$^{-1}$ and evaluate typical time averaged loading factors for mass ($\eta_\mathrm{M}$ $\sim$ 1.0, in good agreement with recent observations) and energy ($\eta_\mathrm{E}$ $\sim$ 0.01). The bulk of the mass of the wind is transported by the warm ($T < 5 \times 10^5$K) phase, while there is a similar amount of energy transported in the warm and the hot phases ($T > 5 \times 10^5$K). We find an average opening angle of 30 degrees for the wind, decreasing with higher altitude above the midplane. The wind mass loading is decreasing (flat) for the warm (hot) phase as a function of the star formation surface rate density $\Sigma_{\rm SFR}$, while the energy loading shows inverted trends with $\Sigma_{\rm SFR}$, decreasing for the warm wind and increasing for the hot wind, although with very shallow slopes. These scalings are in good agreement with previous simulations of resolved wind driving in the multi-phase ISM., Comment: submitted to ApJ, 26 pages, 14 figures, comments welcome

Massive, young stars are the main source of energy that maintains multiphase structure and turbulence in the interstellar medium (ISM), and without this "feedback" the star formation rate (SFR) would be much higher than is observed. Rapid energy loss in the ISM and efficient energy recovery by stellar feedback lead to co-regulation of SFRs and the ISM state. Realistic approaches to this problem should solve the dynamical evolution of the ISM, including star formation, and the input of feedback energy self-consistently and accurately. Here, we present the TIGRESS-NCR numerical framework, in which UV radiation, supernovae, cooling and heating processes, and gravitational collapse are modeled explicitly. We use an adaptive ray tracing method for UV radiation transfer from star clusters represented by sink particles, accounting for attenuation by dust and gas. We solve photon-driven chemical equations to determine the abundances of H (time-dependent) and C/O-bearing species (steady-state), which then set cooling and heating rates self-consistently. Applying these methods, we present high-resolution magnetohydrodynamics simulations of differentially rotating local galactic disks representing typical conditions of nearby star-forming galaxies. We analyze ISM properties and phase distributions and show good agreement with existing multiwavelength galactic observations. We measure midplane pressure components (turbulent, thermal, and magnetic) and the weight, demonstrating that vertical dynamical equilibrium holds. We quantify the ratios of pressure components to the SFR surface density, which we call the feedback yields. The TIGRESS-NCR framework will allow for a wide range of parameter exploration, including low metallicity system., Comment: ApJ accepted

We present an efficient heating/cooling method coupled with chemistry and ultraviolet (UV) radiative transfer, which can be applied to numerical simulations of the interstellar medium (ISM). We follow the time-dependent evolution of hydrogen species (H$_2$, H, H$^+$), assume carbon/oxygen species (C, C$^+$, CO, O, and O$^+$) are in formation-destruction balance given the non-steady hydrogen abundances, and include essential heating/cooling processes needed to capture thermodynamics of all ISM phases. UV radiation from discrete point sources and the diffuse background is followed through adaptive ray tracing and a six-ray approximation, respectively, allowing for H$_2$ self-shielding; cosmic ray (CR) heating and ionization are also included. To validate our methods and demonstrate their application for a range of density, metallicity, and radiation field, we conduct a series of tests, including the equilibrium curves of thermal pressure vs. density, the chemical and thermal structure in photo-dissociation regions, H I-to-H$_2$ transitions, and the expansion of H II regions and radiative supernova remnants. Careful treatment of photochemistry and CR ionization is essential for many aspects of ISM physics, including identifying the thermal pressure at which cold and warm neutral phases co-exist. We caution that many current heating and cooling treatments used in galaxy formation simulations do not reproduce the correct thermal pressure and ionization fraction in the neutral ISM. Our new model is implemented in the MHD code Athena and incorporated in the TIGRESS simulation framework, for use in studying the star-forming ISM in a wide range of environments., Comment: 57 pages, 22 figures; accepted for publication in ApJS

We present a rich, multiwavelength, multiscale database built around the PHANGS-ALMA CO$\,$(2-1) survey and ancillary data. We use this database to present the distributions of molecular cloud populations and sub-galactic environments in 80 PHANGS galaxies, to characterize the relationship between population-averaged cloud properties and host galaxy properties, and to assess key timescales relevant to molecular cloud evolution and star formation. We show that PHANGS probes a wide range of kpc-scale gas, stellar, and star formation rate (SFR) surface densities, as well as orbital velocities and shear. The population-averaged cloud properties in each aperture correlate strongly with both local environmental properties and host galaxy global properties. Leveraging a variable selection analysis, we find that the kpc-scale surface densities of molecular gas and SFR tend to possess the most predictive power for the population-averaged cloud properties. Once their variations are controlled for, galaxy global properties contain little additional information, which implies that the apparent galaxy-to-galaxy variations in cloud populations are likely mediated by kpc-scale environmental conditions. We further estimate a suite of important timescales from our multiwavelength measurements. The cloud-scale free-fall time and turbulence crossing time are ${\sim}5{-}20$ Myr, comparable to previous cloud lifetime estimates. The timescales for orbital motion, shearing, and cloud-cloud collisions are longer, ${\sim}100$ Myr. The molecular gas depletion time is $1{-}3$ Gyr and shows weak to no correlations with the other timescales in our data. We publish our measurements online and expect them to have broad utility to future studies of molecular clouds and star formation., Comment: 32 pages + 6 appendices. AJ in press. Data products available at https://www.canfar.net/storage/list/phangs/RELEASES/Sun_etal_2022 . Associated software package available at https://doi.org/10.5281/zenodo.6584842 . Slides summarizing the key results can be found at https://www.dropbox.com/s/7llnrob9wi0isfk/Sun_et_PHANGS_2022.pptx?dl=0

The star formation rate (SFR) in galactic disks depends on both the quantity of available interstellar medium (ISM) gas and its physical state. Conversely, the ISM's physical state depends on the SFR, because the "feedback" energy and momentum injected by recently-formed massive stars is crucial to offsetting losses from turbulent dissipation and radiative cooling. The ISM's physical state also responds to the gravitational field that confines it, with increased weight driving higher pressure. In a quasi-steady state, it is expected that the mean total pressure of different thermal phases will match each other, that the component pressures and total pressure will satisfy thermal and dynamical equilibrium requirements, and that the SFR will adjust as needed to provide the requisite stellar radiation and supernova feedback. The pressure-regulated, feedback-modulated (PRFM) theory of the star-forming ISM formalizes these ideas, leading to a prediction that the SFR per unit area, Sigma_SFR, will scale nearly linearly with ISM weight W. In terms of large-scale gas surface density Sigma, stellar plus dark matter density rho_sd, and effective ISM velocity dispersion sigma_eff, an observable weight estimator is W~P_DE=pi G Sigma^2/2+Sigma(2G rho_sd)^{1/2} sigma_eff, and this is predicted to match the total midplane pressure P_tot. Using a suite of multiphase magnetohydrodynamic simulations run with the TIGRESS computational framework, we test the principles of the PRFM model and calibrate the total feedback yield Upsilon_tot = P_tot/Sigma_SFR ~ 1000 km/s, as well as its components. We compare results from TIGRESS to theory, previous numerical simulations, and observations, finding excellent agreement., Comment: 34 pages, 15 figures; now accepted to ApJ

The observed star formation rate of the Milky Way can be explained by applying a metallicity-dependent factor to convert CO luminosity to molecular gas mass and a star formation efficiency per free-fall time that depends on the virial parameter of a molecular cloud. These procedures also predict the trend of star formation rate surface density with Galactocentric radius. The efficiency per free-fall time variation with virial parameter plays the major role in bringing theory into agreement with observations for the total star formation rate, while the metallicity dependence of the CO luminosity to mass conversion is most notable in the variation with Galactocentric radius. Application of these changes resolves a factor of over 100 discrepancy between observed and theoretical star formation rates that has been known for nearly 50 years., Comment: In press, Astrophysical Journal Letters 11 pages, 6 figures

We study the propagation of mildly-relativistic cosmic rays (CRs) in multiphase interstellar medium environments with conditions typical of nearby disk galaxies. We employ the techniques developed in Armillotta+21 to post-process three high-resolution TIGRESS magnetohydrodynamic simulations modeling local patches of star-forming galactic disks. Together, the three simulations cover a wide range of gas surface density, gravitational potential, and star formation rate (SFR). Our prescription for CR propagation includes the effects of advection by the background gas, streaming along the magnetic field at the local ion Alfv\'en speed, and diffusion relative to the Alfv\'en waves, with the diffusion coefficient set by the balance between streaming-driven Alfv\'en wave excitation and damping mediated by local gas properties. We find that the combined transport processes are more effective in environments with higher SFR. These environments are characterized by higher-velocity hot outflows (created by clustered supernovae) that rapidly advect CRs away from the galactic plane. As a consequence, the ratio of midplane CR pressure to midplane gas pressures decreases with increasing SFR. We also use the post-processed simulations to make predictions regarding potential dynamical impacts of CRs. The relatively flat CR pressure profiles near the midplane argue that they would not provide significant support against gravity for most of the ISM mass. However, the CR pressure gradients are larger than the other pressure gradients in the extra-planar region (|z|>0.5 kpc), suggesting that CRs may affect the dynamics of galactic fountains and/or winds. The degree of this impact is expected to increase in environments with lower SFR., Comment: Accepted for publication in ApJ

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., Comment: 37 pages, 5 figures; To appear in Protostars and Planets VII; Editors: Shu-ichiro Inutsuka, Yuri Aikawa, Takayuki Muto, Kengo Tomida, and Motohide Tamura

Observations indicate that the star formation rate (SFR) of nuclear rings varies considerably with time and is sometimes asymmetric rather than being uniform across a ring. To understand what controls temporal and spatial distributions of ring star formation, we run semi-global, hydrodynamic simulations of nuclear rings subject to time-varying and/or asymmetric mass inflow rates. These controlled variations in the inflow lead to variations in the star formation, while the ring orbital period ($18\,{\rm Myr}$) and radius ($600\,{\rm pc}$) remain approximately constant. We find that both the mass inflow rate and supernova feedback affect the ring SFR. An oscillating inflow rate with period $\Delta \tau_\text{in}$ and amplitude 20 causes large-amplitude (a factor of $\gtrsim 5$), quasi-periodic variations of the SFR, when $\Delta \tau_\text{in} \gtrsim 50\,{\rm Myr}$. We find that the time-varying ISM weight and midplane pressure track each other closely, establishing an instantaneous vertical equilibrium. The measured time-varying depletion time is consistent with the prediction from self-regulation theory provided the time delay between star formation and supernova feedback is taken into account. The supernova feedback is responsible only for small-amplitude (a factor of $\sim 2$) fluctuations of the SFR with a timescale $\lesssim 40\,{\rm Myr}$. Asymmetry in the inflow rate does not necessarily lead to asymmetric star formation in nuclear rings. Only when the inflow rate from one dust lane is suddenly increased by a large factor, the rings undergo a transient period of lopsided star formation., Comment: 15 pages, 6 figures, accepted for publication in ApJ. Replaced with minor text corrections

Stellar winds contain enough energy to easily disrupt the parent cloud surrounding a nascent star cluster, and for this reason have been considered candidates for regulating star formation. However, direct observations suggest most wind power is lost, and Lancaster21a,b recently proposed that this is due to efficient mixing and cooling processes. Here, we simulate star formation with wind feedback in turbulent, self-gravitating clouds, extending our previous work. Our simulations cover clouds with initial surface density $10^2-10^4$ $M_{\odot} \, {\rm pc}^{-2}$, and show that star formation and residual gas dispersal is complete within 2 - 8 initial cloud free-fall times. The "Efficiently Cooled" model for stellar wind bubble evolution predicts enough energy is lost for the bubbles to become momentum-driven, we find this is satisfied in our simulations. We also find that wind energy losses from turbulent, radiative mixing layers dominate losses by "cloud leakage" over the timescales relevant for star formation. We show that the net star formation efficiency (SFE) in our simulations can be explained by theories that apply wind momentum to disperse cloud gas, allowing for highly inhomogeneous internal cloud structure. For very dense clouds, the SFE is similar to those observed in extreme star-forming environments. Finally, we find that, while self-pollution by wind material is insignificant in cloud conditions with moderate density (only $\lesssim 10^{-4}$ of the stellar mass originated in winds), our simulations with conditions more typical of a super star cluster have star particles that form with as much as 1\% of their mass in wind material., Comment: 20 pages, 5 figures, submitted to ApJL, comments welcome

Cosmic ray transport on galactic scales depends on the detailed properties of the magnetized, multiphase interstellar medium (ISM). In this work, we post-process a high-resolution TIGRESS magnetohydrodynamic simulation modeling a local galactic disk patch with a two-moment fluid algorithm for cosmic ray transport. We consider a variety of prescriptions for the cosmic rays, from a simple purely diffusive formalism with constant scattering coefficient, to a physically-motivated model in which the scattering coefficient is set by critical balance between streaming-driven Alfv\'en wave excitation and damping mediated by local gas properties. We separately focus on cosmic rays with kinetic energies of $\sim 1$ GeV (high-energy) and $\sim 30$~MeV (low-energy), respectively important for ISM dynamics and chemistry. We find that simultaneously accounting for advection, streaming, and diffusion of cosmic rays is crucial for properly modeling their transport. Advection dominates in the high-velocity, low-density, hot phase, while diffusion and streaming are more important in higher density, cooler phases. Our physically-motivated model shows that there is no single diffusivity for cosmic-ray transport: the scattering coefficient varies by four or more orders of magnitude, maximal at density $n_\mathrm{H} \sim 0.01\, \mathrm{cm}^{-3}$. Ion-neutral damping of Alfv\'en waves results in strong diffusion and nearly uniform cosmic ray pressure within most of the mass of the ISM. However, cosmic rays are trapped near the disk midplane by the higher scattering rate in the surrounding lower-density, higher-ionization gas. The transport of high-energy cosmic rays differs from that of low-energy cosmic rays, with less effective diffusion and greater energy losses for the latter., Comment: Accepted for publication in ApJ

Nuclear rings are sites of intense star formation at the centers of barred galaxies. To understand what determines the structure and star formation rate (SFR; $\dot{M}_{\rm SF}$) of nuclear rings, we run semi-global, hydrodynamic simulations of nuclear rings subject to constant mass inflow rates $\dot{M}_{\rm in}$. We adopt the TIGRESS framework of Kim \& Ostriker to handle radiative heating and cooling, star formation, and related supernova (SN) feedback. We find that the SN feedback is never strong enough to destroy the ring or quench star formation everywhere in the ring. Under the constant $\dot{M}_{\rm in}$, the ring star formation is very steady and persistent, with the SFR exhibiting only mild temporal fluctuations. The ring SFR is tightly correlated with the inflow rate as $\dot{M}_{\rm SF}\approx 0.8\dot{M}_{\rm in}$, for a range of $\dot{M}_{\rm in}=0.125-8\,M_\odot\,{\rm yr}^{-1}$. Within the ring, vertical dynamical equilibrium is maintained, with the midplane pressure (powered by SN feedback) balancing the weight of the overlying gas. The SFR surface density is correlated nearly linearly with the midplane pressure, as predicted by the pressure-regulated, feedback-modulated star formation theory. Based on our results, we argue that the ring SFR is causally controlled by $\dot{M}_\text{in}$, while the ring gas mass adapts to the SFR to maintain the vertical dynamical equilibrium under the gravitational field arising from both gas and stars., Comment: 31 pages, 22 figures, accepted for publication in ApJ

In a companion paper, we develop a theory for the evolution of stellar wind driven bubbles in dense, turbulent clouds. This theory proposes that turbulent mixing at a fractal bubble-shell interface leads to highly efficient cooling, in which the vast majority of the input wind energy is radiated away. This energy loss renders the majority of the bubble evolution momentum-driven rather than energy-driven, with expansion velocities and pressures orders of magnitude lower than in the classical Weaver77 solution. In this paper, we validate our theory with three-dimensional, hydrodynamic simulations. We show that extreme cooling is not only possible, but is generic to star formation in turbulent clouds over more than three orders of magnitude in density. We quantify the few free parameters in our theory, and show that the momentum exceeds the wind input rate by only a factor ~ 1.2-4. We verify that the bubble/cloud interface is a fractal with dimension ~ 2.5-2.7. The measured turbulent amplitude (v_t ~ 200-400 km/s) in the hot gas near the interface is shown to be consistent with theoretical requirements for turbulent diffusion to efficiently mix and radiate away most of the wind energy. The fraction of energy remaining after cooling is only 1-\Theta ~ 0.1-0.01, decreasing with time, explaining observations that indicate low hot-gas content and weak dynamical effects of stellar winds., Comment: Accepted for publication in ApJ, 36 pages, 25 figures, but short summary and conclusion. Comments welcome

Winds from massive stars have velocities of 1000 km/s or more, and produce hot, high pressure gas when they shock. We develop a theory for the evolution of bubbles driven by the collective winds from star clusters early in their lifetimes, which involves interaction with the turbulent, dense interstellar medium of the surrounding natal molecular cloud. A key feature is the fractal nature of the hot bubble's surface. The large area of this interface with surrounding denser gas strongly enhances energy losses from the hot interior, enabled by turbulent mixing and subsequent cooling at temperatures T = 10^4-10^5 K where radiation is maximally efficient. Due to the extreme cooling, the bubble radius scales differently (R ~ t^1/2) from the classical Weaver77 solution, and has expansion velocity and momentum lower by factors of 10-10^2 at given R, with pressure lower by factors of 10^2 - 10^3. Our theory explains the weak X-ray emission and low shell expansion velocities of observed sources. We discuss further implications of our theory for observations of the hot bubbles and cooled expanding shells created by stellar winds, and for predictions of feedback-regulated star formation in a range of environments. In a companion paper, we validate our theory with a suite of hydrodynamic simulations., Comment: Accepted for publication in ApJ, comments welcome, 23 pages 4 figures

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

Cosmic rays (CRs) have critical impacts in the multiphase interstellar medium (ISM), driving dynamical motions in low-density plasma and modifying the ionization state, temperature, and chemical composition of higher-density atomic and molecular gas. We present a study of CR propagation between the ionized ISM and a neutral cloud. Using one-dimensional magnetohydrodynamic particle-in-cell simulations which include ion-neutral drag to damp Alfv$\acute{\text{e}}$n waves in the cloud, we self-consistently evolve the kinetic physics of CRs and fluid dynamics of the multiphase gas. By introducing the cloud in our periodic domain, our simulations break translational symmetry and allow the emergence of spatial structure in the CR distribution function. A negative spatial gradient forms across the fully-ionized ISM region while a positive gradient forms across the neutral cloud. We connect our results with CR hydrodynamics formulations by computing the wave-particle scattering rates as predicted by quasilinear, fluid, and Fokker-Planck theory. For momenta where the mean free path is short relative to the box size, we find excellent agreement among all scattering rates. By exploring different cloud sizes and ion-neutral collision rates, we show that our results are robust. Our work provides a first-principles verification of CR hydrodynamics when particles stream down their pressure gradient, and opens a pathway toward comprehensive calibrations of transport coefficients from self-generated Alfv$\acute{\text{e}}$n wave scattering with CRs., Comment: 25 pages, 10 figures, submitted to ApJ

We explore the physics of the gyro-resonant cosmic ray streaming instability (CRSI) including the effects of ion-neutral (IN) damping. This is the main damping mechanism in (partially-ionized) atomic and molecular gas, which are the primary components of the interstellar medium (ISM) by mass. Limitation of CRSI by IN damping is important in setting the amplitude of Alfv\'en waves that scatter cosmic rays and control galactic-scale transport. Our study employs the MHD-PIC hybrid fluid-kinetic numerical technique to follow linear growth as well as post-linear and saturation phases. During the linear phase of the instability -- where simulations and analytical theory are in good agreement -- IN damping prevents wave growth at small and large wavelengths, with the unstable bandwidth lower for higher ion-neutral collision rate $\nu_{\rm in}$. Purely MHD effects during the post-linear phase extend the wave spectrum towards larger $k$. In the saturated state, the cosmic ray distribution evolves toward greater isotropy (lower streaming velocity) by scattering off of Alv\'en waves excited by the instability. In the absence of low-$k$ waves, CRs with sufficiently high momentum are not isotropized. The maximum wave amplitude and rate of isotropization of the distribution function decreases at higher $\nu_{\rm in}$. When the IN damping rate approaches the maximum growth rate of CSRI, wave growth and isotropization is suppressed. Implications of our results for CR transport in partially ionized ISM phases are discussed., Comment: 20 pages, 15 figures, accepted for publication in ApJ

Supernova-driven galactic winds are multiphase streams of gas that are often observed flowing at a range of velocities out of star-forming regions in galaxies. In this study, we use high resolution 3D simulations of multiphase galactic winds modeled with the hydrodynamics code Cholla to investigate the connection between numerical studies and observations. Using a simulated interaction between a hot $T\sim 10^{7}\,\rm K$ supernova-driven wind and a cool $T\sim 10^{4}\,\rm K$ cloud of interstellar material, we create mock observables, including the optical depth and covering fraction of six commonly observed ions (Si II, C II, Si IV, C IV, N V, and O VI) as a function of gas velocity. We compare our mock observables to surveys of galactic winds in the literature, finding good agreement with velocities and profiles of the low ions. We then compute "empirical" values for the optical depth and covering fraction following observational techniques, and compare them to the values calculated directly from the simulation data. We find that the empirically computed values tend to underestimate the "true" value of $\tau$ for ions with high optical depth, and overestimate the "true" value of $\tau$ for ions with low optical depth, relative to the simulated data.

Molecular clouds are supported by turbulence and magnetic fields, but quantifying their influence on cloud lifecycle and star formation efficiency (SFE) remains an open question. We perform radiation MHD simulations of star-forming giant molecular clouds (GMCs) with UV radiation feedback, in which the propagation of UV radiation via ray-tracing is coupled to hydrogen photochemistry. We consider 10 GMC models that vary in either initial virial parameter ($1\le\alpha_{v,0}\le 5$) or dimensionless mass-to-magnetic flux ratio (0.5-8 and $\infty$); the initial mass $10^5M_{\odot}$ and radius 20pc are fixed. Each model is run with five different initial turbulence realizations. In most models, the duration of star formation and the timescale for molecular gas removal (primarily by photoevaporation) are 4-8Myr. Both the final SFE ($\epsilon_*$) and time-averaged SFE per freefall time ($\epsilon_{ff}$) are reduced by strong turbulence and magnetic fields. The median $\epsilon_*$ ranges between 2.1% and 9.5%. The median $\epsilon_{ff}$ ranges between 1.0% and 8.0% and anticorrelates with $\alpha_{v,0}$, in qualitative agreement with previous analytic theory and simulations. However, the time-dependent $\alpha_{v}(t)$ and $\epsilon_{ff,obs}(t)$ based on instantaneous gas properties and cluster luminosity are positively correlated due to rapid evolution, making observational validation of star formation theory difficult. Our median $\epsilon_{ff,obs}(t)\approx$ 2% is similar to observed values. We show that the traditional virial parameter estimates the true gravitational boundedness within a factor of 2 on average, but neglect of magnetic support and velocity anisotropy can sometimes produce large departures. Magnetically subcritical GMCs are unlikely to represent sites of massive star formation given their unrealistic columnar outflows, prolonged lifetime, and low escape fraction of radiation., Comment: 42 pages, 20 figures, 2 tables. Accepted by ApJ, with new additions since v1. For a simulation movie, see https://youtu.be/mjGr4vD6MnA

Supernova (SN) explosions are a major feedback mechanism regulating star formation in galaxies through their momentum input. We review the observations of SNRs in radiative stages in the Milky Way to validate the theoretical results on the momentum/energy injection from a single SN explosion. For seven SNRs where we can observe fast-expanding, atomic radiative shells, we show that the shell momentum inferred from HI 21 cm line observations is in the range of (0.5--4.5)$\times 10^5$ $M_\odot$ km s$^{-1}$. In two SNRs (W44 and IC 443), shocked molecular gas with momentum comparable to that of the atomic SNR shells has been also observed. We compare the momentum and kinetic/thermal energy of these seven SNRs with the results from 1D and 3D numerical simulations. The observation-based momentum and kinetic energy agree well with the expected momentum/energy input from an SN explosion of $\sim 10^{51}$ erg. It is much more difficult to use data/model comparisons of thermal energy to constrain the initial explosion energy, however, due to rapid cooling and complex physics at the hot/cool interface in radiative SNRs. We discuss the observational and theoretical uncertainties of these global parameters and explosion energy estimates for SNRs in complex environments., Comment: 44 pages, 13 figures, accepted for publication in Astrophysical Journal

Galactic outflows have density, temperature, and velocity variations at least as large as that of the multiphase, turbulent interstellar medium (ISM) from which they originate. We have conducted a suite of parsec-resolution numerical simulations using the TIGRESS framework, in which outflows emerge as a consequence of interaction between supernovae (SNe) and the star-forming ISM. The outflowing gas is characterized by two distinct thermal phases, cool (T<10^4 K) and hot (T>10^6 K), with most mass carried by the cool phase and most energy and newly-injected metals carried by the hot phase. Both components have a broad distribution of outflow velocity, and especially for cool gas this implies a varying fraction of escaping material depending on the halo potential. Informed by the TIGRESS results, we develop straightforward analytic formulae for the joint probability density functions (PDFs) of mass, momentum, energy, and metal loading as distributions in outflow velocity and sound speed. The model PDFs have only two parameters, SFR surface density \Sigma_SFR and the metallicity of the ISM, and fully capture the behavior of the original TIGRESS simulation PDFs over \Sigma_SFR~(10^{-4},1)M_sun/kpc^2/yr. Employing PDFs from resolved simulations will enable galaxy formation subgrid model implementations with wind velocity and temperature (as well as total loading factors) that are based on theoretical predictions rather than empirical tuning. This is a critical step to incorporate advances from TIGRESS and other high-resolution simulations in future cosmological hydrodynamics and semi-analytic galaxy formation models. We release a python package to prototype our model and to ease its implementation., Comment: ApJL accepted. For associated python package, see https://twind.readthedocs.io/en/latest/

CO is the most widely used observational tracer of molecular gas. The observable CO luminosity is translated to H_2 mass via a conversion factor, X_CO, which is a source of uncertainty and bias. Despite variations in X_CO, the empirically-determined solar neighborhood value is often applied across different galactic environments. To improve understanding of X_CO, we employ 3D magnetohydrodynamics simulations of the interstellar medium (ISM) in galactic disks with a large range of gas surface densities, allowing for varying metallicity, far-ultraviolet (FUV) radiation, and cosmic ray ionization rate (CRIR). With the TIGRESS simulation framework we model the three-phase ISM with self-consistent star formation and feedback, and post-process outputs with chemistry and radiation transfer to generate synthetic CO(1--0) and (2--1) maps. Our models reproduce the observed CO excitation temperatures, line-widths, and line ratios in nearby disk galaxies. X_CO decreases with increasing metallicity, with a power-law slope of -0.8 for the (1--0) line and -0.5 for the (2--1) line. X_CO also decreases at higher CRIR, and is insensitive to the FUV radiation. As density increases, X_CO first decreases due to increasing excitation temperature, and then increases when the emission is fully saturated. We provide fits between X_CO and observable quantities such as the line ratio, peak antenna temperature, and line brightness, which probe local gas conditions. These fits, which allow for varying beam size, may be used in observations to calibrate out systematic biases. We also provide estimates of the CO-dark H_2 fraction at different gas surface densities, observational sensitivities, and beam sizes., Comment: Accepted by ApJ

Using the PHANGS-ALMA CO (2-1) survey, we characterize molecular gas properties on ${\sim}$100 pc scales across 102,778 independent sightlines in 70 nearby galaxies. This yields the best synthetic view of molecular gas properties on cloud scales across the local star-forming galaxy population obtained to date. Consistent with previous studies, we observe a wide range of molecular gas surface densities (3.4 dex), velocity dispersions (1.7 dex), and turbulent pressures (6.5 dex) across the galaxies in our sample. Under simplifying assumptions about sub-resolution gas structure, the inferred virial parameters suggest that the kinetic energy of the molecular gas typically exceeds its self-gravitational binding energy at ${\sim}$100 pc scales by a modest factor (1.3 on average). We find that the cloud-scale surface density, velocity dispersion, and turbulent pressure (1) increase towards the inner parts of galaxies, (2) are exceptionally high in the centers of barred galaxies (where the gas also appears less gravitationally bound), and (3) are moderately higher in spiral arms than in inter-arm regions. The galaxy-wide averages of these gas properties also correlate with the integrated stellar mass, star formation rate, and offset from the star-forming main sequence of the host galaxies. These correlations persist even when we exclude regions with extraordinary gas properties in galaxy centers, which contribute significantly to the inter-galaxy variations. Our results provide key empirical constraints on the physical link between molecular cloud populations and their galactic environment., Comment: 9 pages + appendices, ApJL in press. Data tables available at https://www.canfar.net/storage/list/phangs/RELEASES/Sun_etal_2020b

Physical and chemical properties of the interstellar medium (ISM) at sub-galactic ($\sim$kpc) scales play an indispensable role in controlling the ability of gas to form stars. As part of the SMAUG (Simulating Multiscale Astrophysics to Understand Galaxies) project, in this paper, we use the TNG50 cosmological simulation to explore the physical parameter space of 8 resolved ISM properties in star-forming regions to constrain the areas of this hyperspace over which most star-forming environments exist. We deconstruct our simulated galaxies spanning a wide range of mass (M$_\star = 10^{7-11}$ M$_\odot$) and redshift ($0 \leq z \leq 3$) into kpc-sized regions, and statistically analyze the gas/stellar surface densities, gas metallicity, vertical stellar velocity dispersion, epicyclic frequency and dark-matter volumetric density representative of each region in the context of their star formation activity and galactic environment (radial galactocentric location). By examining the star formation rate (SFR) weighted distributions of these properties, we show that stars primarily form in two spatially distinct environmental regimes, which are brought about by an underlying bi-component radial SFR surface density profile in galaxies. We examine how the relative prominence of these two regimes depends on host galaxy mass and cosmic time. We also compare our findings with those from integral field spectroscopy observations and achieve a good overall agreement. Further, using dimensionality reduction, we characterise the aforementioned hyperspace to reveal a high-degree of multicollinearity in relationships amongst ISM properties that drive the distribution of star formation at kpc-scales. Based on this, we show that a reduced 3D representation underpinned by a multi-variate radius relationship is sufficient to capture most of the variance in the original 8D space., Comment: Updated intro, minor text and abstract modifications. More information on the SMAUG project and the other first result papers can be found here: https://www.simonsfoundation.org/flatiron/center-for-computational-astrophysics/galaxy-formation/smaug/papersplash1

Large scale outflows in star-forming galaxies are observed to be ubiquitous, and are a key aspect of theoretical modeling of galactic evolution in a cosmological context, the focus of the SMAUG (Simulating Multiscale Astrophysics to Understand Galaxies) project. Gas blown out from galactic disks, similar to gas within galaxies, consists of multiple phases with large contrasts of density, temperature, and other properties. To study multiphase outflows as emergent phenomena, we run a suite of ~pc-resolution local galactic disk simulations using the TIGRESS framework. Explicit modeling of the interstellar medium (ISM), including star formation and self-consistent radiative heating plus supernova feedback, regulates ISM properties and drives the outflow. We investigate the scaling of outflow mass, momentum, energy, and metal loading factors with galactic disk properties, including star formation rate (SFR) surface density (\Sigma_SFR~10^{-4}-1 M_sun/kpc^2/yr), gas surface density (~1-100 M_sun/pc^2), and total midplane pressure (or weight) (~10^3-10^6 k_B cm^{-3} K). The main components of outflowing gas are mass-delivering cool gas (T~10^4 K) and energy/metal-delivering hot gas (T~10^6 K). Cool mass outflow rates measured at outflow launch points (one or two scale heights) are 1-100 times the SFR (decreasing with \Sigma_SFR), although in massive galaxies most mass falls back due to insufficient outflow velocity. The hot galactic outflow carries mass comparable to 10% of the SFR, together with 10-20% of the energy and 30-60% of the metal mass injected by SN feedback. The characteristic outflow velocities of both phases scale very weakly with SFR, as v_out \propto \Sigma_SFR^{0.1~0.2}, consistent with observations. Importantly, our analysis demonstrates that in any physically-motivated cosmological wind model, it is crucial to include at least two distinct thermal wind components., Comment: Accepted for publication in ApJ. More information on the SMAUG project and the other first result papers can be found here: https://www.simonsfoundation.org/flatiron/center-for-computational-astrophysics/galaxy-formation/smaug/papersplash1

It has been hypothesized that photons from young, massive star clusters are responsible for maintaining the ionization of diffuse warm ionized gas seen in both the Milky Way and other disk galaxies. For a theoretical investigation of the warm ionized medium (WIM), it is crucial to solve radiation transfer equations where the ISM and clusters are modeled self-consistently. To this end, we employ a Solar neighborhood model of TIGRESS, a magnetohydrodynamic simulation of the multiphase, star-forming ISM, and post-process the simulation with an adaptive ray tracing method to transfer UV radiation from star clusters. We find that the WIM volume filling factor is highly variable, and sensitive to the rate of ionizing photon production and ISM structure. The mean WIM volume filling factor rises to ~0.15 at |z|~1 kpc. Approximately half of ionizing photons are absorbed by gas and half by dust; the cumulative ionizing photon escape fraction is 1.1%. Our time-averaged synthetic H$\alpha$ line profile matches WHAM observations on the redshifted (outflowing) side, but has insufficient intensity on the blueshifted side. Our simulation matches the Dickey-Lockman neutral density profile well, but only a small fraction of snapshots have high-altitude WIM density consistent with Reynolds Layer estimates. We compute a clumping correction factor C = /sqrt~0.2 that is remarkably constant with distance from the midplane and time; this can be used to improve estimates of ionized gas mass and mean electron density from observed H$\alpha$ surface brightness profiles in edge-on galaxies., Comment: Accepted to ApJ; 31 pages, 17 figures

Spiral arms greatly affect gas flows and star formation in disk galaxies. We use local three-dimensional simulations of the vertically-stratified, self-gravitating, differentially-rotating, interstellar medium (ISM) subject to a stellar spiral potential to study the effects of spiral arms on star formation and formation of arm spurs/feathers. We adopt the TIGRESS framework of Kim & Ostriker (2017) to handle radiative heating and cooling, star formation, and ensuing supernova (SN) feedback. We find that more than 90% of star formation takes place in spiral arms, but the global star formation rate (SFR) in models with spiral arms is enhanced by less than a factor of 2 compared to the no-arm counterpart. This results from a quasi-linear relationship between the SFR surface density Sigma_SFR and the gas surface density Sigma, and supports the picture that spiral arms do not trigger star formation but rather concentrate star-forming regions. Correlated SN feedback produces gaseous spurs/feathers downstream from arms in both magnetized and unmagnetized models. These spurs/feathers are short-lived and have magnetic fields parallel to their length, in contrast to the longer-lived features with perpendicular magnetic fields induced by gravitational instability. SN feedback drives the turbulent component of magnetic fields, with the total magnetic field strength sublinearly proportional to Sigma. The total midplane pressure varies by a factor of ~10 between arm and interarm regions but agrees locally with the total vertical ISM weight, while Sigma_SFR is locally consistent with the prediction of pressure-regulated, feedback-modulated theory., Comment: Accepted for publication in the ApJ

In typical environments of star-forming clouds, converging supersonic turbulence generates shock-compressed regions, and can create strongly-magnetized sheet-like layers. Numerical MHD simulations show that within these post-shock layers, dense filaments and embedded self-gravitating cores form via gathering material along the magnetic field lines. As a result of the preferred-direction mass collection, a velocity gradient perpendicular to the filament major axis is a common feature seen in simulations. We show that this prediction is in good agreement with recent observations from the CARMA Large Area Star Formation Survey (CLASSy), from which we identified several filaments with prominent velocity gradients perpendicular to their major axes. Highlighting a filament from the northwest part of Serpens South, we provide both qualitative and quantitative comparisons between simulation results and observational data. In particular, we show that the dimensionless ratio $C_v \equiv {\Delta v_h}^2/(GM/L)$, where $\Delta v_h$ is half of the observed perpendicular velocity difference across a filament, and $M/L$ is the filament's mass per unit length, can distinguish between filaments formed purely due to turbulent compression and those formed due to gravity-induced accretion. We conclude that the perpendicular velocity gradient observed in the Serpens South northwest filament can be caused by gravity-induced anisotropic accretion of material from a flattened layer. Using synthetic observations of our simulated filaments, we also propose that a density-selection effect may explain observed subfilaments (one filament breaking into two components in velocity space) as reported in Dhabal et al. (2018)., Comment: 11 pages, 7 figures, accepted for publication in MNRAS

A common situation in galactic and intergalactic gas involves cold dense gas in motion relative to hot diffuse gas. Kelvin-Helmholtz instability creates a turbulent mixing layer and populates the intermediate-temperature phase, which often cools rapidly. The energy lost to cooling is balanced by the advection of hot high enthalpy gas into the mixing layer, resulting in growth and acceleration of the cold phase. This process may play a major role in determining the interstellar medium and circumgalactic medium phase structure, and accelerating cold gas in galactic winds and cosmic filaments. Cooling in these mixing layers occurs in a thin corrugated sheet, which we argue has an area with fractal dimension $D=5/2$ and a thickness that adjusts to match the hot phase mixing time to the cooling time. These cooling sheet properties form the basis of a new model for how the cooling rate and hot gas inflow velocity depend on the size $L$, cooling time $t_{\rm cool}$, relative velocity $v_{\rm rel}$, and density contrast $\rho_{\rm cold}/\rho_{\rm hot}$ of the system. Entrainment is expected to be enhanced in environments with short $t_{\rm cool}$, large $v_{\rm rel}$, and large $\rho_{\rm cold}/\rho_{\rm hot}$. Using a large suite of three dimensional hydrodynamic simulations, we demonstrate that this fractal cooling layer model accurately captures the energetics and evolution of turbulent interfaces and can therefore be used as a foundation for understanding multiphase mixing with strong radiative cooling., Comment: 11 pages, 5 figures, submitted to ApJL. Movies can be found here https://dfielding14.github.io/movies/