Search

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 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., Comment: 7 pages, 8 figures; accepted to MNRAS (July 24th 2023)

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

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

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

We present a novel, physically-motivated sub-grid model for HII region feedback within the moving mesh code Arepo, accounting for both the radiation pressure-driven and thermal expansion of the ionised gas surrounding young stellar clusters. We apply this framework to isolated disc galaxy simulations with mass resolutions between $10^3~{\rm M}_\odot$ and $10^5~{\rm M}_\odot$ per gas cell. Each simulation accounts for the self-gravity of the gas, the momentum and thermal energy from supernovae, the injection of mass by stellar winds, and the non-equilibrium chemistry of hydrogen, carbon and oxygen. We reduce the resolution-dependence of our model by grouping those HII regions with overlapping ionisation front radii. The Str\"{o}mgren radii of the grouped HII regions are at best marginally-resolved, so that the injection of purely-thermal energy within these radii has no effect on the interstellar medium. By contrast, the injection of momentum increases the fraction of cold and molecular gas by more than 50 per cent at mass resolutions of $10^3~{\rm M}_\odot$, and decreases its turbulent velocity dispersion by $\sim 10~{\rm kms}^{-1}$. The mass-loading of galactic outflows is decreased by an order of magnitude. The characteristic lifetime of the least-massive molecular clouds ($M/{\rm M}_\odot < 5.6 \times 10^4$) is reduced from $\sim 18$ Myr to $<10$ Myr, indicating that HII region feedback is effective in destroying these clouds. Conversely, the lifetimes of intermediate-mass clouds ($5.6 \times 10^4 < M/{\rm M}_\odot < 5 \times 10^5$) are elongated by $\sim 7$ Myr, likely due to a reduction in supernova clustering. The derived cloud lifetimes span the range from $10$-$40$ Myr, in agreement with observations. All results are independent of whether the momentum is injected from a 'spherical' or a 'blister-type HII region., Comment: Accepted for publication in MNRAS. Comments welcome!

The centre of the Milky Way is the site of several high-energy processes that have strongly impacted the inner regions of our Galaxy. Activity from the super-massive black hole, Sgr A*, and/or stellar feedback from the inner molecular ring expel matter and energy from the disc in the form of a galactic wind. Multiphase gas has been observed within this outflow, from hot highly-ionized, to warm ionized and cool atomic gas. To date, however, there has been no evidence of the cold and dense molecular phase. Here we report the first detection of molecular gas outflowing from the centre of our Galaxy. This cold material is associated with atomic hydrogen clouds travelling in the nuclear wind. The morphology and the kinematics of the molecular gas, resolved on ~1 pc scale, indicate that these clouds are mixing with the warmer medium and are possibly being disrupted. The data also suggest that the mass of molecular gas driven out is not negligible and could impact the rate of star formation in the central regions. The presence of this cold, dense, high-velocity gas is puzzling, as neither Sgr A* at its current level of activity, nor star formation in the inner Galaxy seem viable sources for this material., Comment: Published in the August 19 issue of Nature. This is the authors' version before final edits. Published version is available at http://www.nature.com/articles/s41586-020-2595-z

We use the hydrodynamical simulation of our inner Galaxy presented in Armillotta et al. (2019) to study the gas distribution and kinematics within the CMZ. We use a resolution high enough to capture the gas emitting in dense molecular tracers such as NH3 and HCN, and simulate a time window of 50 Myr, long enough to capture phases during which the CMZ experiences both quiescent and intense star formation. We then post-process the simulated CMZ to calculate its spatially-dependent chemical and thermal state, producing synthetic emission data cubes and maps of both HI and the molecular gas tracers CO, NH3 and HCN. We show that, as viewed from Earth, gas in the CMZ is distributed mainly in two parallel and elongated features extending from positive longitudes and velocities to negative longitudes and velocities. The molecular gas emission within these two streams is not uniform, and it is mostly associated to the region where gas flowing towards the Galactic Center through the dust lanes collides with gas orbiting within the ring. Our simulated data cubes reproduce a number of features found in the observed CMZ. However, some discrepancies emerge when we use our results to interpret the position of individual molecular clouds. Finally, we show that, when the CMZ is near a period of intense star formation, the ring is mostly fragmented as a consequence of supernova feedback, and the bulk of the emission comes from star-forming molecular clouds. This correlation between morphology and star formation rate should be detectable in observations of extragalactic CMZs., Comment: 19 pages, 11 figures, accepted for publication in MNRAS

Galactic outflows produced by stellar feedback are known to be multiphase in nature. Both observations and simulations indicate that the material within several kpc of galactic disk mid-planes consists of warm clouds embedded within a hot wind. A theoretical understanding of the outflow phenomenon, including both winds and fountain flows, requires study of the interactions among thermal phases. We develop a method to quantify these interactions via measurements of mass, momentum, and energy flux exchanges using temporally and spatially averaged quantities and conservation laws. We apply this method to a star-forming ISM MHD simulation based on the TIGRESS framework, for Solar neighbourhood conditions. To evaluate the extent of interactions among the phases, we first examine the validity of the ``ballistic model,'' which predicts trajectories of the warm phase ($5050\,\rm{K}5\times10^5\,\rm{K}$) to the warm phase. The large energy flux from the hot outflow that is transferred to the warm and intermediate phases is quickly radiated away. A simple interaction model implies an effective warm cloud size in the fountain flow of a few 100~pc, showing that warm-hot flux exchange mainly involves a few large clouds rather than many small ones., Comment: 22 pages, Submitted to ApJ

Ongoing surveys are in the process of measuring the chemical abundances in large numbers of stars, with the ultimate goal of reconstructing the formation history of the Milky Way using abundances as tracers. However, interpretation of these data requires that we understand the relationship between stellar distributions in chemical and physical space, i.e., how similar in chemical abundance do we expect a pair of stars to be as a function of the distance between their formation sites. We investigate this question by simulating the gravitational collapse of a turbulent molecular cloud extracted from a galaxy-scale simulation, seeded with chemical inhomogeneities with different initial spatial scales. We follow the collapse from galactic scales down to resolutions scales of $\approx 10^{-3}$ pc, and find that, during this process, turbulence mixes the metal patterns, reducing the abundance scatter initially present in the gas by an amount that depends on the initial scale of inhomogeneity of each metal field. However, we find that regardless of the initial spatial structure of the metals at the onset of collapse, the final stellar abundances are highly correlated on distances below a few pc, and nearly uncorrelated on larger distances. Consequently, the star formation process defines a natural size scale of $\sim 1$ pc for chemically-homogenous star clusters, suggesting that any clusters identified as homogenous in chemical space must have formed within $\sim 1$ pc of one another. However, in order to distinguish different star clusters in chemical space, observations across multiple elements will be required, and the elements that are likely to be most efficient at separating distinct clusters in chemical space are those whose correlation length in the ISM is of order tens of pc, comparable to the sizes of individual molecular clouds., Comment: 15 pages, 10 figures; submitted to MNRAS