Search for: All records

The cold neutral medium (CNM) is an important part of the galactic gas cycle and a precondition for the formation of molecular and star-forming gas, yet its distribution is still not fully understood. In this work, we present extremely high resolution simulations of spiral galaxies with time-dependent chemistry such that we can track the formation of the CNM, its distribution within the galaxy, and its correlation with star formation. We find no strong radial dependence between the CNM fraction and total neutral atomic hydrogen (H i) due to the decreasing interstellar radiation field counterbalancing the decreasing gas column density at larger galactic radii. However, the CNM fraction does increase in spiral arms where the CNM distribution is clumpy, rather than continuous, overlapping more closely with H2. The CNM does not extend out radially as far as H i, and the vertical scale height is smaller in the outer galaxy compared to H i with no flaring. The CNM column density scales with total mid-plane pressure and disappears from the gas phase below values of PT/kB = 1000 K cm−3. We find that the star formation rate density follows a similar scaling law with CNM column density to the total gas Kennicutt–Schmidt law. In the outer galaxy, we produce realistic vertical velocity dispersions in the H i purely from galactic dynamics, but our models do not predict CNM at the extremely large radii observed in H i absorption studies of the Milky Way. We suggest that extended spiral arms might produce isolated clumps of CNM at these radii.

 

We explore the impact of diffusive cosmic rays (CRs) on the evolution of the interstellar medium (ISM) under varying assumptions of supernova explosion environment. In practice, we systematically vary the relative fractions of supernovae (SN) occurring in star-forming high-density gas and those occurring in random locations decoupled from star-forming gas to account for SN from run-away stars or explosions in regions that have been cleared by prior SN, stellar winds, or radiation. We find that in the simple system of a periodic stratified gas layer the ISM structure will evolve to one of two solutions: a ‘peak driving’ state where warm gas is volume filling or a ‘thermal runaway’ state where hot gas is volume filling. CR pressure and transport are important factors that strongly influence the solution state the ISM reaches and have the ability to flip the ISM between solutions. Observable signatures such as gamma-ray emission and H i gas are explored. We find that gamma-ray luminosity from pion decay is largely consistent with observations for a range of model parameters. The thickness of the H i gas layer may be too compact, however, this may be due to a large cold neutral fraction of mid-plane gas. The volume fraction of hot gas evolves to stable states in both solutions, but neither settles to a Milky Way-like configuration, suggesting that additional physics omitted here (e.g. a cosmological circumgalactic medium, radiation transport, or spectrally resolved and spatially varying CR transport) may be required.

 

Molecular hydrogen allows cooling in primordial gas, facilitating its collapse into Population III stars within primordial halos. Lyman–Werner (LW) radiation from these stars can escape the halo and delay further star formation by destroying H2 in other halos. As cosmological simulations show that increasing the background LW field strength increases the average halo mass required for star formation, we perform follow-up simulations of selected halos to investigate the knock-on effects this has on the Population III IMF. We follow 5 halos for each of the J21  = 0, 0.01, and 0.1 LW field strengths, resolving the pre-stellar core density of 10−6 g cm−3 (1018 cm−3) before inserting sink particles and following the fragmentation behaviour for hundreds of years further. We find that the mass accreted onto sinks by the end of the simulations is proportional to the mass within the ∼10−2 pc molecular core, which is not correlated to the initial mass of the halo. As such, the IMFs for masses above the brown dwarf limit show little dependence on the LW strength, although they do show variance in the number of low-mass clumps formed. As the range of background LW field strengths tested here covers the most likely values from literature, we conclude that the IMF for so-called Pop III.2 stars is not significantly different from the initial population of Pop III.1 stars. The primordial IMF therefore likely remains unchanged until the formation of the next generation of Population II stars.

 

Determining how the 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₂conversion factor (α_CO) that impact the assessment of the gas column densities and thus of the SFE. To address such issues, we investigate the dependence ofα_COon the local CO velocity dispersion at 150 pc scales using a new set of dust-basedα_COmeasurements and propose a newα_COprescription 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–Atacama Large Millimeter/submillimeter Array survey. We find increasing SFE toward high-surface-density regions like galaxy centers, while using a constant or metallicity-basedα_COresults 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 three to four times shorter than in nonbarred galaxy centers or the disks. Across the galaxy disks, the depletion time is consistently around 2–3 Gyr, regardless of the choice ofα_COprescription. 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 to an enhanced SFE compared to nonbarred centers or disk regions.

 

The Central Molecular Zone (the central ∼500 pc of the Milky Way) hosts molecular clouds in an extreme environment of strong shear, high gas pressure and density, and complex chemistry. G0.253+0.016, also known as ‘the Brick’, is the densest, most compact, and quiescent of these clouds. High-resolution observations with the Atacama Large Millimetre/submillimetre Array (ALMA) have revealed its complex, hierarchical structure. In this paper we compare the properties of recent hydrodynamical simulations of the Brick to those of the ALMA observations. To facilitate the comparison, we post-process the simulations and create synthetic ALMA maps of molecular line emission from eight molecules. We correlate the line emission maps to each other and to the mass column density and find that HNCO is the best mass tracer of the eight emission lines within the simulations. Additionally, we characterize the spatial structure of the observed and simulated cloud using the density probability distribution function (PDF), spatial power spectrum, fractal dimension, and moments of inertia. While we find good agreement between the observed and simulated data in terms of power spectra and fractal dimensions, there are key differences in the density PDFs and moments of inertia, which we attribute to the omission of magnetic fields in the simulations. This demonstrates that the presence of the Galactic potential can reproduce many cloud properties, but additional physical processes are needed to fully explain the gas structure.

 

Many studies concluded that magnetic fields suppress star formation in molecular clouds and Milky Way like galaxies. However, most of these studies are based on fully developed fields that have reached the saturation level, with little work on investigating how an initial weak primordial field affects star formation in low metallicity environments. In this paper, we investigate the impact of a weak initial field on low metallicity dwarf galaxies. We perform high-resolution arepo simulations of five isolated dwarf galaxies. Two models are hydrodynamical, two start with a primordial magnetic field of 10$^{-6} \, \mu$G and different sub-solar metallicities, and one starts with a saturated field of 10$^{-2} \, \mu$G. All models include a non-equilibrium, time-dependent chemical network that includes the effects of gas shielding from the ambient ultraviolet field. Sink particles form directly from the gravitational collapse of gas and are treated as star-forming clumps that can accrete gas. We vary the ambient uniform far ultraviolet field, and cosmic ray ionization rate between 1 per cent and 10 per cent of solar values. We find that the magnetic field has little impact on the global star formation rate (SFR), which is in tension with some previously published results. We further find that the initial field strength has little impact on the global SFR. We show that an increase in the mass fractions of both molecular hydrogen and cold gas, along with changes in the perpendicular gas velocity dispersion and the magnetic field acting in the weak-field model, overcome the expected suppression in star formation.

 

Abstract The CO-to-H 2 conversion factor ( α CO ) is central to measuring the amount and properties of molecular gas. It is known to vary with environmental conditions, and previous studies have revealed lower α CO in the centers of some barred galaxies on kiloparsec scales. To unveil the physical drivers of such variations, we obtained Atacama Large Millimeter/submillimeter Array bands (3), (6), and (7) observations toward the inner ∼2 kpc of NGC 3627 and NGC 4321 tracing 12 CO, 13 CO, and C 18 O lines on ∼100 pc scales. Our multiline modeling and Bayesian likelihood analysis of these data sets reveal variations of molecular gas density, temperature, optical depth, and velocity dispersion, which are among the key drivers of α CO . The central 300 pc nuclei in both galaxies show strong enhancement of temperature T k ≳ 100 K and density n H 2 > 10 3 cm −3 . Assuming a CO-to-H 2 abundance of 3 × 10 −4 , we derive 4–15 times lower α CO than the Galactic value across our maps, which agrees well with previous kiloparsec-scale measurements. Combining the results with our previous work on NGC 3351, we find a strong correlation of α CO with low- J 12 CO optical depths ( τ CO ), as well as an anticorrelation with T k . The τ CO correlation explains most of the α CO variation in the three galaxy centers, whereas changes in T k influence α CO to second order. Overall, the observed line width and 12 CO/ 13 CO 2–1 line ratio correlate with τ CO variation in these centers, and thus they are useful observational indicators for α CO variation. We also test current simulation-based α CO prescriptions and find a systematic overprediction, which likely originates from the mismatch of gas conditions between our data and the simulations. 

The Central Molecular Zone (CMZ; the central ∼500 pc of the Galaxy) is a kinematically unusual environment relative to the Galactic disc, with high-velocity dispersions and a steep size–linewidth relation of the molecular clouds. In addition, the CMZ region has a significantly lower star formation rate (SFR) than expected by its large amount of dense gas. An important factor in explaining the low SFR is the turbulent state of the star-forming gas, which seems to be dominated by rotational modes. However, the turbulence driving mechanism remains unclear. In this work, we investigate how the Galactic gravitational potential affects the turbulence in CMZ clouds. We focus on the CMZ cloud G0.253+0.016 (‘the Brick’), which is very quiescent and unlikely to be kinematically dominated by stellar feedback. We demonstrate that several kinematic properties of the Brick arise naturally in a cloud-scale hydrodynamics simulation, that takes into account the Galactic gravitational potential. These properties include the line-of-sight velocity distribution, the steepened size–linewidth relation, and the predominantly solenoidal nature of the turbulence. Within the simulation, these properties result from the Galactic shear in combination with the cloud’s gravitational collapse. This is a strong indication that the Galactic gravitational potential plays a crucial role in shaping the CMZ gas kinematics, and is a major contributor to suppressing the SFR, by inducing predominantly solenoidal turbulent modes.

 

Galactic bars can drive cold gas inflows towards the centres of galaxies. The gas transport happens primarily through the so-called bar dust lanes, which connect the galactic disc at kpc scales to the nuclear rings at hundreds of pc scales much like two gigantic galactic rivers. Once in the ring, the gas can fuel star formation activity, galactic outflows, and central supermassive black holes. Measuring the mass inflow rates is therefore important to understanding the mass/energy budget and evolution of galactic nuclei. In this work, we use CO datacubes from the PHANGS-ALMA survey and a simple geometrical method to measure the bar-driven mass inflow rate on to the nuclear ring of the barred galaxy NGC 1097. The method assumes that the gas velocity in the bar lanes is parallel to the lanes in the frame co-rotating with the bar, and allows one to derive the inflow rates from sufficiently sensitive and resolved position–position–velocity diagrams if the bar pattern speed and galaxy orientations are known. We find an inflow rate of $\dot{M}=(3.0 \pm 2.1)\, \rm M_\odot \, yr^{-1}$ averaged over a time span of 40 Myr, which varies by a factor of a few over time-scales of ∼10 Myr. Most of the inflow appears to be consumed by star formation in the ring, which is currently occurring at a star formation rate (SFR) of $\simeq\!1.8\!-\!2 \, \rm M_\odot \, yr^{-1}$, suggesting that the inflow is causally controlling the SFR in the ring as a function of time.

 

Abstract The Galactic bar plays a critical role in the evolution of the Milky Way’s Central Molecular Zone (CMZ), driving gas toward the Galactic Center via gas flows known as dust lanes. To explore the interaction between the CMZ and the dust lanes, we run hydrodynamic simulations in arepo , modeling the potential of the Milky Way’s bar in the absence of gas self-gravity and star formation physics, and we study the flows of mass using Monte Carlo tracer particles. We estimate the efficiency of the inflow via the dust lanes, finding that only about a third (30% ± 12%) of the dust lanes’ mass initially accretes onto the CMZ, while the rest overshoots and accretes later. Given observational estimates of the amount of gas within the Milky Way’s dust lanes, this suggests that the true total inflow rate onto the CMZ is 0.8 ± 0.6 M ⊙ yr −1 . Clouds in this simulated CMZ have sudden peaks in their average density near the apocenter, where they undergo violent collisions with inflowing material. While these clouds tend to counter-rotate due to shear, co-rotating clouds occasionally occur due to the injection of momentum from collisions with inflowing material (∼52% are strongly counter-rotating, and ∼7% are strongly co-rotating of the 44 cloud sample). We investigate the formation and evolution of these clouds, finding that they are fed by many discrete inflow events, providing a consistent source of gas to CMZ clouds even as they collapse and form stars.