skip to main content


Title: A new method for spatially resolving the turbulence-driving mixture in the ISM with application to the Small Magellanic Cloud
ABSTRACT

Turbulence plays a crucial role in shaping the structure of the interstellar medium. The ratio of the three-dimensional density contrast ($\sigma _{\rho /\rho _0}$) to the turbulent sonic Mach number ($\mathcal {M}$) of an isothermal, compressible gas describes the ratio of solenoidal to compressive modes in the turbulent acceleration field of the gas, and is parameterized by the turbulence driving parameter: $b=\sigma _{\rho /\rho _0}/\mathcal {M}$. The turbulence driving parameter ranges from b = 1/3 (purely solenoidal) to b = 1 (purely compressive), with b = 0.38 characterizing the natural mixture (1/3 compressive, 2/3 solenoidal) of the two driving modes. Here, we present a new method for recovering $\sigma _{\rho /\rho _0}$, $\mathcal {M}$, and b, from observations on galactic scales, using a roving kernel to produce maps of these quantities from column density and centroid velocity maps. We apply our method to high-resolution ${\rm H}\,\rm{\small I}$ emission observations of the Small Magellanic Cloud (SMC) from the GASKAP-HI survey. We find that the turbulence driving parameter varies between b ∼ 0.3 and 1.0 within the main body of the SMC, but the median value converges to b ∼ 0.51, suggesting that the turbulence is overall driven more compressively (b > 0.38). We observe no correlation between the b parameter and ${\rm H}\,\rm{\small I}$ or H α intensity, indicating that compressive driving of ${\rm H}\,\rm{\small I}$ turbulence cannot be determined solely by observing ${\rm H}\,\rm{\small I}$ or H α emission density, and that velocity information must also be considered. Further investigation is required to link our findings to potential driving mechanisms such as star-formation feedback, gravitational collapse, or cloud–cloud collisions.

 
more » « less
NSF-PAR ID:
10465615
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
526
Issue:
1
ISSN:
0035-8711
Page Range / eLocation ID:
p. 982-999
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    ABSTRACT We report three-dimensional hydrodynamical simulations of shocks (${\cal M_{\rm shock}}\ge 4$) interacting with fractal multicloud layers. The evolution of shock–multicloud systems consists of four stages: a shock-splitting phase in which reflected and refracted shocks are generated, a compression phase in which the forward shock compresses cloud material, an expansion phase triggered by internal heating and shock re-acceleration, and a mixing phase in which shear instabilities generate turbulence. We compare multicloud layers with narrow ($\sigma _{\rho }=1.9\bar{\rho }$) and wide ($\sigma _{\rho }=5.9\bar{\rho }$) lognormal density distributions characteristic of Mach ≈ 5 supersonic turbulence driven by solenoidal and compressive modes. Our simulations show that outflowing cloud material contains imprints of the density structure of their native environments. The dynamics and disruption of multicloud systems depend on the porosity and the number of cloudlets in the layers. ‘Solenoidal’ layers mix less, generate less turbulence, accelerate faster, and form a more coherent mixed-gas shell than the more porous ‘compressive’ layers. Similarly, multicloud systems with more cloudlets quench mixing via a shielding effect and enhance momentum transfer. Mass loading of diffuse mixed gas is efficient in all models, but direct dense gas entrainment is highly inefficient. Dense gas only survives in compressive clouds, but has low speeds. If normalized with respect to the shock-passage time, the evolution shows invariance for shock Mach numbers ≥10 and different cloud-generating seeds, and slightly weaker scaling for lower Mach numbers and thinner cloud layers. Multicloud systems also have better convergence properties than single-cloud systems, with a resolution of eight cells per cloud radius being sufficient to capture their overall dynamics. 
    more » « less
  2. Abstract

    While it is well known that cosmic rays (CRs) can gain energy from turbulence via second-order Fermi acceleration, how this energy transfer affects the turbulent cascade remains largely unexplored. Here, we show that damping and steepening of the compressive turbulent power spectrum are expected once the damping timetdampρv2/ĖCRECR1becomes comparable to the turbulent cascade time. Magnetohydrodynamic simulations of stirred compressive turbulence in a gas-CR fluid with diffusive CR transport show clear imprints of CR-induced damping, saturating atĖCRϵ˜, whereϵ˜is the turbulent energy input rate. In that case, almost all of the energy in large-scale motions is absorbed by CRs and does not cascade down to grid scale. Through a Hodge–Helmholtz decomposition, we confirm that purely compressive forcing can generate significant solenoidal motions, and we find preferential CR damping of the compressive component in simulations with diffusion and streaming, rendering small-scale turbulence largely solenoidal, with implications for thermal instability and proposed resonant scattering ofE≳ 300 GeV CRs by fast modes. When CR transport is streaming dominated, CRs also damp large-scale motions, with kinetic energy reduced by up to 1 order of magnitude in realisticECREgscenarios, but turbulence (with a reduced amplitude) still cascades down to small scales with the same power spectrum. Such large-scale damping implies that turbulent velocities obtained from the observed velocity dispersion may significantly underestimate turbulent forcing rates, i.e.,ϵ˜ρv3/L.

     
    more » « less
  3. ABSTRACT

    We present new observations of the Mystic Mountains cloud complex in the Carina Nebula using the ALMA Atacama Compact Array (ACA) to quantify the impact of strong UV radiation on the structure and kinematics of the gas. Our Band 6 observations target CO, 13CO, and C18O; we also detect DCN J=3–2 and 13CS J=5–4. A dendrogram analysis reveals that the Mystic Mountains are a coherent structure, with continuous emission over −10.5 km s−1 < v < −2 km s−1. We perform multiple analyses to isolate non-thermal motions in the Mystic Mountains including computing the turbulent driving parameter, b, which indicates whether compressive or solenoidal modes dominate. Each analysis yields values similar to other pillars in Carina that have been observed in a similar way but are subject to an order of magnitude less intense ionizing radiation. We find no clear correlation between the velocity or turbulent structure of the gas and the incident radiation, in contrast to other studies targeting different regions of Carina. This may reflect differences in the initial densities of regions that go on to collapse into pillars and those that still look like clouds or walls in the present day. Pre-existing over-densities that enable pillar formation may also explain why star formation in the pillars appears more evolved (from the presence of jets) than in other heavily irradiated but non-pillar-like regions. High resolution observations of regions subject to an array of incident radiation are required to test this hypothesis.

     
    more » « less
  4. null (Ed.)
    Context. The excitation of the filamentary gas structures surrounding giant elliptical galaxies at the center of cool-core clusters, also known as brightest cluster galaxies (BCGs), is key to our understanding of active galactic nucleus (AGN) feedback, and of the impact of environmental and local effects on star formation. Aims. We investigate the contribution of thermal radiation from the cooling flow surrounding BCGs to the excitation of the filaments. We explore the effects of small levels of extra heating (turbulence), and of metallicity, on the optical and infrared lines. Methods. Using the C LOUDY code, we modeled the photoionization and photodissociation of a slab of gas of optical depth A V  ≤ 30 mag at constant pressure in order to calculate self-consistently all of the gas phases, from ionized gas to molecular gas. The ionizing source is the extreme ultraviolet (EUV) and soft X-ray radiation emitted by the cooling gas. We tested these models comparing their predictions to the rich multi-wavelength observations from optical to submillimeter, now achieved in cool core clusters. Results. Such models of self-irradiated clouds, when reaching sufficiently large A V , lead to a cloud structure with ionized, atomic, and molecular gas phases. These models reproduce most of the multi-wavelength spectra observed in the nebulae surrounding the BCGs, not only the low-ionization nuclear emission region like optical diagnostics, [O  III ] λ 5007 Å/H β , [N  II ] λ 6583 Å/H α , and ([S  II ] λ 6716 Å+[S  II ] λ 6731 Å)/H α , but also the infrared emission lines from the atomic gas. [O  I ] λ 6300 Å/H α , instead, is overestimated across the full parameter space, except for very low A V . The modeled ro-vibrational H 2 lines also match observations, which indicates that near- and mid-infrared H 2 lines are mostly excited by collisions between H 2 molecules and secondary electrons produced naturally inside the cloud by the interaction between the X-rays and the cold gas in the filament. However, there is still some tension between ionized and molecular line tracers (i.e., CO), which requires optimization of the cloud structure and the density of the molecular zone. The limited range of parameters over which predictions match observations allows us to constrain, in spite of degeneracies in the parameter space, the intensity of X-ray radiation bathing filaments, as well as some of their physical properties like A V or the level of turbulent heating rate. Conclusions. The reprocessing of the EUV and X-ray radiation from the plasma cooling is an important powering source of line emission from filaments surrounding BCGs. C LOUDY self-irradiated X-ray excitation models coupled with a small level of turbulent heating manage to simultaneously reproduce a large number of optical-to-infrared line ratios when all the gas phases (from ionized to molecular) are modeled self-consistently. Releasing some of the simplifications of our model, like the constant pressure, or adding the radiation fields from the AGN and stars, as well as a combination of matter- and radiation-bounded cloud distribution, should improve the predictions of line emission from the different gas phases. 
    more » « less
  5. Abstract

    The Large Magellanic Cloud (LMC) has an extensive Hαemission halo that traces an extended, warm ionized component of its interstellar medium. Using the Wisconsin HαMapper telescope, we present the first kinematic Hαsurvey of an extensive region around the LMC, from (,b) = (264.°5, − 45.°5) to (295.°5, − 19.°5), covering +150 ≤vLSR≤ + 390 km s−1. We find that ionized hydrogen exists throughout the galaxy and extends several degrees beyond detected neutral hydrogen emission(logNHI/cm218.3)as traced by 21 cm in current surveys. Using the column density structure of the neutral gas and stellar line-of-sight depths as a guide, we estimate the upper limit mass of the ionized component of the LMC to be roughlyMionized≈ (0.6–1.8) × 109M, which is comparable to the total neutral atomic gas mass in the same region (Mneutral≈ 0.76–0.85 × 109M). Considering only the atomic phases, we findMionized/Mionized+neutral, to be 46%–68% throughout the LMC and its extended halo. Additionally, we find an ionized gas cloud that extends off of the LMC at (,b) ≈ (285°, − 28°) into a region previously identified as the Leading Arm complex. This gas is moving at a similar line-of-sight velocity as the LMC and hasMionized/Mionized+neutral= 13%–51%. This study, combined with previous studies of the SMC and extended structures of the Magellanic Clouds, continues to suggest that warm, ionized gas is as massive and dynamically important as the neutral gas in the Magellanic System.

     
    more » « less