Search for: All records

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

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

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. 

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

Abstract The formation of globular clusters and their relation to the distribution of dark matter have long puzzled astronomers. One of the most recently proposed globular cluster formation channels ties ancient star clusters to the large-scale streaming velocity of baryons relative to dark matter in the early universe. These streaming velocities affect the global infall of baryons into dark matter halos, the high-redshift halo mass function, and the earliest generations of stars. In some cases, streaming velocities may result in dense regions of dark matter-free gas that becomes Jeans unstable, potentially leading to the formation of compact star clusters. We investigate this hypothesis using cosmological hydrodynamical simulations that include a full chemical network and the formation and destruction of H₂, a process crucial for the formation of the first stars. We find that high-density gas in regions with significant streaming velocities is indeed somewhat offset from the centers of dark matter halos, but this offset is typically significantly smaller than the virial radius. Gas outside of dark matter halos never reaches Jeans-unstable densities in our simulations. We postulate that low-level (Z≈ 10⁻³Z_⊙) metal enrichment by Population III supernovae may enable cooling in the extra-virial regions, allowing gas outside of dark matter halos to cool to the cosmic microwave background temperature and become Jeans unstable. Follow-up simulations that include both streaming velocities and metal enrichment by Population III supernovae are needed to understand if streaming velocities provide one path for the formation of globular clusters in the early universe. 

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

ABSTRACT To investigate how molecular clouds react to different environmental conditions at a galactic scale, we present a catalogue of giant molecular clouds (GMCs) resolved down to masses of ∼10 M⊙ from a simulation of the entire disc of an interacting M51-like galaxy and a comparable isolated galaxy. Our model includes time-dependent gas chemistry, sink particles for star formation, and supernova feedback, meaning we are not reliant on star formation recipes based on threshold densities and can follow the physics of the cold molecular phase. We extract GMCs from the simulations and analyse their properties. In the disc of our simulated galaxies, spiral arms seem to act merely as snowplows, gathering gas, and clouds without dramatically affecting their properties. In the centre of the galaxy, on the other hand, environmental conditions lead to larger, more massive clouds. While the galaxy interaction has little effect on cloud masses and sizes, it does promote the formation of counter-rotating clouds. We find that the identified clouds seem to be largely gravitationally unbound at first glance, but a closer analysis of the hierarchical structure of the molecular interstellar medium shows that there is a large range of virial parameters with a smooth transition from unbound to mostly bound for the densest structures. The common observation that clouds appear to be virialized entities may therefore be due to CO bright emission highlighting a specific level in this hierarchical binding sequence. The small fraction of gravitationally bound structures found suggests that low galactic star formation efficiencies may be set by the process of cloud formation and initial collapse. 

ABSTRACT We use hydrodynamical simulations to study the Milky Way’s central molecular zone (CMZ). The simulations include a non-equilibrium chemical network, the gas self-gravity, star formation, and supernova feedback. We resolve the structure of the interstellar medium at sub-parsec resolution while also capturing the interaction between the CMZ and the bar-driven large-scale flow out to $$R\sim 5\, {\rm kpc}$$. Our main findings are as follows: (1) The distinction between inner (R ≲ 120 pc) and outer (120 ≲ R ≲ 450 pc) CMZ that is sometimes proposed in the literature is unnecessary. Instead, the CMZ is best described as single structure, namely a star-forming ring with outer radius R ≃ 200 pc which includes the 1.3° complex and which is directly interacting with the dust lanes that mediate the bar-driven inflow. (2) This accretion can induce a significant tilt of the CMZ out of the plane. A tilted CMZ might provide an alternative explanation to the ∞-shaped structure identified in Herschel data by Molinari et al. (3) The bar in our simulation efficiently drives an inflow from the Galactic disc (R ≃ 3 kpc) down to the CMZ (R ≃ 200 pc) of the order of $$1\rm \, M_\odot \, yr^{-1}$$, consistent with observational determinations. (4) Supernova feedback can drive an inflow from the CMZ inwards towards the circumnuclear disc of the order of $${\sim}0.03\, \rm M_\odot \, yr^{-1}$$. (5) We give a new interpretation for the 3D placement of the 20 and 50 km s−1 clouds, according to which they are close (R ≲ 30 pc) to the Galactic Centre, but are also connected to the larger scale streams at R ≳ 100 pc.