No Arabic abstract
The connection between the pre-stellar core mass function (CMF) and the stellar initial mass function (IMF) lies at the heart of all star formation theories. In this paper, we study the earliest phases of star formation with a series of high-resolution numerical simulations that include the formation of sinks. In particular, we focus on the transition from cores to sinks within a massive molecular filament. We compare the CMF and IMF between magnetized and unmagnetized simulations, and between different resolutions. We find that selecting cores based on their kinematic virial parameter excludes collapsing objects because they host large velocity dispersions. Selecting only the thermally unstable magnetized cores, we observe that their mass-to-flux ratio spans almost two orders of magnitude for a given mass. We also see that, when magnetic fields are included, the CMF peaks at higher core mass values with respect to pure hydrodynamical simulations. Nonetheless, all models produce sink mass functions with a high-mass slope consistent with Salpeter. Finally, we examine the effects of resolution and find that, in isothermal simulations, even models with very high dynamical range fail to converge in the mass function. Our main conclusion is that, although the resulting CMFs and IMFs have similar slopes in all simulations, the cores have slightly different sizes and kinematical properties when a magnetic field is included. However, a core selection based on the mass-to-flux ratio alone is not enough to alter the shape of the CMF, if we do not take thermal stability into account. Finally, we conclude that extreme care should be given to resolution issues when studying sink formation with an isothermal equation of state.
We present high-angular-resolution ALMA (Atacama Large Millimeter Array) images of N$_{2}$H$^{+}$ (1--0) that has been combined with those from the Nobeyama telescope toward OMC-2 and OMC-3 filamentary regions. The filaments (with typical widths of $sim$ 0.1 pc) and dense cores are resolved. The measured 2D velocity gradients of cores are between 1.3 and 16.7 km,s$^{-1}$,pc$^{-1}$, corresponding to a specific angular momentum ($J/M$) between 0.0012 and 0.016 pc,km,s$^{-1}$. With respect to the core size $R$, the specific angular momentum follows a power law $J/M propto R^{1.52~pm~0.14}$. The ratio ($beta$) between the rotational energy and gravitational energy ranges from 0.00041 to 0.094, indicating insignificant support from rotation against gravitational collapse. We further focus on the alignment between the cores rotational axes, which is defined to be perpendicular to the direction of the velocity gradient ($theta_{G}$), and the direction of elongation of filaments ($theta_{f}$) in this massive star-forming region. The distribution of the angle between $theta_{f}$ and $theta_{G}$ was f ound to be random, i.e. the cores rotational axes have no discernible correlation with the elongation of their hosting filament. This implies that, in terms of angular momentum, the cores have evolved to be dynamically independent from their natal filaments.
We investigate the dependence of the peak of the IMF on the physics of the so-called first Larson core, which corresponds to the point where the dust becomes opaque to its own radiation. We perform numerical simulations of collapsing clouds of $1000 M_odot$ for various gas equation of state (eos), paying great attention to the numerical resolution and convergence. The initial conditions of these numerical experiments are varied in the companion paper. We also develop analytical models that we confront to our numerical results. If an isothermal eos is used, we show that the peak of the IMF shifts to lower masses with improved numerical resolution. When an adiabatic eos is employed, numerical convergence is obtained. The peak position varies with the eos and we find that the peak position is about ten times the mass of the first Larson core. By analyzing the stability of non-linear density fluctuations in the vicinity of a point mass and then summing over a reasonable density distribution, we find that tidal forces exert a strong stabilizing effect and likely lead to a preferential mass several times larger than that of the first Larson core. We propose that in a sufficiently massive and cold cloud, the peak of the IMF is determined by the thermodynamics of the high density adiabatic gas as well as the stabilizing influence of tidal forces. The resulting characteristic mass is about ten times the mass of the first Larson core, which altogether leads to a few tenths of solar masses. Since these processes are not related to the large scale physical conditions and to the environment, our results suggest a possible explanation for the apparent universality of the peak of the IMF.
Multi-phase filamentary structures around Brightest Cluster Galaxies are likely a key step of AGN-feedback. We observed molecular gas in 3 cool cluster cores: Centaurus, Abell S1101, and RXJ1539.5 and gathered ALMA and MUSE data for 12 other clusters. Those observations show clumpy, massive and long, 3--25 kpc, molecular filaments, preferentially located around the radio bubbles inflated by the AGN (Active Galactic Nucleus). Two objects show nuclear molecular disks. The optical nebula is certainly tracing the warm envelopes of cold molecular filaments. Surprisingly, the radial profile of the H$alpha$/CO flux ratio is roughly constant for most of the objects, suggesting that (i) between 1.2 to 7 times more cold gas could be present and (ii) local processes must be responsible for the excitation. Projected velocities are between 100--400 km s$^{-1}$, with disturbed kinematics and sometimes coherent gradients. This is likely due to the mixing in projection of several thin unresolved filaments. The velocity fields may be stirred by turbulence induced by bubbles, jets or merger-induced sloshing. Velocity and dispersions are low, below the escape velocity. Cold clouds should eventually fall back and fuel the AGN. We compare the filaments radial extent, r$_{fil}$, with the region where the X-ray gas can become thermally unstable. The filaments are always inside the low-entropy and short cooling time region, where t$_{cool}$/t$_{ff}$<20 (9 of 13 sources). The range t$_{cool}$/t$_{ff}$, 8-23 at r$_{fil}$, is likely due to (i) a more complex gravitational potential affecting the free-fall time (e.g., sloshing, mergers); (ii) the presence of inhomogeneities or uplifted gas in the ICM, affecting the cooling time. For some of the sources, r$_{fil}$ lies where the ratio of the cooling time to the eddy-turnover time, t$_{cool}$/t$_{eddy}$, is approximately unity.
The stellar mass spectrum is an important property of the stellar cluster and a fundamental quantity to understand our Universe. The fragmentation of diffuse molecular cloud into stars is subject to physical processes such as gravity, turbulence, thermal pressure, and magnetic field. The final mass of a star is believed to be a combined outcome of a virially unstable reservoir and subsequent accretion. We aim to clarify the roles of different supporting energies, notably the thermal pressure and the magnetic field, in determining the stellar mass. Following previous studies by Lee & Hennebelle (2018a,b), we perform a series of numerical experiments of stellar cluster formation inside an isolated molecular clump. By changing the effective equation of state (EOS) of the diffuse gas (that is to say gas whose density is below the critical density at which dust becomes opaque to its radiation) and the strength of the magnetic field, we investigate whether any characteristic mass is introduced into the fragmentation processes. The EOS of the diffuse gas, including the bulk temperature and the polytropic index, does not affect significantly the shape of the stellar mass spectrum. The presence of magnetic field slightly modifies the shape of the mass spectrum only when extreme values are applied. This study confirms that the peak of the IMF is primarily determined by the adiabatic high-density end of the EOS that mimics the radiation inside the high-density gas. Furthermore, the shape of the mass spectrum is mostly sensitive to the density PDF, and the magnetic field has likely only a secondary role. In particular, we stress that the Jeans mass at the mean cloud density and at the critical density are not responsible of setting the peak.
We study the core mass function (CMF) of the massive protocluster G286.21+0.17 with the Atacama Large Millimeter/submillimeter Array via 1.3~mm continuum emission at a resolution of 1.0arcsec (2500~au). We have mapped a field of 5.3arcmin$times$5.3arcmin centered on the protocluster clump. We measure the CMF in the central region, exploring various core detection algorithms, which give source numbers ranging from 60 to 125, depending on parameter selection. We estimate completeness corrections due to imperfect flux recovery and core identification via artificial core insertion experiments. For masses $Mgtrsim1:M_odot$, the fiducial dendrogram-identified CMF can be fit with a power law of the form ${rm{d}}N/{rm{d}}{rm{log}}Mpropto{M}^{-alpha}$ with $alpha simeq1.24pm0.17$, slightly shallower than, but still consistent with, the index of the Salpeter stellar initial mass function of 1.35. Clumpfind-identified CMFs are significantly shallower with $alphasimeq0.64pm0.13$. While raw CMFs show a peak near $1:M_odot$, completeness-corrected CMFs are consistent with a single power law extending down to $sim 0.5:M_odot$, with only a tentative indication of a shallowing of the slope around $sim1:M_odot$. We discuss the implications of these results for star and star cluster formation theories.