skip to main content


Title: Turbulent dissipation, CH+ abundance, H2 line luminosities, and polarization in the cold neutral medium
ABSTRACT In the cold neutral medium, high out-of-equilibrium temperatures are created by intermittent dissipation processes, including shocks, viscous heating, and ambipolar diffusion. The high-temperature excursions are thought to explain the enhanced abundance of CH+ observed along diffuse molecular sightlines. Intermittent high temperatures should also have an impact on H2 line luminosities. We carry out simulations of magnetohydrodynamic (MHD) turbulence in molecular clouds including heating and cooling, and post-process them to study H2 line emission and hot-gas chemistry, particularly the formation of CH+. We explore multiple magnetic field strengths and equations of state. We use a new H2 cooling function for $n_{\text{H}}\le 10^5\, {\text{cm}}^{-3}$, $T\le 5000\, {\text{K}}$, and variable H2 fraction. We make two important simplifying assumptions: (i) the H2/H fraction is fixed everywhere and (ii) we exclude from our analysis regions where the ion–neutral drift velocity is calculated to be greater than 5 km s−1. Our models produce H2 emission lines in accord with many observations, although extra excitation mechanisms are required in some clouds. For realistic root-mean-square (rms) magnetic field strengths (≈10 μG) and velocity dispersions, we reproduce observed CH+ abundances. These findings contrast with those of Valdivia et al. (2017) Comparison of predicted dust polarization with observations by Planck suggests that the mean field is ≳5 µG, so that the turbulence is sub-Alfvénic. We recommend future work treating ions and neutrals as separate fluids to more accurately capture the effects of ambipolar diffusion on CH+ abundance.  more » « less
Award ID(s):
1908123 1408723
NSF-PAR ID:
10286268
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
500
Issue:
3
ISSN:
0035-8711
Page Range / eLocation ID:
3290 to 3308
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT

    Magnetic fields are a defining yet enigmatic aspect of the interstellar medium, with their three-dimensional (3D) mapping posing a substantial challenge. In this study, we harness the innovative velocity gradient technique (VGT), underpinned by magnetohydrodynamic turbulence theories, to map the magnetic field structure by applying it to the atomic neutral hydrogen (H i) emission line and the molecular tracer 12CO. We construct the tomography of the magnetic field in the low-mass star-forming region L1688, utilizing two approaches: (1) VGT-H i combined with the Galactic rotational curve, and (2) stellar polarization paired with precise star parallax measurements. Our analysis reveals that the magnetic field orientations deduced from stellar polarization undergo a distinct directional change in the vicinity of L1688, providing evidence that the misalignment between VGT-H i and stellar polarization stems from the influence of the molecular cloud’s magnetic field on the polarization of starlight. When comparing VGT-12CO to stellar polarization and Planck polarization data, we observe that VGT-12CO effectively reconciles the misalignment noted with VGT-H i, showing statistical alignment with Planck polarization measurements. This indicates that VGT-12CO could be integrated with VGT-H i, offering vital insights into the magnetic fields of molecular clouds, thereby enhancing the accuracy of our 3D magnetic field reconstructions.

     
    more » « less
  2. Context. The 1°.3 (G1.3) and 1°.6 (G1.6) cloud complexes in the central molecular zone (CMZ) of our Galaxy have been proposed to possibly reside at the intersection region of the X1 and X2 orbits for several reasons. This includes the detection of co-spatial low- and high-velocity clouds, high velocity dispersion, high fractional molecular abundances of shock-tracing molecules, and kinetic temperatures that are higher than for usual CMZ clouds. Aims. By investigating the morphology and deriving physical properties as well as chemical composition, we want to find the origin of the turbulent gas and, in particular, whether evidence of an interaction between clouds can be identified. Methods. We mapped both cloud complexes in molecular lines in the frequency range from 85 to 117 GHz with the IRAM 30 m telescope. The APEX 12m telescope was used to observe higher frequency transitions between 210 and 475 GHz from selected molecules that are emitted from higher energy levels. We performed non-local thermodynamic equilibrium (non-LTE) modelling of the emission of an ensemble of CH 3 CN lines to derive kinetic temperatures and H 2 volume densities. These were used as starting points for non-LTE modelling of other molecules, for which column densities and abundances were determined and compared with values found for other sources in the CMZ. Results. The kinematic structure of G1.3 reveals an ‘emission bridge’ at intermediate velocities (~150 km s −1 ) connecting low-velocity (~100 km s −1 ) and high-velocity (~180 km s −1 ) gas and an overall fluffy shell-like structure. These may represent observational evidence of cloud-cloud interactions. Low- and high-velocity gas components in G1.6 do not show this type of evidence of an interaction, suggesting that they are spatially separated. We selected three positions in each cloud complex for further analysis. Each position reveals several gas components at various peak velocities and of various line widths. We derived kinetic temperatures of 60–100 K and H 2 volume densities of 10 4 –10 5 cm −3 in both complexes. Molecular abundances relative to H 2 suggest a similar chemistry of the two clouds, which is moreover similar to that of other GC clouds and, especially, agrees well with that of G+0.693 and G−0.11. Conclusions. We conclude that G1.3 may indeed exhibit signs of cloud-cloud interactions. In particular, we propose an interaction of gas that is accreted from the near-side dust lane to the CMZ, with gas pre-existing at this location. Low- and high-velocity components in G1.6 are rather coincidentally observed along the same line of sight. They may be associated with either overshot decelerated gas from the far-side dust line or actual CMZ gas and high-velocity gas moving on a dust lane. These scenarios would be in agreement with numerical simulations. 
    more » « less
  3. Abstract

    The mean plane-of-sky magnetic field strength is traditionally obtained from the combination of polarization and spectroscopic data using the Davis–Chandrasekhar–Fermi (DCF) technique. However, we identify the major problem of the DCF technique to be its disregard of the anisotropic character of MHD turbulence. On the basis of the modern MHD turbulence theory we introduce a new way of obtaining magnetic field strength from observations. Unlike the DCF technique, the new technique uses not the dispersion of the polarization angle and line-of-sight velocities, but increments of these quantities given by the structure functions. To address the variety of astrophysical conditions for which our technique can be applied, we consider turbulence in both media with magnetic pressure higher than the gas pressure, corresponding, e.g., to molecular clouds, and media with gas pressure higher than the magnetic pressure, corresponding to the warm neutral medium. We provide general expressions for arbitrary admixtures of Alfvén, slow, and fast modes in these media and consider in detail particular cases relevant to diffuse media and molecular clouds. We successfully test our results using synthetic observations obtained from MHD turbulence simulations. We demonstrate that our differential measure approach, unlike the DCF technique, can be used to measure the distribution of magnetic field strengths, can provide magnetic field measurements with limited data, and is much more stable in the presence of induced large-scale variations of nonturbulent nature. Furthermore, our study uncovers the deficiencies of earlier DCF research.

     
    more » « less
  4. Context. Magnetic flux emergence from the solar interior has been shown to be a key mechanism for unleashing a wide variety of phenomena. However, there are still open questions concerning the rise of the magnetized plasma through the atmosphere, mainly in the chromosphere, where the plasma departs from local thermodynamic equilibrium (LTE) and is partially ionized. Aims. We aim to investigate the impact of the nonequilibrium (NEQ) ionization and recombination and molecule formation of hydrogen, as well as ambipolar diffusion, on the dynamics and thermodynamics of the flux emergence process. Methods. Using the radiation-magnetohydrodynamic Bifrost code, we performed 2.5D numerical experiments of magnetic flux emergence from the convection zone up to the corona. The experiments include the NEQ ionization and recombination of atomic hydrogen, the NEQ formation and dissociation of H 2 molecules, and the ambipolar diffusion term of the generalized Ohm’s law. Results. Our experiments show that the LTE assumption substantially underestimates the ionization fraction in most of the emerged region, leading to an artificial increase in the ambipolar diffusion and, therefore, in the heating and temperatures as compared to those found when taking the NEQ effects on the hydrogen ion population into account. We see that LTE also overestimates the number density of H 2 molecules within the emerged region, thus mistakenly magnifying the exothermic contribution of the H 2 molecule formation to the thermal energy during the flux emergence process. We find that the ambipolar diffusion does not significantly affect the amount of total unsigned emerged magnetic flux, but it is important in the shocks that cross the emerged region, heating the plasma on characteristic times ranging from 0.1 to 100 s. We also briefly discuss the importance of including elements heavier than hydrogen in the equation of state so as not to overestimate the role of ambipolar diffusion in the atmosphere. 
    more » « less
  5. ABSTRACT Molecular hydrogen is the main constituent of dense molecular clouds, but is expected to also be a dominant constituent in many environments where CO can no longer be seen, the so-called ‘CO-dark molecular gas’. Based on comparisons of ultraviolet spectroscopy of H2 and optical line observations (4300 Å), CH is a prime candidate to trace H2. Since the optical line (and the UV lines at 3143, 3890, and 3878 Å) require bright background sources (and the CH N = 2←1 ground state rotation line at 149 µm requires space-based, or stratospheric, observations), the hyperfine structure transition at 3335 MHz is a potentially important tool for probing the CO-dark molecular gas. However, the excitation of this transition is complicated, and has often been found to be inverted, making column density determinations uncertain. To clarify the potential use of the 3.3-GHz line as a proxy for H2, we have observed the CH 3335-MHz line with the Arecibo 305-m radio telescope along 16 lines of sight towards stars with existing measurements of the 4300-Å line. By comparing the CH column densities from optical and UV absorption lines to the CH radio emission line, we can derive the excitation temperature (Tex) of the 3335-MHz transition. We obtain a wide range of excitation temperatures for nine lines of sight, including some with |Tex| < 5 K. The common assumption that Tex for the 3335-MHz line is always much larger than the background temperature (Tbg) is not always warranted and can lead to significant errors in the value of N(CH). 
    more » « less