skip to main content

Title: How do velocity structure functions trace gas dynamics in simulated molecular clouds?
Context. Supersonic disordered flows accompany the formation and evolution of molecular clouds (MCs). It has been argued that this is turbulence that can support against gravitational collapse and form hierarchical sub-structures. Aims. We examine the time evolution of simulated MCs to investigate: What physical process dominates the driving of turbulent flows? How can these flows be characterised? Are they consistent with uniform turbulence or gravitational collapse? Do the simulated flows agree with observations? Methods. We analysed three MCs that have formed self-consistently within kiloparsec-scale numerical simulations of the interstellar medium (ISM). The simulated ISM evolves under the influence of physical processes including self-gravity, stratification, magnetic fields, supernova-driven turbulence, and radiative heating and cooling. We characterise the flows using velocity structure functions (VSFs) with and without density weighting or a density cutoff, and computed in one or three dimensions. However, we do not include optical depth effects that can hide motions in the densest gas, limiting comparison of our results with observations. Results. In regions with sufficient resolution, the density-weighted VSFs initially appear to follow the expectations for uniform turbulence, with a first-order power-law exponent consistent with Larson’s size-velocity relationship. Supernova blast wave impacts on MCs produce short-lived coherent motions at large scales, increasing the scaling exponents for a crossing time. Gravitational contraction drives small-scale motions, producing scaling coefficients that drop or even turn negative as small scales become dominant. Removing the density weighting eliminates this effect as it emphasises the diffuse ISM. Conclusions. We conclude that two different effects coincidentally reproduce Larson’s size velocity relationship. Initially, uniform turbulence dominates, so the energy cascade produces VSFs that are consistent with Larson’s relationship. Later, contraction dominates and the density-weighted VSFs become much shallower or even inverted, but the relationship of the global average velocity dispersion of the MCs to their radius follows Larson’s relationship, reflecting virial equilibrium or free-fall collapse. The injection of energy by shocks is visible in the VSFs, but decays within a crossing time.  more » « less
Award ID(s):
1815461 1109395
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Astronomy & Astrophysics
Page Range / eLocation ID:
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract The interstellar medium (ISM) is turbulent over vast scales and in various phases. In this paper, we study turbulence with different tracers in four nearby star-forming regions: Orion, Ophiuchus, Perseus, and Taurus. We combine the APOGEE-2 and Gaia surveys to obtain the full six-dimensional measurements of positions and velocities of young stars in these regions. The velocity structure functions (VSFs) of the stars show a universal scaling of turbulence. We also obtain H α gas kinematics in these four regions from the Wisconsin H-Alpha Mapper. The VSFs of the H α are more diverse compared to those of stars. In regions with recent supernova activities, they show characteristics of local energy injections and higher amplitudes compared to the VSFs of stars and of CO from the literature. Such difference in amplitude of the VSFs can be explained by the different energy and momentum transport from supernovae into different phases of the ISM, thus resulting in higher levels of turbulence in the warm ionized phase traced by H α . In regions without recent supernova activities, the VSFs of young stars, H α , and CO are generally consistent, indicating well-coupled turbulence between different phases. Within individual regions, the brighter parts of the H α gas tend to have a higher level of turbulence than the low-emission parts. Our findings support a complex picture of the Milky Way ISM, where turbulence can be driven at different scales and inject energy unevenly into different phases. 
    more » « less
  2. null (Ed.)
    ABSTRACT Pressure balance plays a central role in models of the interstellar medium (ISM), but whether and how pressure balance is realized in a realistic multiphase ISM is not yet well understood. We address this question by using a set of FIRE-2 cosmological zoom-in simulations of Milky Way-mass disc galaxies, in which a multiphase ISM is self-consistently shaped by gravity, cooling, and stellar feedback. We analyse how gravity determines the vertical pressure profile as well as how the total ISM pressure is partitioned between different phases and components (thermal, dispersion/turbulence, and bulk flows). We show that, on average and consistent with previous more idealized simulations, the total ISM pressure balances the weight of the overlying gas. Deviations from vertical pressure balance increase with increasing galactocentric radius and with decreasing averaging scale. The different phases are in rough total pressure equilibrium with one another, but with large deviations from thermal pressure equilibrium owing to kinetic support in the cold and warm phases, which dominate the total pressure near the mid-plane. Bulk flows (e.g. inflows and fountains) are important at a few disc scale heights, while thermal pressure from hot gas dominates at larger heights. Overall, the total mid-plane pressure is well-predicted by the weight of the disc gas and we show that it also scales linearly with the star formation rate surface density (ΣSFR). These results support the notion that the Kennicutt–Schmidt relation arises because ΣSFR and the gas surface density (Σg) are connected via the ISM mid-plane pressure. 
    more » « less
  3. The intracluster medium (ICM) in the centers of galaxy clusters is heavily influenced by the “feedback” from supermassive black holes (SMBHs). Feedback can drive turbulence in the ICM and turbulent dissipation can potentially be an important source of heating. Due to the limited spatial and spectral resolutions of X-ray telescopes, direct observations of turbulence in the hot ICM have been challenging. Recently, we developed a new method to measure turbulence in the ICM using multiphase filaments as tracers. These filaments are ubiquitous in cluster centers and can be observed at very high resolution using optical and radio telescopes. We study the kinematics of the filaments by measuring their velocity structure functions (VSFs) over a wide range of scales in the centers of ∼ 10 galaxy clusters. We find features of the VSFs that correlate with the SMBHs activities, suggesting that SMBHs are the main driver of gas motions in the centers of galaxy clusters. In all systems, the VSF is steeper than the classical Kolmogorov expectation and the slopes vary from system to system. One theoretical explanation is that the VSFs we have measured so far mostly reflect the motion of the driver (jets and bubbles) rather than the cascade of turbulence. We show that in Abell 1795, the VSF of the outer filaments far from the SMBH flattens on small scales to a Kolmogorov slope, suggesting that the cascade is only detectable farther out with the current telescope resolution. The level of turbulent heating computed at small scales is typically an order of magnitude lower than that estimated at the driving scale. Even though SMBH feedback heavily influences the kinematics of the ICM in cluster centers, the level of turbulence it drives is rather low, and turbulent heating can only offset ≲ 10% of the cooling loss, consistent with the findings of numerical simulations. 
    more » « less

    We present the first empirical constraints on the turbulent velocity field of the diffuse circumgalactic medium around four luminous quasi-stellar objects (QSOs) at z ≈ 0.5–1.1. Spatially extended nebulae of ≈50–100 physical kpc in diameter centred on the QSOs are revealed in [O ii] $\lambda \lambda \, 3727,3729$ and/or [O iii] $\lambda \, 5008$ emission lines in integral field spectroscopic observations obtained using Multi-Unit Spectroscopic Explorer on the Very Large Telescope. We measure the second- and third-order velocity structure functions (VSFs) over a range of scales, from ≲5 kpc to ≈20–50 kpc, to quantify the turbulent energy transfer between different scales in these nebulae. While no constraints on the energy injection and dissipation scales can be obtained from the current data, we show that robust constraints on the power-law slope of the VSFs can be determined after accounting for the effects of atmospheric seeing, spatial smoothing, and large-scale bulk flows. Out of the four QSO nebulae studied, one exhibits VSFs in spectacular agreement with the Kolmogorov law, expected for isotropic, homogeneous, and incompressible turbulent flows. The other three fields exhibit a shallower decline in the VSFs from large to small scales. However, with a limited dynamic range in the spatial scales in seeing-limited data, no constraints can be obtained for the VSF slopes of these three nebulae. For the QSO nebula consistent with the Kolmogorov law, we determine a turbulence energy cascade rate of ≈0.2 cm2 s−3. We discuss the implication of the observed VSFs in the context of QSO feeding and feedback in the circumgalactic medium.

    more » « less

    Interaction of three-dimensional magnetic fields, turbulence, and self-gravity in the molecular cloud is crucial in understanding star formation but has not been addressed so far. In this work, we target the low-mass star-forming region L1688 and use the spectral emissions of 12CO, 13CO, C18O, and H i, as well as polarized dust emissions. To obtain the 3D direction of the magnetic field, we employ the novel polarization fraction analysis. In combining with the plane-of-the-sky (POS) magnetic field strength derived from the Davis–Chandrasekhar–Fermi (DCF) method and the new differential measure analysis (DMA) technique, we present the first measurement of L1688’s three-dimensional magnetic field, including its orientation and strength. We find that L1688’s magnetic field has two statistically different inclination angles. The low-intensity tail has an inclination angle ≈55° on average, while that of the central dense clump is ≈30°. We find the global mean value of total magnetic field strength is Btot ≈ $135 \,\mathrm{\mu }{\rm G}$ from DCF and Btot ≈ $75 \,\mathrm{\mu }{\rm G}$ from DMA. We use the velocity gradient technique (VGT) to separate the magnetic fields’ POS orientation associated with L1688 and its foreground/background. The magnetic fields’ orientations are statistically coherent. The probability density function of H2 column density and VGT reveal that L1688 is potentially undergoing gravitational contraction at large scale ≈1.0 pc and gravitational collapse at small scale ≈0.2 pc. The gravitational contraction mainly along the magnetic field resulting in an approximate power-law relation $B_{\rm tot}\propto n_{\rm H}^{1/2}$ when volume density nH is less than approximately 6.0 × 103 cm−3.

    more » « less