No Arabic abstract
We revisit an alternate explanation for the turbulent nature of molecular clouds - namely, that velocity dispersions matching classical predictions of driven turbulence can be generated by the passage of clumpy material through a shock. While previous work suggested this mechanism can reproduce the observed Larson relation between velocity dispersion and size scale ($sigma propto L^{Gamma}$ with $Gamma approx 0.5$), the effects of self-gravity and magnetic fields were not considered. We run a series of smoothed particle magnetohydrodynamics experiments, passing clumpy gas through a shock in the presence of a combination of self-gravity and magnetic fields. We find powerlaw relations between $sigma$ and $L$ throughout, with indices ranging from $Gamma=0.3-1.2$. These results are relatively insensitive to the strength and geometry of magnetic fields, provided that the shock is relatively strong. $Gamma$ is strongly sensitive to the angle between the gas bulk velocity and the shock front, and the shock strength (compared to the gravitational boundness of the pre-shock gas). If the origin of the $sigma-L$ relation is in clumpy shocks, deviations from the standard Larson relation constrain the strength and behaviour of shocks in spiral galaxies.
We expand our study on the dispersion of polarization angles in molecular clouds. We show how the effect of signal integration through the thickness of the cloud as well as across the area subtended by the telescope beam inherent to dust continuum measurements can be incorporated in our analysis to correctly account for its effect on the measured angular dispersion and inferred turbulent to large-scale magnetic field strength ratio. We further show how to evaluate the turbulent magnetic field correlation scale from polarization data of sufficient spatial resolution and high enough spatial sampling rate. We apply our results to the molecular cloud OMC-1, where we find a turbulent correlation length of approximately 16 mpc, a turbulent to large-scale magnetic field strength ratio of approximately 0.5, and a plane-of-the-sky large-scale magnetic field strength of approximately 0.76 mG.
Supernovae from core-collapse of massive stars drive shocks into the molecular clouds from which the stars formed. Such shocks affect future star formation from the molecular clouds, and the fast-moving, dense gas with compressed magnetic fields is associated with enhanced cosmic rays. This paper presents new theoretical modeling, using the Paris-Durham shock model, and new observations, using the Stratospheric Observatory for Infrared Astronomy (SOFIA), of the H$_2$ S(5) pure rotational line from molecular shocks in the supernova remnant IC443. We generate MHD models for non-steady-state shocks driven by the pressure of the IC443 blast wave into gas of densities $10^3$ to $10^5$ cm$^{-3}$. We present the first detailed derivation of the shape of the velocity profile for emission from H$_2$ lines behind such shocks, taking into account the shock age, preshock density, and magnetic field. For preshock densities $10^3$-$10^5$ cm$^{-3}$, the the predicted shifts of line centers, and the line widths, of the H$_2$ lines range from 20-2, and 30-4 km/s, respectively. The a priori models are compared to the observed line profiles, showing that clumps C and G can be explained by shocks into gas with density 10$^3$ to $2times 10^4$ cm$^{-3}$ and strong magnetic fields. For clump B2 (a fainter region near clump B), the H$_2$ spectrum requires a J-type shock into moderate density (~100 cm$^{-3}$) with the gas accelerated to 100 km/s from its pre-shock location. Clump B1 requires both a magnetic-dominated C-type shock (like for clumps C and G) and a J-type shock (like for clump B1) to explain the highest observed velocities. The J-type shocks that produce high-velocity molecules may be locations where the magnetic field is nearly parallel to the shock velocity, which makes it impossible for a C-type shock (with ions and neutrals separated) to form.
Magnetic fields are dynamically important in the diffuse interstellar medium. Understanding how gravitationally bound, star-forming clouds form requires modeling of the fields in a self-consistent, supernova-driven, turbulent, magnetized, stratified disk. We employ the FLASH magnetohydrodynamics code to follow the formation and early evolution of clouds with final masses of 3-8 $times 10^3 M_{odot}$ within such a simulation. We use the codes adaptive mesh refinement capabilities to concentrate numerical resolution in zoom-in regions covering single clouds, allowing us to investigate the detailed dynamics and field structure of individual self-gravitating clouds in a consistent background medium. Our goal is to test the hypothesis that dense clouds are dynamically evolving objects far from magnetohydrostatic equilibrium. We find that the cloud envelopes are magnetically supported with field lines parallel to density gradients and flow velocity, as indicated by the histogram of relative orientations and other statistical measures. In contrast, the dense cores of the clouds are gravitationally dominated, with gravitational energy exceeding internal, kinetic, or magnetic energy and accelerations due to gravity exceeding those due to magnetic or thermal pressure gradients. In these regions field directions vary strongly, with a slight preference towards being perpendicular to density gradients, as shown by three-dimensional histograms of relative orientation.
We expand on the dispersion analysis of polarimetry maps toward applications to interferometry data. We show how the filtering of low-spatial frequencies can be accounted for within the idealized Gaussian turbulence model, initially introduced for single-dish data analysis, to recover reliable estimates for correlation lengths of magnetized turbulence, as well as magnetic field strengths (plane-of-the-sky component) using the Davis-Chandrasekhar-Fermi method. We apply our updated technique to TADPOL/CARMA data obtained on W3(OH), W3 Main, and DR21(OH). For W3(OH) our analysis yields a turbulence correlation length $deltasimeq19$ mpc, a ratio of turbulent-to-total magnetic energy $leftlangle B_{mathrm{t}}^{2}rightrangle /leftlangle B^{2}rightrangle simeq0.58$, and a magnetic field strength $B_{0}sim1.1:mathrm{mG}$; for W3 Main $deltasimeq22$ mpc, $leftlangle B_{mathrm{t}}^{2}rightrangle /leftlangle B^{2}rightrangle simeq0.74$, and $B_{0}sim0.7:mathrm{mG}$; while for DR21(OH) $deltasimeq12$ mpc, $leftlangle B_{mathrm{t}}^{2}rightrangle /leftlangle B^{2}rightrangle simeq0.70$, and $B_{0}sim1.2:mathrm{mG}$.
The Zeeman effect and dust grain alignment are two major methods for probing magnetic fields (B-fields) in molecular clouds, largely motivated by the study of star formation, as the B-field may regulate gravitational contraction and channel turbulence velocity. This review summarizes our observations of B-fields over the past decade, along with our interpretation. Galactic B-fields anchor molecular clouds down to cloud cores with scales around 0.1 pc and densities of 10^4-5 H_2/cc. Within the cores, turbulence can be slightly super-Alfvenic, while the bulk volumes of pa-rental clouds are sub-Alfvenic. The consequences of these largely ordered cloud B-fields on fragmentation and star formation are observed. The above paradigm is very different from the generally accepted theory during the first decade of the century, when cloud turbulence was assumed to be highly super-Alfvenic. Thus, turbulence anisotropy and turbulence-induced ambipolar diffusion are also revisited.