Search

Crucial for star formation is the interplay between gravity and turbulence. The observed cloud virial parameter, $\alpha_{\mathrm{vir}}$, which is the ratio of twice the turbulent kinetic energy to the gravitational energy, is found to vary significantly in different environments, where the scatter among individual star-forming clouds can exceed an order of magnitude. Therefore, a strong dependence of the initial mass function (IMF) on $\alpha_{\mathrm{vir}}$ may challenge the notion of a universal IMF. To determine the role of $\alpha_{\mathrm{vir}}$ on the IMF, we compare the star-particle mass functions obtained in high-resolution magnetohydrodynamical simulations including jet and heating feedback, with $\alpha_{\mathrm{vir}}=0.0625$, $0.125$, and $0.5$. We find that varying $\alpha_{\mathrm{vir}}$ from $\alpha_{\mathrm{vir}}\sim0.5$ to $\alpha_{\mathrm{vir}}<0.1$ shifts the peak of the IMF to lower masses by a factor of $\sim2$ and increases the star formation rate by a similar factor. The dependence of the IMF and star formation rate on $\alpha_{\mathrm{vir}}$ is non-linear, with the dependence subsiding at $\alpha_{\mathrm{vir}}<0.1$. Our study shows a systematic dependence of the IMF on $\alpha_{\mathrm{vir}}$. Yet, it may not be measurable easily in observations, considering the uncertainties, and the relatively weak dependence found in this study., Comment: 16 pages, 12 figures, 2 tables, Published in MNRAS, minor changes

We study the wide-binary eccentricity ($e$) distribution in young star clusters and the role of turbulence in setting the form of the $e$ distribution using magnetohydrodynamical (MHD) simulations of star cluster formation. The simulations incorporate gravity, turbulence, magnetic fields, protostellar heating, and jets/outflows. We find that (1) simulations that employ purely compressive turbulence driving produce binaries with a superthermal $e$ distribution ($\alpha>1$ in $p(e) \propto e^\alpha$), while simulations with purely solenoidal driving or natural mixture of driving modes produce subthermal/thermal distributions ($\alpha \leq$ 1), (2) the $e$ distribution over the full range of binary separations in our simulations is set at the early stages of the star cluster formation process, (3) while binaries (separation of $r_{\mathrm{pair}} \leq 1000\, \mathrm{AU}$) have subthermal to thermal $e$ distributions ($\alpha \sim 0.8$), wide binaries ($r_{\mathrm{pair}} > 1000\, \mathrm{AU}$) have a superthermal distribution ($\alpha \sim 1.8$), and (4) low-mass binary systems (system masses of $M_{\mathrm{sys}} \leq 0.8\, \mathrm{M_\odot}$) have a highly superthermal distribution ($\alpha \sim 2.4$), whereas high-mass systems ($M_{\mathrm{sys}} > 0.8\, \mathrm{M_\odot}$) exhibit a subthermal/thermal distribution ($\alpha \sim 0.8$). The binary eccentricity distribution is often modelled as a thermal distribution. However, our results suggest that the $e$ distribution depends on the range of separation of the sampled binaries, which agrees with the findings from recent Gaia observations. We conclude that the dependence of the $e$ distribution on the binary separation and mass is linked to the binary formation mechanism governed by the turbulent properties of the parent cloud., Comment: 15 pages, 9 figures, 1 table (added additional citation, published in MNRAS)

Turbulence is a critical ingredient for star formation, yet its role for the initial mass function (IMF) is not fully understood. Here we perform magnetohydrodynamical (MHD) simulations of star cluster formation including gravity, turbulence, magnetic fields, stellar heating and outflow feedback to study the influence of the mode of turbulence driving on IMF. We find that simulations that employ purely compressive turbulence driving (COMP) produce a higher fraction of low-mass stars as compared to simulations that use purely solenoidal driving (SOL). The characteristic (median) mass of the sink particle (protostellar) distribution for COMP is shifted to lower masses by a factor of ~ 1.5 compared to SOL. Our simulation IMFs capture the important features of the observed IMF form. We find that turbulence-regulated theories of the IMF match our simulation IMFs reasonably well in the high-mass and low-mass range, but underestimate the number of very low-mass stars, which form towards the later stages of our simulations and stop accreting due to dynamical interactions. Our simulations show that for both COMP and SOL, the multiplicity fraction is an increasing function of the primary mass, although the multiplicity fraction in COMP is higher than that of SOL for any primary mass range. We find that binary mass ratio distribution is independent of the turbulence driving mode. The average specific angular momentum of the sink particles in SOL is a factor of 2 higher than that for COMP. Overall, we conclude that the turbulence driving mode plays a significant role in shaping the IMF., Comment: 26 pages, 16 figures, 4 tables, accepted for publication in MNRAS, minor changes, part of methodology is same as in Mathew & Federrath 2021. arXiv admin note: text overlap with arXiv:2106.06521

The interaction of turbulence, magnetic fields, self-gravity, and stellar feedback within molecular clouds is crucial for understanding star formation. We study the effects of self-gravity and outflow feedback on the properties of the turbulent velocity via the structure function over length scales from $\sim$ 0.01 pc to 2 pc. We analyze a series of three-dimensional, magnetohydrodynamical (MHD) simulations of star cluster formation. We find outflow feedback can change the scaling of velocity fluctuations but still roughly being in between Kolmogorov and Burgers turbulence. We observe that self-gravity and protostellar outflows increase the velocity fluctuations over all length scales. Outflows can amplify the velocity fluctuations by up to a factor of $\sim$7 on scales $\sim$ 0.01 - 0.2 pc and drive turbulence up to a scale of $\sim$ 1 pc. The amplified velocity fluctuations provide more support against gravity and enhance fragmentation on small scales. The self-gravity's effect is more significant on smaller dense clumps and it increases the fraction of the compressive velocity component up to a scale of $\sim$ 0.2 pc. However, outflow feedback drives both solenoidal and compressive modes, but it induces a higher fraction of solenoidal modes relative to compressive modes. Thus, with outflows, the dense core ends up with a slightly higher fraction of solenoidal modes. We find that the compressible fraction is fairly constant with about 1/3 on scales $\sim$ 0.1 - 0.2 pc. The combined effect of enhanced velocity dispersion and reduced compressive fraction contributes to a reduction in the star formation rate., Comment: 11 pages, 9 figures, accepted for publication in MNRAS

We perform a series of three-dimensional, magnetohydrodynamical (MHD) simulations of star cluster formation including gravity, turbulence, magnetic fields, stellar radiative heating and outflow feedback. We observe that the inclusion of protostellar outflows (1) reduces the star formation rate by a factor of $\sim2$, (2) increases fragmentation, and (3) shifts the initial mass function (IMF) to lower masses by a factor of $2.0\pm0.2$, without significantly affecting the overall shape of the IMF. The form of the sink particle (protostellar objects) mass distribution obtained from our simulations matches the observational IMFs reasonably well. We also show that turbulence-based theoretical models of the IMF agree well with our simulation IMF in the high-mass and low-mass regime, but do not predict any brown dwarfs, whereas our simulations produce a considerable number of sub-stellar objects, which are produced by dynamical interactions (ejections). We find that these dynamical interactions also play a key role for the binary separation distribution and stellar kinematics in general. Our numerical model of star cluster formation also reproduces the observed mass dependence of multiplicity. Our multiplicity fraction estimates generally concur with the observational estimates for different spectral types. We further calculate the specific angular momentum of all the sink particles and find that the average value of $1.5 \times 10^{19}\, \mathrm{cm^2\, s^{-1}}$ is consistent with observational data. The specific angular momentum of our sink particles lies in the range typical of protostellar envelopes and binaries. We conclude that the IMF is controlled by a combination of gravity, turbulence, magnetic fields, radiation and outflow feedback., Comment: 21 pages, 20 figures, 2 tables, accepted for publication in MNRAS, included brief discussions on stellar kinematics and binary separation distribution, few other minor changes including additional references

Explaining the initial mass function (IMF) of stars is a long-standing problem in astrophysics. The number of complex mechanisms involved in the process of star cluster formation, such as turbulence, magnetic fields and stellar feedback, make understanding and modeling the IMF a challenging task. In this paper, we aim to assert the importance of stellar heating feedback in the star cluster formation process and its effect on the shape of the IMF. We use an analytical sub-grid model to implement the radiative feedback in fully three-dimensional magnetohydrodynamical (MHD) simulations of star cluster formation, with the ultimate objective of obtaining numerical convergence on the IMF. We compare a set of MHD adaptive-mesh-refinement (AMR) simulations with three different implementations of the heating of the gas: 1) a polytropic equation of state (EOS), 2) a spherically symmetric stellar heating feedback, and 3) our newly developed polar heating model that takes into account the geometry of the accretion disc and the resulting shielding of stellar radiation by dust. For each of the three heating models, we analyse the distribution of stellar masses formed in ten molecular cloud simulations with different realizations of the turbulence to obtain a statistically representative IMF. We conclude that stellar heating feedback has a profound influence on the number of stars formed and plays a crucial role in controlling the IMF. We find that the simulations with the polar heating model achieve the best convergence on the observed IMF., Comment: 10 pages, 5 figures, 1 table; added references, published in MNRAS