Time-dependent simulations of steady C-type shocks

Using a time-dependent multifluid, magnetohydrodynamic code, we calculated the structure of steady perpendicular and oblique C-type shocks in dusty plasmas. We included relevant processes to describe mass transfer between the different fluids, radiative cooling by emission lines and grain charging and studied the effect of single-sized and multiple sized grains on the shock structure. Our models are the first of oblique fast-mode molecular shocks in which such a rigorous treatment of the dust grain dynamics has been combined with a self-consistent calculation of the thermal and ionisation structures including appropriate microphysics. At low densities the grains do not play any significant role in the shock dynamics. At high densities, the ionisation fraction is sufficiently low that dust grains are important charge and current carriers and, thus, determine the shock structure. We find that the magnetic field in the shock front has a significant rotation out of the initial upstream plane. This is most pronounced for single-sized grains and small angles of the shock normal with the magnetic field. Our results are similar to previous studies of steady C-type shocks showing that our method is efficient, rigorous and robust. Unlike the method employed in the previous most detailed treatment of dust in steady oblique fast-mode shocks, ours allows a reliable calculation even when chemical or other conditions deviate from local statistical equilibrium. We are also able to model transient phenomena.

💡 Research Summary

The authors present a comprehensive study of steady C‑type (continuous) shocks in dusty molecular gas using a time‑dependent, multi‑fluid magnetohydrodynamic (MHD) code. The model treats the neutral gas, ions, electrons, and up to several grain fluids as separate components, each governed by its own continuity, momentum, and energy equations. Mass exchange between fluids is driven by cosmic‑ray ionisation, recombination (electron–Mg⁺, dissociative recombination of HCO⁺), and adsorption of ions/electrons onto grain surfaces. Grain charging follows the formalism of Havnes, Hartquist & Pilipp (1987), with the average grain charge determined by the balance of electron‑grain and ion‑grain collision rates.

The electric field is expressed in terms of Pedersen, Hall, and ambipolar resistivities (r‖, r_H, r_AD), allowing the magnetic induction equation to include the full resistance matrix. The resistivities depend strongly on the Hall parameter β_i = α_i|B|/(K_ni ρ_n), which measures the ratio of gyro‑frequency to neutral collision frequency for each charged species. This formulation captures the Hall effect that can rotate the magnetic field out of the upstream plane.

Radiative cooling is implemented with a detailed chemical network based on Pilipp, Hartquist & Havnes (1990) and includes line cooling from O I, CO, H₂ (rotational‑vibrational), and H₂O (rotational). Electron cooling by H₂ excitation is added for temperatures above ~1000 K. The ion temperature is not solved independently; instead it is derived from the neutral temperature and the drift speed using T_i = T_n + (m_n/3k_B)(v_n−v_i)².

Initial conditions are taken directly from the classic Draine, Roberge & Dalgarno (1983) models to enable direct comparison. Two upstream densities are examined: n_H = 10⁴ cm⁻³ (low‑density) and n_H = 10⁶ cm⁻³ (high‑density). The magnetic field strength scales as B ∝ n_H^½, giving B = 10⁻⁴ G and 10⁻³ G respectively, with an Alfvén speed of 2.2 km s⁻¹ in both cases. Shocks propagate at 25 km s⁻¹, corresponding to an Alfvénic Mach number M_A ≈ 11, i.e., a strong fast‑mode shock. The upstream magnetic field lies in the x‑y plane and makes angles of 30°, 45°, 60°, or 90° (perpendicular) with the shock normal.

The simulations evolve the discontinuous initial jump at x = 0 forward in time until a steady profile is reached in the shock frame. A uniform grid of 400 cells spans roughly ±10–6 L_g, where L_g is the grain‑neutral drag length scale, ensuring that both ion‑neutral and grain‑neutral coupling lengths are resolved.

Key results:

1. Low‑density regime (n_H = 10⁴ cm⁻³): The ionisation fraction remains relatively high, so ions and electrons dominate the charge and current transport. Grains are essentially stationary relative to the neutrals. The shock width is about four times larger than in the DRD steady‑state models because the authors compute the ionisation balance self‑consistently rather than imposing a constant ion flux. Consequently, the ion‑neutral drag is weaker and magnetic pressure balances over a larger distance. The maximum neutral temperature in the shock front reaches only ~530 K (vs. ~1200 K in DRD), reflecting the longer cooling time afforded by the broader precursor.

2. High‑density regime (n_H = 10⁶ cm⁻³): The fractional ionisation drops below 10⁻⁷, and grains become the primary charge carriers. In the single‑size grain model (a = 0.4 µm), the Hall conductivity is strongly enhanced, leading to a noticeable rotation of the magnetic field out of the upstream x‑y plane. The rotation is most pronounced for small shock‑field angles (30°) and diminishes as the angle approaches perpendicular. When a distribution of grain sizes is introduced, the Hall effect is partially averaged out, reducing the field rotation.

3. Oblique shocks: For angles of 30°, 45°, and 60°, the transverse magnetic components (B_y, B_z) evolve differently for ions/electrons versus grains, reflecting their distinct drift speeds. The grain drift speed relative to the neutrals is always smaller than the ion drift, because grains experience both Pedersen and Hall forces, whereas ions are dominated by the Pedersen term in the parameter regime studied.

4. Transient and non‑equilibrium chemistry: Because the code integrates the full time‑dependent equations, the authors naturally capture the evolution of species such as O and H₂O, whose abundances can change dramatically over 10⁴–10⁵ yr in the shock precursor. This capability overcomes a major limitation of earlier steady‑state approaches that required the assumption of local chemical equilibrium.

5. Methodological robustness: The authors demonstrate that their scheme converges to the same steady profiles as earlier multi‑fluid studies (Pilipp & Hartquist 1994; Draine et al. 1983) while offering the flexibility to handle arbitrary chemical networks, grain size distributions, and transient phenomena. The explicit treatment of the magnetic field imposes a Courant‑type timestep restriction, but the second‑order Godunov solver for the neutral fluid and the upwind scheme for the charged fluids ensure numerical stability.

In summary, this work provides the first fully self‑consistent, time‑dependent simulations of steady fast‑mode C‑type shocks that include a rigorous treatment of dust dynamics, grain charging, and detailed cooling. It confirms that dust grains dominate the electrodynamics at high densities, that the Hall effect can rotate the magnetic field in oblique shocks, and that non‑equilibrium chemistry can be handled without sacrificing the ability to reach a steady solution. These results have direct implications for interpreting molecular line emission from shocked regions in star‑forming clouds, protostellar outflows, and other astrophysical environments where dust‑laden, weakly ionised plasmas are prevalent.