Star formation is a complex multi-scale phenomenon that is of significant importance for astrophysics in general. Stars and star formation are key pillars in observational astronomy from local star forming regions in the Milky Way up to high-redshift galaxies. From a theoretical perspective, star formation and feedback processes (radiation, winds, and supernovae) play a pivotal role in advancing our understanding of the physical processes at work, both individually and of their interactions. In this review we will give an overview of the main processes that are important for the understanding of star formation. We start with an observationally motivated view on star formation from a global perspective and outline the general paradigm of the life-cycle of molecular clouds, in which star formation is the key process to close the cycle. After that we focus on the thermal and chemical aspects in star forming regions, discuss turbulence and magnetic fields as well as gravitational forces. Finally, we review the most important stellar feedback mechanisms.

This white paper addresses key challenges for the design of next-decade Cosmic Microwave Background (CMB) experiments, and for assessing their capability to extract cosmological information from CMB polarization. We focus here on the challenges posed by foreground emission, CMB lensing, and instrumental systematics to detect the signal that arises from gravitational waves sourced by inflation and parameterized by $r$, at the level of $r \sim 10^{-3}$ or lower, as proposed for future observational efforts. We argue that more accurate and robust analysis and simulation tools are required for these experiments to realize their promise. We are optimistic that the capability to simulate the joint impact of foregrounds, CMB lensing, and systematics can be developed to the level necessary to support the design of a space mission at $r \sim 10^{-4}$ in a few years. We make the case here for supporting such work. Although ground-based efforts present additional challenges (e.g., atmosphere, ground pickup), which are not addressed here, they would also benefit from these improved simulation capabilities.

We show that simulations of magnetohydrodynamic turbulence in the multiphase interstellar medium yield an $E/B$ ratio for polarized emission from Galactic dust in broad agreement with recent $Planck$ measurements. In addition, the $B$-mode spectra display a scale dependence that is consistent with observations over the range of scales resolved in the simulations. The simulations present an opportunity to understand the physical origin of the $E/B$ ratio, and a starting point for more refined models of Galactic emission of use for both current and future CMB experiments.

We present results of high-resolution numerical simulations of compressible 2D turbulence forced at intermediate spatial scales with a solenoidal white-in-time external acceleration. A case with an isothermal equation of state, low energy injection rate, and turbulent Mach number $M\approx0.34$ without energy condensate is studied in detail. Analysis of energy spectra and fluxes shows that the classical dual-cascade picture familiar from the incompressible case is substantially modified by compressibility effects. While the small-scale direct enstrophy cascade remains largely intact, a large-scale energy flux loop forms with the direct acoustic energy cascade compensating for the inverse transfer of solenoidal kinetic energy. At small scales, the direct enstrophy and acoustic energy cascades are fully decoupled at small Mach numbers and hence the corresponding spectral energy slopes comply with theoretical predictions, as expected. At large scales, dispersion of acoustic waves on vortices softens the dilatational velocity spectrum, while the pseudo-sound component of the potential energy associated with coherent vortices steepens the potential energy spectrum.

Three-dimensional, compressible, magnetohydrodynamic turbulence of an isothermal, self-gravitating fluid is analyzed using two-point statistics in the asymptotic limit of large Reynolds numbers (both kinetic and magnetic). Following an alternative formulation proposed by S. Banerjee and S. Galtier (Phys. Rev. E,93, 033120, 2016) and S. Banerjee and S. Galtier (J. Phys. A, Math. and Theor.,50, 015501, 2017), an exact relation has been derived for the total energy transfer. This approach results in a simpler relation expressed entirely in terms of mixed second-order structure functions. The kinetic, thermodynamic, magnetic and gravitational contributions to the energy transfer rate can be easily separated in the present form. By construction, the new formalism includes such additional effects as global rotation, the Hall term in the induction equation, etc. The analysis shows that solid-body rotation cannot alter the energy flux rate of compressible turbulence. However, the contribution of a uniform background magnetic field to the flux is shown to be non-trivial unlike in the incompressible case. Finally, the compressible, turbulent energy flux rate does not vanish completely due to simple alignments, which leads to a zero turbulent energy flux rate in the incompressible case.

Self-organization is a property of dissipative nonlinear processes that are governed by an internal driver and a positive feedback mechanism, which creates regular geometric and/or temporal patterns and decreases the entropy, in contrast to random processes. Here we investigate for the first time a comprehensive number of 16 self-organization processes that operate in planetary physics, solar physics, stellar physics, galactic physics, and cosmology. Self-organizing systems create spontaneous {\sl order out of chaos}, during the evolution from an initially disordered system to an ordered stationary system, via quasi-periodic limit-cycle dynamics, harmonic mechanical resonances, or gyromagnetic resonances. The internal driver can be gravity, rotation, thermal pressure, or acceleration of nonthermal particles, while the positive feedback mechanism is often an instability, such as the magneto-rotational instability, the Rayleigh-B\'enard convection instability, turbulence, vortex attraction, magnetic reconnection, plasma condensation, or loss-cone instability. Physical models of astrophysical self-organization processes involve hydrodynamic, MHD, and N-body formulations of Lotka-Volterra equation systems.

Self-gravitating isothermal supersonic turbulence is analyzed in the asymptotic limit of large Reynolds numbers. Based on the inviscid invariance of total energy, an exact relation is derived for homogeneous, (not necessarily isotropic) turbulence. A modified definition for the two-point energy correlation functions is used to comply with the requirement of detailed energy equipartition in the acoustic limit. In contrast to the previous relations (Galtier and Banerjee, Phys. Rev. Lett., 107, 134501, 2011; Banerjee and Galtier, Phys. Rev. E, 87, 013019, 2013), the current exact relation shows that the pressure dilatation terms plays practically no role in the energy cascade. Both the flux and source terms are written in terms of two-point differences. Sources enter the relation in a form of mixed second-order structure functions. Unlike kinetic and thermodynamic potential energy, gravitational contribution is absent from the flux term. An estimate shows that for the isotropic case, the correlation between density and gravitational acceleration may play an important role in modifying the energy transfer in self-gravitating turbulence. The exact relation is also written in an alternative form in terms of two-point correlation functions, which is then used to describe scale-by-scale energy budget in spectral space.

Large-scale turbulence in fluid layers and other quasi-two-dimensional compressible systems consists of planar vortices and waves. Separately, wave turbulence usually produces a direct energy cascade, while solenoidal planar turbulence transports energy to large scales by an inverse cascade. Here, we consider turbulence at finite Mach numbers when the interaction between acoustic waves and vortices is substantial. We employ solenoidal pumping at intermediate scales and show how both direct and inverse energy cascades are formed starting from the pumping scale. We show that there is an inverse cascade of kinetic energy up to a scale $\ell$, where a typical velocity reaches the speed of sound; this creates shock waves, which provide for a compensating direct cascade. When the system size is less than $\ell$, the steady state contains a system-size pair of long-living condensate vortices connected by a system of shocks. Thus turbulence in fluid layers processes energy via a loop: Most energy first goes to large scales via vortices and is then transported by waves to small-scale dissipation.

We explore the structure and statistics of multiphase, magnetized ISM turbulence in the local Milky Way by means of driven periodic box numerical MHD simulations. Using the higher order-accurate piecewise-parabolic method on a local stencil (PPML), we carry out a small parameter survey varying the mean magnetic field strength and density while fixing the rms velocity to observed values. We quantify numerous characteristics of the transient and steady-state turbulence, including its thermodynamics and phase structure, kinetic and magnetic energy power spectra, structure functions, and distribution functions of density, column density, pressure, and magnetic field strength. The simulations reproduce many observables of the local ISM, including molecular clouds, such as the ratio of turbulent to mean magnetic field at 100 pc scale, the mass and volume fractions of thermally stable HI, the lognormal distribution of column densities, the mass-weighted distribution of thermal pressure, and the linewidth-size relationship for molecular clouds. Our models predict the shape of magnetic field probability density functions (PDFs), which are strongly non-Gaussian, and the relative alignment of magnetic field and density structures. Finally, our models show how the observed low rates of star formation per free-fall time are controlled by the multiphase thermodynamics and large-scale turbulence.

We revisit the origin of Larson's scaling laws describing the structure and kinematics of molecular clouds. Our analysis is based on recent observational measurements and data from a suite of six simulations of the interstellar medium, including effects of self-gravity, turbulence, magnetic field, and multiphase thermodynamics. Simulations of isothermal supersonic turbulence reproduce observed slopes in linewidth-size and mass-size relations. Whether or not self-gravity is included, the linewidth-size relation remains the same. The mass-size relation, instead, substantially flattens below the sonic scale, as prestellar cores start to form. Our multiphase models with magnetic field and domain size 200 pc reproduce both scaling and normalization of the first Larson law. The simulations support a turbulent interpretation of Larson's relations. This interpretation implies that: (i) the slopes of linewidth-size and mass-size correlations are determined by the inertial cascade; (ii) none of the three Larson laws is fundamental; (iii) instead, if one is known, the other two follow from scale invariance of the kinetic energy transfer rate. It does not imply that gravity is dynamically unimportant. The self-similarity of structure established by the turbulence breaks in star-forming clouds due to the development of gravitational instability in the vicinity of the sonic scale. The instability leads to the formation of prestellar cores with the characteristic mass set by the sonic scale. The high-end slope of the core mass function predicted by the scaling relations is consistent with the Salpeter power-law index.

Supersonic turbulence plays an important role in a number of extreme astrophysical and terrestrial environments, yet its understanding remains rudimentary. We use data from a three-dimensional simulation of supersonic isothermal turbulence to reconstruct an exact fourth-order relation derived analytically from the Navier-Stokes equations (Galtier and Banerjee, Phys. Rev. Lett., vol. 107, 2011, p. 134501). Our analysis supports a Kolmogorov-like inertial energy cascade in supersonic turbulence previously discussed on a phenomenological level. We show that two compressible analogues of the four-fifths law exist describing fifth- and fourth-order correlations, but only the fourth-order relation remains `universal' in a wide range of Mach numbers from incompressible to highly compressible regimes. A new approximate relation valid in the strongly supersonic regime is derived and verified. We also briefly discuss the origin of bottleneck bumps in simulations of compressible turbulence.

This paper describes the open-source code Enzo, which uses block-structured adaptive mesh refinement to provide high spatial and temporal resolution for modeling astrophysical fluid flows. The code is Cartesian, can be run in 1, 2, and 3 dimensions, and supports a wide variety of physics including hydrodynamics, ideal and non-ideal magnetohydrodynamics, N-body dynamics (and, more broadly, self-gravity of fluids and particles), primordial gas chemistry, optically-thin radiative cooling of primordial and metal-enriched plasmas (as well as some optically-thick cooling models), radiation transport, cosmological expansion, and models for star formation and feedback in a cosmological context. In addition to explaining the algorithms implemented, we present solutions for a wide range of test problems, demonstrate the code's parallel performance, and discuss the Enzo collaboration's code development methodology.

Comparisons of the integrated thermal pressure support of gas against its gravitational potential energy lead to critical mass scales for gravitational instability such as the Jeans and the Bonnor-Ebert masses, which play an important role in analysis of many physical systems, including the heuristics of numerical simulations. In a strict theoretical sense, however, neither the Jeans nor the Bonnor-Ebert mass are meaningful when applied locally to substructure in a self-gravitating turbulent medium. For this reason, we investigate the local support by thermal pressure, turbulence, and magnetic fields against gravitational compression through an approach that is independent of these concepts. At the centre of our approach is the dynamical equation for the divergence of the velocity field. We carry out a statistical analysis of the source terms of the local compression rate (the negative time derivative of the divergence) for simulations of forced self-gravitating turbulence in periodic boxes with zero, weak, and moderately strong mean magnetic fields (measured by the averages of the magnetic and thermal pressures). We also consider the amplification of the magnetic field energy by shear and by compression. Thereby, we are able to demonstrate that the support against gravity is dominated by thermal pressure fluctuations, although magnetic pressure also yields a significant contribution. The net effect of turbulence in the highly supersonic regime, however, is to enhance compression rather than supporting overdense gas even if the vorticity is very high. This is incommensurate with the support of the highly dynamical substructures in magneto-turbulent fluids being determined by local virial equilibria of volume energies without surface stresses.

Using data from a large-scale three-dimensional simulation of supersonic isothermal turbulence, we have tested the validity of an exact flux relation derived analytically from the Navier--Stokes equation by Falkovich, Fouxon and Oz [2010 New relations for correlation functions in Navier--Stokes turbulence. J. Fluid Mech. 644, 465]. That relation, for compressible barotropic fluids, was derived assuming turbulence generated by a large-scale force. However, compressible turbulence in simulations is usually initialized and maintained by a large-scale acceleration, as in gravity-driven astrophysical flows. We present a new approximate flux relation for isothermal turbulence driven by a large-scale acceleration, and find it in reasonable agreement with the simulation results.

We revisit the origin of Larson's scaling relations, which describe the structure and kinematics of molecular clouds, based on recent observations and large-scale simulations of supersonic turbulence. Using dimensional analysis, we first show that both linewidth-size and mass-size correlations observed on scales 0.1-50 pc can be explained by a simple conceptual theory of compressible turbulence without resorting to the often assumed virial equilibrium or detailed energy balance condition. The scaling laws can be consistently interpreted as a signature of supersonic turbulence with no need to invoke gravity. We then show how self-similarity of structure established by the turbulence breaks in star-forming clouds through development of gravitational instabilities in the vicinity of the sonic scale, l_s~0.1 pc, leading to the formation of prestellar cores.

We examine the effects of self-gravity and magnetic fields on supersonic turbulence in isothermal molecular clouds with high resolution simulations and adaptive mesh refinement. These simulations use large root grids (512^3) to capture turbulence and four levels of refinement to capture high density, for an effective resolution of 8,196^3. Three Mach 9 simulations are performed, two super-Alfv\'enic and one trans-Alfv\'enic. We find that gravity splits the clouds into two populations, one low density turbulent state and one high density collapsing state. The low density state exhibits properties similar to non-self-gravitating in this regime, and we examine the effects of varied magnetic field strength on statistical properties: the density probability distribution function is approximately lognormal; velocity power spectral slopes decrease with field strength; alignment between velocity and magnetic field increases with field; the magnetic field probability distribution can be fit to a stretched exponential. The high density state is characterized by self-similar spheres; the density PDF is a power-law; collapse rate decreases with increasing mean field; density power spectra have positive slopes, P({\rho},k) \propto k; thermal-to-magnetic pressure ratios are unity for all simulations; dynamic-to-magnetic pressure ratios are larger than unity for all simulations; magnetic field distribution is a power-law. The high Alfv\'en Mach numbers in collapsing regions explain recent observations of magnetic influence decreasing with density. We also find that the high density state is found in filaments formed by converging flows, consistent with recent Herschel observations. Possible modifications to existing star formation theories are explored.

We study clustering of inertial particles in turbulent flows and discuss its applications to dust particles in protoplanetary disks. Using numerical simulations, we compute the radial distribution function (RDF), which measures the probability of finding particle pairs at given distances, and the probability density function of the particle concentration. The clustering statistics depend on the Stokes number, $St$, defined as the ratio of the particle friction timescale, $\tau_{\rm p} $, to the Kolmogorov timescale in the flow. In the dissipation range, the clustering intensity strongly peaks at $St \simeq 1$, and the RDF for $St \sim 1$ shows a fast power-law increase toward small scales, suggesting that turbulent clustering may considerably enhance the particle collision rate. Clustering at inertial-range scales is of particular interest to the problem of planetesimal formation. At these scales, the strongest clustering is from particles with $\tau_{\rm p}$ in the inertial range. Clustering of these particles occurs primarily around a scale where the eddy turnover time is $\sim\tau_{\rm p}$. Particles of different sizes tend to cluster at different locations, leading to flat RDFs between different particles at small scales. In the presence of multiple particle sizes, the overall clustering strength decreases as the particle size distribution broadens. We discuss particle clustering in recent models for planetesimal formation. We point out that, in the model based on turbulent clustering of chondrule-size particles, the probability of finding strong clusters that can seed planetesimals may have been significantly overestimated. We discuss various clustering mechanisms in simulations of planetesimal formation by gravitational collapse of dense clumps of meter-size particles, in particular the contribution from turbulent clustering due to the limited numerical resolution.

We employ simulations of supersonic super-Alfvenic turbulence decay as a benchmark test problem to assess and compare the performance of nine astrophysical MHD methods actively used to model star formation. The set of nine codes includes: ENZO, FLASH, KT-MHD, LL-MHD, PLUTO, PPML, RAMSES, STAGGER, and ZEUS. We present a comprehensive set of statistical measures designed to quantify the effects of numerical dissipation in these MHD solvers. We compare power spectra for basic fields to determine the effective spectral bandwidth of the methods and rank them based on their relative effective Reynolds numbers. We also compare numerical dissipation for solenoidal and dilatational velocity components to check for possible impacts of the numerics on small-scale density statistics. Finally, we discuss convergence of various characteristics for the turbulence decay test and impacts of various components of numerical schemes on the accuracy of solutions. We show that the best performing codes employ a consistently high order of accuracy for spatial reconstruction of the evolved fields, transverse gradient interpolation, conservation law update step, and Lorentz force computation. The best results are achieved with divergence-free evolution of the magnetic field using the constrained transport method, and using little to no explicit artificial viscosity. Codes which fall short in one or more of these areas are still useful, but they must compensate higher numerical dissipation with higher numerical resolution. This paper is the largest, most comprehensive MHD code comparison on an application-like test problem to date. We hope this work will help developers improve their numerical algorithms while helping users to make informed choices in picking optimal applications for their specific astrophysical problems.

We use deep adaptive mesh refinement simulations of isothermal self-gravitating supersonic turbulence to study the imprints of gravity on the mass density distribution in molecular clouds. The simulations show that the density distribution in self-gravitating clouds develops an extended power-law tail at high densities on top of the usual lognormal. We associate the origin of the tail with self-similar collapse solutions and predict the power index values in the range from -7/4 to -3/2 that agree with both simulations and observations of star-forming molecular clouds.

We provide a brief overview of recent advances and outstanding issues in simulations of interstellar turbulence, including isothermal models for interior structure of molecular clouds and larger-scale multiphase models designed to simulate the formation of molecular clouds. We show how self-organization in highly compressible magnetized turbulence in the multiphase ISM can be exploited in simple numerical models to generate realistic initial conditions for star formation.

We use three-dimensional numerical simulations to study self-organization in supersonic turbulence in molecular clouds. Our numerical experiments describe decaying and driven turbulent flows with an isothermal equation of state, sonic Mach numbers from 2 to 10, and various degrees of magnetization. We focus on properties of the velocity field and, specifically, on the level of its potential (dilatational) component as a function of turbulent Mach number, magnetic field strength, and scale. We show how extreme choices of either purely solenoidal or purely potential forcing can reduce the extent of the inertial range in the context of periodic box models for molecular cloud turbulence. We suggest an optimized forcing to maximize the effective Reynolds number in numerical models.

Is the turbulence in cluster-forming regions internally driven by stellar outflows or the consequence of a large-scale turbulent cascade? We address this question by studying the turbulent energy spectrum in NGC 1333. Using synthetic 13CO maps computed with a snapshot of a supersonic turbulence simulation, we show that the VCS method of Lazarian and Pogosyan provides an accurate estimate of the turbulent energy spectrum. We then apply this method to the 13CO map of NGC 1333 from the COMPLETE database. We find the turbulent energy spectrum is a power law, E(k) k^-beta, in the range of scales 0.06 pc < ell < 1.5 pc, with slope beta=1.85\pm 0.04. The estimated energy injection scale of stellar outflows in NGC 1333 is ell_inj 0.3 pc, well resolved by the observations. There is no evidence of the flattening of the energy spectrum above the scale ell_inj predicted by outflow-driven simulations and analytical models. The power spectrum of integrated intensity is also a nearly perfect power law in the range of scales 0.16 pc < ell < 7.9 pc, with no feature above ell_inj. We conclude that the observed turbulence in NGC 1333 does not appear to be driven primarily by stellar outflows.

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.

We have implemented an Adaptive Mesh Refinement criterion explicitly designed to increase spatial resolution around discontinuities in the velocity field in ENZO cosmological simulations. With this technique, shocks and turbulent eddies developed during the hierarchical assembly of galaxy clusters are followed with unprecedented spatial resolution, even at large distances from the clusters center. By measuring the spectral properties of the gas velocity field, its time evolution and the properties of shocks for a reference galaxy cluster, we investigate the connection between accretion processes and the onset of chaotic motions in the simulated Inter Galactic Medium over a wide range of scales

Stable, accurate, divergence-free simulation of magnetized supersonic turbulence is a severe test of numerical MHD schemes and has been surprisingly difficult to achieve due to the range of flow conditions present. Here we present a new, higher order-accurate, low dissipation numerical method which requires no additional dissipation or local "fixes" for stable execution. We describe PPML, a local stencil variant of the popular PPM algorithm for solving the equations of compressible ideal magnetohydrodynamics. The principal difference between PPML and PPM is that cell interface states are evolved rather that reconstructed at every timestep, resulting in a compact stencil. Interface states are evolved using Riemann invariants containing all transverse derivative information. The conservation laws are updated in an unsplit fashion, making the scheme fully multidimensional. Divergence-free evolution of the magnetic field is maintained using the higher order-accurate constrained transport technique of Gardiner and Stone. The accuracy and stability of the scheme is documented against a bank of standard test problems drawn from the literature. The method is applied to numerical simulation of supersonic MHD turbulence, which is important for many problems in astrophysics, including star formation in dark molecular clouds. PPML accurately reproduces in three-dimensions a transition to turbulence in highly compressible isothermal gas in a molecular cloud model. The low dissipation and wide spectral bandwidth of this method make it an ideal candidate for direct turbulence simulations.

We use three-dimensional simulations to study the statistics of supersonic turbulence in molecular clouds. Our numerical experiments describe driven turbulent flows with an isothermal equation of state, Mach numbers around 10, and various degrees of magnetization. We first support the so-called 1/3-rule of Kritsuk et al. 2007 with our new data from a larger 2048^3 simulation. We then attempt to extend the 1/3-rule to supersonic MHD turbulence and get encouraging preliminary results based on a set of 512^3 simulations. Our results suggest an interesting new approach to tackle universal scaling relations and intermittency in supersonic MHD 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.

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.

We use large-scale three-dimensional simulations of supersonic Euler turbulence to study the physics of a highly compressible cascade. Our numerical experiments describe non-magnetized driven turbulent flows with an isothermal equation of state and an rms Mach number of 6. We find that the inertial range velocity scaling deviates strongly from the incompressible Kolmogorov laws. We propose an extension of Kolmogorov's K41 phenomenology that takes into account compressibility by mixing the velocity and density statistics and preserves the K41 scaling of the density-weighted velocity v=rho^{1/3}u. We show that low-order statistics of 'v' are 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 third-order structure function of 'v' scales linearly with separation. We 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.

The Padoan and Nordlund model of the stellar initial mass function (IMF) is derived from low order statistics of supersonic turbulence, neglecting gravity (e.g. gravitational fragmentation, accretion and merging). In this work the predictions of that model are tested using the largest numerical experiments of supersonic hydrodynamic (HD) and magneto-hydrodynamic (MHD) turbulence to date (~1000^3 computational zones) and three different codes (Enzo, Zeus and the Stagger Code). The model predicts a power law distribution for large masses, related to the turbulence energy power spectrum slope, and the shock jump conditions. This power law mass distribution is confirmed by the numerical experiments. The model also predicts a sharp difference between the HD and MHD regimes, which is recovered in the experiments as well, implying that the magnetic field, even below energy equipartition on the large scale, is a crucial component of the process of turbulent fragmentation. These results suggest that the stellar IMF of primordial stars may differ from that in later epochs of star formation, due to differences in both gas temperature and magnetic field strength. In particular, we find that the IMF of primordial stars born in turbulent clouds may be narrowly peaked around a mass of order 10 solar masses, as long as the column density of such clouds is not much in excess of 10^22 cm^-2.

We test a method of estimating the power spectrum of turbulence in molecular clouds based on the comparison of power spectra of integrated intensity maps and single-velocity-channel maps, suggested by Lazarian and Pogosyan. We use synthetic 13CO data from non-LTE radiative transfer calculations based on density and velocity fields of a simulation of supersonic hydrodynamic turbulence. We find that the method yields the correct power spectrum with good accuracy. We then apply the method to the Five College Radio Astronomy Observatory 13CO map of the Perseus region, from the COMPLETE website. We find a power law power spectrum with slope beta=1.81+-0.10. The values of beta as a function of velocity resolution are also confirmed using the lower resolution map of the same region obtained with the AT&T Bell Laboratories antenna. Because of its small uncertainty, this result provides a useful constraint for numerical codes used to simulate molecular cloud turbulence.

We present the results of three-dimensional simulations of supersonic Euler turbulence with grid resolutions up to 1024^3 points. Our numerical experiments describe nonmagnetized driven turbulent flows with an isothermal equation of state and an rms Mach number of 6. We demonstrate that the inertial range scaling properties of turbulence in this strongly compressible regime deviate substantially from a Kolmogorov-like behavior previously recovered for mildly compressible transonic flows.

Maps of estimated dust column density in molecular clouds are usually assumed to reliably trace the total gas column density structure. In this work we present results showing a clear discrepancy between the dust and the gas distribution in the Taurus molecular cloud complex. We compute the power spectrum of a 2MASS extinction map of the Taurus region and find it is much shallower than the power spectrum of a 13CO map of the same region previously analyzed. This discrepancy may be explained as the effect of grain growth on the grain extinction efficiency. However, this would require a wide range of maximum grain sizes, which is ruled out based on constraints from the extinction curve and the available grain models. We show that major effects due to CO formation and depletion are also ruled out. Our result may therefore suggest the existence of intrinsic spatial fluctuations of the dust to gas ratio, with amplitude increasing toward smaller scales. Preliminary results of numerical simulations of trajectories of inertial particles in turbulent flows illustrate how the process of clustering of dust grains by the cloud turbulence may lead to observable effects. However, these results cannot be directly applied to large scale supersonic and magnetized turbulence at present.

We performed a series of three-dimensional numerical simulations of supersonic homogeneous Euler turbulence with adaptive mesh refinement (AMR) and effective grid resolution up to 1024^3 zones. Our experiments describe non-magnetized driven supersonic turbulent flows with an isothermal equation of state. Mesh refinement on shocks and shear is implemented to cover dynamically important structures with the highest resolution subgrids and calibrated to match the turbulence statistics obtained from the equivalent uniform grid simulations. We found that at a level of resolution slightly below 512^3, when a sufficient integral/dissipation scale separation is first achieved, the fraction of the box volume covered by the AMR subgrids first becomes smaller than unity. At the higher AMR levels subgrids start covering smaller and smaller fractions of the whole volume, which scale with the Reynolds number as Re^{-1/4}. We demonstrate the consistency of this scaling with a hypothesis that the most dynamically important structures in intermittent supersonic turbulence are strong shocks with a fractal dimension of two. We show that turbulence statistics derived from AMR simulations and simulations performed on uniform grids agree surprisingly well, even though only a fraction of the volume is covered by AMR subgrids. Based on these results, we discuss the signature of dissipative structures in the statistical properties of supersonic turbulence and their role in overall flow dynamics.

Accretion rates of order 10^-8 M_\odot/yr are observed in young protostars of approximately a solar mass with evidence of circumstellar disks. The accretion rate is significantly lower for protostars of smaller mass, approximately proportional to the second power of the stellar mass, \dot{M}_accr\propto M^2. The traditional view is that the observed accretion is the consequence of the angular momentum transport in isolated protostellar disks, controlled by disk turbulence or self--gravity. However, these processes are not well understood and the observed protostellar accretion, a fundamental aspect of star formation, remains an unsolved problem. In this letter we propose the protostellar accretion rate is controlled by accretion from the large scale gas distribution in the parent cloud, not by the isolated disk evolution. Describing this process as Bondi--Hoyle accretion, we obtain accretion rates comparable to the observed ones. We also reproduce the observed dependence of the accretion rate on the protostellar mass. These results are based on realistic values of the ambient gas density and velocity, as inferred from numerical simulations of star formation in self--gravitating turbulent clouds.

The origin of brown dwarfs (BDs) is an important component of the theory of star formation, because BDs are approximately as numerous as solar mass stars. It has been suggested that BDs originate from the gravitational fragmentation of protostellar disks, a very different mechanism from the formation of hydrogen burning stars. We propose that BDs are instead formed by the process of turbulent fragmentation, like more massive stars. In numerical simulations of turbulence and star formation we find that gravitationally unstable density peaks of BD mass are commonly formed by the turbulent flow. These density peaks collapse into BD mass objects with circumstellar disks, like more massive protostars. We rely on numerical experiments with very large resolution, achieved with adaptive mesh refinement (AMR). The turbulence simulation presented here is the first AMR turbulence experiment ever attempted and achieves an effective resolution of 1024^3 computational zones. The star formation simulation achieves an effective resolution of (10^6)^3 computational zones, from a cloud size of 5pc to protostellar disks resolved down to 1AU.

We employ a generalization of the She & Leveque model to study velocity scaling relations based on our simulations of thermal instability--induced turbulence. Being a by-product of the interstellar phase transition, such multiphase turbulence tends to be more intermittent than compressible isothermal turbulence. Due to radiative cooling, which promotes nonlinear instabilities in supersonic flows, the Hausdorff dimension of the most singular dissipative structures, D, can be as high as 2.3, while in supersonic isothermal turbulence D is limited by shock dissipation to D < 2. We also show that single-phase velocity statistics carry only incomplete information on the turbulent cascade in a multiphase medium. We briefly discuss the possible implications of these results on the hierarchical structure of molecular clouds and on star formation.

In this paper we introduce Enzo, a 3D MPI-parallel Eulerian block-structured adaptive mesh refinement cosmology code. Enzo is designed to simulate cosmological structure formation, but can also be used to simulate a wide range of astrophysical situations. Enzo solves dark matter N-body dynamics using the particle-mesh technique. The Poisson equation is solved using a combination of fast fourier transform (on a periodic root grid) and multigrid techniques (on non-periodic subgrids). Euler's equations of hydrodynamics are solved using a modified version of the piecewise parabolic method. Several additional physics packages are implemented in the code, including several varieties of radiative cooling, a metagalactic ultraviolet background, and prescriptions for star formation and feedback. We also show results illustrating properties of the adaptive mesh portion of the code. Enzo is publicly available and can be downloaded at http://cosmos.ucsd.edu/enzo/ .

We use three-dimensional hydrodynamic numerical simulations to study phase transformations occurring in a clumpy interstellar gas exposed to time-dependent volumetric heating. To mimic conditions in the Galactic interstellar medium, we take a numerical model of a turbulent multiphase medium from Kritsuk & Norman (2002) computed in a periodic box with mean density n_0=0.25 cm^-3 and mean pressure P_0/k~10^3.4 K cm^-3. A second model with n_0=1 cm^-3 is also considered. Variations of the heating rate on a timescale of 1-10 Myr applied thereafter cause pressure variations in the gas and shifting of the thermal equilibrium curve in the phase plane. This stimulates mass transfer between the gas phases via thermal instability, converting 5-10% of the thermal energy into kinetic energy of gas motions. The experiments demonstrate that recurrent substantial heating episodes can maintain turbulence at this level. Possible applications to the interstellar gas heated by variable far-ultraviolet background radiation produced by short-living massive stars are discussed.

We study the dynamics of phase transitions in the interstellar medium by means of three-dimensional hydrodynamic numerical simulations. We use a realistic cooling function and generic nonequilibrium initial conditions to follow the formation history of a multiphase medium in detail in the absence of gravity. We outline a number of qualitatively distinct stages of this process, including a linear isobaric evolution, transition to an isochoric regime, formation of filaments and voids (also known as "thermal" pancakes), the development and decay of supersonic turbulence, an approach to pressure equilibrium, and final relaxation of the multiphase medium. We find that 1%-2% of the initial thermal energy is converted into gas motions in one cooling time. The velocity field then randomizes into turbulence that decays on a dynamical timescale E_k ~ t^-n, 1 < n < 2. While not all initial conditions yield a stable two-phase medium, we examine such a case in detail. We find that the two phases are well mixed with the cold clouds possessing a fine-grained structure near our numerical resolution limit. The amount of gas in the intermediate unstable phase roughly tracks the rms turbulent Mach number, peaking at 25% when M_rms ~ 8, decreasing to 11% when M_rms ~ 0.4.

We use hydrodynamic simulations with adaptive grid refinement to study the dependence of hot gas flows in X-ray luminous giant elliptical galaxies on the efficiency of heat supply to the gas. We consider a number of potential heating mechanisms including Type Ia supernovae and sporadic nuclear activity of a central supermassive black hole. As a starting point for this research we use an equilibrium hydrostatic recycling model (Kritsuk 1996). We show that a compact cooling inflow develops, if the heating is slightly insufficient to counterbalance radiative cooling of the hot gas in the central few kiloparsecs. An excessive heating in the centre, instead, drives a convectively unstable outflow. We model the onset of the instability and a quasi-steady convective regime in the core of the galaxy in two-dimensions assuming axial symmetry. Provided the power of net energy supply in the core is not too high, the convection remains subsonic. The convective pattern is dominated by buoyancy driven large-scale mushroom-like structures. Unlike in the case of a cooling inflow, the X-ray surface brightness of an (on average) isentropic convective core does not display a sharp maximum at the centre. A hybrid model, which combines a subsonic peripheral cooling inflow with an inner convective core, appears to be stable. We also discuss observational implications of these results.

Using numerical techniques we studied the global stability of cooling flows in giant elliptical galaxies. As an initial equilibrium state we choose the hydrostatic gas recycling model (Kritsuk 1996). Non-equilibrium radiative cooling, stellar mass loss, heating by type Ia supernovae, distributed mass deposition, and thermal conductivity are included. Although the recycling model reproduces the basic X-ray observables, it appears to be unstable with respect to the development of inflow or outflow. In spherically symmetry the inflows are subject to a central cooling catastrophe, while the outflows saturate in a form of a subsonic galactic wind. Two-dimensional axisymmetric random velocity perturbations of the equilibrium model trigger the onset of a cooling catastrophe, which develops in an essentially non-spherical way. The simulations show a patchy pattern of mass deposition and the formation of hollow gas jets, which penetrate through the outflow down to the galaxy core. The X-ray observables of such a hybrid gas flow mimic those of the equilibrium recycling model, but the gas temperature exhibits a central depression. The mass deposition rate M_dot consists of two contributions of similar size: (i) a hydrostatic one resembling that of the equilibrium model, and (ii) a dynamical one which is related to the jets and is more concentrated to the centre. For a model galaxy, like NGC 4472, our 2D simulations predict M_dot = 2 M_sun/yr within the cooling radius for the advanced non-linear stage of the instability. We discuss the implications of these results to H_alpha nebulae and star formation in cooling flow galaxies and emphasize the need for high-resolution 3D simulations.

A two-component isothermal equilibrium model is applied to reproduce basic structural properties of dynamically hot stellar systems immersed in their massive dark haloes. The origin of the fundamental plane relation for giant ellipticals is naturally explained as a consequence of dynamical equilibrium in the context of the model. The existence of two galactic families displaying different behaviour in the luminosity--surface-brightness diagram is shown to be a result of a smooth transition from dwarfs, dominated by dark matter near the centre, to giants dominated by the luminous stellar component. The comparison of empirical scaling relations with model predictions suggests that probably a unique dissipative process was operating during the violent stage of development of stellar systems in the dark haloes, and the depth of the potential well controlled the observed luminosity of the resulting galaxies. The interpretation also provides some restrictions on the properties of dark haloes implied by the fundamental scaling laws.

An equilibrium model is presented for a hydrostatic isothermal hot gas distribution in a gravitational well of a giant elliptical galaxy immersed in a massive dark halo. The self-consistently determined gravitational potential of the system provides the optical surface brightness distribution of the galaxy, matching the de Vaucouleurs R^1/4-law, if the sound speed in the gas, stellar velocity dispersion, and velocity dispersion of dark matter are related so that optical and X-ray brightness profiles match each other. The thermal equilibrium of the gas is described in the framework of a one-zone model, which incorporates radiative cooling, stellar mass loss, supernova heating, and mass sink due to local condensational instability. The sequence of saddle-node and saddle-connection bifurcations provides a first-order phase transition, which separates the stable hot phase in the cooling medium for a reasonably high efficiency of condensation and precludes catastrophic cooling of the gas. The bifurcations occur naturally, due to the shape of the radiative cooling function, and this holds for a wide range of gas metallicities. The one-zone model implies a scaling relation for the equilibrium stellar and gas densities, which allows a hydrostatic, thermally stable distribution for the hot isothermal coronal gas. The resulting model for gas recycling reproduces the basic optical and X-ray properties of relatively isolated giant elliptical galaxies with hot haloes. Implications for cD galaxies in X-ray dominant clusters and groups are briefly discussed in the general context of galaxy formation and evolution.