# Density Statistics of Compressible MHD Turbulence

###### Abstract

Density is the turbulence statistics that is most readily available from observations. Different regimes of turbulence correspond to different density spectra. For instance, the viscosity-damped regime of MHD turbulence relevant, for instance, to partially ionized gas, can be characterized by shallow and very anisotropic spectrum of density. This spectrum can result in substantial variations of the column densities. Addressing MHD turbulence in the regime when viscosity is not important over the inertial range, we demonstrate with numerical simulations that it is possible to reproduce both the observed Kolmogorov spectrum of density fluctuations observed in ionized gas by measuring scintillations and more shallow spectra that are obtained through the emission measurements. We show that in supersonic turbulence the high density peaks dominate shallow isotropic spectrum, while the small-scale underlying turbulence that fills most of the volume has the Kolmogorov spectrum and demonstrates scale-dependent anisotropy. The limitations of the spectrum in studying turbulence induce searches of alternative statistics. We demonstrate that a measure called ”bispectrum” may be a promising tool. Unlike spectrum, the bispectrum preserves the information about wave phases.

^{1}

^{1}footnotetext: Astronomy Department, University of Wisconsin, Madison, WI 53706

^{2}

^{2}footnotetext: Astronomical Observatory, Jagiellonian University, Kraków, Poland

## 1. Density Fluctuations as a Probe of MHD Turbulence

The paradigm of interstellar medium has undergone substantial changes recently. Instead of quiescent medium with hanging and slowly evolving clouds a turbulent picture emerged (see review by Ballesteros-Parredes et al. 2006, McKee & Ostriker 2007 and ref. therein). With magnetic field being dynamically important and dominating the gas pressure in molecular clouds, this calls for studies of compressible magnetohydrodynamic (MHD) turbulence.

Key ideas in describing MHD turbulence can be traced back to Iroshnikov (1963) and Kraichnan (1965) work and the classical work that followed (see Montgometry & Turner 1981, Higdon 1984, Montgomery, Brown & Matthaeus 1987). A more recent model by Goldreich & Sridhar (1995, henceforth GS95) has been successfully tested by numerical simulations (Cho & Vishniac 2000, Maron & Goldreich 2001, Cho, Lazarian & Vishniac 2002, Cho & Lazarian 2002, 2003, henceforth CL03).

Density statistics is the easiest to infer from observations. In comparison, velocity spectra require elaborate techniques to be used (see Lazarian 2006 and ref. therein). How informative can be density fluctuations for understanding of MHD turbulence is the subject of the present review. However, we should mention that for many important astrophysical applications, e.g. for interstellar chemistry, for star formation, for propagation of radiation etc. the density fluctuations themselves are crucially important.

Subsonic compressible MHD is rather well studied topic today. It is suggestive that there may be an analogy between the subsonic MHD turbulence and its incompressible counterpart, namely, GS95 model. Therefore the correspondence between the the two revealed in CL03 is expected.

It could be easily seen, that in the low-beta case density is perturbed mainly due to the slow mode (CL03). Slow modes are sheared by Alfvén turbulence, therefore they exhibit Kolmogorov scaling and GS95 anisotropy for low Mach numbers. However, for high Mach numbers we expect shocks to develop. Density will be perturbed mainly by those shocks.

One can also approach the problem from the point of view of underlying hydrodynamic equations. It is well known that there is a multiplicative symmetry with respect to density in the ideal flow equations for an isothermal fluid (see e.g. Passot & Vazquez-Semadeni, 1998). This assume that if there is some stochastic process disturbing the density it should be a multiplicative process with respect to density, rather than additive, and the distribution for density values should be lognormal, rather that normal. The aforementioned work shows that for 1D numerical simulations of high-Mach hydrodynamics the distribution is approximately lognormal, having power-law tails in case of .

In MHD the above described symmetry is broken by the magnetic tension. However, numerical studies in Beresnyak, Lazarian & Cho (2005, henceforth BLC05) show an approximate correspondence of the PDFs obtained with MHD simulations to the lognormal scaling.

With a high sonic Mach we expect a considerable amount of shocks arise. In a sub-Alfvénic case, however, we expect oblique shocks be disrupted by Alfvénic shearing, and, as most of the shocks are generated randomly by driving, almost all of them will be sheared to smaller shocks. The evolution of the weak shocks will be again governed by the sonic speed, and structures from shearing as in low Mach case should arise.

We also note that shearing will not affect probability density function (PDF) of the density, but have to affect its spectra and structure function (SF) scaling. In other words, we deal with two distinct physical processes, one of which, random multiplication or division of density in presence of shocks, affect PDF, while the other, Alfvénic shearing has to affect anisotropy and scaling of the structure function of the density. The numerical studies that we describe below confirm the theoretical considerations above.

## 2. Viscosity-Damped MHD Turbulence as a Special Case

Turbulence can be viewed as a cascade of energy from a large injection energy scale to dissipation at a smaller scale. The latter is being established by equating the rate of turbulent energy transfer to the rate of energy damping arising, for instance, from viscosity. Naively, one does not expect to see any turbulent phenomena below such a scale.

Such reasoning may not be true in the presence of magnetic field, however. Consider magnetized fluid with viscosity much larger than magnetic diffusivity , which is the case of a high magnetic Prantl number fluid. The partially ionized gas can serve as an example of such a fluid up to the scales of ion-neutral decoupling (see a more rigorous treatment in Lazarian, Vishniac & Cho 2004, henceforth LVC04). Fully ionized plasma is a more controversial example. For instance, It is well known that for plasma the diffusivities are different along and perpendicular to magnetic field lines. Therefore, the plasma the Prandtl number is huge if we use the parallel diffusivity . A treatment of the fully ionized plasma as a high Prandtl number medium is advocated in Schechochihin et al. (2004, henceforth SCTMM).

The turbulent cascade in a fluid with isotropic proceeds up to a scale at which the cascading rate, which for the Kolmogorov turbulence, i.e. , is determined by the eddy turnover rate gets equal to the damping rate . Assuming that the energy is injected at the scale and the injection velocity is , the damping scale is , where is the Reynolds number .The corresponding scale varies from pc for Warm Neutral Medium to for molecular clouds (LVC04). It is evident, however, that magnetic fields at the scale at which hydrodynamic cascade would stop are still sheered by eddies at the larger scales. This should result in creating magnetic structures at scales (LVC04). Note, that magnetic field in this regime is not a passive scalar, but important dynamically. Thus, GS95 turbulence range an inertial range in the partially ionized gas from the injection scale to , but below is the domain of the viscosity-damped turbulence, where plasma/gas compressed within current sheets can contribute to the observed extreme scattering events and the formation of the Small Ionized or Neutral Structures (SINS) (Lazarian 2007). Most compressions arises in the situations when the magnetic fields exhibits perfect reversals.

Cho, Lazarian & Vishniac (2002) reported a new viscosity-damped regime of turbulence using incompressible MHD simulations (see an example of more recent simulations in Fig. 1). In this regime, unlike hydro turbulence, motions, indeed, do not stop at the viscosity damping scale, but magnetic fluctuations protrude to smaller scales. Interestingly enough, these magnetic fluctuations drive small amplitude velocity fluctuations at scales . CL03 confirmed these results with compressible simulations and speculated that these small scale magnetic fluctuations can compress ambient gas to produce ubiquitous tiny structures observed in ISM (see also Fig. 1b). According to the model in LVC04, while the spectrum of averaged over volume magnetic fluctuations scales as , the pressure within intermittent magnetic structures increases with the decrease of the scale as , while the filling factor , the latter being consistent with numerical simulations.

In Fig. 1 we show the results of our high resolution compressible MHD simulations that exhibit strong density fluctuations at the scales below the one at which hydrodynamic turbulence would be damped. According to LVC04 model, the viscosity-damped regime is ubiquitous in turbulent partially ionized gas. Some of its discussed consequences, such as an intermittent resumption of the turbulence the fluid of ions as magnetic fluctuations reach the ion-neutral decoupling scale, are important for radio scintillations.

## 3. Simulations 3D statistics and Simulated Observations

### 3.1. Spectrum

The power spectrum (PS) of density fluctuations is an important property of a compressible flow. In some cases, the spectrum of density can be derived analytically. For nearly turbulent motions in the presence of a strong magnetic field, the spectrum of density scales similarly to the pressure, i.e. if we consider the polytropic equation of state (Biskamp 2003). In weakly magnetized nearly incompressible MHD turbulence, however, velocities convect density fluctuations passively inducing the spectrum (Montgomery et al. 1987). In supersonic flows, these relations are not valid anymore because of shocks accumulating matter into the local and highly dense structures. Due to the high contrast of density, the linear relation is no longer valid, and the spectrum of density cannot be related to pressure so straightforwardly. In addition, the strong asymmetry of density fluctuations suggests the need to analyze the logarithm of density instead of density itself.

In Figure 2 we present the PS of fluctuations of density for models with different . As expected, we note a strong growth of the amplitude of density fluctuations with the sonic Mach number at all scales. This behavior is observed both in sub-Alfvénic as well as in super-Alfvénic turbulence (see Fig. 2). In Table 1 we calculate the spectral index of density and the logarithm of density within the inertial ranges estimated from the PS of velocity. The width of the inertial range is shown by the range of solid lines with slopes and in all spectra plots. It is estimated to be within . In Table 1 we also show the errors of estimation which combine the error of the fitting of the spectral index at each time snapshot and the standard deviation of variance of in time. The slopes of the density spectra do not change significantly with for subsonic experiments and correspond to analytical estimations (about , which is slightly less than , for turbulence with and about , which is slightly more than , for weakly magnetized turbulence with ). Such an agreement confirms the validity of the theoretical approximations. Those, nevertheless, do not cover the entire parameter space. While the fluid motions become supersonic, they strongly influence the density structure, making the small-scale structures more pronounced, which implies flattening of the spectra of density fluctuations (see values for in Table 1) (see also BLC05). Although, the spectral indices are clearly different for sub- and super-Alfvénic turbulence, their errors are relatively large (up to ). These errors have been calculated by summing the maximum value of the uncertainty of fits for individual spectra at each time snapshot and the standard deviation of time variation of spectral indices. The uncertainty of fit contributes the most to the total error bar (about 60%-70%).

0.23 | -2.3 | -2.3 | 0.26 | -1.7 | -1.7 |
---|---|---|---|---|---|

0.33 | -2.2 | -2.2 | 0.36 | -1.7 | -1.7 |

0.68 | -2.0 | -2.1 | 0.74 | -1.6 | -1.6 |

2.20 | -1.3 | -2.0 | 2.34 | -1.2 | -1.6 |

7.0 | -0.5 | -1.7 | 7.1 | -0.6 | -1.5 |

The values have been estimated within the inertial range for models with (left) and (right). Errors of spectral indices combine the errors of estimation at each time snapshot and the standard deviation of variance in time. Errors for sonic Mach numbers are the standard deviation of their variance in time calculated over the period starting from to the last available snapshot.

### 3.2. Anisotropies

Another question is the anisotropy of density and the logarithm of density structures. For subsonic turbulence it is natural to assume that the density anisotropies will mimic velocity anisotropies in GS95 picture. This was confirmed in CL03, who, however, observed that for supersonic turbulence the contours of density isocorrelation get round, corresponding to isotropy. BLC05, however, showed that anisotropies restore the GS95 form if instead of density one studies the logarithm of density. This is due to the suppression of the influence of the high density peaks, which arise from shocks. It is these peaks that mask the anisotropy of weaker, but more widely spread density fluctuations.

In Figure LABEL:fig:anisotropy we show lines that mark the corresponding separation lengths for the second-order SFs parallel and perpendicular to the local mean magnetic field.^{1}^{1}1The local mean magnetic field was computed using the procedure of smoothing by a 3D Gaussian profile with the width equal to the separation length. Because the volume of smoothing grows with the separation length , the direction of the local mean magnetic field might change with at an arbitrary point. This is an extension of the procedures employed in Cho, Lazarian & Vishniac (2002). In the case of subAlfvénic turbulence, the degree of anisotropy for density is very difficult to estimate due to the high dispersion of points. However, rough estimates suggest more isotropic density structures, because the points extend along the line . For models with the points in Figure LABEL:fig:anisotropy have lower dispersion, and the anisotropy is more like the type from GS95, i.e. . In both Alfvénic regimes, the anisotropy of density does not change significantly with . Plots for the logarithm of density show more smooth relations between parallel and perpendicular SFs. The dispersion of points is very small. Moreover, we note the change of anisotropy with the scale. Lower values of SFs correspond to lower values of the separation length (small-scale structures), so we might note that the logarithm of density structures are more isotropic than the GS95 model at small scales, but the anisotropy grows a bit larger than the GS95 prediction at larger scales. This difference is somewhat larger in the case of models with stronger external magnetic field (compare plots in the left and right columns of Figure LABEL:fig:anisotropy), which may signify their dependence on the strength of . The anisotropy of structures is marginally dependent on the sonic Mach number, similar to the density structures. All these observations allow us to confirm the previous studies (see BLC05) that suggest that the anisotropy depends not only on the scale but on .

### 3.3. Bispectrum

Attempts to use multipoint statistics are a more traditional way to remove the constraints that the use of two point statistics, e.g. power spectra entails. Unfortunately, very high quality data is needed to obtain the multipoint statistics. Among multipoint statistics, bispectrum (see Scoccimarro 1997) seems the most promising. This is partially because it has been successfully used in the studies of the Large Scale Structure of the Universe.

Bispectrum is a Fourier transform of the three point correlation function and if the power spectrum is defined as

(1) |

where is the Dirac delta function that is zero apart from the case when , the bispectrum is

(2) |

The bispectra of density shown in Fig. 4 give information as to how shocks and magnetic fields effect turbulence. It has been shown by Kowal, Lazarian & Beresnyak (2007) and Beresnyak et al. (2005) that in supersonic turbulence shocks produce compressed density and a shallower spectrum. Looking at Fig. 4 it is clear that these shocks play a crucial role in the correlation of modes. The subsonic cases show little correlation between any points except the case of . The compressed densities from supersonic turbulence are critical for correlations between frequencies since waves get closer together and therefore have a much higher interaction rate. We find an good agreement with Kowal, Lazarian & Beresnyak (2007) in that the density structures are affected by the presence of magnetic field. When there is a weak magnetic field present it is clear from Figure 4 that the system lacks the stronger correlations that are characteristic of the sub-Alfv́enic models. For super-Alfv́enic simulations correlations in frequency are not as readily made due to large dispersion of density structure. The bispectrum of supersonic hydro models are similar to the super-Alfv́enic, supersonic cases. A measure related to bispectrum is bicohence (see Koronovskii & Hramov 2002).

## 4. Summary

Above we discussed the statistics of density in compressible MHD turbulence. We analyzed spectra, structure functions for experiments with different sonic and Alfvén Mach numbers. Our results are as follows:

The viscosity-damped regime of MHD turbulence relevant to partially ionized gas can be characterized by shallow and very anisotropic spectrum of density. This spectrum can result in large variations of the column densities.

The amplitude of density fluctuations strongly depends on both in weakly and strongly magnetized turbulent plasmas.

The flattening observed in density spectra is due to the contribution of the highly dense small-scale structures generated in the supersonic turbulence.

Fluctuations of the logarithm of density are much more regular than those of density. The logarithm of density exhibits the GS95 scalings and anisotropies.

Bispectra show strong correlations for supersonic models than for subsonic ones. It suggests the importance of shocks in mode correlations.

Acknowledgement of the support of the NSF Center for Magnetic Self-Organization in Laboratory and Astrophysical Plasma.

## References

- (1) Ballesteros-Paredes, J., Klessen, R., Mac Low, M. & Vasquez-Semadeni, E. 2006, in “Protostars and Planets V”, University of Arizona Press, Tucson, p.63
- Beresnyak et al. (2005) Beresnyak, A., Lazarian, A. & Cho, J., 2005, ApJ, 624, L93, BLC05
- Biskamp (2003) Biskamp, D., 2003, Magnetohydrodynamic Turbulence, Cambridge University Press
- (4) Cho, J. & Lazarian A. 2002, Phys. Rev. Lett., 88, number 24, 5001-(1-4)
- (5) Cho, J. & Lazarian, A. 2003, MNRAS, 345, 325-339
- Cho, Lazarian & Vishniac (2002) Cho, J., Lazarian, A. & Vishniac, E.T. 2002, ApJ, 564, 291
- Cho & Vishniac (2000) Cho, J. & Vishniac, E. T., 2000, ApJ, 539, 273
- Elmegreen & Scalo (2004) Elmegreen, B. & Scalo, J. 2004, ARA&A, 42, 211
- Goldreich & Sridhar (1995) Goldreich, P. & Sridhar, H. 1995, ApJ 438, 763
- (10) Higdon, J.C. 1984, ApJ, 285, 109
- (11) Iroshnikov, P.S. 1963, AZh, 40, 742
- Kowal, Lazarian & Beresnyak (2007) Kowal, G., Lazarian, A. & Beresnyak, A., 2007, ApJ, 658, 423
- (13) Koronovskii, A. & Hramov, A. 2002, Plasma Phys. Rep. 28, 666
- (14) Kraichnan, R. 1965, Phys. Fluids, 8, 1385
- (15) Lazarian, A. 2007, in ”SINS-Small Ionized and Neutral Strucutres in Diffuse ISM”, eds. M. Haverkorn & W.M. Goss, ASP 365, 324
- (16) Lazarian, A. 2006, AIP 874, 301-315
- Lazarian, Vishniac,& Cho (2004) Lazarian, A., Vishniac, E., & Cho, J. 2004, ApJ, 603, 180
- Lazarian & Beresnyak (2005) Lazarian, A. & Beresnyak, A., 2005, in The Magnetized Plasma in Galaxy Evolution, Kraków, Jagiellonian University Press, p. 56
- (19) Maron, J. & Goldreich, P. 2001, ApJ, 554, 1175
- (20) McKee, C. & Ostriker, E. 2007, ARA&A, 45, 565
- (21) Montgomrey D.C. & Turner L. 1981, Phys. Fluids, 24, 825
- Montgomery et al. (1987) Montgomery, D., Brown, M. R., Matthaeus, W. H., 1987, JGR, 92, 282