Understanding intermittent dissipation in the interstellar medium is integral to understanding the chemical abundances and line luminosities observed therein. In the cold neutral medium, continuous heating processes due to cosmic rays and photoelectric emission from dust grains determine a unique equilibrium temperature varying with density. Excursions from this are due to intermittent 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 sight-lines. Intermittent high temperatures should also have an impact on H$_2$ line luminosities. We carry out simulations of MHD turbulence in molecular clouds including heating and cooling, and post-process them to study H$_2$ 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 H$_2$ cooling function for $n_{\rm H} \leq 10^5\,{\rm cm}^{-3}$, $T\leq 5000\,{\rm K}$, and variable H$_2$ fraction. Our models produce H$_2$ emission lines in accord with many observations, although extra excitation mechanisms are required in some clouds. For realistic r.m.s. magnetic field strengths ($\approx 10$ $\mu$G) and velocity dispersions, we reproduce observed CH$^+$ abundances. Comparison of predicted dust polarization with observations by ${\it Planck}$ suggests that the mean field $\gtrsim 5\,\mu$G, so that the turbulence is sub-Alfv\'enic. We recommend future work treating ions and neutrals as separate fluids to more accurately capture the effects of ambipolar diffusion on CH$^+$ abundance.

We study the non-linear evolution of the acoustic 'Resonant Drag Instability' (RDI) using numerical simulations. The acoustic RDI is excited in a dust-gas mixture when dust grains stream through gas, interacting with sound waves to cause a linear instability. We study this process in a periodic box by accelerating neutral dust with an external driving force. The instability grows as predicted by linear theory, eventually breaking into turbulence and saturating. As in linear theory, the non-linear behavior is characterized by three regimes - high, intermediate, and low wavenumbers - the boundary between which is determined by the dust-gas coupling strength and the dust-to-gas mass ratio. The high and intermediate wavenumber regimes behave similarly to one another, with large dust-to-gas ratio fluctuations while the gas remains largely incompressible. The saturated state is highly anisotropic: dust is concentrated in filaments, jets, or plumes along the direction of acceleration, with turbulent vortex-like structures rapidly forming and dissipating in the perpendicular directions. The low-wavenumber regime exhibits large fluctuations in gas and dust density, but the dust and gas remain more strongly coupled in coherent 'fronts' perpendicular to the acceleration. These behaviors are qualitatively different from those of dust 'passively' driven by external hydrodynamic turbulence, with no back-reaction force from dust onto gas. The virulent nature of these instabilities has interesting implications for dust-driven winds in a variety of astrophysical systems, including around cool-stars, in dusty torii around active-galactic-nuclei, and in and around giant molecular clouds.