skip to main content

Title: Ubiquitous velocity fluctuations throughout the molecular interstellar medium
The density structure of the interstellar medium determines where stars form and release energy, momentum and heavy elements, driving galaxy evolution1-4. Density variations are seeded and amplified by gas motion, but the exact nature of this motion is unknown across spatial scales and galactic environments5. Although dense star-forming gas probably emerges from a combination of instabilities6,7, convergent flows8 and turbulence9, establishing the precise origin is challenging because it requires gas motion to be quantified over many orders of magnitude in spatial scale. Here we measure10-12 the motion of molecular gas in the Milky Way and in nearby galaxy NGC 4321, assembling observations that span a spatial dynamic range 10-1-103 pc. We detect ubiquitous velocity fluctuations across all spatial scales and galactic environments. Statistical analysis of these fluctuations indicates how star-forming gas is assembled. We discover oscillatory gas flows with wavelengths ranging from 0.3-400 pc. These flows are coupled to regularly spaced density enhancements that probably form via gravitational instabilities13,14. We also identify stochastic and scale-free velocity and density fluctuations, consistent with the structure generated in turbulent flows9. Our results demonstrate that the structure of the interstellar medium cannot be considered in isolation. Instead, its formation and evolution are controlled by more » nested, interdependent flows of matter covering many orders of magnitude in spatial scale. « less
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more » ; ; ; ; ; ; « less
Award ID(s):
Publication Date:
Journal Name:
Nature Astronomy
Sponsoring Org:
National Science Foundation
More Like this
  1. Context. Molecular filaments and hubs have received special attention recently thanks to new studies showing their key role in star formation. While the (column) density and velocity structures of both filaments and hubs have been carefully studied, their magnetic field (B-field) properties have yet to be characterized. Consequently, the role of B-fields in the formation and evolution of hub-filament systems is not well constrained. Aims. We aim to understand the role of the B-field and its interplay with turbulence and gravity in the dynamical evolution of the NGC 6334 filament network that harbours cluster-forming hubs and high-mass star formation. Methods. We present new observations of the dust polarized emission at 850 μ m toward the 2 pc × 10 pc map of NGC 6334 at a spatial resolution of 0.09 pc obtained with the James Clerk Maxwell Telescope (JCMT) as part of the B-field In STar-forming Region Observations (BISTRO) survey. We study the distribution and dispersion of the polarized intensity ( PI ), the polarization fraction ( PF ), and the plane-of-the-sky B-field angle ( χ B_POS ) toward the whole region, along the 10 pc-long ridge and along the sub-filaments connected to the ridge and the hubs. We derivedmore »the power spectra of the intensity and χ B POS along the ridge crest and compared them with the results obtained from simulated filaments. Results. The observations span ~3 orders of magnitude in Stokes I and PI and ~2 orders of magnitude in PF (from ~0.2 to ~ 20%). A large scatter in PI and PF is observed for a given value of I . Our analyses show a complex B-field structure when observed over the whole region (~ 10 pc); however, at smaller scales (~1 pc), χ B POS varies coherently along the crests of the filament network. The observed power spectrum of χ B POS can be well represented with a power law function with a slope of − 1.33 ± 0.23, which is ~20% shallower than that of I . We find that this result is compatible with the properties of simulated filaments and may indicate the physical processes at play in the formation and evolution of star-forming filaments. Along the sub-filaments, χ B POS rotates frombeing mostly perpendicular or randomly oriented with respect to the crests to mostly parallel as the sub-filaments merge with the ridge and hubs. This variation of the B-field structure along the sub-filaments may be tracing local velocity flows of infalling matter in the ridge and hubs. Our analysis also suggests a variation in the energy balance along the crests of these sub-filaments, from magnetically critical or supercritical at their far ends to magnetically subcritical near the ridge and hubs. We also detect an increase in PF toward the high-column density ( N H 2 ≳ 10 23  cm −2 ) star cluster-forming hubs. These latter large PF values may be explained by the increase in grain alignment efficiency due to stellar radiation from the newborn stars, combined with an ordered B-field structure. Conclusions. These observational results reveal for the first time the characteristics of the small-scale (down to ~ 0.1 pc) B-field structure of a 10 pc-long hub-filament system. Our analyses show variations in the polarization properties along the sub-filaments that may be tracing the evolution of their physical properties during their interaction with the ridge and hubs. We also detect an impact of feedback from young high-mass stars on the local B-field structure and the polarization properties, which could put constraints on possible models for dust grain alignment and provide important hints as to the interplay between the star formation activity and interstellar B-fields.« less
  2. The spatial decorrelation of dense molecular gas and young stars observed on ≲ 1 kiloparsec scales in nearby galaxies indicates rapid dispersal of star-forming regions by stellar feedback. We explore the sensitivity of this decorrelation to different processes controlling the structure of the interstellar medium, the abundance of molecular gas, star formation, and feedback in a suite of simulations of an isolated dwarf galaxy with structural properties similar to NGC300 that self-consistently model radiative transfer and molecular chemistry. Our fiducial simulation reproduces the magnitude of decorrelation and its scale dependence measured in NGC300, and we show that this agreement is due to different aspects of feedback, including H2 dissociation, gas heating by the locally variable UV field, early mechanical feedback, and supernovae. In particular, early radiative and mechanical feedback affect the correlation on ≲100 pc scales, while supernovae play a significant role on ≳100 pc scales. The correlation is also sensitive to the choice of the local star formation efficiency per freefall time, eps_ff, which provides a strong observational constraint on eps_ff when the global star formation rate is independent of its value. Finally, we explicitly show that the degree of correlation between the peaks of molecular gas and starmore »formation density is directly related to the distribution of the lifetimes of star-forming regions.« less
  3. Context. The interaction between active galactic nuclei (AGNs) and their host galaxies is scarcely resolved. Narrow-line Seyfert 1 (NLS1) galaxies are believed to represent AGN at early stages of their evolution and to allow one to observe feeding and feedback processes at high black hole accretion rates. Aims. We aim to constrain the properties of the ionised gas outflow in Mrk 1044, a nearby super-Eddington accreting NLS1. Based on the outflow energetics and the associated timescales, we estimate the outflow’s future impact on the ongoing host galaxy star formation on different spatial scales. Methods. We applied a spectroastrometric analysis to observations of Mrk 1044’s nucleus obtained with the adaptive-optics-assisted narrow field mode of the VLT/MUSE instrument. This allowed us to map two ionised gas outflows traced by [O  III ], which have velocities of −560 ± 20 km s −1 and −144 ± 5 km s −1 . Furthermore, we used an archival spectrum from the Space Telescope Imaging Spectrograph on HST to identify two Ly- α absorbing components that escape from the centre with approximately twice the velocity of the ionised gas components. Results. Both [O  III ] outflows are spatially unresolved and located close to the AGN (< 1 pc). They havemore »gas densities higher than 10 5 cm −3 , which implies that the BPT diagnostic cannot be used to constrain the underlying ionisation mechanism. We explore whether an expanding shell model can describe the velocity structure of Mrk 1044’s multi-phase outflow. In the ionised gas emission, an additional outflowing component, which is spatially resolved, is present. It has a velocity of −211 ± 22 km s −1 and a projected size of 4.6 ± 0.6 pc. Our kinematic analysis suggests that significant turbulence is present in the interstellar medium around the nucleus, which may lead to a condensation rain, potentially explaining the efficient feeding of Mrk 1044’s AGN. Within the innermost 0.5″ (160 pc), we detect modest star formation hidden by the beam-smeared emission from the outflow. Conclusions. We estimate that the multi-phase outflow was launched < 10 4 yr ago. Together with the star formation in the vicinity of the nucleus, this suggests that Mrk 1044’s AGN phase started only recently. The outflow carries enough mass and energy to impact the host galaxy star formation on different spatial scales, highlighting the complexity of the AGN feeding and feedback cycle in its early stages.« less

    Formation of supermassive black holes (BHs) remains a theoretical challenge. In many models, especially beginning from stellar relic ‘seeds,’ this requires sustained super-Eddington accretion. While studies have shown BHs can violate the Eddington limit on accretion disc scales given sufficient ‘fuelling’ from larger scales, what remains unclear is whether or not BHs can actually capture sufficient gas from their surrounding interstellar medium (ISM). We explore this in a suite of multiphysics high-resolution simulations of BH growth in magnetized, star-forming dense gas complexes including dynamical stellar feedback from radiation, stellar mass-loss, and supernovae, exploring populations of seeds with masses $\sim 1\!-\!10^{4}\, \mathrm{M}_{\odot }$. In this initial study, we neglect feedback from the BHs: so this sets a strong upper limit to the accretion rates seeds can sustain. We show that stellar feedback plays a key role. Complexes with gravitational pressure/surface density below $\sim 10^{3}\, \mathrm{M}_{\odot }\, {\rm pc^{-2}}$ are disrupted with low star formation efficiencies so provide poor environments for BH growth. But in denser cloud complexes, early stellar feedback does not rapidly destroy the clouds but does generate strong shocks and dense clumps, allowing $\sim 1{{\ \rm per\ cent}}$ of randomly initialized seeds to encounter a dense clump withmore »low relative velocity and produce runaway, hyper-Eddington accretion (growing by orders of magnitude). Remarkably, mass growth under these conditions is almost independent of initial BH mass, allowing rapid intermediate-mass black hole (IMBH) formation even for stellar-mass seeds. This defines a necessary (but perhaps not sufficient) set of criteria for runaway BH growth: we provide analytic estimates for the probability of runaway growth under different ISM conditions.

    « less
  5. Abstract

    Understanding the interplay of stellar feedback and turbulence in the interstellar medium (ISM) is essential to modeling the evolution of galaxies. To determine the timescales over which stellar feedback drives turbulence in the ISM, we performed a spatially resolved, multiwavelength study of the nearby star-forming dwarf galaxy UGC 4305. As indicators of turbulence on local scales (400 pc), we utilized ionized gas velocity dispersion derived from IFU Hαobservations and atomic gas velocity dispersion and energy surface densities derived from Hisynthesis observations with the Very Large Array. These indicators of turbulence were tested against star formation histories over the past 560 Myr derived from color–magnitude diagrams using Spearman’s rank correlation coefficient. The strongest correlation identified at the 400 pc scale is between measures of Hiturbulence and star formation 70–140 Myr ago. We repeated our analysis of UGC 4305's current turbulence and past star formation activity on multiple physical scales (∼560 and 800 pc) to determine whether there are indications of changes in the correlation timescale with changes to the physical scale. No notable correlations were found at larger physical scales, emphasizing the importance of analyzing star formation-driven turbulence as a local phenomenon.