We address the problem of the origin of massive stars, namely the origin, path and timescale of the mass flows that create them. Based on extensive numerical simulations, we propose a scenario where massive stars are assembled by large-scale, converging, inertial flows that naturally occur in supersonic turbulence. We refer to this scenario of massive-star formation as the {\it Inertial-Inflow Model}. This model stems directly from the idea that the mass distribution of stars is primarily the result of {\it turbulent fragmentation}. Under this hypothesis, the statistical properties of the turbulence determine the formation timescale and mass of prestellar cores, posing definite constraints on the formation mechanism of massive stars. We quantify such constraints by the analysis of a simulation of supernova-driven turbulence in a 250-pc region of the interstellar medium, describing the formation of hundreds of massive stars over a time of approximately 30 Myr. Due to the large size of our statistical sample, we can say with full confidence that {\it massive stars in general do not form from the collapse of massive cores, nor from competitive accretion}, as both models are incompatible with the numerical results. We also compute synthetic continuum observables in Herschel and ALMA bands. We find that, depending on the distance of the observed regions, estimates of core mass based on commonly-used methods may exceed the actual core masses by up to two orders of magnitude, and that there is essentially no correlation between estimated and real core masses.

We have investigated the formation and kinematics of sub-mm continuum cores in the Orion A molecular cloud. A comparison between sub-mm continuum and near infrared extinction shows a continuum core detection threshold of $A_V\sim$ 5-10 mag. The threshold is similar to the star formation extinction threshold of $A_V\sim$ 7 mag proposed by recent work, suggesting a universal star formation extinction threshold among clouds within 500 pc to the Sun. A comparison between the Orion A cloud and a massive infrared dark cloud G28.37+0.07 indicates that Orion A produces more dense gas within the extinction range 15 mag $\lesssim A_V \lesssim$ 60 mag. Using data from the CARMA-NRO Orion Survey, we find that dense cores in the integral-shaped filament (ISF) show sub-sonic core-to-envelope velocity dispersion that is significantly less than the local envelope line dispersion, similar to what has been found in nearby clouds. Dynamical analysis indicates that the cores are bound to the ISF. An oscillatory core-to-envelope motion is detected along the ISF. Its origin is to be further explored.

The Planck Catalogue of Galactic Cold Clumps (PGCC) contains over 13000 sources detected based on their cold dust signature. They are believed to consist of a mixture of quiescent, pre-stellar, and already star-forming objects. We extracted PGCC-type objects from cloud simulations and examined their physical and polarisation properties. The comparison with the PGCC catalogue helps to characterise the PGCC sample and provides valuable tests for numerical simulations of interstellar medium. We used several MHD snapshots to define the density field of our models. Sub-millimetre images of the surface brightness and polarisation were obtained with radiative transfer calculations. We examined the statistics of synthetic cold clump catalogues and examined the variations of the clump polarisation fraction p. The clump sizes, aspect ratios, and temperatures in the synthetic catalogue are similar to the PGCC. The fluxes and column densities are smaller by a factor of a few. Rather than with an increased dust opacity, this could be explained by increasing the average column density of the models by a factor of two to three, close to N(H2)= 10^22 cm-2. When the line of sight is parallel to the mean magnetic field, the polarisation fraction tends to increase towards the clump centres, contrary to observations. When the field is perpendicular, the polarisation fraction tends to decrease towards the clumps, but the drop in $p$ is small (e.g. from p~8% to p~7%). Magnetic field geometry reduces the polarisation fraction in the simulated clumps by only \Delta p~1% on average. The larger drop seen towards the actual PGCC clumps suggests some loss of grain alignment in the dense medium, such as predicted by the radiative torque mechanism. The statistical study is not able to quantify dust opacity changes at the scale of the PGCC clumps.

The Vela-Puppis region is known to host the Vela OB2 association as well as several young clusters featuring OB and pre-main sequence stars. Several spatial and kinematic subgroups have been identified in the recent years. By grouping stars based on their positions and velocity we can address the quetion of the dynamical history of the region and the mechanisms that drove stellar formation. The Gaia DR2 astrometry and photometry enables us to characterise the 3D spatial and 3D kinematic distribution of young stars and to estimate the ages of the identified components. We use an unsupervised classification method to group stars based on their proper motions and parallax. We perform a study of the expansion rates of the different identified groups from 3D velocities, and from corrected tangential velocities. We make use of theoretical isochrones to estimate ages. The young stars can be separated into seven main groups of different ages and kinematical distribution. All groups are found to be expanding, although the expansion is mostly not isotropic. The size of the region, the age substructure, and the anistropic expansion rates are compatible with a prolonged period of star formation in a turbulent molecular cloud, and that the dispersion of the stars cannot be explained by gas expulsion alone.

We study density fluctuations in supersonic turbulence using both theoretical methods and numerical simulations. A theoretical formulation is developed for the probability distribution function (PDF) of the density at steady state, connecting it to the conditional statistics of the velocity divergence. Two sets of numerical simulations are carried out, using either a Riemann solver to evolve the Euler equations or a finite-difference method to evolve the Navier-Stokes (N-S) equations. After confirming the validity of our theoretical formulation with the N-S simulations, we examine the effects of dynamical processes on the PDF, showing that the nonlinear term in the divergence equation amplifies the right tail of the PDF and reduces the left one, the pressure term reduces both the right and left tails, and the viscosity term, counter-intuitively, broadens the right tail of the PDF. Despite the inaccuracy of the velocity divergence from the Riemann runs, as found in our previous work, we show that the density PDF from the Riemann runs is consistent with that from the N-S runs. Taking advantage of their much higher effective resolution, we then use the Riemann runs to study the dependence of the PDF on the Mach number, $\mathcal{M}$, up to $\mathcal{M}\sim30$. The PDF width, $\sigma_{s}$, follows the relation $\sigma_{s}^2 = \ln (1+b^2 {\mathcal M}^2)$, with $b\approx0.38$. However, the PDF exhibits a negative skewness that increases with increasing $\mathcal{M}$, so much of the growth of $\sigma_{s}$ is accounted for by the growth of the left PDF tail, while the growth of the right tail tends to saturate. Thus, the usual prescription that combines a lognormal shape with the standard variance-Mach number relation greatly overestimates the right PDF tail at large $\mathcal{M}$, which may have a significant impact on theoretical models of star formation.

We examine the accuracy of spatial derivatives computed from numerical simulations of supersonic turbulence. Two sets of simulations, carried out using a finite-volume code that evolves the hydrodynamic equations with an approximate Riemann solver and a finite-difference code that solves the Navier-Stokes equations, are tested against a number of criteria based on the continuity equation, including exact results at statistically steady state. We find that the spatial derivatives in the Navier-Stokes runs are accurate and satisfy all the criteria. In particular, they satisfy our exact results that the conditional mean velocity divergence, $\langle \nabla \cdot {\bs u}|s\rangle$, where $s$ is the logarithm of density, and the conditional mean of the advection of $s$, $\langle {\bs u} \cdot \nabla s|s\rangle$, vanish at steady state for all density values, $s$. On the other hand, the Riemann solver simulations fail all the tests that require accurate evaluation of spatial derivatives, resulting in apparent violation of the continuity equation, even if the solver enforces mass conservation. In particular, analysis of the Riemann simulations may lead to the incorrect conclusion that the $\pdv$ work tends to preferentially convert kinetic energy into thermal energy, inconsistent with the exact result that the energy exchange by $\pdv$ work is symmetric in barotropic supersonic turbulence at steady state. The inaccuracy of spatial derivatives is a general problem in the post-processing of simulations of supersonic turbulence with Riemann solvers. Solutions from such simulations must be used with caution in post-processing studies concerning the spatial gradients.

We present an initial overview of the filamentary structure in the Orion A molecular cloud utilizing a high angular and velocity resolution C$^{18}$O(1-0) emission map that was recently produced as part of the CARMA-NRO Orion Survey. The main goal of this study is to build a credible method to study varying widths of filaments which has previously been linked to star formation in molecular clouds. Due to the diverse star forming activities taking place throughout its $\sim$20 pc length, together with its proximity of 388 pc, the Orion A molecular cloud provides an excellent laboratory for such an experiment to be carried out with high resolution and high sensitivity. Using the widely-known structure identification algorithm, DisPerSE, on a 3-dimensional (PPV) C$^{18}$O cube, we identified 625 relatively short (the longest being 1.74 pc) filaments over the entire cloud. We study the distribution of filament widths using FilChaP, a python package that we have developed and made publicly available. We find that the filaments identified in a 2 square degree PPV cube do not overlap spatially, except for the complex OMC-4 region that shows distinct velocity components along the line of sight. The filament widths vary between 0.02 and 0.3 pc depending on the amount of substructure that a filament possesses. The more substructure a filament has, the larger is its width. We also find that despite this variation, the filament width shows no anticorrelation with the central column density which is in agreement with previous Herschel observations.

Context. Polarized dust emission can be used as a probe for magnetic fields at different scales. In particular, the analysis of the relative orientation between the density structures and the magnetic field projected on the plane of the sky can provide information on the role of magnetic fields in shaping the structure of molecular clouds. Aims. The relative orientation has been investigated using different methods. Our goal is to explicitly compare two of them: the Rolling Hough Transform (RHT) and the gradient technique. Methods. We have applied the RHT method and the gradient technique to two regions, identified in surface brightness maps at 353 GHz generated via MHD simulations post-processed using radiative transfer modelling. For both methods we have derived the relative orientation between the magnetic field and the density structures, to which we have applied two different statistics, the histogram of relative orientation (HRO) statistic and the projected Rayleigh statistic (PRS). Results. The samples of pixels selected by each method are substantially different. When the methods are applied to the same pixel selection, the results in terms of relative orientation as a function of column density are consistent between each other. When each method is applied to its own pixel selection, the differences are apparent. In Region 1, the RHT method roughly reproduces the observed trend of the relative orientation becoming more perpendicular for increasing column density, while the gradient method, applied at the same resolution as RHT, gives the opposite trend, with the relative orientation moving towards a more parallel alignment. In Region 2, the situation is reversed. The inconsistent results, which are due to the different pixel selections operated by the methods and to the intrinsic differences between them, provide complementary valuable information. Abridged.

The probabilistic approach to turbulence is applied to investigate density fluctuations in supersonic turbulence. We derive kinetic equations for the probability distribution function (PDF) of the logarithm of the density field, $s$, in compressible turbulence in two forms: a first-order partial differential equation involving the average divergence conditioned on the flow density, $\langle \nabla \cdot {\bs u} | s\rangle$, and a Fokker-Planck equation with the drift and diffusion coefficients equal to $-\langle {\bs u} \cdot \nabla s | s\rangle$ and $\langle {\bs u} \cdot \nabla s | s\rangle$, respectively. Assuming statistical homogeneity only, the detailed balance at steady state leads to two exact results, $\langle \nabla \cdot {\bs u} | s \rangle =0$, and $\langle {\bs u} \cdot \nabla s | s\rangle=0$. The former indicates a balance of the flow divergence over all expanding and contracting regions at each given density. The exact results provide an objective criterion to judge the accuracy of numerical codes with respect to the density statistics in supersonic turbulence. We also present a method to estimate the effective numerical diffusion as a function of the flow density and discuss its effects on the shape of the density PDF.

We present the first results from a new, high resolution, $^{12}$CO(1-0), $^{13}$CO(1-0), and C$^{18}$O(1-0) molecular line survey of the Orion A cloud, hereafter referred to as the CARMA-NRO Orion Survey. CARMA observations have been combined with single-dish data from the Nobeyama 45m telescope to provide extended images at about 0.01 pc resolution, with a dynamic range of approximately 1200 in spatial scale. Here we describe the practical details of the data combination in uv space, including flux scale matching, the conversion of single dish data to visibilities, and joint deconvolution of single dish and interferometric data. A $\Delta$-variance analysis indicates that no artifacts are caused by combining data from the two instruments. Initial analysis of the data cubes, including moment maps, average spectra, channel maps, position-velocity diagrams, excitation temperature, column density, and line ratio maps provides evidence of complex and interesting structures such as filaments, bipolar outflows, shells, bubbles, and photo-eroded pillars. The implications for star formation processes are profound and follow-up scientific studies by the CARMA-NRO Orion team are now underway. We plan to make all the data products described here generally accessible; some are already available at https://dataverse.harvard.edu/dataverse/CARMA-NRO-Orion

We address the turbulent fragmentation scenario for the origin of the stellar initial mass function (IMF), using a large set of numerical simulations of randomly driven supersonic MHD turbulence. The turbulent fragmentation model successfully predicts the main features of the observed stellar IMF assuming an isothermal equation of state without any stellar feedback. As a test of the model, we focus on the case of a magnetized isothermal gas, neglecting stellar feedback, while pursuing a large dynamic range in both space and timescales covering the full spectrum of stellar masses from brown dwarfs to massive stars. Our simulations represent a generic 4 pc region within a typical Galactic molecular cloud, with a mass of 3000 Msun and an rms velocity 10 times the isothermal sound speed and 5 times the average Alfven velocity, in agreement with observations. We achieve a maximum resolution of 50 au and a maximum duration of star formation of 4.0 Myr, forming up to a thousand sink particles whose mass distribution closely matches the observed stellar IMF. A large set of medium-size simulations is used to test the sink particle algorithm, while larger simulations are used to test the numerical convergence of the IMF and the dependence of the IMF turnover on physical parameters predicted by the turbulent fragmentation model. We find a clear trend toward numerical convergence and strong support for the model predictions, including the initial time evolution of the IMF. We conclude that the physics of isothermal MHD turbulence is sufficient to explain the origin of the IMF.

The magnetic field of molecular clouds (MCs) plays an important role in the process of star formation: it determins the statistical properties of supersonic turbulence that controls the fragmentation of MCs, controls the angular momentum transport during the protostellar collapse, and affects the stability of circumstellar disks. In this work, we focus on the problem of the determination of the magnetic field strength. We review the idea that the MC turbulence is super-Alfv\'{e}nic, and we argue that MCs are bound to be born super-Alfv\'{e}nic. We show that this scenario is supported by results from a recent simulation of supernova-driven turbulence on a scale of 250 pc, where the turbulent cascade is resolved on a wide range of scales, including the interior of MCs.

We compute the star formation rate (SFR) in molecular clouds (MCs) that originate {\it ab initio} in a new, higher-resolution simulation of supernova-driven turbulence. Because of the large number of well-resolved clouds with self-consistent boundary and initial conditions, we obtain a large range of cloud physical parameters with realistic statistical distributions, an unprecedented sample of star-forming regions to test SFR models and to interpret observational surveys. We confirm the dependence of the SFR per free-fall time, $SFR_{\rm ff}$, on the virial parameter, $\alpha_{\rm vir}$, found in previous simulations, and compare a revised version of our turbulent fragmentation model with the numerical results. The dependences on Mach number, ${\cal M}$, gas to magnetic pressure ratio, $\beta$, and compressive to solenoidal power ratio, $\chi$ at fixed $\alpha_{\rm vir}$ are not well constrained, because of random scatter due to time and cloud-to-cloud variations in $SFR_{\rm ff}$. We find that $SFR_{\rm ff}$ in MCs can take any value in the range $0 \le SFR_{\rm ff} \lesssim 0.2$, and its probability distribution peaks at a value $SFR_{\rm ff}\approx 0.025$, consistent with observations. The values of $SFR_{\rm ff}$ and the scatter in the $SFR_{\rm ff}$--$\alpha_{\rm vir}$ relation are consistent with recent measurements in nearby MCs and in clouds near the Galactic center. Although not explicitly modeled by the theory, the scatter is consistent with the physical assumptions of our revised model and may also result in part from a lack of statistical equilibrium of the turbulence, due to the transient nature of MCs.

The compressibility of molecular cloud (MC) turbulence plays a crucial role in star formation models, because it controls the amplitude and distribution of density fluctuations. The relation between the compressive ratio (the ratio of powers in compressive and solenoidal motions) and the statistics of turbulence has been previously studied systematically only in idealized simulations with random external forces. In this work, we analyze a simulation of large-scale turbulence (250 pc) driven by supernova (SN) explosions that has been shown to yield realistic MC properties. We demonstrate that SN driving results in MC turbulence with a broad lognormal distribution of the compressive ratio, with a mean value $\approx 0.3$, lower than the equilibrium value of $\approx 0.5$ found in the inertial range of isothermal simulations with random solenoidal driving. We also find that the compressibility of the turbulence is not noticeably affected by gravity, nor are the mean cloud radial (expansion or contraction) and solid-body rotation velocities. Furthermore, the clouds follow a general relation between the rms density and the rms Mach number similar to that of supersonic isothermal turbulence, though with a large scatter, and their average gas density PDF is described well by a lognormal distribution, with the addition of a high-density power-law tail when self-gravity is included.

We present a comparison of molecular clouds (MCs) from a simulation of supernova-driven interstellar medium (ISM) turbulence with real MCs from the Outer Galaxy Survey. The radiative transfer calculations to compute synthetic CO spectra are carried out assuming the CO relative abundance depends only on gas density, according to four different models. Synthetic MCs are selected above a threshold brightness temperature value, $T_{\rm B,min}=1.4$ K, of the $J=1-0$ $^{12}$CO line, generating 16 synthetic catalogs (four different spatial resolutions and four CO abundance models), each containing up to several thousands MCs. The comparison with the observations focuses on the mass and size distributions and on the velocity-size and mass-size Larson relations. The mass and size distributions are found to be consistent with the observations, with no significant variations with spatial resolution or chemical model, except in the case of the unrealistic model with constant CO abundance. The velocity-size relation is slightly too steep for some of the models, while the mass-size relation is a bit too shallow for all models only at a spatial resolution $dx\approx 1$ pc. The normalizations of the Larson relations show a clear dependence on spatial resolution, for both the synthetic and the real MCs. The comparison of the velocity-size normalization suggests that the SN rate in the Perseus arm is approximately 70\% or less of the rate adopted in the simulation. Overall, the realistic properties of the synthetic clouds confirm that supernova-driven turbulence can explain the origin and dynamics of MCs.

Turbulence is ubiquitous in molecular clouds (MCs), but its origin is still unclear because MCs are usually assumed to live longer than the turbulence dissipation time. Interstellar medium (ISM) turbulence is likely driven by SN explosions, but it has never been demonstrated that SN explosions can establish and maintain a turbulent cascade inside MCs consistent with the observations. In this work, we carry out a simulation of SN-driven turbulence in a volume of (250 pc)$^3$, specifically designed to test if SN driving alone can be responsible for the observed turbulence inside MCs. We find that SN driving establishes a velocity scaling consistent with the usual scaling laws of supersonic turbulence, suggesting that previous idealized simulations of MC turbulence, driven with a random, large-scale volume force, were correctly adopted as appropriate models for MC turbulence, despite the artificial driving. We also find that the same scaling laws extend to the interior of MCs, and that the velocity-size relation of the MCs selected from our simulation is consistent with that of MCs from the Outer-Galaxy Survey, the largest MC sample available. The mass-size relation and the mass and size probability distributions also compare successfully with those of the Outer Galaxy Survey. Finally, we show that MC turbulence is super-Alfv\'{e}nic with respect to both the mean and rms magnetic-field strength. We conclude that MC structure and dynamics are the natural result of SN-driven turbulence.

Context: Understanding how protostars accrete their mass is a central question of star formation. One aspect of this is trying to understand whether the time evolution of accretion rates in deeply embedded objects is best characterised by a smooth decline from early to late stages or by intermittent bursts of high accretion. Aims: We create synthetic observations of deeply embedded protostars in a large numerical simulation of a molecular cloud, which are compared directly to real observations. The goal is to compare episodic accretion events in the simulation to observations and to test the methodology used for analysing the observations. Methods: Simple freeze-out and sublimation chemistry is added to the simulation, and synthetic C$^{18}$O line cubes are created for a large number of simulated protostars. The spatial extent of C$^{18}$O is measured for the simulated protostars and compared directly to a sample of 16 deeply embedded protostars observed with the Submillimeter Array. If CO is distributed over a larger area than predicted based on the protostellar luminosity, it may indicate that the luminosity has been higher in the past and that CO is still in the process of refreezing. Results: Approximately 1% of the protostars in the simulation show extended C$^{18}$O emission, as opposed to approximately 50% in the observations, indicating that the magnitude and frequency of episodic accretion events in the simulation is too low relative to observations. The protostellar accretion rates in the simulation are primarily modulated by infall from the larger scales of the molecular cloud, and do not include any disk physics. The discrepancy between simulation and observations is taken as support for the necessity of disks, even in deeply embedded objects, to produce episodic accretion events of sufficient frequency and amplitude.

Old Globular Clusters (GCs) in the Milky Way have ages of about 13 Gyr, placing their formation time in the reionization epoch. We propose a novel scenario for the formation of these systems based on the merger of two or more atomic cooling halos at high-redshift (z>6). First generation stars are formed as an intense burst in the center of a minihalo that grows above the threshold for hydrogen cooling (halo mass M_h~10^8 Msun) by undergoing a major merger within its cooling timescale (~150 Myr). Subsequent minor mergers and sustained gas infall bring new supply of pristine gas at the halo center, creating conditions that can trigger new episodes of star formation. The dark-matter halo around the GC is then stripped during assembly of the host galaxy halo. Minihalo merging is efficient only in a short redshift window, set by the LCDM parameters, allowing us to make a strong prediction on the age distribution for old GCs. From cosmological simulations we derive an average merging redshift <z>=9 and narrow distribution Dz=2, implying average GC age <t_age>=13.0+/-0.2 Gyr including ~0.2 Gyr of star formation delay. Qualitatively, our scenario reproduces other general old GC properties (characteristic masses and number of objects, metallicity versus galactocentric radius anticorrelation, radial distribution), but unlike age, these generally depend on details of baryonic physics. In addition to improved age measurements, direct validation of the model at z~10 may be within reach of ultradeep gravitationally lensed observations with the James Webb Space Telescope.

Motivated by its importance for modeling dust particle growth in protoplanetary disks, we study turbulence-induced collision statistics of inertial particles as a function of the particle friction time, tau_p. We show that turbulent clustering significantly enhances the collision rate for particles of similar sizes with tau_p corresponding to the inertial range of the flow. If the friction time, tau_p,h, of the larger particle is in the inertial range, the collision kernel per unit cross section increases with increasing friction time, tau_p,l, of the smaller particle, and reaches the maximum at tau_p,l = tau_p,h, where the clustering effect peaks. This feature is not captured by the commonly-used kernel formula, which neglects the effect of clustering. We argue that turbulent clustering helps alleviate the bouncing barrier problem for planetesimal formation. We also investigate the collision velocity statistics using a collision-rate weighting factor to account for higher collision frequency for particle pairs with larger relative velocity. For tau_p,h in the inertial range, the rms relative velocity with collision-rate weighting is found to be invariant with tau_p,l and scales with tau_p,h roughly as ~ tau_p,h^(1/2). The weighting factor favors collisions with larger relative velocity, and including it leads to more destructive and less sticking collisions. We compare two collision kernel formulations based on spherical and cylindrical geometries. The two formulations give consistent results for the collision rate and the collision-rate weighted statistics, except that the spherical formulation predicts more head-on collisions than the cylindrical formulation.

We investigate the role of mass infall in the formation and evolution of protostars. To avoid ad hoc initial and boundary conditions, we consider the infall resulting self-consistently from modeling the formation of stellar clusters in turbulent molecular clouds. We show that infall rates in turbulent clouds are comparable to accretion rates inferred from protostellar luminosities or measured in pre-main-sequence stars. They should not be neglected in modeling the luminosity of protostars and the evolution of disks, even after the embedded protostellar phase. We find large variations of infall rates from protostar to protostar, and large fluctuations during the evolution of individuals protostars. In most cases, the infall rate is initially of order 10$^{-5}$\msun\ yr$^{-1}$, and may either decay rapidly in the formation of low-mass stars, or remain relatively large when more massive stars are formed. The simulation reproduces well the observed characteristic values and scatter of protostellar luminosities and matches the observed protostellar luminosity function. The luminosity problem is therefore solved once realistic protostellar infall histories are accounted for, with no need for extreme accretion episodes. These results are based on a simulation of randomly-driven magneto-hydrodynamic turbulence on a scale of 4pc, including self-gravity, adaptive-mesh refinement to a resolution of 50AU, and accreting sink particles. The simulation yields a low star formation rate, consistent with the observations, and a mass distribution of sink particles consistent with the observed stellar initial mass function during the whole duration of the simulation, forming nearly 1,300 sink particles over 3.2 Myr.

Motivated by its important role in the collisional growth of dust particles in protoplanetary disks, we investigate the probability distribution function (PDF) of the relative velocity of inertial particles suspended in turbulent flows. Using the simulation from our previous work, we compute the relative velocity PDF as a function of the friction timescales, tau_p1 and tau_p2, of two particles of arbitrary sizes. The friction time of particles included in the simulation ranges from 0.1 tau_eta to 54T_L, with tau_eta and T_L the Kolmogorov time and the Lagrangian correlation time of the flow, respectively. The relative velocity PDF is generically non-Gaussian, exhibiting fat tails. For a fixed value of tau_p1, the PDF is the fattest for equal-size particles (tau_p2~tau_p1), and becomes thinner at both tau_p2<tau_ p1 and tau_p2>tau_p1. Defining f as the friction time ratio of the smaller particle to the larger one, we find that, at a given f in 1/2<f<1, the PDF fatness first increases with the friction time, tau_ph, of the larger particle, peaks at tau_ph~tau_eta, and then decreases as tau_ph increases further. For 0<f<1/4, the PDF shape becomes continuously thinner with increasing tau_ph. The PDF is nearly Gaussian only if tau_ph is sufficiently large (>>T_L). These features are successfully explained by the Pan & Padoan model. Using our simulation data and some simplifying assumptions, we estimated the fractions of collisions resulting in sticking, bouncing, and fragmentation as a function of the dust size in protoplanerary disks, and argued that accounting for non-Gaussianity of the collision velocity may help further alleviate the bouncing barrier problem.

We extend our earlier work on turbulence-induced relative velocity between equal-size particles (Pan and Padoan, Paper I) to particles of arbitrarily different sizes. The Pan and Padoan (PP10) model shows that the relative velocity between different particles has two contributions, named the generalized shear and acceleration terms, respectively. The generalized shear term represents the particles' memory of the spatial flow velocity difference across the particle distance in the past, while the acceleration term is associated with the temporal flow velocity difference on individual particle trajectories. Using the simulation of Paper I, we compute the root-mean-square relative velocity, <w^2>^1/2, as a function of the friction times, tau_p1 and tau_p2, of the two particles, and show that the PP10 prediction is in satisfactory agreement with the data, confirming its physical picture. For a given tau_p1 below the Lagrangian correlation time of the flow, T_L, <w^2>^1/2 as a function of tau_p2 shows a dip at tau_p2~tau_p1, indicating tighter velocity correlation between similar particles. Defining a ratio f=tau_pl/tau_ph, with tau_pl and tau_ph the friction times of the smaller and larger particles, we find that <w^2>^1/2 increases with decreasing f due to the generalized acceleration contribution, which dominates at f<1/4. At a fixed f, our model predicts that <w^2>^1/2 scales as tau_ph^1/2 for tau_p,h in the inertial range of the flow, stays roughly constant for T_L <tau_ph < T_L/f, and finally decreases as tau_ph^-1/2 for tau_ph>>T_L/f. The acceleration term is independent of the particle distance, r, and thus reduces the r-dependence of <w^2>^1/2 in the bidisperse case.

We review recent advances in the analytical and numerical modeling of the star formation rate in molecular clouds and discuss the available observational constraints. We focus on molecular clouds as the fundamental star formation sites, rather than on the larger-scale processes that form the clouds and set their properties. Molecular clouds are shaped into a complex filamentary structure by supersonic turbulence, with only a small fraction of the cloud mass channeled into collapsing protostars over a free-fall time of the system. In recent years, the physics of supersonic turbulence has been widely explored with computer simulations, leading to statistical models of this fragmentation process, and to the prediction of the star formation rate as a function of fundamental physical parameters of molecular clouds, such as the virial parameter, the rms Mach number, the compressive fraction of the turbulence driver, and the ratio of gas to magnetic pressure. Infrared space telescopes, as well as ground-based observatories have provided unprecedented probes of the filamentary structure of molecular clouds and the location of forming stars within them.

We use the adaptive mesh refinement code RAMSES to model the formation of protoplanetary disks in realistic star formation environments. The resolution scales over up to 29 powers of two ($\sim$ 9 orders of magnitude) covering a range from outer scales of 40 pc to inner scales of 0.015 AU. The accretion rate from a 1.5 solar mass envelope peaks near $10^{-4}$ $\mspy$ about 6 kyr after sink particle formation and then decays approximately exponentially, reaching $10^{-6}$ $\mspy$ in 100 kyr. The models suggest universal scalings of physical properties with radius during the main accretion phase, with kinetic and / or magnetic energy in approximate balance with gravitational energy. Efficient accretion is made possible by the braking action of the magnetic field, which nevertheless allows a near-Keplerian disk to grow to a 100 AU size. The magnetic field strength ranges from more than 10 G at 0.1 AU to less than 1 mG at 100 AU, and drives a time dependent bipolar outflow, with a collimated jet and a broader disk wind.

We study the relative velocity of inertial particles suspended in turbulent flows and discuss implications for dust particle collisions in protoplanetary disks. We simulate a weakly compressible turbulent flow, evolving 14 particle species with friction timescale, tau_p, covering the entire range of scales in the flow. The particle Stokes numbers, St, measuring the ratio of tau_p to the Kolmogorov timescale, are in the range from ~0.1 to ~800. Using simulation results, we show that the model by Pan & Padoan (PP10) gives satisfactory predictions for the rms relative velocity between identical particles. The probability distribution function (PDF) of the relative velocity is found to be highly non-Gaussian. The PDF tails are well described by a 4/3 stretched exponential function for particles with tau_p ~ 1-2 T_L, where T_L is the Lagrangian correlation timescale, consistent with a prediction based on PP10. The PDF approaches Gaussian only for very large particles with tau_p >~ 54 T_L. We split particle pairs at given distances into two types with low and high relative speeds, referred to as continuous and caustic types, respectively, and compute their contributions to the collision kernel. Although amplified by the effect of clustering, the continuous contribution vanishes in the limit of infinitesimal particle distance, where the caustic contribution dominates. The caustic kernel per unit cross section rises rapidly as St increases toward ~1, reaches a maximum at tau_p ~ 2 T_L, and decreases as tau_p^{-1/2} for tau_p >> T_L.

We show that supersonic MHD turbulence yields a star formation rate (SFR) as low as observed in molecular clouds (MCs), for characteristic values of the free-fall time divided by the dynamical time, $t_{\rm ff}/t_{\rm dyn}$, the alfv\'{e}nic Mach number, ${\cal M}_{\rm a}$, and the sonic Mach number, ${\cal M}_{\rm s}$. Using a very large set of deep adaptive-mesh-refinement simulations, we quantify the dependence of the SFR per free-fall time, $\epsilon_{\rm ff}$, on the above parameters. Our main results are: i) $\epsilon_{\rm ff}$ decreases exponentially with increasing $t_{\rm ff}/t_{\rm dyn}$, but is insensitive to changes in ${\cal M}_{\rm s}$, for constant values of $t_{\rm ff}/t_{\rm dyn}$ and ${\cal M}_{\rm a}$. ii) Decreasing values of ${\cal M}_{\rm a}$ (stronger magnetic fields) reduce $\epsilon_{\rm ff}$, but only to a point, beyond which $\epsilon_{\rm ff}$ increases with a further decrease of ${\cal M}_{\rm a}$. iii) For values of ${\cal M}_{\rm a}$ characteristic of star-forming regions, $\epsilon_{\rm ff}$ varies with ${\cal M}_{\rm a}$ by less than a factor of two. We propose a simple star-formation law, based on the empirical fit to the minimum $\epsilon_{\rm ff}$, and depending only on $t_{\rm ff}/t_{\rm dyn}$: $\epsilon_{\rm ff} \approx \epsilon_{\rm wind} \exp(-1.6 \,t_{\rm ff}/t_{\rm dyn})$. Because it only depends on the mean gas density and rms velocity, this law is straightforward to implement in simulations and analytical models of galaxy formation and evolution.

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.

The observed similarities between the mass function of prestellar cores (CMF) and the stellar initial mass function (IMF) have led to the suggestion that the IMF is already largely determined in the gas phase. However, theoretical arguments show that the CMF may differ significantly from the IMF. In this Letter, we study the relation between the CMF and the IMF, as predicted by the IMF model of Padoan and Nordlund. We show that 1) the observed mass of prestellar cores is on average a few times smaller than that of the stellar systems they generate; 2) the CMF rises monotonically with decreasing mass, with a noticeable change in slope at approximately 3-5 solar masses, depending on mean density; 3) the selection of cores with masses larger than half their Bonnor-Ebert mass yields a CMF approximately consistent with the system IMF, rescaled in mass by the same factor as our model IMF, and therefore suitable to estimate the local efficiency of star formation, and to study the dependence of the IMF peak on cloud properties; 4) only one in five pre-brown-dwarf core candidates is a true progenitor to a brown dwarf.

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 study the reliability of mass estimates obtained for molecular cloud cores using sub-millimetre and infrared dust emission. We use magnetohydrodynamic simulations and radiative transfer to produce synthetic observations with spatial resolution and noise levels typical of Herschel surveys. We estimate dust colour temperatures using different pairs of intensities, calculate column densities and compare the estimated masses with the true values. We compare these results to the case when all five Herschel wavelengths are available. We investigate the effects of spatial variations of dust properties and the influence of embedded heating sources. Wrong assumptions of dust opacity and its spectral index beta can cause significant systematic errors in mass estimates. These are mainly multiplicative and leave the slope of the mass spectrum intact, unless cores with very high optical depth are included. Temperature variations bias colour temperature estimates and, in quiescent cores with optical depths higher than for normal stable cores, masses can be underestimated by up to one order of magnitude. When heated by internal radiation sources the observations recover the true mass spectra. The shape, although not the position, of the mass spectrum is reliable against observational errors and biases introduced in the analysis. This changes only if the cores have optical depths much higher than expected for basic hydrostatic equilibrium conditions. Observations underestimate the value of beta whenever there are temperature variations along the line of sight. A bias can also be observed when the true beta varies with wavelength. Internal heating sources produce an inverse correlation between colour temperature and beta that may be difficult to separate from any intrinsic beta(T) relation of the dust grains. This suggests caution when interpreting the observed mass spectra and the spectral indices.

We report results from a cosmological simulation with non-equilibrium chemistry of 21 species, including H2, HD, and LiH molecular cooling. Starting from cosmological initial conditions, we focus on the evolution of the central 1.8 Kpc region of a 3 x 10^7 Msun halo. The crossing of a few 10^6 Msun halos and the gas accretion through larger scale filaments generate a turbulent environment within this region. Due to the short cooling time caused by the non-equilibrium formation of H2, the supersonic turbulence results in a very fragmented mass distribution, where dense, gravitationally unstable clumps emerge from a complex network of dense filaments. At z=10.87, we find approximately 25 well defined, gravitationally unstable clumps, with masses of 4 x 10^3-9 x 10^5 Msun, temperatures of approximately 300K, and cooling times much shorter than the free-fall time. Only the initial phase of the collapse of individual clumps is spatially resolved in the simulation. Depending on the density reached in the collapse, the estimated average Bonnor-Ebert masses are in the range 200-800 Msun. We speculate that each clump may further fragment into a cluster of stars with a characteristic mass in the neighborhood of 50 Msun. This process at z ~ 11 may represent the dominant mode of Pop. III star formation, causing a rapid chemical enrichment of the protogalactic environment.

This work presents a new physical model of the star formation rate (SFR), verified with an unprecedented set of large numerical simulations of driven, supersonic, self-gravitating, magneto-hydrodynamic (MHD) turbulence, where collapsing cores are captured with accreting sink particles. The model depends on the relative importance of gravitational, turbulent, magnetic, and thermal energies, expressed through the virial parameter, alpha_vir, the rms sonic Mach number, M_S,0, and the ratio of mean gas pressure to mean magnetic pressure, beta_0. The SFR is predicted to decrease with increasing alpha_vir (stronger turbulence relative to gravity), to increase with increasing M_S,0 (for constant values of alpha_vir), and to depend weakly on beta_0 for values typical of star forming regions (M_S,0 ~ 4-20 and beta_0 ~ 1-20). In the unrealistic limit of beta_0 -> infinity, that is in the complete absence of a magnetic field, the SFR increases approximately by a factor of three, which shows the importance of magnetic fields in the star formation process, even when they are relatively weak (super-Alfvenic turbulence). In this non-magnetized limit, our definition of the critical density for star formation has the same dependence on alpha_vir, and almost the same dependence on M_S,0, as in the model of Krumholz and McKee, although our physical derivation does not rely on the concepts of local turbulent pressure and sonic scale. However, our model predicts a different dependence of the SFR on alpha_vir and M_S,0 than the model of Krumholz and McKee. The star-formation simulations used to test the model result in an approximately constant SFR, after an initial transient phase. Both the value of the SFR and its dependence on the virial parameter found in the simulations are shown to agree very well with the theoretical predictions.

In this work, we present the mass and magnetic distributions found in a recent Adaptive Mesh Refinement (AMR) MHD simulation of supersonic, \sa, self gravitating turbulence. Powerlaw tails are found in both volume density and magnetic field probability density functions, with $P(\rho) \propto \rho^{-1.67}$ and $P(B)\propto B^{-2.74}$. A power law is also found between magnetic field strength and density, with $B\propto \rho^{0.48}$, throughout the collapsing gas. The mass distribution of gravitationally bound cores is shown to be in excellent agreement with recent observation of prestellar cores. The mass to flux distribution of cores is also found to be in excellent agreement with recent Zeeman splitting measurements.

We present a model for the relative velocity of inertial particles in turbulent flows. Our general formulation shows that the relative velocity has contributions from two terms, referred to as the generalized acceleration and generalized shear terms, because they reduce to the well known acceleration and shear terms in the Saffman-Turner limit. The generalized shear term represents particles' memory of the flow velocity difference along their trajectories and depends on the inertial particle pair dispersion backward in time. The importance of this backward dispersion in determining the particle relative velocity is emphasized. We find that our model with a two-phase separation behavior, an early ballistic phase and a later tracer-like phase, as found by recent simulations for the forward (in time) dispersion of inertial particle pairs, gives good fits to the measured relative speeds from simulations at low Reynolds numbers. In the monodisperse case with identical particles, the generalized acceleration term vanishes and the relative velocity is determined by the generalized shear term. At large Reynolds numbers, our model gives a $St^{1/2}$ dependence of the relative velocity on the Stokes number $St$ in the inertial range for both the ballistic behavior and the Richardson separation law. This leads to the same inertial-range scaling for the two-phase separation that well fits the simulation results. Our calculations for the bidisperse case show that, with the friction timescale of one particle fixed, the relative speed as a function of the other particle's friction time has a dip when the two timescales are similar. We find that the primary contribution at the dip is from the generalized shear term, while the generalized acceleration term is dominant for particles of very different sizes.

(abridged) We develop a model which describes the coevolution of the mass function of dense cores and of the IMF in a protocluster clump. In the model, cores injected in the clump evolve under the effect of gas accretion. Accretion onto the cores follows a time-dependent accretion rate that describes accretion in a turbulent medium. Once the accretion timescales of cores exceed their contraction timescales, they are turned into stars. We include the effect of feedback by the newly formed massive stars through their stellar winds. A fraction of the wind's energy is assumed to counter gravity and disperse the gas from the protocluster and as a consequence, quench further star formation. The latter effect sets the final IMF of the cluster. We apply our model to a clump that is expected to resemble the progenitor clump of the Orion Nebula Cluster (ONC). Our model is able to reproduce both the shape and normalization of the ONC's IMF and the mass function of dense cores in Orion. The complex features of the ONC's IMF,i.e., a shallow slope in the mass range ~0.3-2.5 Msol,a steeper slope in the mass range ~2.5-12 Msol, and a nearly flat tail at the high mass end are reproduced. The model predicts a 'rapid' star formation process with an age spread for the stars of 2.3 10^5 yr which is consistent with the fact that 80% of the ONC's stars have ages of <=0.3 Myr. The model predicts a primordial mass segregation with the most massive stars being born in the region between 2-4 times the core radius of the cluster. In parallel, the model also reproduces, simultaneously, the mass function of dense cores in Orion. We study the effects of varying the model parameters on the resulting IMF and show that the IMF of stellar clusters is expected to show significant variations, provided variations in the clumps and cores properties exist.

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.

Recent measurements of the Zeeman effect in dark-cloud cores provide important tests for theories of cloud dynamics and prestellar core formation. In this Letter we report results of simulated Zeeman measurements, based on radiative transfer calculations through a snapshot of a simulation of supersonic and super-Alfv\'enic turbulence. We have previously shown that the same simulation yields a relative mass-to-flux ratio (core versus envelope) in agreement with the observations (and in contradiction with the ambipolar-drift model of core formation). Here we show that the mass-to-flux and turbulent-to-magnetic-energy ratios in the simulated cores agree with observed values as well. The mean magnetic field strength in the simulation is very low, \bar{B}=0.34 \muG, presumably lower than the mean field in molecular clouds. Nonetheless, high magnetic field values are found in dense cores, in agreement with the observations (the rms field, amplified by the turbulence, is B_{rms}=3.05 \muG). We conclude that a strong large-scale mean magnetic field is not required by Zeeman effect measurements to date, although it is not ruled out by this work.

Polarization carries information about the magnetic fields in interstellar clouds. The observations of polarized dust emission are used to study the role of magnetic fields in the evolution of molecular clouds and the initial phases of star-formation. We study the grain alignment with realistic simulations, assuming the radiative torques to be the main mechanism that spins the grains up. The aim is to study the efficiency of the grain alignment as a function of cloud position and to study the observable consequences of these spatial variations. Our results are based on the analysis of model clouds derived from MHD simulations. The continuum radiative transfer problem is solved with Monte Carlo methods to estimate the 3D distribution of dust emission and the radiation field strength affecting the grain alignment. We also examine the effect of grain growth in cores. We are able to reproduce the results of Cho & Lazarian using their assumptions. However, the anisotropy factor even in the 1D case is lower than their assumption of $\gamma = 0.7$, and thus we get less efficient radiative torques. Compared with our previous paper, the polarization degree vs. intensity relation is steeper because of less efficient grain alignment within dense cores. Without grain growth, the magnetic field of the cores is poorly recovered above a few $A_{\rm V}$. If grain size is doubled in the cores, the polarization of dust emission can trace the magnetic field lines possibly up to $A_{\rm V} \sim 10$ magnitudes. However, many of the prestellar cores may be too young for grain coagulation to play a major role. The inclusion of direction dependent radiative torque efficiency weakens the alignment. Even with doubled grain size, we would not expect to probe the magnetic field past a few magnitudes in $A_{\rm V}$.

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.

This White Paper, submitted to the National Academy of Sciences' Astro2010 Decadal Review Committee, focuses on 2 central themes in the study of young brown dwarfs -- their formation mechanism and disk characteristics -- which are of direct relevance to fundamental questions of stellar and planetary origins and properties.

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.

To evaluate the effect of turbulent heating in the thermal balance of interstellar clouds, we develop an extension of the log-Poisson intermittency model to supersonic turbulence. The model depends on a parameter, d, interpreted as the dimension of the most dissipative structures. By comparing the model with the probability distribution of the turbulent dissipation rate in a simulation of supersonic and super-Alfvenic turbulence, we find a best-fit value of d=1.64. We apply this intermittency model to the computation of the mass-weighted probability distribution of the gas temperature of molecular clouds, high-mass star-forming cores, and cold diffuse HI clouds. Our main results are: i) The mean gas temperature in molecular clouds can be explained as the effect of turbulent heating alone, while cosmic ray heating may dominate only in regions where the turbulent heating is low; ii) The mean gas temperature in high-mass star-forming cores with typical FWHM of ~6 km/s (corresponding to a 1D rms velocity of 2.5 km/s) may be completely controlled by turbulent heating, which predicts a mean value of approximately 36 K, two to three times larger than the mean gas temperature in the absence of turbulent heating; iii) The intermittency of the turbulent heating can generate enough hot regions in cold diffuse HI clouds to explain the observed CH+ abundance, if the rms velocity on a scale of 1 pc is at least 3 km/s, in agreement with previous results based on incompressible turbulence. Because of its importance in the thermal balance of molecular clouds and high-mass star-forming cores, the process of turbulent heating may be central in setting the characteristic stellar mass and in regulating molecular chemical reactions.

We present synthetic OH Zeeman splitting measurements of a super-Alfvenic model of molecular clouds. We select dense cores from synthetic 13CO maps computed from the largest simulation to date of supersonic and super-Alfvenic turbulence. The synthetic Zeeman splitting measurements in the cores yield a relation between the magnetic field strength, B, and the column density, N, in good agreement with the observations. The large scatter in B at a fixed value of N is partly due to intrinsic variations in the magnetic field strength from core to core. We also compute the relative mass-to-flux ratio between the center of the cores and their envelopes, ${\cal R}_{\mu}$, and show that super-Alfvenic turbulence produces a significant scatter also in ${\cal R}_{\mu}$, including negative values (field reversal between core center and envelope). We find ${\cal R}_{\mu}<1$ for 70% of the cores, and ${\cal R}_{\mu}<0$ for 12%. Of the cores with $|B_{\rm LOS}|>10$ \muG, 81% have ${\cal R}_{\mu}<1$. These predictions of the super-Alfvenic model are in stark contrast to the ambipolar drift model of core formation, where only ${\cal R}_{\mu}>1$ is allowed.

With current near-infrared (NIR) instruments the near-infrared light scattered from interstellar clouds can be mapped over large areas. The surface brightness carries information on the line-of-sight dust column density. Therefore, scattered light could provide an important tool to study mass distribution in quiescent interstellar clouds at a high, even sub-arcsecond resolution. We wish to confirm the assumption that light scattering dominates the surface brightness in all NIR bands. Furthermore, we want to show that scattered light can be used for an accurate estimation of dust column densities in clouds with Av in the range 1-15mag. We have obtained NIR images of a quiescent filament in the Corona Australis molecular cloud. The observations provide maps of diffuse surface brightness in J, H, and Ks bands. Using the assumption that signal is caused by scattered light we convert surface brightness data into a map of dust column density. The same observations provide colour excesses for a large number of background stars. These data are used to derive an extinction map of the cloud. The two, largely independent tracers of the cloud structure are compared. Results. In regions below Av=15m both diffuse surface brightness and background stars lead to similar column density estimates. The existing differences can be explained as a result of normal observational errors and bias in the sampling of extinctions provided by the background stars. There is no indication that thermal dust emission would have a significant contribution even in the Ks band. The results show that, below Av=15mag, scattered light does provide a reliable way to map cloud structure. Compared with the use of background stars it can also in practice provide a significantly higher spatial resolution.

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.

Our aim is to study the polarization of thermal dust emission to see if the alignment of grain by radiative torques could explain the observed relation between the degree of polarization and the intensity in dense cores. Predictions are made for polarimetry observations with the Planck satellite. Our results are based on model clouds derived from MHD simulations of magnetized turbulent flows, while the continuum radiative transfer problem is solved with Monte Carlo methods in order to estimate the three-dimensional distribution of dust emission and the radiation field strength affecting the grain alignment. The influence of grain alignment efficiency is examined in the calculated polarization maps. We are able to reproduce the P/I-relation with the grain alignment by radiative torques. The decrease in intrinsic polarization and total emission means that sub-mm polarimetry carries only little information about the magnetic fields in dense cores with high visual extinction. The interpretation of the observations will be further complicated by the unknown magnetic field geometry and the fact that what is observed as individual cores may, in fact, be a superposition of several density enhancements. According to our calculations, Planck will be able to map dust polarization reliably when A_V > 2 mag at spatial resolution of 15'.

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.

With current wide-field near-infrared (NIR) instruments the scattered light in the near-infrared can be mapped over large areas. Below A_V ~ 10 the surface brightness is directly proportional to the column density, and at slightly higher column densities the saturation of the intensity values can be corrected using the ratios of the intensity in different NIR bands. NIR scattered light provides a promising new method for the mapping of quiescent interstellar clouds. We develop a method to convert the observed near-infrared surface brightness into estimates of the column density. We study and quantify the effect that different error sources could have on the accuracy of such estimates. We also propose to reduce systematic errors by combining surface brightness data with extinction measurements derived from the near-infrared colour excess of background stars. Our study is based on a set of three-dimensional magnetohydrodynamic turbulence simulations. Maps of near-infrared scattered light are obtained with radiative transfer calculations, and the maps are converted back into column density estimates using the proposed method. The results are compared with the true column densities. Extinction measurements are simulated using the same turbulence simulations, and are used as a complementary column density tracer. We find that NIR intensities can be converted into a reliable estimate of the column density in regions with A_V up to almost 20mag.

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.

Developing a theory of low-mass star formation ($\sim 0.1$ to 3~M$_{\odot}$) remains one of the most elusive and important goals of theoretical astrophysics. The star-formation process is the outcome of the complex dynamics of interstellar gas involving non-linear interactions of turbulence, gravity, magnetic field and radiation. The evolution of protostellar condensations, from the moment they are assembled by turbulent flows to the time they reach stellar densities, spans an enormous range of scales, resulting in a major computational challenge for simulations. Since the previous Protostars and Planets conference, dramatic advances in the development of new numerical algorithmic techniques have been successfully implemented on large scale parallel supercomputers. Among such techniques, Adaptive Mesh Refinement and Smooth Particle Hydrodynamics have provided frameworks to simulate the process of low-mass star formation with a very large dynamic range. It is now feasible to explore the turbulent fragmentation of molecular clouds and the gravitational collapse of cores into stars self-consistently within the same calculation. The increased sophistication of these powerful methods comes with substantial caveats associated with the use of the techniques and the interpretation of the numerical results. In this review, we examine what has been accomplished in the field and present a critique of both numerical methods and scientific results. We stress that computational simulations should obey the available observational constraints and demonstrate numerical convergence. Failing this, results of numerical simulations do not advance our understanding of low-mass star formation.

We discuss the possibility of mapping interstellar clouds at unprecedentedly high spatial resolution by means of near-IR imaging of their scattered light. We calculate the scattering of the interstellar radiation field by a cloud model obtained from the simulation of a supersonic turbulent flow. Synthetic maps of scattered light are computed in the J, H and K bands and are found to allow an accurate estimate of column density, in the range of visual extinction between 1 and 20 magnitudes. We provide a formalism to convert the intensity of scattered light at these near-IR bands into a total gas column density. We also show that this new method of mapping interstellar clouds is within the capability of existing near-IR facilities, which can achieve a spatial resolution of up to ~ 0.1 arcsec. This opens new perspectives in the study of interstellar dust and gas structure on very small scales. The validity of the method has been recently demonstrated by the extraordinary images of the Perseus region obtained by Foster & Goodman (2005), which motivated this investigation.

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.

Low-energy cosmic rays (CRs) are confined by self-generated MHD waves in the mostly neutral ISM. We show that the CR transport equation can be expressed as a continuity equation for the CR number density involving an effective convection velocity. Assuming balance between wave growth and ion-neutral damping, this equation gives a steady-state condition n_cr ~ n_i^1/2 up to a critical density for free streaming. This relation naturally accounts for the heretofore unexplained difference in CR ionization rates derived for dense diffuse clouds (McCall et al. 2003) and dark clouds, and predicts large spatial variations in the CR heating rate and pressure.

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.

We propose to interpret the stellar IMF as a property of the turbulence in the star--forming gas. Gravitationally unstable density enhancements in the turbulent flow collapse and form stars. Their mass distribution can be derived analytically from the power spectrum of the turbulent flow and the isothermal shock jump conditions in the magnetized gas. For a power spectrum index \beta=1.74, consistent with Larson's velocity dispersion--size relation as well as with new numerical and analytic results on supersonic turbulence, we obtain a power law mass distribution of dense cores with a slope equal to 3/(4-\beta)=1.33, consistent with the slope of Salpeter's stellar IMF. Below one solar mass, the mass distribution flattens and turns around at a fraction of a solar mass, as observed for the stellar IMF in a number of stellar clusters, because only the densest cores are gravitationally unstable. The mass distribution at low masses is determined by the Log--Normal distribution of the gas density. The intermittent nature of this distribution is responsible for the generation of a significant number of collapsing cores of brown dwarf mass.

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 present describe a new computer code that solves the radiative transfer problem on multi-resolution grids. If the cloud model is from an MHD simulation on a regular cartesian grid, criteria based for example on local density or velocity gradients are used to refine the grid by dividing selected cells into sub-cells. Division can be repeated hierarchically. Alternatively, if the cloud model is from an MHD simulations with adaptive mesh refinement, the same multi-resolution grid used for the MHD simulation is adopted in the radiative transfer calculations. High discretization is often needed only in a small fraction of the total volume. This makes it possible to simulate spectral line maps with good accuracy, also minimizing the total number of cells and the computational cost (time and memory). Multi-resolution models are compared with models on regular grids. In the case of moderate optical depths (tau ~ a few) an accuracy of 10% can be reached with multi-resolution models where only 10% of the cells of the full grid are used. For optically thick species (tau~100), the same accuracy is achieved using 15% of the cells. The relation between accuracy and number of cells is not found to be significantly different in the two MHD models we have studied. The new code is used to study differences between LTE and non-LTE spectra and between isothermal and non-isothermal cloud models. We find significant differences in line ratios and individual spectral line profiles of the isothermal and LTE models relative to the more realistic non-isothermal case. The slope of the power spectrum of integrated intensity is instead very similar in all models.

Hundreds of brown dwarfs (BDs) have been discovered in the last few years in stellar clusters and among field stars. BDs are almost as numerous as hydrogen burning stars and so a theory of star formation should also explain their origin. The ``mystery'' of the origin of BDs is that their mass is two orders of magnitude smaller than the average Jeans' mass in star--forming clouds, and yet they are so common. In this work we investigate the possibility that gravitationally unstable protostellar cores of BD mass are formed directly by the process of turbulent fragmentation. Supersonic turbulence in molecular clouds generates a complex density field with a very large density contrast. As a result, a fraction of BD mass cores formed by the turbulent flow are dense enough to be gravitationally unstable. We find that with density, temperature and rms Mach number typical of cluster--forming regions, turbulent fragmentation can account for the observed BD abundance.

The determination of the fractional electron abundance, Xe, in protostellar cores relies on observations of molecules, such as DCO+, H13CO+ and CO, and on chemical models to interpret their abundance. Studies of protostellar cores have revealed significant variations of Xe from core to core within a range 10^-8<Xe<10^-6. The physical origin of these large variations in Xe is not well understood, unless unlikely variations in the cosmic ray ionization rate or ad hoc values of metal depletion are assumed. In this work we explore other potential causes of these variations in Xe, including core age, extinction and density. We compute numerically the intensity of the radiation field within a density distribution generated by supersonic turbulence. Taking into account the lines of sight in all directions, the effective visual extinction in dense regions is found to be always much lower than the extinction derived from the column density along a fixed line of sight. Dense cores with volume and column densities comparable to observed protostellar cores have relatively low mass-averaged visual extinction, 2mag <= A_V <= 5mag, such that photo-ionization can sometimes be as important as cosmic ray ionization. Chemical models, including gas-grain chemistry and time dependent gas depletion and desorption, are computed for values of visual extinction in the range 2mag <= A_V <= 6mag, and for a hydrogen gas density of 10^4cm^-3$, typical of protostellar cores. The models presented here can reproduce some of the observed variations of ion abundance from core to core as the combined effect of visual extinction and age variations. The range of electron abundances predicted by the models is relatively insenstive to density over 10^4 to 10^6 cm^{-3}.

Supersonic turbulent flows of magnetized gas are believed to play an important role in the dynamics of star-forming clouds in galaxies. Understanding statistical properties of such flows is crucial for developing a theory of star formation. In this letter we propose a unified approach for obtaining the velocity scaling in compressible and super--Alfv\'{e}nic turbulence, valid for arbitrary sonic Mach number, \ms. We demonstrate with numerical simulations that the scaling can be described with the She--L\'{e}v\^{e}que formalism, where only one parameter, interpreted as the Hausdorff dimension of the most intense dissipative structures, needs to be varied as a function of \ms. Our results thus provide a method for obtaining the velocity scaling in interstellar clouds once their Mach numbers have been inferred from observations.

We present a new set of stellar interior and synthesis models for predicting the integrated emission from stellar populations in star clusters and galaxies of arbitrary age and metallicity. This work differs from existing spectral synthesis codes in a number of important ways, namely (1) the incorporation of new stellar evolutionary tracks, with sufficient resolution in mass to sample rapid stages of stellar evolution; (2) a physically consistent treatment of evolution in the HR diagram, including the approach to the main sequence and the effects of mass loss on the giant and horizontal-branch phases. Unlike several existing models, ours yield consistent ages when used to date a coeval stellar population from a wide range of spectral features and colour indexes. We rigorously discuss degeneracies in the age-metallicity plane and show that inclusion of spectral features blueward of 4500 AA, suffices to break any remaining degeneracy and that with moderate S/N spectra (10 per 20AA, resolution element) age and metallicity are not degenerate. We also study sources of systematic errors in deriving the age of a single stellar population and conclude that they are not larger than 10-15%. We illustrate the use of single stellar populations by predicting the colors of primordial proto-galaxies and show that one can first find them and then deduce the form of the IMF for the early generation of stars in the universe. Finally, we provide accurate analytic fitting formulas for ultra fast computation of colors of single stellar populations. The models can be found at http://www.physics.upenn.edu/~raulj

The magnetic field strength in molecular clouds is a fundamental quantity for theories of star formation. It is estimated by Zeeman splitting measurements in a few dense molecular cores, but its volume--averaged value within large molecular clouds (over several parsecs) is still uncertain. In this work we provide a new method to constrain the average magnetic field strength in molecular clouds. We compare the power spectrum of gas density of molecular clouds with that of two $350^3$ numerical simulations of supersonic MHD turbulence. The numerical simulation with approximate equipartition of kinetic and magnetic energies (model A) yields the column density power spectrum $P(k)\propto k^{-2.25\pm 0.01}$, the super--Alfv\'{e}nic simulation (model B) $P(k)\propto k^{-2.71\pm 0.01}$. The column density power spectrum of the Perseus, Taurus and Rosetta molecular cloud complexes is found to be well approximated by a power law, $P_{\rm o}(k)\propto k^{-a}$, with $a=2.74\pm 0.07$, $2.74\pm 0.08$ and $2.76\pm 0.08$ respectively. We conclude that the observations are consistent with the presence of super--Alfv\'{e}nic turbulence in molecular clouds (model B) while model A is inconsistent (more than 99% confidence) with the observations.

The photoelectric heating is believed to be the main heating mechanism in cool HI clouds. The heating rate can be estimated through observations of the [CII] line emission, since this is the main coolant in regions where the photoelectric effect dominates the heating. Comparison of the [CII] emission with the far-infrared (FIR) emission allows to constrain the efficiency of the photoelectric heating, using model calculations that take into account the strength of the radiation field. Recent [CII] observations carried out with the ISO satellite have made this study possible. In this work we study the correlation between FUV absorption and FIR emission using three-dimensional models. The density distributions are obtained with numerical simulations of compressible magneto-hydrodynamic turbulence, with rms sonic Mach numbers 0.6<M<10. The FIR intensities are solved with detailed radiative transfer calculations. The [CII] line radiation is estimated assuming the [CII] line cooling equals the FUV absorption multiplied by the efficiency of the photoelectric heating, epsilon. The average ratio between the predicted [CII] and FIR emissions is found to be remarkably constant between different models, implying that the derived values of epsilon should not depend on the rms Mach number. The comparison with empirical data from translucent, high latitude clouds yields an estimate of the photoelectric heating efficiency of 2.9 10^-2. This value confirms previous theoretical predictions. Our models show that most of the scatter in the observed [CII] and FIR intensities can be understood as a result of the highly fragmented density field in turbulent HI clouds. The scatter can be reproduced with models with supersonic turbulence, while subsonic turbulence fails to generate the observed scatter.

We compute the spectral correlation function (SCF) of 13CO J=1--0 maps of molecular cloud complexes. The SCF is a power law over approximately an order of magnitude in spatial separation in every map. The power law slope of the SCF, alpha, its normalization, S_0(1pc), and the spectral line width averaged over the whole map, sigma_v, are computed for all the observational maps. The values of alpha, S_0(1pc) and sigma_v are combined to obtain empirical correlations to be used as tests for theoretical models of molecular clouds. Synthetic spectral maps are computed from different theoretical models, including solutions of the magneto-hydrodynamic (MHD) equations with different values of the rms Mach number of the flow and stochastic models with different power spectra of the velocity field. In order to compute the radiative transfer from the MHD models it is necessary to assign the models a physical scale and a physical density. When these assignments are made according to Larson type relations the best fit to the observational correlations is obtained. Unphysical stochastic models are instead ruled out by the empirical correlations. MHD models with equipartition of magnetic and kinetic energy of turbulence do not reproduce the observational data when their average magnetic field is oriented approximately parallel to the line of sight.

We present a numerical code for continuum radiative transfer that is based on the idea of a `library' describing the relation between the intensity of the local radiation field and the resulting dust emission. With this information and local intensities at a few reference wavelengths, the radiative transfer equation can be integrated through the source and an approximation of the emission spectrum is obtained. Tests with small models for which the radiative transfer problem can be solved directly show that with our method one can easily obtain an accuracy of a few per cent. We show spectra computed from three-dimensional MHD simulations containing up to 128^3 cells. The models represent starless, inhomogeneous interstellar clouds embedded in the normal interstellar radiation field. The intensity ratios between IRAS bands show large variations that follow the filamentary density distribution. The power law index of the spatial power spectrum of the column density map is -2.8. In infrared maps temperature variations increase the power at high spatial frequencies, and the index varies between -2.5 and -2.7. Assuming constant dust properties throughout the cloud, the IRAS ratio <I_60/I_100> decreases in densest cores only by a factor of ~4 compared with the value in diffuse medium. In observations the ratio can decrease twice as much even in optically thinner clouds. This requires that most of the small grains are removed in these regions, and possibly a modification of the properties of large grains.

Supersonic turbulence fragments the interstellar medium into dense sheets, filaments, cores and large low density voids. The turbulence is driven on large scales, probably predominantly by supernovae. The scaling properties of supersonic turbulence are well described by a new analytical theory. The distribution of proto-stellar core masses depends primarily on the power spectrum of the turbulent flow, and on the jump conditions for near isothermal shocks in a magnetized gas. For the predicted velocity power spectrum index beta=1.74, consistent with results of numerical experiments of supersonic turbulence as well as with Larson's velocity-size relation, one obtains by scaling arguments a power law mass distribution of dense cores with a slope equal to 3/(4-beta) = 1.33, consistent with the slope of the Salpeter stellar initial mass function (IMF). Results from numerical simulations confirm this scaling. Both the analytical model for the stellar IMF and its numerical estimate show that turbulent fragmentation can also explain the origin of brown dwarfs.

We compute the structure function scaling of a 2MASS extinction map of the Taurus molecular cloud complex. The scaling exponents of the structure functions of the extinction map follow the Boldyrev's velocity structure function scaling of supersonic turbulence. This confirms our previous result based on a spectral map of 13CO J=1-0 covering the same region and suggests that supersonic turbulence is important in the fragmentation of this star--forming cloud.

We compute the structure function scaling of the integrated intensity images of two J=1-0 13CO maps of Taurus and Perseus. The scaling exponents of the structure functions follow the velocity scaling of supersonic turbulence, suggesting that turbulence plays an important role in the fragmentation of cold interstellar clouds. The data also allows to verify the validity of the two basic assumptions of the hierarchical symmetry model, originally proposed for the derivation of the velocity structure function scaling. This shows that the same hierarchical symmetry holds for the projected density field of cold interstellar clouds.

The morphology and kinematics of molecular clouds (MCs) are best explained as the consequence of super--sonic turbulence. Super--sonic turbulence fragments MCs into dense sheets, filaments and cores and large low density ``voids'', via the action of highly radiative shocks. We refer to this process as "turbulent fragmentation". In this work we derive the mass distribution of gravitationally unstable cores generated by the process of turbulent fragmentation. The mass distribution above one solar mass depends primarily on the power spectrum of the turbulent flow and on the jump conditions for isothermal shocks in a magnetized gas. For a power spectrum index \beta=-1.74, consistent with Larson's velocity dispersion--size relation as well as with new numerical and analytic results on super--sonic turbulence, we obtain a power law mass distribution of dense cores with a slope equal to 3/(4-\beta) = 1.33, consistent with the slope of the stellar IMF. Below one solar mass, the mass distribution flattens and turns around at a fraction of a solar mass, as observed for the stellar IMF in a number of stellar clusters, because only the densest cores are gravitationally unstable. The mass distribution at low masses is determined by the probability distribution of the gas density, which is known to be approximately Log--Normal for an isothermal turbulent gas. The intermittent nature of the turbulent density distribution is thus responsible for the existence of a significant number of small collapsing cores, even of sub--stellar mass. Since turbulent fragmentation is unavoidable in super--sonically turbulent molecular clouds, and given the success of the present model in predicting the observed shape of the stellar IMF, we conclude that turbulent fragmentation is essential to the origin of the stellar IMF.

The interstellar medium (ISM) provides a unique laboratory for highly supersonic, driven hydrodynamics turbulence. We present a theory of such turbulence, confirm it by numerical simulations, and use the results to explain observational properties of interstellar molecular clouds, the regions where stars are born.

We present a direct numerical and analytical study of driven supersonic MHD turbulence that is believed to govern the dynamics of star-forming molecular clouds. We describe statistical properties of the turbulence by measuring the velocity difference structure functions up to the fifth order. In particular, the velocity power spectrum in the inertial range is found to be close to E(k) \~ k^{-1.74}, and the velocity difference scales as <|\Delta u|>~ L^{0.42}. The results agree well with the Kolmogorov--Burgers analytical model suggested for supersonic turbulence in [astro-ph/0108300]. We then generalize the model to more realistic, fractal structure of molecular clouds, and show that depending on the fractal dimension of a given molecular cloud, the theoretical value for the velocity spectrum spans the interval [-1.74 ... -1.89], while the corresponding window for the velocity difference scaling exponent is [0.42 ... 0.78].

We have computed line emission cooling rates for the main cooling species in models of interstellar molecular clouds. The models are based on numerical simulations of super-sonic magneto-hydrodynamic (MHD) turbulence. Non-LTE radiative transfer calculations have been performed to properly account for the complex density and velocity structures in the MHD simulations. Three models are used. Two of the models are based on MHD simulations with different magnetic field strength and the third includes the computation of self-gravity (in the super-Alfvenic regime of turbulence). The density and velocity fields in the simulations are determined self-consistently by the dynamics of super-sonic turbulence. The models are intended to represent molecular clouds with linear size L~6 pc and mean density <n>~300 cm^-3, with the density exceeding 10^4 cm^-3 in the densest cores. We present 12CO, 13CO, C18O, O2, OI, CI and H2O cooling rates in isothermal clouds with kinetic temperatures 10-80K. Analytical approximations are derived for the cooling rates. The inhomogeneity of the models reduces photon trapping and enhances the cooling in the densest parts of the clouds. Compared with earlier models the cooling rates are less affected by optical depth effects and are therefore higher. The main effects comes, however, from the density variation since cooling efficiency increases with density. This is very important for the cooling of the clouds as a whole since most cooling is provided by gas with density above the average.

We model the polarized thermal dust emission from protostellar cores that are assembled by super-sonic turbulent flows in molecular clouds. Self-gravitating cores are selected from a three dimensional simulation of super-sonic and super-Alfvenic magneto-hydrodynamic (MHD) turbulence. The polarization is computed in two ways. In model A it is assumed that dust properties and grain alignment efficiency are uniform; in model B it is assumed that grains are not aligned at visual extinction larger than 3 mag. The main results of this work are: i) Values of the degree of polarization P between 1 and 10% are typical, despite the super-Alfvenic nature of the turbulence; ii) A steep decrease of P with increasing values of the sub-mm dust continuum intensity I is always found in self--gravitating cores selected from the MHD simulations, if grains are not aligned above a certain value of visual extinction (model B); iii) The same behavior is hard to reproduce if grains are aligned independently of visual extinction (model A); iv) The Chandrasekhar-Fermi formula, corrected by a factor f=0.4, provides an approximate estimate of the average magnetic field strength in the cores. Sub-mm dust continuum polarization maps of quiescent protostellar cores and Bok globules always show a decrease in P with increasing value of I consistent with the predictions of our model B. We therefore conclude that sub-mm polarization maps of quiescent cores do not map the magnetic field inside the cores at visual extinction larger than approximately 3 mag. There is no inconsistency between the results from optical and near-IR polarized absorption of background stars, and the observed polarization of sub-mm dust continuum from quiescent cores. In both cases, grains at large visual extinction appear to be virtually unaligned.

We propose a new method to measure and map the gas scale height of nearby disk galaxies. This method is applied successfully to the Australia Telescope Compact Array interferometric HI survey of the Large Magellanic Cloud (LMC); it could also be applied to a significant number of nearby disk galaxies, thanks to the next generation of interferometric facilities, such as the extended VLA and CARMA. The method consists of computing the Spectral Correlation Function (SCF) for a spectral-line map of a face-on galaxy. The SCF quantifies the correlation between spectra at different map positions as a function of their separation, and is sensitive to the properties of both the gas mass distribution and the gas velocity field. It is likely that spatial correlation properties of the gas density and velocity fields in a galactic disk are sensitive to the value of the scale height of the gas disk. A scale-free turbulent cascade is unlikely to extend to scales much larger than the disk scale height, as the disk dynamics on those larger scales should be dominated by two dimensional motions. We find a clear feature in the SCF of the LMC HI disk, on the scale of approximately 180 pc, which we identify as the disk scale height. We are also tentatively able to map variations of the scale height over the disk.

Super-sonic turbulence fragments molecular clouds (MC) into a very complex density field with density contrasts of several orders of magnitude. A fraction of the gas is locked into dense and gravitationally bound cores, which collapse as proto-stars. This process can be studied with numerical simulations of super-sonic self-gravitating turbulence. In this work, we use numerical simulations of magneto-hydrodynamic (MHD), super-sonic, super-Alfv\'{e}nic and self-gravitating turbulence to compute the mass distribution of collapsing proto-stellar cores, which are selected as local density maxima. We find that the mass distribution of collapsing cores is consistent with the stellar initial mass function (IMF), suggesting that super-sonic turbulence may be responsible for the generation of the IMF. To support this conclusion we also show that the physical properties of the numerically selected cores are in agreement with the properties of observed NH$_3$ cores and that their magnetic field strength is consistent with Zeeman splitting measurements. In turbulent MCs, star formation occurs via the gravitational collapse of super-critical cores, formed by the turbulent flow, sub-critical cores being irrelevant for the process of star formation.

The fragmentation of molecular clouds (MC) into protostellar cores is a central aspect of the process of star formation. Because of the turbulent nature of super-sonic motions in MCs, it has been suggested that dense structures such as filaments and clumps are formed by shocks in a turbulent flow. In this work we present strong evidence in favor of the turbulent origin of the fragmentation of MCs. The most generic result of turbulent fragmentation is that dense post shock gas traces a gas component with a smaller velocity dispersion than lower density gas, since shocks correspond to regions of converging flows, where the kinetic energy of the turbulent motion is dissipated. Using synthetic maps of spectra of molecular transitions, computed from the results of numerical simulations of super--sonic turbulence, we show that the dependence of velocity dispersion on gas density generates an observable relation between the rms velocity centroid and the integrated intensity (column density), Sigma(V_0)-I, which is indeed found in the observational data. The comparison between the theoretical model (maps of synthetic 13CO spectra), with 13CO maps from the Perseus, Rosette and Taurus MC complexes, shows excellent agreement in the Sigma(V_0)-I relation. The Sigma(V_0)-I relation of different observational maps with the same total rms velocity are remarkably similar, which is a strong indication of their origin from a very general property of the fluid equations, such as the turbulent fragmentation process.

The effects of noise and sampling on the ``Spectral Correlation Function'' (SCF) introduced by Rosolowsky et al. 1999 are studied using observational data, numerical simulations of magneto-hydrodynamic turbulence, and simple models of Gaussian spectral line profiles. The most significant innovations of this paper are: i) the normalization of the SCF based on an analytic model for the effect of noise; ii) the computation of the SCF as a function of the spatial lag between spectra within a map. By computing of the dependence of the SCF on the spatial lag, S_o(dr), we have been able to conclude that: 1) S_o(dr) is a power law, with slope in the range of scales l_i-l_o. 2) The correlation outer scale, l_o, is determined by the size of the map, and no evidence for a true departure from self-similarity on large scales has been found. 3) The correlation inner scale, l_i, is a true estimate of the smallest self-similar scale in a map. 4) The spectral slope in a given region, is independent of velocity resolution (above a minimum resolution threshold), spatial resolution, and average spectrum quality. 5) Molecular transitions which trace higher gas density yield larger values of the spectral slope than transitions tracing lower gas density. 6) Nyquist sampling, bad pixels in detector arrays, and reference sharing data acquisition need to be taken into account for a correct determination of the SCF at lag dr=1. The value of the spectral slope, however, can be computed correctly without a detailed knowledge of observational procedures.

The high resolution and sensitivity provided by proposed Atacama Large Millimeter Array will reveal new small scale structures in many sources, e.g. star forming regions. Such inhomogeneities may not have been considered in the analysis of past observations but they will be essential to the understanding of future data. Radiative transfer methods are needed to interpret the observations and the presence of complicated source structures and small scale inhomogeneities requires 3D modelling. We will describe studies of molecular line emission we have made using models of inhomogeneous molecular clouds. These are based on MHD models and include e.g. thermal balance calculations for molecular clouds.

The first stars to form in protogalaxies must have primordial chemical composition. We refer to a protogalaxy that is forming stars of primordial composition, or very low metallicity ($Z\le 0.01 Z_{\odot}$), as ``Primordial Proto-Galaxy'' (PPG). PPGs contain little or no dust, and therefore their spectral energy distribution can be modelled from the rest-frame ultraviolet to the infrared without accounting for dust extinction and emission. We present the results of computing the photometric properties of high redshift PPGs at near-d-to-mid infrared wavelengths, that will soon be available with the new generation of infrared space telescopes, such as SIRTF and the NGST. We show that: i) PPGs at very high redshift (5<z<10) should be easily selected from deep near/mid IR surveys with a colour-colour criterion; ii) PPGs at redshift 5<z<10 can be detected at 8 $\mu$m with the NGST, if they have a constant star formation rate of at least 100 M$_{\odot}/$yr; iii) once the redshift of a PPG photometric candidate is determined, its near-to-mid infrared colors should provide strong constraints on the stellar IMF at zero or very low metallicity.

We have computed the evolving ultraviolet-millimeter spectral energy distribution (SED) produced by proto-galactic starbursts at high redshift, incorporating the chemical evolution of the interstellar medium in a consistent manner. Dust extinction is calculated in a novel way, that is not based on empirical calibrations of extinction curves, but rather on the lifetime of molecular clouds which delays the emergence of each successive generation of stars at ultraviolet wavelengths by typically 15 Myr. The predicted rest-frame far-infrared to millimeter-wave emission includes the calculation of molecular emission-line luminosities ($^{12}$CO and O$_{2}$ among other molecules) consistent with the evolving chemical abundances. Here we present details of this new model, along with the results of comparing its predictions with several high-redshift observables, namely: the ultraviolet SEDs of Lyman-limit galaxies, the high-redshift radiogalaxies 4C41.17 and 8C1435, the SCUBA sub-millimeter survey of the Hubble Deep Field (HDF), and the SEDs of intermediate-redshift elliptical galaxies. With our new reddening method, we are able to fit the spectrum of the Lyman-limit galaxy 1512-cB58, and we find an extinction of about 1.9 magnitudes at 1600 {\AA}. The model also predicts that most Lyman-limit galaxies should have a value of $\alpha$ inside that range, as it is observed. The 850 $\mu$m flux density of a typical Lyman-limit galaxy is expected to be only $\simeq 0.5$ mJy, and therefore the optical counterparts of the most luminous sub-mm sources in the HDF (or any other currently feasible sub-mm survey) are unlikely to be Lyman-break galaxies.

Although thermal pressure is unimportant dynamically in most molecular gas, the temperature is an important diagnostic of dynamical processes and physical conditions. This is the first of two papers on thermal equilibrium in molecular clouds. We present calculations of frictional heating by ion-neutral (or ambipolar) drift in three-dimensional simulations of turbulent, magnetized molecular clouds. We show that ambipolar drift heating is a strong function of position in a turbulent cloud, and its average value can be significantly larger than the average cosmic ray heating rate. The volume averaged heating rate per unit volume due to ambipolar drift, H_AD ~ |JxB|^2 ~ B^4/L_B^2, is found to depend on the rms Alfvenic Mach number, M_A, and on the average field strength, as H_AD ~ M_A^2<|B|>^4. This implies that the typical scale of variation of the magnetic field, L_B, is inversely proportional to M_A, which we also demonstrate.

In this work we explore the effects of adopting initial mass functions (IMFs) variable in time on the chemical evolution of the Galaxy. In order to do that we adopt a chemical evolution model which assumes two main infall episodes for the formation of the Galaxy and which proved to be successful in reproducing the majority of the observational constraints, at least for the case of a constant IMF. Different variable IMFs are tested with this model, all assuming that massive stars are preferentially formed in ambients of low metallicity. This implies that massive stars are formed preferentially at early times and at large galactocentric distances. Our numerical results have shown that all the variable IMFs proposed so far are unable to reproduce all the relevant observational constraints in the Galaxy and that a constant IMF still reproduces better the observations. In particular, variable IMFs of the kind explored here are unable to reproduce the observed abundance gradients even when allowing for changes in other chemical evolution model parameters as, for instance, the star formation rate. As a consequence of this we conclude that the G-dwarf metallicity distribution is best explained by infall with a large timescale and a constant IMF, since it is possible to find variable IMFs of the kind studied here, reproducing the G-dwarf metallicity but this worsen the agreement with other observational constraints.

We compare the statistical properties of J=1-0 13CO spectra observed in the Perseus Molecular Cloud with synthetic J=1-0 13CO spectra, computed solving the non-LTE radiative transfer problem for a model cloud obtained as solutions of the three dimensional magneto-hydrodynamic (MHD) equations. The model cloud is a randomly forced super-Alfvenic and highly super-sonic turbulent isothermal flow. The purpose of the present work is to test if idealized turbulent flows, without self-gravity, stellar radiation, stellar outflows, or any other effect of star formation, are inconsistent or not with statistical properties of star forming molecular clouds. We present several statistical results that demonstrate remarkable similarity between real data and the synthetic cloud. Statistical properties of molecular clouds like Perseus are appropriately described by random super-sonic and super-Alfvenic MHD flows. Although the description of gravity and stellar radiation are essential to understand the formation of single protostars and the effects of star formation in the cloud dynamics, the overall description of the cloud and of the initial conditions for star formation can apparently be provided on intermediate scales without accounting for gravity, stellar radiation, and a detailed modeling of stellar outflows. We also show that the relation between equivalent line width and integrated antenna temperature indicates the presence of a relatively strong magnetic field in the core B1, in agreement with Zeeman splitting measurements.

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.

The question of the shape of the density PDF for supersonic turbulence is addressed, using both analytical and numerical methods. For isothermal supersonic turbulence, the PDF is Log-Normal, with a width that scales approximately linearly with the Mach number. For a polytropic equation of state, with an effective gamma smaller than one, the PDF becomes skewed and reminiscent of (but not equal to) a power law on the high density side.

We investigate in detail the hypothesis that low surface brightness galaxies (LSB) differ from ordinary galaxies simply because they form in halos with large spin parameters. We compute star formation rates using the Schmidt law, assuming the same gas infall dependence on surface density as used in models of the Milky Way. We build stellar population models, predicting colours, spectra, and chemical abundances. We compare our predictions with observed values of metallicity and colours for LSB galaxies and find excellent agreement with all observables. In particular, integrated colours, colour gradients, surface brightness and metallicity match very well to the observed values of LSBs for models with ages larger than 7 Gyr and high values ($\lambda > 0.05$) for the spin parameter of the halos. We also compute the global star formation rate (SFR) in the Universe due to LSBs and show that it has a flatter evolution with redshift than the corresponding SFR for normal discs. We furthermore compare the evolution in redshift of $[Zn/H]$ for our models to those observed in Damped Lyman $\alpha$ systems by \scite{Pettini+97} and show that Damped Lyman $\alpha$ systems abundances are consistent with the predicted abundances at different radii for LSBs. Finally, we show how the required late redshift of collapse of the halo may constrain the power spectrum of fluctuations.

The dynamics of molecular clouds is characterized by supersonic random motions in the presence of a magnetic field. We study this situation using numerical solutions of the three-dimensional compressible magneto-hydrodynamic (MHD) equations in a regime of highly supersonic random motions. The non-LTE radiative transfer calculations are performed through the complex density and velocity fields obtained as solutions of the MHD equations, and more than 5x10^5 synthetic molecular spectra are obtained. We use a numerical flow without gravity or external forcing. The flow is super-Alfvenic and corresponds to model A of Padoan and Nordlund (1997). Synthetic data consist of sets of 90x90 synthetic spectra with 60 velocity channels, in five molecular transitions: J=1-0 and J=2-1 for 12CO and 13CO, and J=1-0 for CS. Though we do not consider the effects of stellar radiation, gravity, or mechanical energy input from discrete sources, our models do contain the basic physics of magneto-fluid dynamics and non-LTE radiation transfer and are therefore more realistic than previous calculations. As a result, these synthetic maps and spectra bear a remarkable resemblance to the corresponding observations of real clouds.

We compare the statistical properties of J=1-0 13CO spectra observed in the Perseus Molecular Cloud with synthetic J=1-0 13CO spectra from a 5 pc model cloud. The synthetic spectra are computed solving the non-LTE radiative transfer problem for a model cloud obtained as solutions of the 3-D magneto-hydrodynamic (MHD) equations in both the highly supersonic and super-Alfvenic regimes of random flows. We present several statistical results that demonstrate remarkable similarity between real data and the synthetic cloud. The three-dimensional structure and dynamics of molecular clouds like Perseus are appropriately described by random supersonic and super-Alfvenic MHD flows. Although the description of gravity and stellar radiation are essential to understand the formation of single protostars and the effects of star formation in the cloud dynamics, the overall description of the cloud and of the initial conditions for star formation can apparently be described on intermediate scales without accounting for gravity and stellar radiation.

Gravitational lensing studies suggest that the Universe may contain a population of dark galaxies; we investigate this intriguing possibility and propose a mechanism to explain their nature. In this mechanism a dark galaxy is formed with a low density disk in a dark halo of high spin parameter; such galaxies can have surface densities below the critical Toomre value for instabilities to develop, and following Kennicutt's work we expect these galaxies to have low star formation rates. The only stellar component of the galaxies is a halo system, formed during the collapse of the proto-galactic cloud. We compute synthetic stellar population models and show that, at a redshift $z=0.5$, such galaxies have apparent magnitudes $B \simeq 28, R \simeq 26$ and $I \simeq 25$, and could be unveiled by deep searches with the Hubble Space Telescope. Dark galaxies have an initial short blue phase and then become essentially invisible, therefore they may account for the blue population of galaxies at high redshift. We find a strong mass-dependence in the fraction of dark galaxies, and predict that spiral galaxies will not be found in halos with masses less than about $10^9 M_\odot$, if $\Omega=1$. Above about $10^{12}M_\odot$, all halos can produce luminous disks. The mass-dependence of the galaxy-formation efficiency introduces the possibility of `spin bias' -- luminous galaxies being associated preferentially with strongly-clustered high-mass halos. A further prediction is that the slope of the faint-end luminosity function for galaxies will be flatter than the associated halo mass function.

In this work we use models of molecular clouds (MC), and non-LTE radiative transfer calculations, to obtain a theoretical calibration of the relation between LTE 13CO column density and true column density in MCs. The cloud models consist of 3 dimensional grids of density and velocity fields obtained as solutions of the compressible magneto-hydrodynamic equations in a 128x128x128 periodic grid in both the supersonic and super-Alfvenic regimes. Due to the random nature of the velocity field and the presence of shocks, the density spans a continuous range of values covering over 5-6 orders of magnitude (from ~0.1 to ~10^5 cm^-3). As a result, the LTE column density can be calibrated over 3 orders of magnitude. We find that LTE column density of molecular clouds typically underestimates the mean 13CO true column density by a factor ranging from 1.3 to 7. These results imply that the standard LTE methods for the derivation of column densities from CO data systematically underestimate the true values independent of other major sources of uncertainty such as the relative abundance of CO.

Supersonic random motions are observed in dark clouds and are traditionally interpreted as Alfv\'{e}n 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: beta ~1 and beta << 1, where beta is the ratio of gas pressure and magnetic pressure: beta=P_g/P_m. Our results, combined with observational tests, show that the model with beta ~ 1 is consistent with the observed properties of molecular clouds, while the model with beta << 1 has several properties that are in conflict with the observations. We also 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. Magnetically dominated cores form early in the evolution, while later on the intermittency in the density field wins out.

We propose that the stellar initial mass function (IMF) is universal in the sense that its functional form arises as a consequence of the statistics of random supersonic flows. A model is developed for the origin of the stellar IMF, that contains a dependence on the average physical parameters (temperature, density, velocity dispersion) of the large scale site of star formation. The model is based on recent numerical experiments of highly supersonic random flows that have a strong observational counterpart. It is shown that a Miller-Scalo like IMF is naturally produced by the model for the typical physical conditions in molecular clouds. A more ``massive'' IMF in star bursts is also predicted.

We present new age and distance determinations for the Galactic Globular Clusters M55 and M5, using the luminosity function method (Jimenez & Padoan 1996, Padoan & Jimenez 1997). We find an age of $11.8 \pm 1.5$ Gyr for M55 and $11.1 \pm 0.7$ Gyr for M5. This confirms previous results (Jimenez et al. 1996, Sandquist et al. 1996) and allows to conclude that the oldest stars in the Universe are not older than 14 Gyr. We also find $m-M=14.13 \pm 0.11$ for M55, and $m-M=14.49 \pm 0.06$ for M5. These values agree with the ones obtained using the tip of the red giant branch (Jimenez et al. 1996) and the sub-dwarf fitting method (Sandquist et al. 1996).

We reconsider the problem of the age of the stellar discs of late-type Low Surface Brightness (LSB) galaxies by making use of a new IMF recently derived from numerical fluid dynamical simulations (Padoan, Nordlund and Jones, 1997). While a Miller-Scalo IMF cannot adequately describe the photometric properties of LSBs, when we apply the new Padoan et al. (1997) IMF to a simple exponential disc model with parameters appropriate to LSBs, we get excellent fits of the colors and color gradients. We then conclude that: a) The star formation history of LSB disc galaxies can be described by an initial burst of a few times $10^7$ yr followed by a quiescent period with only sporadic star formation; b) LSBs' discs are not young. The age of the LSB disc galaxies for which colors have been measured are all larger than about 10 Gyr$.

Using a new physical model for star formation (Padoan 1995) we have tested the possibility that globular clusters (GCs) are formed from primordial mass fluctuations, whose mass scale ($10^8$ - $10^9$ M$_{\odot}$) is selected out of a CDM spectrum by the mechanism of non-equilibrium formation of $H_2$. We show that such clouds are able to convert about 0.003 of their total mass into a bound system (GC) and about 0.02 into halo stars. The metal enriched gas is dispersed away from the GC by supernova explosions and forms the galactic disk. These mass ratios between GCs, halo and disk depend on the predicted IMF which is a consequence of the universal statistics of fluid turbulence. They also depend on the ratio of baryonic over non-baryonic mass ,$X_b$, and are comparable with the values observed in typical spiral galaxies for $X_b \approx 0.1-0.2$. The computed mass and radius for a GC ( $5\times 10^5$ M$_{\odot}$ and 30 pc) are in good agreement with the average values in the Galaxy. The model predicts an exponential cut off in the stellar IMF below 0.1 M$_{\odot}$ in GCs and 0.6 M$_{\odot}$ in the halo. The quite massive star formation in primordial clouds leads to a large number of supernovae and to a high blue luminosity during the first two Gyr of the life of every galaxy.

We extend our previous method to determine globular cluster ages using the luminosity function (Jimenez \& Padoan 1996). We show that the luminosity function depends on both age and distance modulus and that it is possible to distinguish between the two. This method provides at the same time independent determinations of distance and age of a GC by simply counting the number of stars found inside specified luminosity bins. The main uncertainties in other traditional methods for determining GCs ages are absent (e.g. mixing length, color-$T_{\rm eff}$ calibration, morphology of the color-magnitude diagram ). The distance modulus is the biggest uncertainty in determining the age of GCs. Here we show that the age can be determined with small uncertainty for any value of distance modulus using the LF and that the LF allows a determination of the distance modulus itself. This is explained by the fact that the luminosity function is affected by a change in distance-modulus in a way that is different from its time evolution. If GC stellar counts with statistical errors not larger than $3\%$ are available, the age can be determined with an uncertainty of about 0.4 Gyr (independent of distance modulus, mixing length and color calibration) and the distance modulus with an uncertainty of about 0.04 mag.

Lada et al. (1994) have described a method for studying the distribution of dust in dark clouds using infrared imaging surveys. In particular they show that the method provides some information about the structure of the gas (dust) on scales smaller than their resolution. In the present work we clarify the nature of the information provided by their method. We show that: 1) the 3D density field of the gas is well described by a Log-Normal distribution down to very small scales; 2) the power spectrum and the standard deviation of the 3D density field can be constrained; 3) the origin of such a structure of the density field is likely to be the supersonic turbulence in the gas. In fact we find a qualitative and quantitative agreement between the predictions based on recent numerical simulations of supersonic turbulence (Nordlund and Padoan 1996; Padoan, Nordlund and Jones 1996) and the constraints given by the infrared dust extinction measurements.

We present a new method to compute stellar ages in Globular Clusters (GC) that is ten times more precise than the traditional isochrone fitting procedure. The method relies on accurate stellar evolutionary tracks and on photometry for GCs complete down to the main sequence, and it is based on counting number of stars in two different regions of the CMD: the red giant branch and the main-sequence. We have applied this method to the globular cluster M68 and found an age of 16.4$\pm$0.2 Gyr for $(m-M)_V=15.3$. This new method reduces the error associated to the uncertainty in the distance modulus by a factor of two, the error due to the choice of the value for the mixing length parameter to almost zero and the error due to the colour-$T_{\rm eff}$ transformation to zero.

The role played by velocity fields in the fragmentation of a cold medium and in the formation of protostars is studied. The velocity field is modeled with a compressible turbulent flow. A supersonic turbulent velocity field can fragment the medium into clumps of mass smaller than a local Jeans' mass, and therefore stabilize the medium against the formation of protostars. Based on this idea, the protostar formation efficiency and the protostar mass distribution are determined as functions of the following ambient parameters: average density $n_{0}$, average temperature $T_{0}$, r.m.s. turbulent velocity $\sigma_{v,0}$ (or its Mach number ${\cal{M}}_{t}$), postshock cooling time (e.g. chemistry). The main results are: i) the protostar mass distribution and its dependence on the ambient parameters are quantified; ii) the characteristic protostar mass is $M_{J,cl}\propto n_{0}^{-1/2}T_{0}^{2} \sigma_{v,0}^{-1}$; iii) the protostar formation efficiency $e$ is higher for larger mean density, larger mean temperature, lower velocity dispersion on a given scale and longer postshock cooling time (e.g. lower metallicity): $e\propto n_{0}^{{3/2}(\beta-1)}T_{0}^{\beta -1} \sigma_{v,0}^{-5(\beta-1)}L_{0}^{3(\beta-1)}$, where $\beta>1$ is the exponent of the clump mass distribution; iv) the efficiency is quite sensitive to the ambient parameters and therefore to the dynamical evolution of the star forming system.