Showing 1–50 of 116 results for author: Padoan, P

Protoplanetary Disks from Pre-Main Sequence Bondi-Hoyle Accretion

Authors: Paolo Padoan, Liubin Pan, Veli-Matti Pelkonen, Troels Haugboelle, AAke Nordlund

Abstract: Protoplanetary disks are routinely described as finite mass reservoirs left over by the gravitational collapse of the protostar, an assumption that strongly constrains both disk evolution and planet formation models. We propose a different scenario where protoplanetary disks of pre-main sequence stars are assembled primarily by Bondi-Hoyle accretion from the parent gas cloud. We demonstrate that B… ▽ More Protoplanetary disks are routinely described as finite mass reservoirs left over by the gravitational collapse of the protostar, an assumption that strongly constrains both disk evolution and planet formation models. We propose a different scenario where protoplanetary disks of pre-main sequence stars are assembled primarily by Bondi-Hoyle accretion from the parent gas cloud. We demonstrate that Bondi-Hoyle accretion can supply not only the mass, but also the angular momentum necessary to explain the observed size of protoplanetary disks, and we predict the dependence of the disk specific angular momentum on the stellar mass. Our results are based on an analytical derivation of the scaling of the angular momentum in a turbulent flow, which we also confirm with a numerical simulation of supersonic turbulence. This new scenario for disk formation and evolution may alleviate a number of observational problems as well as compel major revisions of disk and planet formation models. △ Less

Submitted 12 May, 2024; originally announced May 2024.

Comments: Submitted to Nature Astronomy

Origin and Evolution of Angular Momentum of Class II Disks

Authors: Veli-Matti Pelkonen, Paolo Padoan, Mika Juvela, Troels Haugbølle, Åke Nordlund

Abstract: Context. While Class II stars have already accreted most of their mass, the continued inflow of fresh material via Bondi-Hoyle accretion acts as an additional mass reservoir for their circumstellar disks. This may explain the observed accretion rates of pre-main- sequence (PMS) stars, as well as observational inconsistencies in the mass and angular momentum balance of their disks. Aims. Using a ne… ▽ More Context. While Class II stars have already accreted most of their mass, the continued inflow of fresh material via Bondi-Hoyle accretion acts as an additional mass reservoir for their circumstellar disks. This may explain the observed accretion rates of pre-main- sequence (PMS) stars, as well as observational inconsistencies in the mass and angular momentum balance of their disks. Aims. Using a new simulation that reproduces the stellar IMF, we want to quantify the role of Bondi-Hoyle accretion in the formation of Class II disks, as well address the prospect of its observational detection with the James Webb Space Telescope (JWST). Methods. We study the mass and angular momentum of the accreting gas using passively-advected tracer particles in the simulation, and we carry out radiative transfer calculations of near-infrared scattering to generate synthetic JWST observations of Bondi-Hoyle trails of PMS stars. Results. Gas accreting on Class II PMS stars approximately 1 Myr after their formation has enough mass and angular momentum to strongly affect the evolution of the preexisting disks. The accreted angular momentum is large enough to also explain the observed size of Class II disks. The orientation of the angular momentum vector can differ significantly from that of the previously accreted gas, which may result in a significant disk warping or misalignment. We also predict that JWST observations of Class II stars will be able to detect Bondi-Hoyle trails with a 90% success rate with only 2 min exposure time, if stars with accretion rates \dot{M} > 5e-10 Msol/yr and luminosity of L > 0.5 Lsol are selected. △ Less

Submitted 10 May, 2024; originally announced May 2024.

Comments: 12 pages, 16 figures, submitted to A&A

Updated kinematics of the Radcliffe Wave: non-synchronous, dipole-like vertical oscillations

Authors: Zhi-Kai Zhu, Min Fang, Zu-Jia Lu, Junzhi Wang, Guang-Xing Li, Shiyu Zhang, Veli-Matti Pelkonen, Paolo Padoan, En-Wei Liang

Abstract: The kinematic structure of the Radcliffe Wave (RW) is crucial for understanding its origin and evolution. In this work, we present an accurate measurement of the vertical velocity $V_Z$ by where the radial velocity (RV) measures are taken into consideration. This is achieved in two ways. First, the velocities are measured towards Young Stellar Objects (YSOs), using their RV and proper motion measu… ▽ More The kinematic structure of the Radcliffe Wave (RW) is crucial for understanding its origin and evolution. In this work, we present an accurate measurement of the vertical velocity $V_Z$ by where the radial velocity (RV) measures are taken into consideration. This is achieved in two ways. First, the velocities are measured towards Young Stellar Objects (YSOs), using their RV and proper motion measurements from APOGEE-2 and Gaia DR3. Second, we combine RV measurements toward clouds with proper motion measurements of associated YSOs to determine the vertical velocities of the clouds. The results reveal that the oscillations in $V_Z$ are not synchronous with the vertical coordinate. The difference is caused by a combination of the effect of the radial velocity which we include in this paper, and the difference in models. By supplementing our analysis with additional young star samples, we find a consistent dipole pattern in $V_Z$. The fact that no significant amplitude differences are found among the analyzed samples indicates that there is no apparent age gradient within the dipole. We propose that RW evolves at a relatively slow rate. The fact that it will take a much longer time for RW to complete a full period compared to the cloud lifetimes challenges its classification as a traditional "wave". This age discrepancy should explain the phase difference, and non-synchronous oscillation found in kinematic studies. △ Less

Submitted 23 February, 2024; originally announced February 2024.

Comments: 13 pages, 9 figures, submitted on 6 Feb 2024

Evaporation Ages: a New Dating Method for Young Star Clusters

Authors: V. -M. Pelkonen, N. Miret-Roig, P. Padoan

Abstract: The ages of young star clusters are fundamental clocks to constrain the formation and evolution of pre-main-sequence stars and their protoplanetary disks and exoplanets. However, dating methods for very young clusters often disagree, casting doubts on the accuracy of the derived ages. We propose a new method to derive the kinematic age of star clusters based on the evaporation ages of their stars.… ▽ More The ages of young star clusters are fundamental clocks to constrain the formation and evolution of pre-main-sequence stars and their protoplanetary disks and exoplanets. However, dating methods for very young clusters often disagree, casting doubts on the accuracy of the derived ages. We propose a new method to derive the kinematic age of star clusters based on the evaporation ages of their stars. The method is validated and calibrated using hundreds of clusters identified in a supernova-driven simulation of the interstellar medium forming stars for approximately 40 Myr within a 250 pc region. We demonstrate that the clusters' evaporation-age uncertainty can be as small as about 10% for clusters with a large enough number of evaporated stars and small but realistic observational errors. We have obtained evaporation ages for a pilot sample of 10 clusters, finding a good agreement with their published isochronal ages. The evaporation ages will provide important constraints for modeling the pre-main-sequence evolution of low-mass stars, as well as to investigate the star-formation and gas-evaporation history of young clusters. These ages can be more accurate than isochronal ages for very young clusters, for which observations and models are more uncertain. △ Less

Submitted 16 February, 2024; v1 submitted 14 November, 2023; originally announced November 2023.

Comments: 13 pages, 11 figures, 2 tables, Submitted to A&A on Nov 14th 2023, Accepted 9 January 2024

Will ALMA Reveal the True Core Mass Function of Protoclusters?

Authors: Paolo Padoan, Veli-Matti Pelkonen, Mika Juvela, Troels Haugbølle, Åke Nordlund

Abstract: Characterizing prestellar cores in star-forming regions is an important step towards the validation of theoretical models of star formation. Thanks to their sub-arcsecond resolution, ALMA observations can potentially provide samples of prestellar cores up to distances of a few kpc, where regions of massive star formation can be targeted. However, the extraction of real cores from dust-continuum ob… ▽ More Characterizing prestellar cores in star-forming regions is an important step towards the validation of theoretical models of star formation. Thanks to their sub-arcsecond resolution, ALMA observations can potentially provide samples of prestellar cores up to distances of a few kpc, where regions of massive star formation can be targeted. However, the extraction of real cores from dust-continuum observations of turbulent star-forming clouds is affected by complex projection effects. In this work, we study the problem of core extraction both in the idealized case of column-density maps and in the more realistic case of synthetic 1.3\,mm ALMA observations. The analysis is carried out on 12 regions of high column density from our 250 pc simulation. We find that derived core masses are highly unreliable, with only {\em a weak correlation between the masses of cores selected in the synthetic ALMA maps and those of the corresponding three-dimensional cores}. The fraction of real three-dimensional cores detected in the synthetic maps increases monotonically with mass and remains always below 50\%. Above $\sim 1\,M_{\odot}$, the core mass function derived from the column-density maps is steeper than that of the three-dimensional cores, while the core mass function from the synthetic ALMA maps has a slope closer to that of the real three-dimensional cores. Because of the mass uncertainties, proper guidance from realistic simulations is essential if ALMA observations of protoclusters at kpc distances are to be used to test star-formation models. △ Less

Submitted 4 May, 2023; v1 submitted 18 January, 2023; originally announced January 2023.

Comments: 20 pages, 21 figures, MNRAS, in press

CARMA-NRO Orion Survey: unbiased survey of dense cores and core mass functions in Orion A

Authors: Hideaki Takemura, Fumitaka Nakamura, Héctor G. Arce, Nicola Schneider, Volker Ossenkopf-Okada, Shuo Kong, Shun Ishii, Kazuhito Dobashi, Tomomi Shimoikura, Patricio Sanhueza, Takashi Tsukagoshi, Paolo Padoan, Ralf S. Klessen, Paul. F. Goldsmith, Blakesley Burkhart, Dariusz C. Lis Álvaro Sánchez-Monge, Yoshito Shimajiri, Ryohei Kawabe

Abstract: 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 asso… ▽ More 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. △ Less

Submitted 21 November, 2022; v1 submitted 18 November, 2022; originally announced November 2022.

Comments: 58 pages, 33 figures, 21 tables, accepted by ApJS

The Dynamical State of Massive Clumps

Authors: Zu-Jia Lu, Veli-Matti Pelkonen, Mika Juvela, Paolo Padoan, Troels Haugbølle, Åke Nordlund

Abstract: The dynamical state of massive clumps is key to our understanding of the formation of massive stars. In this work, we study the kinematic properties of massive clumps using synthetic observations. We have previously compiled a very large catalog of synthetic dust-continuum compact sources from our 250 pc, SN-driven, star formation simulation. Here, we compute synthetic $\rm N_{2}H^{+}$ line profil… ▽ More The dynamical state of massive clumps is key to our understanding of the formation of massive stars. In this work, we study the kinematic properties of massive clumps using synthetic observations. We have previously compiled a very large catalog of synthetic dust-continuum compact sources from our 250 pc, SN-driven, star formation simulation. Here, we compute synthetic $\rm N_{2}H^{+}$ line profiles for a subsample of those sources and compare their properties with the observations and with those of the corresponding three-dimensional (3D) clumps in the simulation. We find that the velocity dispersion of the sources estimated from the $\rm N_{2}H^{+}$ line is a good estimate of that of the 3D clumps, although its correlation with the source size is weaker than the velocity-size correlation of the 3D clumps. The relation between the mass of the 3D clumps, $M_{\rm main}$, and that of the corresponding synthetic sources, $M_{\rm SED}$, has a large scatter and a slope of 0.5, $M_{\rm main} \propto M_{\rm SED}^{0.5}$, due to uncertainties arising from the observational band-merging procedure and from projection effects along the line of sight. As a result, the virial parameters of the 3D clumps are not correlated with the clump masses, even if a negative correlation is found for the compact sources, and the virial parameter of the most massive sources may significantly underestimate that of the associated clumps. △ Less

Submitted 16 November, 2021; originally announced November 2021.

Comments: 19 pages, 12 figures, 1 table. Accepted for publication in MNRAS

Physical Properties and Real Nature of Massive Clumps in the Galaxy

Authors: Zu-Jia Lu, Veli-Matti Pelkonen, Mika Juvela, Paolo Padoan, Troels Haugbølle, Åke Nordlund

Abstract: 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… ▽ More 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 catalog of compact sources following the exact same procedure as for the Hi-GAL compact source catalog. 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 three-dimensional 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. △ Less

Submitted 19 August, 2021; originally announced August 2021.

Comments: 17 pages, 15 figures, submitted to MNRAS

The Core Mass Function in the Orion Nebula Cluster Region: What Determines the Final Stellar Masses?

Authors: Hideaki Takemura, Fumitaka Nakamura, Shuo Kong, Héctor G. Arce, John M. Carpenter, Volker Ossenkopf-Okada, Ralf Klessen, Patricio Sanhueza, Yoshito Shimajiri, Takashi Tsukagoshi, Ryohei Kawabe, Shun Ishii, Kazuhito Dobashi, Tomomi Shimoikura, Paul F. Goldsmith, Álvaro Sánchez-Monge, Jens Kauffmann, Thushara Pillai, Paolo Padoan, Adam Ginsberg, Rowan J. Smith, John Bally, Steve Mairs, Jaime E. Pineda, Dariusz C. Lis , et al. (7 additional authors not shown)

Abstract: Applying dendrogram analysis to the CARMA-NRO C$^{18}$O ($J$=1--0) data having an angular resolution of $\sim$ 8", we identified 692 dense cores in the Orion Nebula Cluster (ONC) 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… ▽ More Applying dendrogram analysis to the CARMA-NRO C$^{18}$O ($J$=1--0) data having an angular resolution of $\sim$ 8", we identified 692 dense cores in the Orion Nebula Cluster (ONC) 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_\odot$, consistent with previous studies. Our CMF has a peak at a subsolar mass of $\sim$ 0.1 $M_\odot$, 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 a 100\% efficiency. Thus, the mass accretion from the surroundings may play a crucial role in determining the final stellar masses of stars. △ Less

Submitted 25 February, 2021; originally announced March 2021.

Comments: 13 pages, 5 figures, accepted by ApJL

From the CMF to the IMF: Beyond the Core-Collapse Model

Authors: V. -M. Pelkonen, P. Padoan, T. Haugbølle, Å. Nordlund

Abstract: Observations have indicated that the prestellar 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 prestellar core, as postulated by the core-collapse model,… ▽ More Observations have indicated that the prestellar 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 prestellar 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 conditions 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%. 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. △ Less

Submitted 14 March, 2021; v1 submitted 5 August, 2020; originally announced August 2020.

Comments: 23 pages, 29 figures. Link to supplementary material and full Table 1: http://www.erda.dk/vgrid/core-mass-function/ . Submitted to MNRAS

The Effect of Supernovae on the Turbulence and Dispersal of Molecular Clouds

Authors: Zu-Jia Lu, Veli-Matti Pelkonen, Paolo Padoan, Liubin Pan, Troels Haugbølle, Åke Nordlund

Abstract: 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 super… ▽ More 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. △ Less

Submitted 18 July, 2020; originally announced July 2020.

Comments: 33 pages, 23 figures, submitted to ApJ

The Origin of Massive Stars: The Inertial-Inflow Model

Authors: Paolo Padoan, Liubin Pan, Mika Juvela, Troels Haugbølle, Åke Nordlund

Abstract: 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 "Inertial-Inflow Model".… ▽ More 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 "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 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 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. △ Less

Submitted 29 July, 2020; v1 submitted 11 November, 2019; originally announced November 2019.

Comments: 35 pages, 30 figures, ApJ in press

The CARMA-NRO Orion Survey: Core Emergence and Kinematics in the Orion A Cloud

Authors: Shuo Kong, Héctor G. Arce, Anneila I. Sargent, Steve Mairs, Ralf S. Klessen, John Bally, Paolo Padoan, Rowan J. Smith, María José Maureira, John M. Carpenter, Adam Ginsburg, Amelia M. Stutz, Paul Goldsmith, Stefan Meingast, Peregrine McGehee, Álvaro Sánchez-Monge, Sümeyye Suri, Jaime E. Pineda, João Alves, Jesse R. Feddersen, Jens Kauffmann, Peter Schilke

Abstract: 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 extinctio… ▽ More 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. △ Less

Submitted 13 August, 2019; originally announced August 2019.

Comments: 21 pages, 11 figures, accepted by ApJ, comments welcome shuo.kong@yale.edu

Synthetic observations of dust emission and polarisation of Galactic cold clumps

Authors: Mika Juvela, Paolo Padoan, Isabelle Ristorcelli, Veli-Matti Pelkonen

Abstract: 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… ▽ More 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 Δ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. △ Less

Submitted 9 August, 2019; originally announced August 2019.

Comments: Accepted to A&A (20 pages, 26 figures)

The Probability Distribution of Density Fluctuations in Supersonic Turbulence

Authors: Liubin Pan, Paolo Padoan, Åke Nordlund

Abstract: 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 E… ▽ More 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, $σ_{s}$, follows the relation $σ_{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 $σ_{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. △ Less

Submitted 2 May, 2019; originally announced May 2019.

Comments: Submitted to ApJ

Inaccuracy of Spatial Derivatives in Simulations of Supersonic Turbulence

Authors: Liubin Pan, Paolo Padoan, Åke Nordlund

Abstract: 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, includin… ▽ More 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. △ Less

Submitted 31 January, 2019; originally announced February 2019.

Comments: Submitted to ApJ

The CARMA-NRO Orion Survey: The filamentary structure as seen in C$^{18}$O emission

Authors: S. T. Suri, A. Sanchez-Monge, P. Schilke, S. D. Clarke, R. J. Smith, V. Ossenkopf-Okada, R. Klessen, P. Padoan, P. Goldsmith, H. G. Arce, J. Bally, J. M. Carpenter, A. Ginsburg, D. Johnstone, J. Kauffmann, S. Kong, D. C. Lis, S. Mairs, T. Pillai, J. E. Pineda, A. Duarte-Cabral

Abstract: 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.… ▽ More 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. △ Less

Submitted 1 January, 2019; originally announced January 2019.

Comments: 19 pages, 20 figures. Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 623, A142 (2019)

Expanding associations in the Vela-Puppis region: 3D structure and kinematics of the young population

Authors: T. Cantat-Gaudin, C. Jordi, N. J. Wright, J. J. Armstrong, A. Vallenari, L. Balaguer-Núñez, P. Ramos, D. Bossini, P. Padoan, V. M. Pelkonen, M. Mapelli, R. D. Jeffries

Abstract: 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 Ga… ▽ More 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. △ Less

Submitted 7 May, 2019; v1 submitted 19 December, 2018; originally announced December 2018.

Comments: accepted A&A

Journal ref: A&A 626, A17 (2019)

Dust polarisation studies on MHD simulations of molecular clouds: methods comparison for the relative orientations analysis

Authors: Elisabetta R. Micelotta, Mika Juvela, Paolo Padoan, Isabelle Ristorcelli, Dana Alina, Johanna Malinen

Abstract: 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… ▽ More 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. △ Less

Submitted 23 October, 2018; originally announced October 2018.

Comments: 15 pages, 9 figures, submitted to A&A, comments welcomed

Journal ref: A&A 647, A121 (2021)

Detailed Balance and Exact Results for Density Fluctuations in Supersonic Turbulence

Authors: Liubin Pan, Paolo Padoan, Åke Nordlund

Abstract: 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,… ▽ More 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. △ Less

Submitted 28 September, 2018; v1 submitted 24 August, 2018; originally announced August 2018.

Comments: Accepted for publication in ApJL

The CARMA-NRO Orion Survey

Authors: Shuo Kong, Héctor G. Arce, Jesse R. Feddersen, John M. Carpenter, Fumitaka Nakamura, Yoshito Shimajiri, Andrea Isella, Volker Ossenkopf-Okada, Anneila I. Sargent, Álvaro Sánchez-Monge, Sümeyye T. Suri, Jens Kauffmann, Thushara Pillai, Jaime E. Pineda, Jin Koda, John Bally, Dariusz C. Lis, Paolo Padoan, Ralf Klessen, Steve Mairs, Alyssa Goodman, Paul Goldsmith, Peregrine McGehee, Peter Schilke, Peter J. Teuben , et al. (13 additional authors not shown)

Abstract: 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… ▽ More 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 $Δ$-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 △ Less

Submitted 30 March, 2018; originally announced March 2018.

Comments: 46 pages, 28 figures, 2 tables, accepted by ApJS

The magnetic field of molecular clouds

Authors: Paolo Padoan

Abstract: 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 ma… ▽ More 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énic, and we argue that MCs are bound to be born super-Alfvé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. △ Less

Submitted 9 December, 2017; originally announced December 2017.

Comments: 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)

The stellar IMF from Isothermal MHD Turbulence

Authors: Troels Haugbølle, Paolo Padoan, Aake Nordlund

Abstract: 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 focu… ▽ More 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. △ Less

Submitted 31 January, 2018; v1 submitted 4 September, 2017; originally announced September 2017.

Comments: 25 pages, 21 figures, Accepted by ApJ

Supernova Driving. IV. The Star Formation Rate of Molecular Clouds

Authors: Paolo Padoan, Troels Haugbølle, Åke Nordlund, Søren Frimann

Abstract: 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 sta… ▽ More 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, $α_{\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, $β$, and compressive to solenoidal power ratio, $χ$ at fixed $α_{\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}$--$α_{\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. △ Less

Submitted 31 March, 2017; v1 submitted 23 February, 2017; originally announced February 2017.

Comments: ApJ, in press

Supernova Driving. III. Synthetic Molecular Cloud Observations

Authors: Paolo Padoan, Mika Juvela, Liubin Pan, Troels Haugbølle, Åke Nordlund

Abstract: 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 br… ▽ More 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. △ Less

Submitted 12 May, 2016; originally announced May 2016.

Comments: accepted for publication in ApJ

Protostellar accretion traced with chemistry: Comparing synthetic C18O maps of embedded protostars to real observations

Authors: Søren Frimann, Jes K. Jørgensen, Paolo Padoan, Troels Haugbølle

Abstract: 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… ▽ More 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. △ Less

Submitted 19 December, 2015; v1 submitted 1 December, 2015; originally announced December 2015.

Comments: Accepted for publication in A&A, 11 pages, 8 figures; v2 contains minor updates to the language

Journal ref: A&A 587, A60 (2016)

Supernova Driving. II. Compressive Ratio in Molecular-Cloud Turbulence

Authors: Liubin Pan, Paolo Padoan, Troels Haugbolle, Aake Nordlund

Abstract: 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 e… ▽ More 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. △ Less

Submitted 19 May, 2016; v1 submitted 15 October, 2015; originally announced October 2015.

Comments: Accepted by ApJ

Supernova Driving. I. The Origin of Molecular Cloud Turbulence

Authors: Paolo Padoan, Liubin Pan, Troels Haugboelle, Ake Nordlund

Abstract: 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,… ▽ More 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é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. △ Less

Submitted 2 March, 2016; v1 submitted 15 September, 2015; originally announced September 2015.

Comments: ApJ, in press

The relative and absolute ages of old globular clusters in the LCDM framework

Authors: Michele Trenti, Paolo Padoan, Raul Jimenez

Abstract: 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 coo… ▽ More 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. △ Less

Submitted 7 July, 2015; v1 submitted 9 February, 2015; originally announced February 2015.

Comments: ApJL accepted, minor content changes to highlight the robust predictions for GC ages; for an animated version of Fig. 1 (minihalo mergers movie) see http://www.ph.unimelb.edu.au/~mtrenti/trenti_padoan_jimenez_2015_fig1_GC_movie.gif

Journal ref: The Astrophysical Journal Letters. 808, L35, 2015

Turbulence-Induced Relative Velocity of Dust particles IV: the Collision Kernel

Authors: Liubin Pan, Paolo Padoan

Abstract: 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, t… ▽ More 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. △ Less

Submitted 8 October, 2014; originally announced October 2014.

Comments: 20 pages, 9 figures, accepted to ApJ

Infall-Driven Protostellar Accretion and the Solution to the Luminosity Problem

Authors: Paolo Padoan, Troels Haugbølle, Åke Nordlund

Abstract: 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 p… ▽ More 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. △ Less

Submitted 7 October, 2014; v1 submitted 5 July, 2014; originally announced July 2014.

Comments: 21 pages, 16 figures, accepted for publication in ApJ

Turbulence-Induced Relative Velocity of Dust Particles III: The Probability Distribution

Authors: Liubin Pan, Paolo Padoan, John Scalo

Abstract: 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… ▽ More 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. △ Less

Submitted 3 July, 2014; originally announced July 2014.

Comments: 25 pages, 19 figures, accepted by ApJ

Turbulence-Induced Relative Velocity of Dust Particles II: The Bidisperse Case

Authors: Liubin Pan, Paolo Padoan, John Scalo

Abstract: 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… ▽ More 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. △ Less

Submitted 3 July, 2014; v1 submitted 15 March, 2014; originally announced March 2014.

Comments: 23 pages, 12 figures, Accepted to ApJ

The Star Formation Rate of Molecular Clouds

Authors: Paolo Padoan, Christoph Federrath, Gilles Chabrier, Neal J. Evans II, Doug Johnstone, Jes K. Jørgensen, Christopher F. McKee, Åke Nordlund

Abstract: 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 supers… ▽ More 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. △ Less

Submitted 18 December, 2013; originally announced December 2013.

Comments: 24 pages, 5 figures, Accepted for publication as a chapter in Protostars and Planets VI, University of Arizona Press (2014), eds. H. Beuther, R. S. Klessen, C. P. Dullemond, Th. Henning

Zooming in on the Formation of Protoplanetary Disks

Authors: A. Nordlund, T. Haugboelle, M. Kuffmeier, P. Padoan, A. Vasileiades

Abstract: 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 particl… ▽ More 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. △ Less

Submitted 17 September, 2013; v1 submitted 9 September, 2013; originally announced September 2013.

Comments: 4 pages, 3 figures, slightly edited / corrected version of contribution to IAU Symposium 299, Victoria

Turbulence-Induced Relative Velocity of Dust Particles I: Identical Particles

Authors: Liubin Pan, Paolo Padoan

Abstract: 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… ▽ More 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. △ Less

Submitted 9 August, 2013; v1 submitted 1 May, 2013; originally announced May 2013.

Comments: 41 pages, major revisions, accepted to ApJ

A Simple Law of Star Formation

Authors: Paolo Padoan, Troels Haugboelle, AAke Nordlund

Abstract: 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é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 t… ▽ More 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é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, $ε_{\rm ff}$, on the above parameters. Our main results are: i) $ε_{\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 $ε_{\rm ff}$, but only to a point, beyond which $ε_{\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, $ε_{\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 $ε_{\rm ff}$, and depending only on $t_{\rm ff}/t_{\rm dyn}$: $ε_{\rm ff} \approx ε_{\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. △ Less

Submitted 27 September, 2012; v1 submitted 18 August, 2012; originally announced August 2012.

Comments: ApJ Letters - in press

The Two States of Star Forming Clouds

Authors: David Collins, Alexei G. Kritsuk, Paolo Padoan, Hui Li, Hao Xu, Sergey D. Ustyugov, Michael L. Norman

Abstract: 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énic… ▽ More 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énic and one trans-Alfvénic. 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(ρ,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én 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. △ Less

Submitted 12 February, 2012; originally announced February 2012.

Comments: 19 pages, 20 figures

The observable prestellar phase of the IMF

Authors: Paolo Padoan, AAke Nordlund

Abstract: 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 Pa… ▽ More 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. △ Less

Submitted 13 October, 2011; v1 submitted 11 August, 2011; originally announced August 2011.

Comments: ApJ Letters, accepted

Turbulent Clustering of Protoplanetary Dust and Planetesimal Formation

Authors: Liubin Pan, Paolo Padoan, John Scalo, Alexei G. Kritsuk, Michael L. Norman

Abstract: 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 St… ▽ More 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, $τ_{\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 $τ_{\rm p}$ in the inertial range. Clustering of these particles occurs primarily around a scale where the eddy turnover time is $\simτ_{\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. △ Less

Submitted 18 June, 2011; originally announced June 2011.

Comments: 24 pages, 13 figures, accepted by ApJ

Comparing Numerical Methods for Isothermal Magnetized Supersonic Turbulence

Authors: Alexei G. Kritsuk, Aake Nordlund, David Collins, Paolo Padoan, Michael L. Norman, Tom Abel, Robi Banerjee, Christoph Federrath, Mario Flock, Dongwook Lee, Pak Shing Li, Wolf-Christian Mueller, Romain Teyssier, Sergey D. Ustyugov, Christian Vogel, Hao Xu

Abstract: 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 n… ▽ More 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. △ Less

Submitted 7 June, 2011; v1 submitted 28 March, 2011; originally announced March 2011.

Comments: 17 pages, 5 color figures, revised version to appear in ApJ, 735, July 2011

Report number: NSF-KITP-11-079

Journal ref: Astrophysical Journal, 737, 13 (2011)

Pop. III stars from turbulent fragmentation at redshift ~ 11

Authors: Joaquin Prieto, Paolo Padoan, Raul Jimenez, Leopoldo Infante

Abstract: 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 t… ▽ More 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. △ Less

Submitted 26 January, 2011; originally announced January 2011.

Journal ref: ApJL, 731, 38 (2011)

Accuracy of core mass estimates in simulated observations of dust emission

Authors: J. Malinen, M. Juvela, D. C. Collins, T. Lunttila, P. Padoan

Abstract: 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 c… ▽ More 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. △ Less

Submitted 18 April, 2011; v1 submitted 23 September, 2010; originally announced September 2010.

Comments: Revised version, 17 pages, 17 figures, submitted to A&A

Journal ref: A&A, 530 (2011), A101

Mass and Magnetic distributions in Self Gravitating Super Alfvenic Turbulence with AMR

Authors: David C. Collins, Paolo Padoan, Michael L. Norman, Hao Xu

Abstract: 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(ρ) \propto ρ^{-1.67}$ and $P(B)\propto B^{-2.74}$. A power law is also found between magnetic field strength and density… ▽ More 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(ρ) \propto ρ^{-1.67}$ and $P(B)\propto B^{-2.74}$. A power law is also found between magnetic field strength and density, with $B\propto ρ^{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. △ Less

Submitted 19 August, 2010; v1 submitted 13 August, 2010; originally announced August 2010.

Comments: 9 pages, 10 figures (3 color). Submitted to the Astrophysical Journal

Relative Velocity of Inertial Particles in Turbulent Flows

Authors: Liubin Pan, Paolo Padoan

Abstract: 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' memo… ▽ More 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. △ Less

Submitted 13 May, 2010; originally announced May 2010.

Comments: Accepted by J. Fluid. Mech. 35 pages

Self-organization in Turbulent Molecular Clouds: Compressional versus Solenoidal Modes

Authors: Alexei G. Kritsuk, Sergey D. Ustyugov, Michael L. Norman, Paolo Padoan

Abstract: 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 (dil… ▽ More 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. △ Less

Submitted 3 December, 2009; originally announced December 2009.

Comments: 7 pages, 1 figure; ASTRONUM-2009 "Numerical Modeling of Space Plasma Flows", Chamonix, France, July 2009, ASP Conf. Proc., in press

Journal ref: ASP Conf.Ser.429:15-21,2010

The Power Spectrum of Turbulence in NGC 1333: Outflows or Large-Scale Driving?

Authors: Paolo Padoan, Mika Juvela, Alexei G. Kritsuk, Michael L. Norman

Abstract: 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… ▽ More 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. △ Less

Submitted 18 November, 2009; v1 submitted 7 October, 2009; originally announced October 2009.

Comments: Submitted to APJ Letters on September 22, 2009 - Accepted on November 18, 2009

Journal ref: Astrophys.J.707:L153-L157,2009

The IMF of stellar clusters: effects of accretion and feedback

Authors: Sami Dib, Mohsen Shadmehri, Paolo Padoan, G. Maheswar, D. K. Ojha, Fazeleh Khajenabi

Abstract: (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… ▽ More (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. △ Less

Submitted 30 January, 2010; v1 submitted 31 August, 2009; originally announced August 2009.

Comments: accepted to MNRAS. 22 pages, 17 figures. 2 new sections added to better present the model parameters and compare to previous work. Some figures are now better presented. Main conclusions unchanged

Simulating Supersonic Turbulence in Magnetized Molecular Clouds

Authors: Alexei G. Kritsuk, Sergey D. Ustyugov, Michael L. Norman, Paolo Padoan

Abstract: 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 fu… ▽ More 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. △ Less

Submitted 3 August, 2009; originally announced August 2009.

Comments: 10 pages, 3 figures, to appear in the proceedings of the DOE/SciDAC 2009 conference

Journal ref: Journal of Physics: Conf. Ser., 180, 012020 (2009)

The Super-Alfvénic Model of Molecular Clouds: Predictions for Mass-to-Flux and Turbulent-to-Magnetic Energy Ratios

Authors: Tuomas Lunttila, Paolo Padoan, Mika Juvela, Åke Nordlund

Abstract: 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énic turbulence. We have previously shown that the same simulation yields a rela… ▽ More 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énic 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. △ Less

Submitted 24 July, 2009; v1 submitted 3 July, 2009; originally announced July 2009.

Comments: 7/24/09: 6 pages, 7 figures, accepted to ApJL. Minor additions to the previous version. 7/3/09: 5 pages, 7 figures, submitted to ApJL

Journal ref: Astrophys.J.702:L37-L41,2009