A Turbulent Model for the Interstellar Medium. II. Magnetic Fields and Rotation

We present results from two-dimensional numerical simulations of a supersonic turbulent flow in the plane of the galactic disk, incorporating shear, thresholded and discrete star formation (SF), self-gravity, rotation and magnetic fields. A test of the model in the linear regime supports the results of the linear theory of Elmegreen (1991). In the fully nonlinear turbulent regime, while some results of the linear theory persist, new effects also emerge. Some exclusively nonlinear effects are: a) Even though there is no dynamo in 2D, the simulations are able to maintain or increase their net magnetic energy in the presence of a seed uniform azimuthal component. b) A well-defined power-law magnetic spectrum and an inverse magnetic cascade are observed in the simulations, indicating full MHD turbulence. Thus, magnetic field energy is generated in regions of SF and cascades up to the largest scales. c) The field has a slight but noticeable tendency to be aligned with density features. d) The magnetic field prevents HII regions from expanding freely, as in the recent results of Slavin \& Cox (1993). e) A tendency to exhibit {\it less} filamentary structures at stronger values of the uniform component of the magnetic field is present in several magnetic runs. f) For fiducial values of the parameters, the flow in general appears to be in rough equipartition between magnetic and kinetic energy. There is no clear domination of either the magnetic or the inertial forces. g) A median value of the magnetic field strength within clouds is $\sim 12\mu$G, while for the intercloud medium a value of $\sim 3\mu$G is found. Maximum contrasts of up to a factor of $\sim 10$ are observed.Comment: Plain TeX file, 25 pages. Gzipped, tarred set of Tex file plus 17 figures and 3 tables (Postscript) available at ftp://kepler.astroscu.unam.mx/incoming/enro/papers/mhdgturb.tar.g

Highly Compressible MHD Turbulence and Gravitational Collapse

We investigate the properties of highly compressible turbulence and its ability to produce self-gravitating structures. The compressibility is parameterized by an effective polytropic exponent gama-eff. In the limit of small gama-eff, the density jump at shocks is shown to be of the order of e^{M^2}, and the production of vorticity by the nonlinear terms appears to be negligible. In the presence of self-gravity, we suggest that turbulence can produce bound structures for gama-eff < 2(1-1/n), where 'n' is the typical dimensionality of the turbulent compressions. We show, by means of numerical simulations, that, for sufficiently small gama-eff, small-scale turbulent density fluctuations eventually collapse even though the medium is globally stable. This result is preserved in the presence of a magnetic field for supercritical mass-to-flux ratios.Comment: 4 pages, 3 postscript figures. Latex, uses aipproc.sty Contribution to the Conference Proc. of the 7th Annual Astrophysics Conference in Maryland, STAR FORMATION, NEAR AND FAR, eds. Stephen S. Holt and Lee G. Mund

Density probability distribution in one-dimensional polytropic gas dynamics

We discuss the generation and statistics of the density fluctuations in highly compressible polytropic turbulence, based on a simple model and one-dimensional numerical simulations. Observing that density structures tend to form in a hierarchical manner, we assume that density fluctuations follow a random multiplicative process. When the polytropic exponent $\gamma$ is equal to unity, the local Mach number is independent of the density, and our assumption leads us to expect that the probability density function (PDF) of the density field is a lognormal. This isothermal case is found to be singular, with a dispersion $\sigma_s^2$ which scales like the square turbulent Mach number $\tilde M^2$, where $s\equiv \ln \rho$ and $\rho$ is the fluid density. This leads to much higher fluctuations than those due to shock jump relations. Extrapolating the model to the case $\gamma \not =1$, we find that, as the Mach number becomes large, the density PDF is expected to asymptotically approach a power-law regime, at high densities when $\gamma<1$, and at low densities when $\gamma>1$. This effect can be traced back to the fact that the pressure term in the momentum equation varies exponentially with $s$, thus opposing the growth of fluctuations on one side of the PDF, while being negligible on the other side. This also causes the dispersion $\sigma_s^2$ to grow more slowly than $\tilde M^2$ when $\gamma\not=1$. In view of these results, we suggest that Burgers flow is a singular case not approached by the high-$\tilde M$ limit, with a PDF that develops power laws on both sides.Comment: 9 pages + 12 postscript figures. Submitted to Phys. Rev.

Macroscopic and Local Magnetic Moments in Si-doped CuGeO3_3 with Neutron and μ\muSR Studies

The temperature-concentration phase diagram of the Si-doped spin-Peierls compound CuGeO$_{3}$ is investigated by means of neutron scattering and muon spin rotation spectroscopy in order to determine the microscopic distribution of the magnetic and lattice dimerised regions as a function of doping. The analysis of the zero-field muon spectra has confirmed the spatial inhomogeneity of the staggered magnetisation that characterises the antiferromagnetic superlattice peaks observed with neutrons. In addition, the variation of the macroscopic order parameter with doping can be understood by considering the evolution of the local magnetic moment as well as of the various regions contributing to the muon signal

Influence of Cooling-Induced Compressibility on the Structure of Turbulent Flows and Gravitational Collapse

We investigate the properties of highly compressible turbulence, the compressibility arising from a small effective polytropic exponent $\gamma_e$ due to cooling. In the limit of small $\gamma_e$, the density jump at shocks is shown to be of the order of $e^{M^2}$. Without self-gravity, the density structures arising in the moderately compressible case consist mostly of patches separated by shocks and behaving like waves, while in the highly compressible case clearly defined long-lived object-like clouds emerge. When the forcing in the momentum equation is purely compressible, the rotational energy decays monotonically in time, indicating that the vortex-stretching term is not efficient in transferring energy to rotational modes. This property may be at the origin of the low amount of rotation found in interstellar clouds. Vorticity production is found to rely heavily on the presence of additional terms in the equations. In the presence of self-gravity, we suggest that turbulence can produce bound structures for $\gamma_e < 2(1-1/n)$, where $n$ is the typical dimensionality of the turbulent compressions. We support this result by means of numerical simulations in which, for sufficiently small $\gamma_e$, small-scale turbulent density fluctuations eventually collapse even though the medium is globally stable. This result is preserved in the presence of a magnetic field for supercritical mass-to-flux ratios. At larger polytropic exponents, turbulence alone is not capable of producing bound structures, and collapse can only occur when the medium is globally unstable. This mechanism is a plausible candidate for the differentiation between primordial and present-day stellar-cluster formation and for the low efficiency of star formation.Comment: 20 pages, 12 Postscript figures. Uses aas2pp4.sty. Accepted in Ap

A Holistic Scenario of Turbulent Molecular Cloud Evolution and Control of the Star Formation Efficiency. First Tests

We compile a holistic scenario for molecular cloud (MC) evolution and control of the star formation efficiency (SFE), and present a first set of numerical tests of it. A {\it lossy} compressible cascade can generate density fluctuations and further turbulence at small scales from large-scale motions, implying that the turbulence in MCs may originate from the compressions that form them. Below a {\it sonic} scale \ls, turbulence cannot induce any further subfragmentation, nor be a dominant support agent against gravity. Since progressively smaller density peaks contain progressively smaller fractions of the mass, we expect the SFE to decrease with decreasing \ls, at least when the cloud is globally supported by turbulence. Our numerical experiments confirm this prediction. We also find that the collapsed mass fraction in the simulations always saturates below 100% efficiency. This may be due to the decreased mean density of the leftover interclump medium, which in real clouds (not confined to a box) should then be more easily dispersed, marking the ``death'' of the cloud. We identify two different functional dependences (``modes'') of the SFE on \ls, which roughly correspond to globally supported and unsupported cases. Globally supported runs with most of the turbulent energy at the largest scales have similar SFEs to those of unsupported runs, providing numerical evidence of the dual role of turbulence, whereby large-scale turbulent modes induce collapse at smaller scales. We tentatively suggest that these modes may correspond to the clustered and isolated modes of star formation, although here they are seen to form part of a continuum rather than being separate modes. Finally, we compare with previous proposals that the relevant parameter is the energy injection scale.Comment: 6 pages, 3 figures. Uses emulateapj. Accepted in ApJ Letter

Do Lognormal Column-Density Distributions in Molecular Clouds Imply Supersonic Turbulence?

Recent observations of column densities in molecular clouds find lognormal distributions with power-law high-density tails. These results are often interpreted as indications that supersonic turbulence dominates the dynamics of the observed clouds. We calculate and present the column-density distributions of three clouds, modeled with very different techniques, none of which is dominated by supersonic turbulence. The first star-forming cloud is simulated using smoothed particle hydrodynamics (SPH); in this case gravity, opposed only by thermal-pressure forces, drives the evolution. The second cloud is magnetically subcritical with subsonic turbulence, simulated using nonideal MHD; in this case the evolution is due to gravitationally-driven ambipolar diffusion. The third cloud is isothermal, self-gravitating, and has a smooth density distribution analytically approximated with a uniform inner region and an r^-2 profile at larger radii. We show that in all three cases the column-density distributions are lognormal. Power-law tails develop only at late times (or, in the case of the smooth analytic profile, for strongly centrally concentrated configurations), when gravity dominates all opposing forces. It therefore follows that lognormal column-density distributions are generic features of diverse model clouds, and should not be interpreted as being a consequence of supersonic turbulence.Comment: 6 pages, 6 figures, accepted for publication in MNRA

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