-
Inefficient star formation in high Mach number environments II. Numerical simulations and comparison with analytical models
Authors:
Noé Brucy,
Patrick Hennebelle,
Tine Colman,
Ralf S. Klessen,
Corentin Le Yhuelic
Abstract:
Predicting the star formation rate (SFR) in galaxies is crucial to understand their evolution and morphology. To do so requires a fine understanding of how dense structures of gas are created and collapse. In that, turbulence and gravity play a major role. Within the gravo-turbulent framework, we assume that turbulence shapes the ISM, creating density fluctuations that, if gravitationally unstable…
▽ More
Predicting the star formation rate (SFR) in galaxies is crucial to understand their evolution and morphology. To do so requires a fine understanding of how dense structures of gas are created and collapse. In that, turbulence and gravity play a major role. Within the gravo-turbulent framework, we assume that turbulence shapes the ISM, creating density fluctuations that, if gravitationally unstable, will collapse and form stars. The goal of this work is to quantify how different regimes of turbulence, characterized by the strength and compressibility of the driving, shape the density field. We are interested in the outcome in terms of SFR and how it compares with existing analytical models for the SFR. We run a series of hydrodynamical simulations of turbulent gas. The simulations are first conducted without gravity, so that the density and velocity are shaped by the turbulence driving. Gravity is then switched on, and the SFR is measured and compared with analytical models. The physics included in these simulations is very close to the one assumed in the classical gravo-turbulent SFR analytical models, which makes the comparison straightforward. We found that the existing analytical models convincingly agree with simulations at low Mach number, but we measure a much lower SFR in the simulation with a high Mach number. We develop, in a companion paper, an updated physically-motivated SFR model that reproduces well the inefficient high Mach regime of the simulations. Our work demonstrates that accurate estimations of the turbulent-driven replenishment time of dense structures and the dense gas spatial distribution are necessary to correctly predict the SFR in the high Mach regime. The inefficient high-Mach regime is a possible explanation for the low SFR found in dense and turbulent environments such as the centers of our Milky Way and other galaxies.
△ Less
Submitted 25 July, 2024; v1 submitted 26 April, 2024;
originally announced April 2024.
-
Inefficient star formation in high Mach number environments I. The turbulent support analytical model
Authors:
Patrick Hennebelle,
Noé Brucy,
Tine Colman
Abstract:
The star formation rate (SFR), the number of stars formed per unit of time, is a fundamental quantity in the evolution of the Universe. While turbulence is believed to play a crucial role in setting the SFR, the exact mechanism remains unclear. Turbulence promotes star formation by compressing the gas, but also slows it down by stabilizing the gas against gravity. Most widely-used analytical model…
▽ More
The star formation rate (SFR), the number of stars formed per unit of time, is a fundamental quantity in the evolution of the Universe. While turbulence is believed to play a crucial role in setting the SFR, the exact mechanism remains unclear. Turbulence promotes star formation by compressing the gas, but also slows it down by stabilizing the gas against gravity. Most widely-used analytical models rely on questionable assumptions, including: $i)$ integrating over the density PDF, a one-point statistical description that ignores spatial correlation, $ii)$ selecting self-gravitating gas based on a density threshold that often ignores turbulent dispersion, $iii)$ assuming the freefall time as the timescale for estimating SFR without considering the need to rejuvenate the density PDF, $iv)$ assuming the density PDF to be lognormal. Improving upon the only existing model that incorporates the spatial correlation of the density field, we present a new analytical model. We calculate the time needed to rejuvenate density fluctuations of a given density and spatial scale, revealing that it is generally much longer than the freefall time, rendering the latter inappropriate for use. We make specific predictions regarding the role of the Mach number, $ M $, and the driving scale of turbulence divided by the mean Jeans length. At low to moderate Mach numbers, turbulence does not reduce and may even slightly promote star formation by broadening the PDF. However, at higher Mach numbers, most density fluctuations are stabilized by turbulent dispersion, leading to a steep drop in the SFR as the Mach number increases.
△ Less
Submitted 25 July, 2024; v1 submitted 26 April, 2024;
originally announced April 2024.
-
Cloud properties across spatial scales in simulations of the interstellar medium
Authors:
Tine Colman,
Noé Brucy,
Philipp Girichidis,
Simon C. O Glover,
Milena Benedettini,
Juan D. Soler,
Robin G. Tress,
Alessio Traficante,
Patrick Hennebelle,
Ralf S. Klessen,
Sergio Molinari,
Marc-Antoine Miville-Deschênes
Abstract:
Molecular clouds (MC) are structures of dense gas in the interstellar medium (ISM), that extend from ten to a few hundred parsecs and form the main gas reservoir available for star formation. Hydrodynamical simulations of varying complexity are a promising way to investigate MC evolution and their properties. However, each simulation typically has a limited range in resolution and different cloud…
▽ More
Molecular clouds (MC) are structures of dense gas in the interstellar medium (ISM), that extend from ten to a few hundred parsecs and form the main gas reservoir available for star formation. Hydrodynamical simulations of varying complexity are a promising way to investigate MC evolution and their properties. However, each simulation typically has a limited range in resolution and different cloud extraction algorithms are used, which complicates the comparison between simulations. In this work, we aim to extract clouds from different simulations covering a wide range of spatial scales. We compare their properties, such as size, shape, mass, internal velocity dispersion and virial state. We apply the Hop cloud detection algorithm on (M)HD numerical simulations of stratified ISM boxes and isolated galactic disk simulations that were produced using Flash Ramses and Arepo We find that the extracted clouds are complex in shape ranging from round objects to complex filamentary networks in all setups. Despite the wide range of scales, resolution, and sub-grid physics, we observe surprisingly robust trends in the investigated metrics. The mass spectrum matches in the overlap between simulations without rescaling and with a high-mass slope of $\mathrm{d} N/\mathrm{d}\ln M\propto-1$ in accordance with theoretical predictions. The internal velocity dispersion scales with the size of the cloud as $σ\propto R^{0.75}$ for large clouds ($R\gtrsim3\,\mathrm{pc}$). For small clouds we find larger sigma compared to the power-law scaling, as seen in observations, which is due to supernova-driven turbulence. Almost all clouds are gravitationally unbound with the virial parameter scaling as $α_\mathrm{vir}\propto M^{-0.4}$, which is slightly flatter compared to observed scaling, but in agreement given the large scatter.
△ Less
Submitted 1 March, 2024;
originally announced March 2024.
-
Stellar Feedback in the Star Formation-Gas Density Relation: Comparison between Simulations and Observations
Authors:
Paolo Suin,
Annie Zavagno,
Tine Colman,
Patrick Hennebelle,
Antoine Verliat,
Delphine Russeil
Abstract:
Context. The impact of stellar feedback on the Kennicutt-Schmidt law (KS law), which relates star formation rate (SFR) to surface gas density, is a topic of ongoing debate. The interpretation of individual cloud observations is challenging due to the various processes at play simultaneously and inherent biases. Therefore, a numerical investigation is necessary to understand the role of stellar fee…
▽ More
Context. The impact of stellar feedback on the Kennicutt-Schmidt law (KS law), which relates star formation rate (SFR) to surface gas density, is a topic of ongoing debate. The interpretation of individual cloud observations is challenging due to the various processes at play simultaneously and inherent biases. Therefore, a numerical investigation is necessary to understand the role of stellar feedback and identify observable signatures.
Aims. We investigate the role of stellar feedback on the KS law, aiming to identify distinct signatures that can be observed and analysed.
Methods. We analyse MHD numerical simulations of a $10^4\,M_{\odot}$ cloud evolving under different feedback prescriptions. The set of simulations contains four types of feedback: with only protostellar jets, with ionising radiation from massive stars $(>8\,M_{\odot})$, with both of them and without any stellar feedback. To compare these simulations with the existing observational results, we analyse their evolution by adopting the same techniques applied in observational studies. Then, we simulate how the same analyses would change if the data were affected by typical observational biases.
Conclusions. The presence of stellar feedback strongly influences the KS relation and the star formation efficiency per free-fall time ($ε_\mathrm{ff}$). Its impact is primarily governed by its influence on the cloud's structure. Although the $ε_\mathrm{ff}$ measured in our clouds results to be higher than what is usually observed in real clouds, upon applying prescriptions to mimic observational biases we recover good agreement with the expected values. Therefore, we can infer that observations tend to underestimate the total SFR. Moreover, this likely indicates that the physics included in our simulations is sufficient to reproduce the basic mechanisms contributing to set $ε_\mathrm{ff}$.
△ Less
Submitted 30 November, 2023;
originally announced November 2023.
-
Synthetic populations of protoplanetary disks. Impact of magnetic fields and radiative transfer
Authors:
U. Lebreuilly,
P. Hennebelle,
T. Colman,
A. Maury,
N. -D. Tung,
L. Testi,
R. Klessen,
S. Molinari,
B. Commerçon,
M. González,
E. Pacetti,
A. Somigliana,
G. Rosotti
Abstract:
Protostellar disks are the product of angular momentum conservation during the protostellar collapse. Understanding their formation is crucial because they are the birthplace of planets and because their formation is tightly related to star formation. Unfortunately, the initial properties of Class 0 disks and their evolution are still poorly constrained observationally and theoretically. We aim to…
▽ More
Protostellar disks are the product of angular momentum conservation during the protostellar collapse. Understanding their formation is crucial because they are the birthplace of planets and because their formation is tightly related to star formation. Unfortunately, the initial properties of Class 0 disks and their evolution are still poorly constrained observationally and theoretically. We aim to better understand the mechanisms that set the statistics of disk properties as well as to study their formation in massive protostellar clumps. We also want to provide the community with synthetic disk populations to better interpret young disk observations. We use the ramses code to model star and disk formation in massive protostellar clumps with MHD including the effect of ambipolar diffusion and RT including the stellar radiative feedback. Those simulations, resolved up to the astronomical unit scale, allow to investigate the formation of disk populations. Magnetic fields play a crucial role in disk formation. A weaker initial field leads to larger and massive disks and weakens the stellar radiative feedback by increasing fragmentation. We find that ambipolar diffusion impacts disk and star formation and leads to very different disk magnetic properties. The stellar radiative feedback also have a strong influence, increasing the temperature and reducing fragmentation. Comparing our disk populations with observations reveals that our models with a mass-to-flux ratio of 10 seems to better reproduce observed disk sizes. This also sheds light on a tension between models and observations for the disk masses. The clump properties and physical modeling impact disk populations significantly. The tension between observations and models for disk mass estimates is critical to solve with synthetic observations in future years, in particular for our comprehension of planet formation.
△ Less
Submitted 30 October, 2023;
originally announced October 2023.
-
A comparison of the Milky Way's recent star formation revealed by dust thermal emission and high-mass stars
Authors:
J. D. Soler,
E. Zari,
D. Elia,
S. Molinari,
C. Mininni,
E. Schisano,
A. Traficante,
R. S. Klessen,
S. C. O. Glover,
P. Hennebelle,
T. Colman,
N. Frankel,
T. Wenger
Abstract:
We present a comparison of the Milky Way's star formation rate (SFR) surface density ($Σ_{\rm SFR}$) obtained with two independent state-of-the-art observational methods. The first method infers $Σ_{\rm SFR}$ from observations of the dust thermal emission from interstellar dust grains in far-infrared wavelengths registered in the Herschel infrared Galactic Plane Survey (Hi-GAL). The second method…
▽ More
We present a comparison of the Milky Way's star formation rate (SFR) surface density ($Σ_{\rm SFR}$) obtained with two independent state-of-the-art observational methods. The first method infers $Σ_{\rm SFR}$ from observations of the dust thermal emission from interstellar dust grains in far-infrared wavelengths registered in the Herschel infrared Galactic Plane Survey (Hi-GAL). The second method determines $Σ_{\rm SFR}$ by modeling the current population of O-, B-, and A-type stars in a 6 kpc $\times$ 6 kpc area around the Sun. We find an agreement between the two methods within a factor of two for the mean SFRs and the SFR surface density profiles. Given the broad differences between the observational techniques and the independent assumptions in the methods for computing the SFRs, this agreement constitutes a significant advance in our understanding of the star formation of our Galaxy and implies that the local SFR has been roughly constant over the past 10\,Myr.
△ Less
Submitted 18 September, 2023; v1 submitted 2 August, 2023;
originally announced August 2023.
-
Large-scale turbulent driving regulates star formation in high-redshift gas-rich galaxies II: Influence of the magnetic field and the turbulent compressive fraction
Authors:
Noé Brucy,
Patrick Hennebelle,
Tine Colman,
Simon Iteanu
Abstract:
The observed star formation rate (SFR) in galaxies is well below what it should be if gravitational collapse alone were at play. It has recently been shown that one candidate that might regulate star formation, the feedback from massive stars, is suitable only if the mean column density at the kiloparsec scale is lower than $\approx 20 M_\odot\cdot\mathrm{pc}^{-2}$. On the other hand, intense larg…
▽ More
The observed star formation rate (SFR) in galaxies is well below what it should be if gravitational collapse alone were at play. It has recently been shown that one candidate that might regulate star formation, the feedback from massive stars, is suitable only if the mean column density at the kiloparsec scale is lower than $\approx 20 M_\odot\cdot\mathrm{pc}^{-2}$. On the other hand, intense large-scale turbulent driving might slow down star formation in high-density environments to values that are compatible with observations. In this work, we explore the effect of the nature and strength of the turbulent driving, as well as the effect of the magnetic field. We performed a large series of feedback-regulated numerical simulations of the interstellar medium in which bidimensional large-scale turbulent driving was also applied. We determined the driving intensity needed to reproduce the Schmidt-Kennicutt relation for several gas column densities, magnetization, and driving compressibility. We confirm that in the absence of turbulent forcing and even with a substantial magnetic field, the SFR is too high, particularly at a high column density, compared to the Schmidt-Kennicutt relation. We find that the SFR outcome strongly depends on the initial magnetic field and on the compressibility of the turbulent driving. As a consequence, a higher magnetic field in high column density environment may lower the energy necessary to sustain a turbulence that is sufficiently intense to regulate star formation. Stellar feedback does not seem to be sufficient to regulate star formation in gas-rich galaxies where large-scale turbulent driving may be needed. The sources of this large-scale turbulence as well as its characteristics, such as its intensity, compressibility, and anisotropy, need to be understood and quantified.
△ Less
Submitted 6 June, 2023; v1 submitted 29 May, 2023;
originally announced May 2023.
-
A panoptic view of the Taurus molecular cloud I. The cloud dynamics revealed by gas emission and 3D dust
Authors:
J. D. Soler,
C. Zucker,
J. E. G. Peek,
M. Heyer,
P. F. Goldsmith,
S. C. O. Glover,
S. Molinari,
R. S. Klessen,
P. Hennebelle,
L. Testi,
T. Colman,
M. Benedettini,
D. Elia,
C. Mininni,
S. Pezzuto,
E. Schisano,
A. Traficante
Abstract:
We present a study of the three-dimensional (3D) distribution of interstellar dust derived from stellar extinction observations toward the Taurus molecular cloud (MC) and its relation with the neutral atomic hydrogen (HI) emission at 21 cm wavelength and the carbon monoxide $^{12}$CO and $^{13}$CO emission in the $J=1\rightarrow0$ transition. We used the histogram of oriented gradients (HOG) metho…
▽ More
We present a study of the three-dimensional (3D) distribution of interstellar dust derived from stellar extinction observations toward the Taurus molecular cloud (MC) and its relation with the neutral atomic hydrogen (HI) emission at 21 cm wavelength and the carbon monoxide $^{12}$CO and $^{13}$CO emission in the $J=1\rightarrow0$ transition. We used the histogram of oriented gradients (HOG) method to match the morphology in a 3D reconstruction of the dust density (3D dust) and the distribution of the gas tracers' emission. The result of the HOG analysis is a map of the relationship between the distances and radial velocities. The HOG comparison between the 3D dust and the HI emission indicates a morphological match at the distance of Taurus but an anti-correlation between the dust density and the HI emission, which uncovers a significant amount of cold HI within the Taurus MC. The HOG between the 3D dust and $^{12}$CO reveals a pattern in radial velocities and distances that is consistent with converging motions of the gas in the Taurus MC, with the near side of the cloud moving at higher velocities and the far side moving at lower velocities. This convergence of flows is likely triggered by the large-scale gas compression caused by the interaction of the Local Bubble and the Per-Tau shell, with Taurus lying at the intersection of the two bubble surfaces.
△ Less
Submitted 22 May, 2023;
originally announced May 2023.
-
The Star Formation Rate of the Milky Way as seen by Herschel
Authors:
D. Elia,
S. Molinari,
E. Schisano,
J. D. Soler,
M. Merello,
D. Russeil,
M. Veneziani,
A. Zavagno,
A. Noriega-Crespo,
L. Olmi,
M. Benedettini,
P. Hennebelle,
R. S. Klessen,
S. Leurini,
R. Paladini,
S. Pezzuto,
A. Traficante,
D. J. Eden,
P. G. Martin,
M. Sormani,
A. Coletta,
T. Colman,
R. Plume,
Y. Maruccia,
C. Mininni
, et al. (1 additional authors not shown)
Abstract:
We present a new derivation of the Milky Way's current star formation rate (SFR) based on the data of the Hi-GAL Galactic plane survey. We estimate the distribution of the SFR across the Galactic plane from the star-forming clumps identified in the Hi-GAL survey and calculate the total SFR from the sum of their contributions. The estimate of the global SFR amounts to $2.0 \pm 0.7$~M$_{\odot}$~yr…
▽ More
We present a new derivation of the Milky Way's current star formation rate (SFR) based on the data of the Hi-GAL Galactic plane survey. We estimate the distribution of the SFR across the Galactic plane from the star-forming clumps identified in the Hi-GAL survey and calculate the total SFR from the sum of their contributions. The estimate of the global SFR amounts to $2.0 \pm 0.7$~M$_{\odot}$~yr$^{-1}$, of which $1.7 \pm 0.6$~M$_{\odot}$~yr$^{-1}$ coming from clumps with reliable heliocentric distance assignment. This value is in general agreement with estimates found in the literature of last decades. The profile of SFR density averaged in Galactocentric rings is found to be qualitatively similar to others previously computed, with a peak corresponding to the Central Molecular Zone and another one around Galactocentric radius $R_\mathrm{gal} \sim 5$~kpc, followed by an exponential decrease as $\log(Σ_\mathrm{SFR}/[\mathrm{M}_\odot~\mathrm{yr}^{-1}~\mathrm{kpc}^{-2}])=a\,R_\mathrm{gal}/[\mathrm{kpc}]+b $, with $a=-0.28 \pm 0.01$. In this regard, the fraction of SFR produced within and outside the Solar circle is 84\% and 16\%, respectively; the fraction corresponding to the far outer Galaxy ($R_\mathrm{gal} > 13.5$~kpc) is only 1\%. We also find that, for $R_\mathrm{gal}>3$~kpc, our data follow a power law as a function of density, similarly to the Kennicutt-Schmidt relation. Finally, we compare the distribution of the SFR density across the face-on Galactic plane and those of median parameters, such as temperature, luminosity/mass ratio and bolometric temperature, describing the evolutionary stage of Hi-GAL clumps. We found no clear correlation between the SFR and the clump evolutionary stage.
△ Less
Submitted 10 November, 2022;
originally announced November 2022.
-
How magnetic field and stellar radiative feedback influences the collapse and the stellar mass spectrum of a massive star forming clump
Authors:
Patrick Hennebelle,
Ugo Lebreuilly,
Tine Colman,
Davide Elia,
Gary Fuller,
Silvia Leurini,
Thomas Nony,
Eugenio Schisano,
Juan D. Soler,
Alessio Traficante,
Ralf S. Klessen,
Sergio Molinari,
Leonardo Testi
Abstract:
In spite of decades of theoretical efforts, the physical origin of the stellar initial mass function (IMF) is still debated. We aim at understanding the influence of various physical processes such as radiative stellar feedback, magnetic field and non-ideal magneto-hydrodynamics on the IMF. We present a series of numerical simulations of collapsing 1000 M$_\odot$ clumps taking into account radiati…
▽ More
In spite of decades of theoretical efforts, the physical origin of the stellar initial mass function (IMF) is still debated. We aim at understanding the influence of various physical processes such as radiative stellar feedback, magnetic field and non-ideal magneto-hydrodynamics on the IMF. We present a series of numerical simulations of collapsing 1000 M$_\odot$ clumps taking into account radiative feedback and magnetic field with spatial resolution down to 1 AU. Both ideal and non-ideal MHD runs are performed and various radiative feedback efficiencies are considered. We also develop analytical models that we confront to the numerical results. The sum of the luminosities produced by the stars in the calculations is computed and it compares well with the bolometric luminosities reported in observations of massive star forming clumps. The temperatures, velocities and densities are also found to be in good agreement with recent observations. The stellar mass spectrum inferred for the simulations is, generally speaking, not strictly universal and in particular varies with magnetic intensity. It is also influenced by the choice of the radiative feedback efficiency. In all simulations, a sharp drop in the stellar distribution is found at about $M_{min} \simeq$ 0.1 M$_\odot$, which is likely a consequence of the adiabatic behaviour induced by dust opacities at high densities. As a consequence, when the combination of magnetic and thermal support is not too large, the mass distribution presents a peak located at 0.3-0.5 M$_\odot$. When magnetic and thermal support are large, the mass distribution is better described by a plateau, i.e. $d N / d \log M \propto M^{-Γ}$, $Γ\simeq 0$. Abridged
△ Less
Submitted 22 October, 2022;
originally announced October 2022.
-
The signature of large scale turbulence driving on the structure of the interstellar medium
Authors:
Tine Colman,
Jean-François Robitaille,
Patrick Hennebelle,
Marc-Antoine Miville-Deschênes,
Noé Brucy,
Ralf S. Klessen,
Simon C. O. Glover,
Juan D. Soler,
Davide Elia,
Alessio Traficante,
Sergio Molinari,
Leonardo Testi
Abstract:
The mechanisms that maintain turbulence in the interstellar medium (ISM) are still not identified. This work investigates how we can distinguish between two fundamental driving mechanisms: the accumulated effect of stellar feedback versus the energy injection from Galactic scales. We perform a series of numerical simulations describing a stratified star forming ISM subject to self-consistent stell…
▽ More
The mechanisms that maintain turbulence in the interstellar medium (ISM) are still not identified. This work investigates how we can distinguish between two fundamental driving mechanisms: the accumulated effect of stellar feedback versus the energy injection from Galactic scales. We perform a series of numerical simulations describing a stratified star forming ISM subject to self-consistent stellar feedback. Large scale external turbulent driving of various intensities is added to mimic galactic driving mechanisms. We analyse the resulting column density maps with a technique called Multi-scale non-Gaussian segmentation that separates the coherent structures and the Gaussian background. This effectively discriminates between the various simulations and is a promising method to understand the ISM structure. In particular the power spectrum of the coherent structures flattens above 60 pc when turbulence is driven only by stellar feedback. When large-scale driving is applied, the turn-over shifts to larger scales. A systematic comparison with the Large Magellanic Cloud (LMC) is then performed. Only 1 out of 25 regions has a coherent power spectrum which is consistent with the feedback-only simulation. A detailed study of the turn-over scale leads us to conclude that regular stellar feedback is not enough to explain the observed ISM structure on scales larger than 60 pc. Extreme feedback in the form of supergiant shells likely plays an important role but cannot explain all the regions of the LMC. If we assume ISM structure is generated by turbulence, another large scale driving mechanism is needed to explain the entirety of the observations.
△ Less
Submitted 1 June, 2022;
originally announced June 2022.
-
Protoplanetary disk birth in massive star forming clumps: the essential role of the magnetic field
Authors:
Ugo Lebreuilly,
Patrick Hennebelle,
Tine Colman,
Benoît Commerçon,
Ralf Klessen,
Anaëlle Maury,
Sergio Molinari,
Leonardo Testi
Abstract:
Protoplanetary disks form through angular momentum conservation in collapsing dense cores. In this work, we perform the first simulations with a maximal resolution down to the astronomical unit (au) of protoplanetary disk formation, through the collapse of 1000 solar mass clumps, treating self-consistently both non-ideal magnetohydrodynamics with ambipolar diffusion as well as radiative transfer i…
▽ More
Protoplanetary disks form through angular momentum conservation in collapsing dense cores. In this work, we perform the first simulations with a maximal resolution down to the astronomical unit (au) of protoplanetary disk formation, through the collapse of 1000 solar mass clumps, treating self-consistently both non-ideal magnetohydrodynamics with ambipolar diffusion as well as radiative transfer in the flux-limited diffusion approximation including stellar feedback. Using the adaptive mesh-refinement code RAMSES, we investigate the influence of the magnetic field on the disks properties with three models. We show that, without magnetic fields, a population dominated by large disks is formed, which is not consistent with Class 0 disk properties as estimated from observations. The inclusion of magnetic field leads, through magnetic braking, to a very different evolution. When it is included, small < 50 au disks represent about half the population. In addition, about ~ 70% of the stars have no disk in this case which suggests that our resolution is still insufficient to preserve the smaller disks. With ambipolar diffusion, the proportion of small disks is also prominent and we report a flat mass distribution around 0.01-0.1 solar mass and a typical disk-to-star mass ratios of ~0.01-0.1. This work shows that the magnetic field and its evolution plays a prominent role in setting the initial properties of disk populations.
△ Less
Submitted 18 July, 2021;
originally announced July 2021.
-
On the origin of the peak of the stellar initial mass function: exploring the tidal screening theory
Authors:
Tine Colman,
Romain Teyssier
Abstract:
Classical theories for the stellar initial mass function (IMF) predict a peak mass which scales with the properties of the molecular cloud. In this work, we explore a new theory proposed by Lee & Hennebelle (2018). The idea is that the tidal field around first Larson cores prevents the formation of other collapsing clumps within a certain radius. The protostar can then freely accrete the gas withi…
▽ More
Classical theories for the stellar initial mass function (IMF) predict a peak mass which scales with the properties of the molecular cloud. In this work, we explore a new theory proposed by Lee & Hennebelle (2018). The idea is that the tidal field around first Larson cores prevents the formation of other collapsing clumps within a certain radius. The protostar can then freely accrete the gas within this radius. This leads to a peak mass of roughly $10 \, M_{\mathrm{1LC}}$, independent of the parent cloud properties. Using simple analytical arguments, we derive a collapse condition for clumps located close to a protostar. We then study the tidal field and the corresponding collapse condition using a series of numerical simulations. We find that the tidal field around protostars is indeed strong enough to prevent clumps from collapsing unless they have high enough densities. For each newly formed protostar, we determine the region in which tidal screening is dominant. We call this the tidal bubble. The mass within this bubble is our estimate for the final mass of the star. Using this formalism, we are able to construct a very good prediction for the final IMF in our simulations. Not only do we correctly predict the peak, but we are also able to reproduced the high and low mass end of the IMF. We conclude that tidal forces are important in determining the final mass of a star and might be the dominant effect in setting the peak mass of the IMF.
△ Less
Submitted 17 November, 2019;
originally announced November 2019.