Search

Systematic surveys of massive clumps have been carried out to study the conditions leading to the formation of massive stars. These clumps are typically at large distances and unresolved, so their physical properties cannot be reliably derived from the observations alone. Numerical simulations are needed to interpret the observations. To this end, we generate synthetic Herschel observations using our large-scale star-formation simulation, where massive stars explode as supernovae driving the interstellar-medium turbulence. From the synthetic observations, we compile a catalogue of compact sources following the exact same procedure as for the Hi-GAL compact source catalogue. We show that the sources from the simulation have observational properties with statistical distributions consistent with the observations. By relating the compact sources from the synthetic observations to their 3D counterparts in the simulation, we find that the synthetic observations overestimate the clump masses by about an order of magnitude on average due to line-of-sight projection, and projection effects are likely to be even worse for Hi-GAL Inner Galaxy sources. We also find that a large fraction of sources classified as protostellar are likely to be starless, and propose a new method to partially discriminate between true and false protostellar sources., ZJL acknowledges financial support from China Scholarship Council (CSC) under grant no. 201606660003. PP and VMP acknowledge support by the Spanish MINECO under project AYA2017-88754-P. Computing resources for this work were provided by the NASA High-End Computing (HEC) Program through the NASA Advanced Supercomputing (NAS) Division at Ames Research Center. We acknowledge PRACE for awarding us access to Curie at GENCI@CEA, France. Storage and computing resources at the University of Copenhagen HPC centre, funded in part by Villum Fonden (VKR023406), were used to carry out part of the data analysis.

The mass distribution of dense cores is a potential key to understand the process of star formation. Applying dendrogram analysis to the CARMA-NRO Orion C$^{18}$O ($J$=1--0) data, we identify 2342 dense cores, about 22 \% of which have virial ratios smaller than 2, and can be classified as gravitationally bound cores. The derived core mass function (CMF) for bound starless cores which are not associate with protostars has a slope similar to Salpeter's initial mass function (IMF) for the mass range above 1 $M_\odot$, with a peak at $\sim$ 0.1 $M_\odot$. We divide the cloud into four parts based on the declination, OMC-1/2/3, OMC-4/5, L1641N/V380 Ori, and L1641C, and derive the CMFs in these regions. We find that starless cores with masses greater than 10 $M_\odot$ exist only in OMC-1/2/3, whereas the CMFs in OMC-4/5, L1641N, and L1641C are truncated at around 5--10 $M_\odot$. From the number ratio of bound starless cores and Class II objects in each subregion, the lifetime of bound starless cores is estimated to be 5--30 free-fall times, consistent with previous studies for other regions. In addition, we discuss core growth by mass accretion from the surrounding cloud material to explain the coincidence of peak masses between IMFs and CMFs. The mass accretion rate required for doubling the core mass within a core lifetime is larger than that of Bondi-Hoyle accretion by a factor of order 2. This implies that more dynamical accretion processes are required to grow cores., Comment: 58 pages, 33 figures, 21 tables, accepted by ApJS

Observations have indicated that the pre-stellar core mass function (CMF) is similar to the stellar initial mass function (IMF), except for an offset towards larger masses. This has led to the idea that there is a one-to-one relation between cores and stars, such that the whole stellar mass reservoir is contained in a gravitationally bound pre-stellar core, as postulated by the core-collapse model, and assumed in recent theoretical models of the stellar IMF. We test the validity of this assumption by comparing the final mass of stars with the mass of their progenitor cores in a high-resolution star formation simulation that generates a realistic IMF under physical condition characteristic of observed molecular clouds. Using a definition of bound cores similar to previous works we obtain a CMF that converges with increasing numerical resolution. We find that the CMF and the IMF are closely related in a statistical sense only; for any individual star there is only a weak correlation between the progenitor core mass and the final stellar mass. In particular, for high-mass stars only a small fraction of the final stellar mass comes from the progenitor core, and even for low-mass stars the fraction is highly variable, with a median fraction of only about 50 per cent. We conclude that the core-collapse scenario and related models for the origin of the IMF are incomplete. We also show that competitive accretion is not a viable alternative., We are grateful to the anonymous referee for a detailed and useful report. PP and VMP acknowledge support by the Spanish MINECO under project AYA2017-88754-P, and financial support from the State Agency for Research of the Spanish Ministry of Science and Innovation through the ‘Unit of Excellence María de Maeztu 2020-2023’ award to the Institute of Cosmos Sciences (CEX2019-000918-M). The research leading to these results has received funding from the Independent Research Fund Denmark through grant No. DFF 8021-00350B (TH). We acknowledge PRACE for awarding us access to Curie at GENCI@CEA, France. Storage and computing resources at the University of Copenhagen HPC centre, funded in part by Villum Fonden (VKR023406), were used to carry out part of the data analysis.

Takemura, H., et al., Applying dendrogram analysis to the CARMA-NRO C18O (J = 1-0) data having an angular resolution of ∼8″, we identified 692 dense cores in the Orion Nebula Cluster region. Using this core sample, we compare the core and initial stellar mass functions in the same area to quantify the step from cores to stars. About 22% of the identified cores are gravitationally bound. The derived core mass function (CMF) for starless cores has a slope similar to Salpeter's stellar initial mass function (IMF) for the mass range above 1 M o˙, consistent with previous studies. Our CMF has a peak at a subsolar mass of ∼0.1 M o˙, which is comparable to the peak mass of the IMF derived in the same area. We also find that the current star formation rate is consistent with the picture in which stars are born only from self-gravitating starless cores. However, the cores must gain additional gas from the surroundings to reproduce the current IMF (e.g., its slope and peak mass), because the core mass cannot be accreted onto the star with 100% efficiency. Thus, the mass accretion from the surroundings may play a crucial role in determining the final stellar masses of stars., P.S. was partially supported by a Grant-in-Aid for Scientific Research (KAKENHI Number 18H01259) of Japan Society for the Promotion of Science (JSPS). R.S.K. acknowledges financial support from the DFG via the collaborative research center (SFB 881, Project-ID 138713538) “The Milky Way System” (subprojects A1, B1, B2, and B8). He is also thankful for subsidies from the Heidelberg Cluster of Excellence STRUCTURES in the framework of Germanyʼs Excellence Strategy (grant EXC-2181/1-390900948) and for funding from the European Research Council (ERC) via the ERC Synergy Grant ECOGAL (grant 855130). P.P. acknowledges support by the Spanish MINECO under project AYA2017- 88754-P, and financial support from the State Agency for Research of the Spanish Ministry of Science and Innovation through the “Unit of Excellence María de Maeztu 2020-2023” award to the Institute of Cosmos Sciences (CEX2019-000918-M). V.O., A.S.M., and P.S. were supported by the Collaborative Research Centre 956, subprojects C1, A6, and C3, funded by the Deutsche Forschungsgemeinschaft (DFG), project ID 184018867. T.G.S.P. gratefully acknowledges support by the National Science Foundation under grant No. AST-2009842. J.E.P. acknowledges the support by the Max Planck Society. We thank the anonymous referee for many useful comments that have improved the presentation

While the importance of supernova feedback in galaxies is well established, its role on the scale of molecular clouds is still debated. In this work, we focus on the impact of supernovae on individual clouds, using a high-resolution magneto-hydrodynamic simulation of a region of 250 pc where we resolve the formation of individual massive stars. The supernova feedback is implemented with real supernovae that are the natural evolution of the resolved massive stars, so their position and timing are self-consistent. We select a large sample of molecular clouds from the simulation to investigate the supernova energy injection and the resulting properties of molecular clouds. We find that molecular clouds have a lifetime of a few dynamical times, less then half of them contract to the point of becoming gravitationally bound, and the dispersal time of bound clouds, of order one dynamical time, is a factor of two shorter than that of unbound clouds. We stress the importance of internal supernovae, that is massive stars that explode inside their parent cloud, in setting the cloud dispersal time, and their huge overdensity compared to models where the supernovae are randomly distributed. We also quantify the energy injection efficiency of supernovae as a function of supernova distance to the clouds. We conclude that intermittent driving by supernovae can maintain molecular-cloud turbulence and may be the main process of cloud dispersal. The role of supernovae in the evolution of molecular clouds cannot be fully accounted for without a self-consistent implementation of their feedback., Comment: 33 pages, 23 figures, submitted to ApJ

We a dress 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 inertial-inflow model. This model stems directly from the idea that the mass distribution of stars is primarily the result of turbulent fragmentation. Under this hypothesis, the statistical properties of 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 analyzing 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 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 the 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., PP acknowledges support by the Spanish MINECO under project AYA2017-88754-P. The work of ˚AN was supported by grant 1323-00199B from the Danish Council for Independent Research (DFF), and by the Centre for Star and Planet Formation, which is funded by the Danish National Research Foundation (DNRF97). LP acknowledges financial support from NSFC under grant No. 11973098. MJ acknowledges the support of the Academy of Finland Grant No. 285769. Computing resources for this work were provided by the NASA High-End Computing (HEC) Program through the NASA Advanced Supercomputing (NAS) Division at Ames Research Center. We thankfully acknowledge the computer resources at MareNostrum and the technical support provided by Barcelona Supercomputing Center (AECT-2018-3-0019). We acknowledge PRACE for awarding us access to Curie at GENCI@CEA, France. Storage and computing resources at the University of Copenhagen HPC centre, funded in part by Villum Fonden (VKR023406), were used to carry out part of the data analysis. We thank the anonymous referee for a very useful report and the following colleagues for reading the first draft and providing comments: Philippe André, Javier Ballesteros-Paredes, Henrik Beuther, Gilles Chabrier, Lee Hartmann, Patrick Hennebelle, Alessio Traficante, Enrique Vazquez-Semadeni. Special thanks to Patricio Sanhueza and Delphine Russeil for providing the electronic tables to compute the observational core mass functions shown in Figure 26, and to Alessio Traficante for many instructive discussions on the observations of infall rates in massive clumps, and for providing the electronic tables to produce the plots in Figures 29 and 30.

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., Comment: Accepted to A&A (20 pages, 26 figures)

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 (N–S) 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 N–S runs are accurate and satisfy all the criteria. In particular, they satisfy our exact results that, at steady state, the average of the velocity divergence conditioned on the flow density and the conditional average of the advection of density both vanish at all density levels. 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 work tends to preferentially convert kinetic energy into thermal energy, which is inconsistent with the exact result that the energy exchange by 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 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., Comment: Submitted to ApJ

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., Comment: 19 pages, 20 figures. Accepted for publication in Astronomy and Astrophysics

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., 21 pages, 11 figures, accepted by ApJ, comments welcome shuo.kong@yale.edu

Context. The all-sky survey from the Planck space telescope has revealed that thermal emission from Galactic dust is polarized on scales ranging from the whole sky down to the inner regions of molecular clouds. Polarized dust emission can therefore be used as a probe for magnetic fields on 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 where star formation takes place. Aims. The orientation of the magnetic field with respect to the density structures has been investigated using different methods. The goal of this paper is to explicitly compare two of these: the Rolling Hough Transform (RHT) and the gradient technique (GRAD). Methods. We generated synthetic surface brightness maps at 353 GHz (850 μm) via magnetohydrodynamic simulations. We applied RHT and GRAD to two morphologically different regions identified in our maps. Region 1 is dominated by a dense and thick filamentary structure with some branches, while Region 2 includes a thinner filament with denser knots immersed in a more tenuous medium. Both methods derive the relative orientation between the magnetic field and the density structures, to which we applied two statistics, the histogram of relative orientation and the projected Rayleigh statistic, to quantify the variations of the relative orientation as a function of column density. Results. Both methods find areas with significant signal, and these areas are substantially different. In terms of relative orientations, in all our considered cases the predominant orientation of the density structures is perpendicular to the direction of the magnetic field. When the methods are applied to the same selected areas the results are consistent with each other in Region 2 but show some noticeable differences in Region 1. In Region 1, RHT globally finds the relative orientation becoming more perpendicular for increasing column density, while GRAD, applied at the same resolution as RHT, gives the opposite trend. These disparities are caused by the intrinsic differences in the methods and in the structures that they select. Conclusions. Our results indicate that the interpretation of the relative orientation between the magnetic field and density structures should take into account the specificity of the methods used to determine such orientation. The combined use of complementary techniques such as RHT and GRAD provides more complete information, which can be advantageously used to better understand the physical mechanisms operating in magnetized molecular clouds.

With this research note, we are releasing the CARMA-NRO Orion Survey data first presented in Kong et al., enhanced with additional coverage of the L1641-C region to the south of the integral-shaped filament. We are including position–position–velocity cubes for the molecular lines 12CO(1–0), 13CO(1–0), and C18O(1–0). The original paper includes details of the data reduction and final sensitivity. The mapped region now spans about 2.°5 along the Orion A cloud in the north–south direction, providing an unprecedented overview of the extended molecular gas. The associated data cubes are publicly available on the Harvard Dataverse 10.7910/DVN/6Q26PN (Version 3). Kong et al. and this research note should be cited when using these data.

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, Comment: 46 pages, 28 figures, 2 tables, accepted by ApJS

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., 25 pages, 21 figures, Accepted by ApJ

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., Comment: ApJ, in press

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., Comment: Accepted by Contrib. Astron. Obs. Skalnate Pleso (Proceedings of the conference Stars with a stable magnetic field: from pre-main sequence to compact remnants, Brno, 2017, editors: E. Paunzen, M. Netopil, V. Petit)

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., accepted for publication in ApJ

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., Comment: ApJ, in press

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., Accepted by ApJ

Supersonic magneto-hydrodynamic (MHD) turbulence in molecular clouds (MCs) plays an important role in the process of star formation. The effect of the turbulence on the cloud fragmentation process depends on the magnetic field strength. In this work we discuss the idea that the turbulence is super-Alfvénic, at least with respect to the cloud mean magnetic field. We argue that MCs are likely to be born super-Alfvénic. We then support this scenario based on a recent simulation of the large-scale warm interstellar medium turbulence. Using small-scale isothermal MHD turbulence simulation, we also show that MCs may remain super-Alfvénic even with respect to their rms magnetic field strength, amplified by the turbulence. Finally, we briefly discuss the comparison with the observations, suggesting that super-Alfvénic turbulence successfully reproduces the Zeeman measurements of the magnetic field strength in dense MC clouds.

This work presents a new physical model of the star formation rate (SFR), tested with a large set of 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, αvir, the rms sonic Mach number, S,0, and the ratio of mean gas pressure to mean magnetic pressure, β0. The SFR is predicted to decrease with increasing αvir (stronger turbulence relative to gravity), and to depend weakly on S,0 and β0, for values typical of star forming regions (S,0≈4-20 and β0≈1-20). 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 agree very well with the theoretical predictions.

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}$., Comment: 12 pages, 15 figures, submitted to A&A 19.12.2008; 09.01.2009: Corrected the name of Juvela; 24.04.2009: revised, added content, 13 pages, 16 figures; 18.06.2009: Language edited, print version

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., Accepted by ApJ, 15 pages, 7 figures

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.

The predictions of the Padoan and Nordlund IMF model are tested using the largest simulations of supersonic hydrodynamic (HD) and magneto-hydrodynamic (MHD) turbulence to date (~10003computational zones). The striking difference between the HD and MHD regimes, predicted by the model, is recovered.

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., Comment: 5 pages, 5 figures, revised version

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., 4 pages, 1 figure included, to appear in ApJ Letters (Corrections after ApJ proofs)

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 discuss a model for driven supersonic, super-Alfvenic MHD turbulence that is believed to govern the structure of molecular clouds. Such turbulence is highly intermittent; we describe its statistical properties by obtaining scaling of velocity-difference structure functions. This scaling was analytically predicted in Boldyrev (2002), confirmed in numerical simulations by Boldyrev et al. (2002), and discovered in observations by Padoan et al. (2003).

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., submitted to ApJ

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., submitted to ApJ

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 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., Comment: 12 pages, 10 figures. Revised version, accepted to A&A

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.

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 ~ 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]., Comment: 17 pages, 6 figures included

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., 20 pages, 9 figures, accepted to ApJ

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., 21 pages, 16 figures, accepted for publication in ApJ

PP is supported by the FP7-PEOPLE- 2010-RG grant PIRG07-GA-2010- 261359. Simulations by PP were carried out on the NASA/Ames Pleiades supercomputer, and under the PRACE project pra50751 running on SuperMUC at the LRZ (project ID pr86li). CF thanks for support from the Australian Research Council for a Discovery Projects Fellowship (Grant DP110102191). NJE was supported by NSF Grant AST-1109116 to the University of Texas at Austin. The research of CFM is supported in part by NSF grant AST-1211729 and NASA grant NNX13AB84G. DJ is supported by the National Research Council of Canada and by a Natural Sciences and Engineering Research Council of Canada (NSERC) Discovery Grant. JKJ is supported by a Lundbeck Foundation Junior Group Leader Fellowship. Research at Centre for Star and Planet Formation was funded by the Danish National Research Foundation and the University of Copenhagens Programme of Excellence. Supercomputing time at Leibniz Rechenzentrum (PRACE projects pr86li, pr89mu, and project pr32lo), at Forschungszentrum J¨ulich (project hhd20), and at DeIC/KU in Copenhagen are gratefully acknowledged.

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 supersonic magnetohydrodynamic (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 (one model is super-Alfvenic, while the other has equipartition of magnetic and kinetic energy). The third model 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 supersonic 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 104 cm-3 in the densest cores. We present 12CO, 13CO, C18O, O2, O I, C I, and H2O cooling rates in isothermal clouds with kinetic temperatures 10-80 K. 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. The main effects come, 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.

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^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., Comment: 37 pages, 9 figures included

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., Comment: Submitted to ApJ

In this work we use models of molecular clouds and non-LTE radiative transfer calculations to compare the column densities of molecular clouds with their LTE 13CO column densities. The cloud models consist of three-dimensional grids of density and velocity fields obtained as solutions of the compressible magnetohydrodynamic equations in a 1283 periodic grid in both the supersonic and super-Alfvenic regimes. Because of the random nature of the velocity field and the presence of shocks, the densities span a continuous range of values covering about 6 orders of magnitude (from ~0.1 to ~105 cm-3). As a result, the LTE column densities can be calculated over 3 orders of magnitude. We find that LTE column densities of molecular clouds typically underestimate the mean 13CO true column densities by factors 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.

In this work we explore the effects on the chemical evolution of the Galaxy of adopting initial mass functions (IMFs) variable in time. To this end, we adopt a chemical evolution model that assumes two main infall episodes for the formation of the Galaxy that has 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 none of the variable IMFs proposed so far are able to reproduce all the relevant observational constraints in the Galaxy, and that a constant IMF still better reproduces the observations. In particular, variable IMFs of the kind explored here are unable to reproduce the observed abundance gradients, even allowing for changes in other chemical evolution model parameters, such as, for instance, the star formation rate. As a consequence, 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 worsens the agreement with other observational constraints.

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., 10 figures, 2 tables included

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., Comment: 33 pages, 12 figures included, 5 jpeg figures not included (fig1a, fig1b, fig3, fig4 fig5), submitted to ApJ

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 reconsider the problem of the age of the stellar disks of late-type low surface brightness galaxies (LSBs) by making use of a new IMF recently derived from numerical fluid dynamics simulations (Padoan, Nordlund, & Jones 1997), and a new synthetic stellar population code, based on Jimenez & MacDonald (1997) evolutionary tracks and Kurucz atmospheric models (Kurucz 1992). We find that the disks of LSBs are not necessarily formed very recently. Their colors seem to indicate that the disks of LSBs started to form stars at least 7 Gyr ago, and more likely about 9 Gyr ago, and therefore contrary to what has been claimed in the literature disks of LSBs are not young.