Simulating Supersonic Turbulence in Magnetized Molecular Clouds

We present results of large-scale three-dimensional simulations of weakly magnetized supersonic turbulence at grid resolutions up to 1024^3 cells. Our numerical experiments are carried out with the Piecewise Parabolic Method on a Local Stencil and assume an isothermal equation of state. The turbulence is driven by a large-scale isotropic solenoidal force in a periodic computational domain and fully develops in a few flow crossing times. We then evolve the flow for a number of flow crossing times and analyze various statistical properties of the saturated turbulent state. We show that the energy transfer rate in the inertial range of scales is surprisingly close to a constant, indicating that Kolmogorov's phenomenology for incompressible turbulence can be extended to magnetized supersonic flows. We also discuss numerical dissipation effects and convergence of different turbulence diagnostics as grid resolution refines from 256^3 to 1024^3 cells.Comment: 10 pages, 3 figures, to appear in the proceedings of the DOE/SciDAC 2009 conferenc

Dissipative Structures in Supersonic Turbulence

We show that density-weighted moments of the dissipation rate, $\epsilon_l$, averaged over a scale $l$, in supersonic turbulence can be successfully explained by the She and L\'ev\^eque model [Phys. Rev. Lett. {\bf 72}, 336 (1994)]. A general method is developed to measure the two parameters of the model, $\gamma$ and $d$, based directly on their physical interpretations as the scaling exponent of the dissipation rate in the most intermittent structures ($\gamma$) and the dimension of the structures ($d$). We find that the best-fit parameters ($\gamma=0.71$ and $d=1.90$) derived from the $\epsilon_l$ scalings in a simulation of supersonic turbulence at Mach 6 agree with their direct measurements, confirming the validity of the model in supersonic turbulence.Comment: 4 pages, 3 figures, accepted by Phys. Rev. Let

Temperature, composition and age of the Kara Sea Shelf sediments in the area of the Marre-Sale Geocryological Station

The paper presents results of the study of the uppermost 20 m-thick layer of the near-Yamal shelf bottom sediments, penetrated in May 2014 by two VSEGINGEO boreholes equipped with LРС loggers, with an aim of the temperature regime dynamics monitoring in the nearshore bottom sediments, both for the research purposes and in as much as the data add value to the forthcoming hydrocarbon resource development on the Russian continental shelf. On the basis of the temperature variation observations during three summer months of 2014, it has been established that marine silty clays and aleurites composing the bottom sediment section, represent relict frozen deposits subjected to cryogenic metamorphism in the subaerial exposure environment. Diatom assemblages occurring in aleurite and clayey deposits consist exclusively of the marine extinct species typical of the Early Eocene Pyxilla gracilis diatom zone. A modern marine sublittoral diatom assemblage is found inhabiting the sands of the upper part of the onshore borehole section

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

Turbulent Mixing in the Interstellar Medium -- an application for Lagrangian Tracer Particles

We use 3-dimensional numerical simulations of self-gravitating compressible turbulent gas in combination with Lagrangian tracer particles to investigate the mixing process of molecular hydrogen (H2) in interstellar clouds. Tracer particles are used to represent shock-compressed dense gas, which is associated with H2. We deposit tracer particles in regions of density contrast in excess of ten times the mean density. Following their trajectories and using probability distribution functions, we find an upper limit for the mixing timescale of H2, which is of order 0.3 Myr. This is significantly smaller than the lifetime of molecular clouds, which demonstrates the importance of the turbulent mixing of H2 as a preliminary stage to star formation.Comment: 10 pages, 5 figures, conference proceedings "Turbulent Mixing and Beyond 2007

Astrophysical turbulence modeling

The role of turbulence in various astrophysical settings is reviewed. Among the differences to laboratory and atmospheric turbulence we highlight the ubiquitous presence of magnetic fields that are generally produced and maintained by dynamo action. The extreme temperature and density contrasts and stratifications are emphasized in connection with turbulence in the interstellar medium and in stars with outer convection zones, respectively. In many cases turbulence plays an essential role in facilitating enhanced transport of mass, momentum, energy, and magnetic fields in terms of the corresponding coarse-grained mean fields. Those transport properties are usually strongly modified by anisotropies and often completely new effects emerge in such a description that have no correspondence in terms of the original (non coarse-grained) fields.Comment: 88 pages, 26 figures, published in Reports on Progress in Physic

The Iron Discrepancy in Elliptical Galaxies after ASCA

We present estimates for the iron content of the stellar and diffused components of elliptical galaxies, as derived respectively from integrated optical spectra and from ASCA X-ray observations. A macroscopic discrepancy emerges between the expected iron abundances in the hot interstellar medium (ISM) and what is indicated by the X-ray observations, especially when allowance is made for the current iron enrichment by Type Ia supernovae. This strong discrepancy, that in some extreme instances may be as large as a factor of $\sim 20$, calls into question our current understanding of supernova enrichment and chemical evolution of galaxies. We discuss several astrophysical implications of the inferred low iron abundances in the ISM, including the chemical evolution of galaxies and cluster of galaxies, the evolution of gas flows in ellipticals, and the heating of the intracluster medium. Some of the consequences appear hard to accept, and in the attempt to avoid some of these difficulties we explore ways of hiding or diluting iron in the ISM of ellipticals. None of these possibilities appears astrophysically plausible, and we alternatively rise the question of the reliability of iron-L line diagonostic tools. Various thin plasma emission models are shown to give iron abundances that may differ significantly, especially at low temperatures (kT \lsim 1 keV). From a collection of ASCA and other X-ray observatory data, it is shown that current thin plasma codes tend to give very low iron abundances when the temperature of the objects is below $\sim 1$ keV. We conclude that -- besides rethinking the chemical evolution of galaxies -- one should also consider the possibility that existing thin plasma models may incorporate inaccurate atomic physics for the ions responsible for the iron-L complex.Comment: 39 pages, TeX file, 5 figures, Accepted for publication in the Ap

The <i>Herschel</i> view of the massive star-forming region NGC 6334

Aims: Fundamental to any theory of high-mass star formation are gravity and turbulence. Their relative importance, which probably changes during cloud evolution, is not known. By investigating the spatial and density structure of the high-mass star-forming complex NGC 6334 we aim to disentangle the contributions of turbulence and gravity. Methods: We used Herschel PACS and SPIRE imaging observations from the HOBYS key programme at wavelengths of 160, 250, 350, and 500 μm to construct dust temperature and column density maps. Using probability distribution functions (PDFs) of the column density determined for the whole complex and for four distinct sub-regions (distinguished on the basis of differences in the column density, temperature, and radiation field), we characterize the density structure of the complex. We investigate the spatial structure using the Δ-variance, which probes the relative amount of structure on different size scales and traces possible energy injection mechanisms into the molecular cloud. Results: The Δ-variance analysis suggests that the significant scales of a few parsec that were found are caused by energy injection due to expanding HII regions, which are numerous, and by the lengths of filaments seen everywhere in the complex. The column density PDFs have a lognormal shape at low densities and a clearly defined power law at high densities for all sub-regions whose slope is linked to the exponent α of an equivalent spherical density distribution. In particular with α = 2.37, the central sub-region is largly dominated by gravity, caused by individual collapsing dense cores and global collapse of a larger region. The collapse is faster than free-fall (which would lead only to α = 2) and thus requires a more dynamic scenario (external compression, flows). The column density PDFs suggest that the different sub-regions are at different evolutionary stages, especially the central sub-region, which seems to be in a more evolved stage

Interstellar MHD Turbulence and Star Formation

This chapter reviews the nature of turbulence in the Galactic interstellar medium (ISM) and its connections to the star formation (SF) process. The ISM is turbulent, magnetized, self-gravitating, and is subject to heating and cooling processes that control its thermodynamic behavior. The turbulence in the warm and hot ionized components of the ISM appears to be trans- or subsonic, and thus to behave nearly incompressibly. However, the neutral warm and cold components are highly compressible, as a consequence of both thermal instability in the atomic gas and of moderately-to-strongly supersonic motions in the roughly isothermal cold atomic and molecular components. Within this context, we discuss: i) the production and statistical distribution of turbulent density fluctuations in both isothermal and polytropic media; ii) the nature of the clumps produced by thermal instability, noting that, contrary to classical ideas, they in general accrete mass from their environment; iii) the density-magnetic field correlation (or lack thereof) in turbulent density fluctuations, as a consequence of the superposition of the different wave modes in the turbulent flow; iv) the evolution of the mass-to-magnetic flux ratio (MFR) in density fluctuations as they are built up by dynamic compressions; v) the formation of cold, dense clouds aided by thermal instability; vi) the expectation that star-forming molecular clouds are likely to be undergoing global gravitational contraction, rather than being near equilibrium, and vii) the regulation of the star formation rate (SFR) in such gravitationally contracting clouds by stellar feedback which, rather than keeping the clouds from collapsing, evaporates and diperses them while they collapse.Comment: 43 pages. Invited chapter for the book "Magnetic Fields in Diffuse Media", edited by Elisabete de Gouveia dal Pino and Alex Lazarian. Revised as per referee's recommendation

Ionization compression impact on dense gas distribution and star formation: probability density functions around H II regions as seen by <i>Herschel</i>

Aims. Ionization feedback should impact the probability distribution function (PDF) of the column density of cold dust around the ionized gas. We aim to quantify this effect and discuss its potential link to the core and initial mass function (CMF/IMF). Methods. We used Herschel column density maps of several regions observed within the HOBYS key program in a systematic way: M 16, the Rosette and Vela C molecular clouds, and the RCW 120 H II region. We computed the PDFs in concentric disks around the main ionizing sources, determined their properties, and discuss the effect of ionization pressure on the distribution of the column density. Results. We fitted the column density PDFs of all clouds with two lognormal distributions, since they present a "double-peak" or an enlarged shape in the PDF. Our interpretation is that the lowest part of the column density distribution describes the turbulent molecular gas, while the second peak corresponds to a compression zone induced by the expansion of the ionized gas into the turbulent molecular cloud. Such a double peak is not visible for all clouds associated with ionization fronts, but it depends on the relative importance of ionization pressure and turbulent ram pressure. A power-law tail is present for higher column densities, which are generally ascribed to the effect of gravity. The condensations at the edge of the ionized gas have a steep compressed radial profile, sometimes recognizable in the flattening of the power-law tail. This could lead to an unambiguous criterion that is able to disentangle triggered star formation from pre-existing star formation. Conclusions. In the context of the gravo-turbulent scenario for the origin of the CMF/IMF, the double-peaked or enlarged shape of the PDF may affect the formation of objects at both the low-mass and the high-mass ends of the CMF/IMF. In particular, a broader PDF is required by the gravo-turbulent scenario to fit the IMF properly with a reasonable initial Mach number for the molecular cloud. Since other physical processes (e.g., the equation of state and the variations among the core properties) have already been said to broaden the PDF, the relative importance of the different effects remains an open question