Source author record

Paolo Padoan

Paolo Padoan appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.GA astro-ph astro-ph.EP astro-ph.SR physics.flu-dyn astro-ph.CO nlin.CD

Catalog footprint

What is connected

36works

7topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2020arXiv

The Origin of Massive Stars: The Inertial-Inflow Model

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.

preprint2016arXiv

Supernova Driving. I. The Origin of Molecular Cloud Turbulence

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.

preprint2016arXiv

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

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.

preprint2016arXiv

Supernova Driving. III. Synthetic Molecular Cloud Observations

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.

preprint2015arXiv

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

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.

preprint2015arXiv

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

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.

preprint2014arXiv

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

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.

preprint2014arXiv

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

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.

preprint2014arXiv

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

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.

preprint2014arXiv

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

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.

preprint2013arXiv

The Star Formation Rate of Molecular Clouds

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.

preprint2013arXiv

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

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.

preprint2012arXiv

A Simple Law of Star Formation

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.

preprint2012arXiv

The Two States of Star Forming Clouds

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.

preprint2011arXiv

Comparing Numerical Methods for Isothermal Magnetized Supersonic Turbulence

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.

preprint2011arXiv

Pop. III stars from turbulent fragmentation at redshift ~ 11

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.

preprint2011arXiv

The observable prestellar phase of the IMF

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.

preprint2011arXiv

Turbulent Clustering of Protoplanetary Dust and Planetesimal Formation

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.

preprint2010arXiv

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

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.

preprint2010arXiv

Relative Velocity of Inertial Particles in Turbulent Flows

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.

preprint2010arXiv

The IMF of stellar clusters: effects of accretion and feedback

(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.

preprint2010arXiv

The Star Formation Rate of Supersonic MHD Turbulence

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

preprint2009arXiv

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

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.

preprint2009arXiv

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

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.

preprint2009arXiv

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

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.

preprint2008arXiv

The Super-Alfvénic Model of Molecular Clouds: Predictions for Zeeman Splitting Measurements

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

preprint2008arXiv

The Temperature of Interstellar Clouds from Turbulent Heating

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

preprint2007arXiv

Two regimes of Turbulent Fragmentation and the stellar IMF from Primordial to Present Day Star Formation

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

preprint2005arXiv

Adaptive Mesh Refinement for Supersonic Molecular Cloud Turbulence

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

preprint2004arXiv

The Stellar IMF as a Property of Turbulence

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

preprint1999arXiv

Super-Sonic Turbulence in the Perseus Molecular Cloud

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

preprint1997arXiv

The Universality of the Stellar IMF

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

preprint1996arXiv

A new method to determine Globular Cluster ages

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

preprint1996arXiv

Ages of globular clusters: breaking the age-distance degeneracy with the luminosity function

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

preprint1996arXiv

On star formation in primordial protoglobular clouds

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

preprint1995arXiv

Supersonic Turbulence and the Fragmentation of a Cold Medium

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

Paolo Padoan

What is connected

Connect this record

See the researcher in context

Building this map preview

36 published item(s)

The Origin of Massive Stars: The Inertial-Inflow Model

Supernova Driving. I. The Origin of Molecular Cloud Turbulence

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

Supernova Driving. III. Synthetic Molecular Cloud Observations

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

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

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

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

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

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

The Star Formation Rate of Molecular Clouds

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

A Simple Law of Star Formation

The Two States of Star Forming Clouds

Comparing Numerical Methods for Isothermal Magnetized Supersonic Turbulence

Pop. III stars from turbulent fragmentation at redshift ~ 11

The observable prestellar phase of the IMF

Turbulent Clustering of Protoplanetary Dust and Planetesimal Formation

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

Relative Velocity of Inertial Particles in Turbulent Flows

The IMF of stellar clusters: effects of accretion and feedback

The Star Formation Rate of Supersonic MHD Turbulence

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

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

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

The Super-Alfvénic Model of Molecular Clouds: Predictions for Zeeman Splitting Measurements

The Temperature of Interstellar Clouds from Turbulent Heating

Two regimes of Turbulent Fragmentation and the stellar IMF from Primordial to Present Day Star Formation

Adaptive Mesh Refinement for Supersonic Molecular Cloud Turbulence

The Stellar IMF as a Property of Turbulence

Super-Sonic Turbulence in the Perseus Molecular Cloud

The Universality of the Stellar IMF

A new method to determine Globular Cluster ages

Ages of globular clusters: breaking the age-distance degeneracy with the luminosity function

On star formation in primordial protoglobular clouds

Supersonic Turbulence and the Fragmentation of a Cold Medium