No Arabic abstract
We present a simulation of isothermal supersonic (rms Mach number $mathcal{M}_{rm rms} sim 3$) turbulent gas with inertial particles (dust) and self-gravity in statistical steady-state, which we compare with a corresponding simulation without self-gravity. The former is in steady state, but close to gravitationally unstable, since we match the scale of the simulation box with Jeans wavelength, which provides the strongest influence of gravity on the dynamics of gas and dust without causing irreversible gravitational collapses. We find that self-gravity does not cause any significant increase in clustering of particles, regardless of particle size, but heavy particles show elevated mean velocities in the presence of self-gravity. The speed distributions are significantly shifted to higher values compared to simulations without self-gravity, but maintains the same shape.
We investigate the clustering and dynamics of nano-sized particles (nano-dust) in high-resolution ($1024^3$) simulations of compressible isothermal hydrodynamic turbulence. It is well-established that large grains will decouple from a turbulent gas flow, while small grains will tend to trace the motion of the gas. We demonstrate that nano-sized grains may cluster in a turbulent flow (fractal small-scale clustering), which increases the local grain density by at least a factor of a few. In combination with the fact that nano-dust grains may be abundant in general, and the increased interaction rate due to turbulent motions, aggregation involving nano dust may have a rather high probability. Small-scale clustering will also affect extinction properties. As an example we present an extinction model based on silicates, graphite and metallic iron, assuming strong clustering of grain sizes in the nanometre range, could explain the extreme and rapidly varying ultraviolet extinction in the host of GRB 140506A.
We conduct numerical experiments to determine the density probability distribution function (PDF) produced in supersonic, isothermal, self-gravitating turbulence of the sort that is ubiquitous in star-forming molecular clouds. Our experiments cover a wide range of turbulent Mach number and virial parameter, allowing us for the first time to determine how the PDF responds as these parameters vary, and we introduce a new diagnostic, the dimensionless star formation efficiency versus density ($epsilon_{rm ff}(s)$) curve, which provides a sensitive diagnostic of the PDF shape and dynamics. We show that the PDF follows a universal functional form consisting of a log-normal at low density with two distinct power law tails at higher density; the first of these represents the onset of self-gravitation, and the second reflects the onset of rotational support. Once the star formation efficiency reaches a few percent, the PDF becomes statistically steady, with no evidence for secular time-evolution at star formation efficiencies from about five to 20 percent. We show that both the Mach number and the virial parameter influence the characteristic densities at which the log-normal gives way to the first power-law, and the first to the second, and we extend (for the former) and develop (for the latter) simple theoretical models for the relationship between these density thresholds and the global properties of the turbulent medium.
The properties of supersonic isothermal turbulence influence a variety of astrophysical phenomena, including the structure and evolution of star forming clouds. This work presents a simple model for the structure of dense regions in turbulence in which the density distribution behind isothermal shocks originates from rough hydrostatic balance between the pressure gradient behind the shock and its deceleration from ram pressure applied by the background fluid. Using simulations of supersonic isothermal turbulence and idealized waves moving through a background medium, we show that the structural properties of dense, shocked regions broadly agree with our analytical model. Our work provides a new conceptual picture for describing the dense regions, which complements theoretical efforts to understand the bulk statistical properties of turbulence and attempts to model the more complex features of star forming clouds like magnetic fields, self-gravity, or radiative properties.
The interstellar turbulence is magnetized and thus anisotropic. The anisotropy of turbulent magnetic fields and velocities is imprinted in the related observables, rotation measures (RMs), and velocity centroids (VCs). This anisotropy provides valuable information on both the direction and strength of the magnetic field. However, its measurement is difficult especially in highly supersonic turbulence in cold interstellar phases due to the distortions by isotropic density fluctuations. By using 3D simulations of supersonic and sub-Alfvenic magnetohydrodynamic(MHD) turbulence, we find that the problem can be alleviated when we selectively sample the volume-filling low-density regions in supersonic MHD turbulence. Our results show that in these low-density regions, the anisotropy of RM and VC fluctuations depends on the Alfvenic Mach number as $rm M_A^{-4/3}$. This anisotropy-$rm M_A$ relation is theoretically expected for sub-Alfv enic MHD turbulence and confirmed by our synthetic observations of $^{12}$CO emission. It provides a new method for measuring the plane-of-the-sky magnetic fields in cold interstellar phases.
The rich structure that we observe in molecular clouds is due to the interplay between strong magnetic fields and supersonic (turbulent) velocity fluctuations. The velocity fluctuations interact with the magnetic field, causing it too to fluctuate. Using numerical simulations, we explore the nature of such magnetic field fluctuations, $vec{delta B}$, over a wide range of turbulent Mach numbers, $mathcal{M} = 2 - 20$ (i.e., from weak to strong compressibility), and Alfven Mach numbers, $mathcal{M}_{text{A}0} = 0.1 - 100$ (i.e., from strong to weak magnetic mean fields, $B_0$). We derive a compressible quasi-static fluctuation model from the magnetohydrodynamical (MHD) equations and show that velocity gradients parallel to the mean magnetic field give rise to compressible modes in sub-Alfvenic flows, which prevents the flow from becoming two-dimensional, as is the case in incompressible MHD turbulence. We then generalise an analytical model for the magnitude of the magnetic fluctuations to include $mathcal{M}$, and find $|vec{delta B}| = delta B = c_ssqrt{pirho_0}mathcal{M}mathcal{M}_{text{A}0}$, where $c_s$ is the sound speed and $rho_0$ is the mean density of gas. This new relation fits well in the strong $B$-field regime. We go on to study the anisotropy between the perpendicular ($ B_{perp}$) and parallel ($ B_{parallel}$) fluctuations and the mean-normalised fluctuations, which we find follow universal scaling relations, invariant of $mathcal{M}$. We provide a detailed analysis of the morphology for the $delta B_{perp}$ and $delta B_{parallel}$ probability density functions and find that eddies aligned with $B_0$ cause parallel fluctuations that reduce $B_{parallel}$ in the most anisotropic simulations. We discuss broadly the implications of our fluctuation models for magnetised gases in the interstellar medium.