Search for: All records

The interstellar medium is threaded by a hierarchy of filaments from large scales (∼100 pc) to small scales (∼0.1 pc). The masses and lengths of these nested structures may reveal important constraints for cloud formation and evolution, but it is difficult to investigate from an evolutionary perspective using single observations. In this work, we extract simulated molecular clouds from the ‘Cloud Factory’ galactic-scale ISM suite in combination with 3D Monte Carlo radiative transfer code polaris to investigate how filamentary structure evolves over time. We produce synthetic dust continuum observations in three regions with a series of snapshots and use the filfinder algorithm to identify filaments in the dust derived column density maps. When the synthetic filaments mass and length are plotted on an mass–length (M–L) plot, we see a scaling relation of L ∝ M0.45 similar to that seen in observations, and find that the filaments are thermally supercritical. Projection effects systematically affect the masses and lengths measured for the filaments, and are particularly severe in crowded regions. In the filament M–L diagram we identify three main evolutionary mechanisms: accretion, segmentation, and dispersal. In particular we find that the filaments typically evolve from smaller to larger masses in the observational M–L plane, indicating the dominant role of accretion in filament evolution. Moreover, we find a potential correlation between line mass and filament growth rate. Once filaments are actively star forming they then segment into smaller sections, or are dispersed by internal or external forces.

 

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.

 

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.

 

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. 

Abstract We present an update to the framework called Simulator of Galaxy Millimeter/submillimeter Emission ( sígame ). sígame derives line emission in the far-infrared (FIR) for galaxies in particle-based cosmological hydrodynamics simulations by applying radiative transfer and physics recipes via a postprocessing step after completion of the simulation. In this version, a new technique is developed to model higher gas densities by parameterizing the probability distribution function (PDF) of the gas density in higher-resolution simulations run with the pseudo-Lagrangian, Voronoi mesh code arepo . The parameterized PDFs are used as a look-up table, and reach higher densities than in previous work. sígame v3 is tested on redshift z = 0 galaxies drawn from the simba cosmological simulation for eight FIR emission lines tracing vastly different phases of the interstellar medium. This version of sígame includes dust radiative transfer with S kirt and high-resolution photoionization models with C loudy , the latter sampled according to the density PDF of the arepo simulations to augment the densities in the cosmological simulation. The quartile distributions of the predicted line luminosities overlap with the observed range for nearby galaxies of similar star formation rate (SFR) for all but two emission lines: [O i ]63 and CO(3–2), which are overestimated by median factors of 1.3 and 1.0 dex, respectively, compared to the observed line–SFR relation of mixed-type galaxies. We attribute the remaining disagreement with observations to the lack of precise attenuation of the interstellar light on sub-grid scales (≲200 pc) and differences in sample selection. 

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 Star formation primarily occurs in filaments where magnetic fields are expected to be dynamically important. The largest and densest filaments trace the spiral structure within galaxies. Over a dozen of these dense (∼10 4 cm −3 ) and long (>10 pc) filaments have been found within the Milky Way, and they are often referred to as “bones.” Until now, none of these bones has had its magnetic field resolved and mapped in its entirety. We introduce the SOFIA legacy project FIELDMAPS which has begun mapping ∼10 of these Milky Way bones using the HAWC+ instrument at 214 μ m and 18.″2 resolution. Here we present a first result from this survey on the ∼60 pc long bone G47. Contrary to some studies of dense filaments in the Galactic plane, we find that the magnetic field is often not perpendicular to the spine (i.e., the center line of the bone). Fields tend to be perpendicular in the densest areas of active star formation and more parallel or random in other areas. The average field is neither parallel nor perpendicular to the Galactic plane or the bone. The magnetic field strengths along the spine typically vary from ∼20 to ∼100 μ G. Magnetic fields tend to be strong enough to suppress collapse along much of the bone, but for areas that are most active in star formation, the fields are notably less able to resist gravitational collapse.