Sub-Alfvenic Non-Ideal MHD Turbulence Simulations with Ambipolar Diffusion: I. Turbulence Statistics

Most numerical investigations on the role of magnetic fields in turbulent molecular clouds (MCs) are based on ideal magneto-hydrodynamics (MHD). However, MCs are weakly ionized, so that the time scale required for the magnetic field to diffuse through the neutral component of the plasma by ambipolar diffusion (AD) can be comparable to the dynamical time scale. We have performed a series of 256^3 and 512^3 simulations on supersonic but sub-Alfvenic turbulent systems with AD using the Heavy-Ion Approximation developed in Li, McKee, & Klein (2006). Our calculations are based on the assumption that the number of ions is conserved, but we show that these results approximately apply to the case of time-dependent ionization in molecular clouds as well. Convergence studies allow us to determine the optimal value of the ionization mass fraction when using the heavy-ion approximation for low Mach number, sub-Alfvenic turbulent systems. We find that ambipolar diffusion steepens the velocity and magnetic power spectra compared to the ideal MHD case. Changes in the density PDF, total magnetic energy, and ionization fraction are determined as a function of the AD Reynolds number. The power spectra for the neutral gas properties of a strongly magnetized medium with a low AD Reynolds number are similar to those for a weakly magnetized medium; in particular, the power spectrum of the neutral velocity is close to that for Burgers turbulence.Comment: 37 pages, 11 figures, 4 table

Gravitational Collapse in Turbulent Molecular Clouds. I. Gasdynamical Turbulence

Observed molecular clouds often appear to have very low star formation efficiencies and lifetimes an order of magnitude longer than their free-fall times. Their support is attributed to the random supersonic motions observed in them. We study the support of molecular clouds against gravitational collapse by supersonic, gas dynamical turbulence using direct numerical simulation. Computations with two different algorithms are compared: a particle-based, Lagrangian method (SPH), and a grid-based, Eulerian, second-order method (ZEUS). The effects of both algorithm and resolution can be studied with this method. We find that, under typical molecular cloud conditions, global collapse can indeed be prevented, but density enhancements caused by strong shocks nevertheless become gravitationally unstable and collapse into dense cores and, presumably, stars. The occurance and efficiency of local collapse decreases as the driving wave length decreases and the driving strength increases. It appears that local collapse can only be prevented entirely with unrealistically short wave length driving, but observed core formation rates can be reproduced with more realistic driving. At high collapse rates, cores are formed on short time scales in coherent structures with high efficiency, while at low collapse rates they are scattered randomly throughout the region and exhibit considerable age spread. We suggest that this naturally explains the observed distinction between isolated and clustered star formation.Comment: Minor revisions in response to referee, thirteen figures, accepted to Astrophys.

The Statistics of Supersonic Isothermal Turbulence

We present results of large-scale three-dimensional simulations of supersonic Euler turbulence with the piecewise parabolic method and multiple grid resolutions up to 2048^3 points. Our numerical experiments describe non-magnetized driven turbulent flows with an isothermal equation of state and an rms Mach number of 6. We discuss numerical resolution issues and demonstrate convergence, in a statistical sense, of the inertial range dynamics in simulations on grids larger than 512^3 points. The simulations allowed us to measure the absolute velocity scaling exponents for the first time. The inertial range velocity scaling in this strongly compressible regime deviates substantially from the incompressible Kolmogorov laws. The slope of the velocity power spectrum, for instance, is -1.95 compared to -5/3 in the incompressible case. The exponent of the third-order velocity structure function is 1.28, while in incompressible turbulence it is known to be unity. We propose a natural extension of Kolmogorov's phenomenology that takes into account compressibility by mixing the velocity and density statistics and preserves the Kolmogorov scaling of the power spectrum and structure functions of the density-weighted velocity v=\rho^{1/3}u. The low-order statistics of v appear to be invariant with respect to changes in the Mach number. For instance, at Mach 6 the slope of the power spectrum of v is -1.69, and the exponent of the third-order structure function of v is unity. We also directly measure the mass dimension of the "fractal" density distribution in the inertial subrange, D_m = 2.4, which is similar to the observed fractal dimension of molecular clouds and agrees well with the cascade phenomenology.Comment: 15 pages, 19 figures, ApJ v665, n2, 200

Statistics of Dissipation and Enstrophy Induced by a Set of Burgers Vortices

Dissipation and enstropy statistics are calculated for an ensemble of modified Burgers vortices in equilibrium under uniform straining. Different best-fit, finite-range scaling exponents are found for locally-averaged dissipation and enstrophy, in agreement with existing numerical simulations and experiments. However, the ratios of dissipation and enstropy moments supported by axisymmetric vortices of any profile are finite. Therefore the asymptotic scaling exponents for dissipation and enstrophy induced by such vortices are equal in the limit of infinite Reynolds number.Comment: Revtex (4 pages) with 4 postscript figures included via psfi

Turbulent dissipation in the ISM: the coexistence of forced and decaying regimes and implications for galaxy formation and evolution

We discuss the dissipation of turbulent kinetic energy Ek in the global ISM by means of 2-D, MHD, non-isothermal simulations in the presence of model radiative heating and cooling. We argue that dissipation in 2D is representative of that in three dimensions as long as it is dominated by shocks rather than by a turbulent cascade. Energy is injected at a few isolated sites in space, over relatively small scales, and over short time periods. This leads to the coexistence of forced and decaying regimes in the same flow. We find that the ISM-like flow dissipates its turbulent energy rapidly. In simulations with forcing, the input parameters are the radius l_f of the forcing region, the total kinetic energy e_k each source deposits into the flow, and the rate of formation of those regions, sfr_OB. The global dissipation time t_d depends mainly on l_f. In terms of measurable properties of the ISM, t_d >= Sigma_g u_rms^2/(e_k sfr_OB), where Sigma_g is the average gas surface density and u_rms is the rms velocity dispersion. For the solar neighborhood, t_d >= 1.5x10^7 yr. The global dissipation time is consistently smaller than the crossing time of the largest energy-containing scales. In decaying simulations, Ek decreases with time as t^-n, where n~0.8-0.9. This suggests a decay with distance d as Ek\propto d^{-2n/(2-n)} in the mixed forced+decaying case. If applicable to the vertical direction, our results support models of galaxy evolution in which stellar energy injection provides significant support for the gas disk thickness, but not models of galaxy formation in which this energy injection is supposed to reheat an intra-halo medium at distances of up to 10-20 times the optical galaxy size, as the dissipation occurs on distances comparable to the disk height.Comment: 23 pages, including figures. To appear in ApJ. Abstract abridge

The Fractal Density Structure in Supersonic Isothermal Turbulence: Solenoidal versus Compressive Energy Injection

In a systematic study, we compare the density statistics in high resolution numerical experiments of supersonic isothermal turbulence, driven by the usually adopted solenoidal (divergence-free) forcing and by compressive (curl-free) forcing. We find that for the same rms Mach number, compressive forcing produces much stronger density enhancements and larger voids compared to solenoidal forcing. Consequently, the Fourier spectra of density fluctuations are significantly steeper. This result is confirmed using the Delta-variance analysis, which yields power-law exponents beta~3.4 for compressive forcing and beta~2.8 for solenoidal forcing. We obtain fractal dimension estimates from the density spectra and Delta-variance scaling, and by using the box counting, mass size and perimeter area methods applied to the volumetric data, projections and slices of our turbulent density fields. Our results suggest that compressive forcing yields fractal dimensions significantly smaller compared to solenoidal forcing. However, the actual values depend sensitively on the adopted method, with the most reliable estimates based on the Delta-variance, or equivalently, on Fourier spectra. Using these methods, we obtain D~2.3 for compressive and D~2.6 for solenoidal forcing, which is within the range of fractal dimension estimates inferred from observations (D~2.0-2.7). The velocity dispersion to size relations for both solenoidal and compressive forcing obtained from velocity spectra follow a power law with exponents in the range 0.4-0.5, in good agreement with previous studies.Comment: 17 pages, 11 figures, ApJ in press, minor changes to language, simulation movies available at http://www.ita.uni-heidelberg.de/~chfeder/videos.shtml?lang=e

A Super-Alfvenic Model of Dark Clouds

Supersonic random motions are observed in dark clouds and are traditionally interpreted as Alfven waves, but the possibility that these motions are super-Alfvenic has not been ruled out. In this work we report the results of numerical experiments in two opposite regimes; M_a ~ 1 and M_a >> 1, where M_a is the initial Alfvenic Mach number --the ratio of the rms velocity to the Alfven speed. Our results show that models with M_a >> 1 are consistent with the observed properties of molecular clouds that we have tested --statistics of extinction measurements, Zeeman splitting measurements of magnetic field strength, line width versus integrated antenna temperature of molecular emission line spectra, statistical B-n relation, and scatter in that relation-- while models with M_a ~ 1 have properties that are in conflict with the observations. We find that both the density and the magnetic field in molecular clouds may be very intermittent. The statistical distributions of magnetic field and gas density are related by a power law, with an index that decreases with time in experiments with decaying turbulence. After about one dynamical time it stabilizes at B ~ n^{0.4}. Magnetically dominated cores form early in the evolution, while later on the intermittency in the density field wins out, and also cores with weak field can be generated, by mass accretion along magnetic field lines.Comment: 10 figures, 2 tables include

Nonlinear theory of mirror instability near threshold

An asymptotic model based on a reductive perturbative expansion of the drift kinetic and the Maxwell equations is used to demonstrate that, near the instability threshold, the nonlinear dynamics of mirror modes in a magnetized plasma with anisotropic ion temperatures involves a subcritical bifurcation,leading to the formation of small-scale structures with amplitudes comparable with the ambient magnetic field

Vorticity alignment results for the three-dimensional Euler and Navier-Stokes equations

We address the problem in Navier-Stokes isotropic turbulence of why the vorticity accumulates on thin sets such as quasi-one-dimensional tubes and quasi-two-dimensional sheets. Taking our motivation from the work of Ashurst, Kerstein, Kerr and Gibson, who observed that the vorticity vector {\boldmath\omega} aligns with the intermediate eigenvector of the strain matrix $S$, we study this problem in the context of both the three-dimensional Euler and Navier-Stokes equations using the variables \alpha = \hat{{\boldmath\xi}}\cdot S\hat{{\boldmath\xi}} and {\boldmath\chi} = \hat{{\boldmath\xi}}\times S\hat{{\boldmath\xi}} where \hat{{\boldmath\xi}} = {\boldmath\omega}/\omega. This introduces the dynamic angle $\phi (x,t) = \arctan(\frac{\chi}{\alpha})$, which lies between {\boldmath\omega} and S{\boldmath\omega}. For the Euler equations a closed set of differential equations for $\alpha$ and {\boldmath\chi} is derived in terms of the Hessian matrix of the pressure $P = \{p_{,ij}\}$. For the Navier-Stokes equations, the Burgers vortex and shear layer solutions turn out to be the Lagrangian fixed point solutions of the equivalent (\alpha,{\boldmath\chi}) equations with a corresponding angle $\phi = 0$. Under certain assumptions for more general flows it is shown that there is an attracting fixed point of the (\alpha,\bchi) equations which corresponds to positive vortex stretching and for which the cosine of the corresponding angle is close to unity. This indicates that near alignment is an attracting state of the system and is consistent with the formation of Burgers-like structures.Comment: To appear in Nonlinearity Nov. 199

Turbulent Control of the Star Formation Efficiency

Supersonic turbulence plays a dual role in molecular clouds: On one hand, it contributes to the global support of the clouds, while on the other it promotes the formation of small-scale density fluctuations, identifiable with clumps and cores. Within these, the local Jeans length \Ljc is reduced, and collapse ensues if \Ljc becomes smaller than the clump size and the magnetic support is insufficient (i.e., the core is ``magnetically supercritical''); otherwise, the clumps do not collapse and are expected to re-expand and disperse on a few free-fall times. This case may correspond to a fraction of the observed starless cores. The star formation efficiency (SFE, the fraction of the cloud's mass that ends up in collapsed objects) is smaller than unity because the mass contained in collapsing clumps is smaller than the total cloud mass. However, in non-magnetic numerical simulations with realistic Mach numbers and turbulence driving scales, the SFE is still larger than observational estimates. The presence of a magnetic field, even if magnetically supercritical, appears to further reduce the SFE, but by reducing the probability of core formation rather than by delaying the collapse of individual cores, as was formerly thought. Precise quantification of these effects as a function of global cloud parameters is still needed.Comment: Invited review for the conference "IMF@50: the Initial Mass Function 50 Years Later", to be published by Kluwer Academic Publishers, eds. E. Corbelli, F. Palla, and H. Zinnecke