Understanding the evolution of self-gravitating, isothermal, magnetized gas is crucial for star formation, as these physical processes have been postulated to set the initial mass function (IMF). We present a suite of isothermal magnetohydrodynamic (MHD) simulations using the GIZMO code, that resolve the formation of individual stars in giant molecular clouds (GMCs), spanning a range of Mach numbers found in observed GMCs. As in past works, the mean and median stellar masses are sensitive to numerical resolution, because they are sensitive to low-mass stars that contribute a vanishing fraction of the overall stellar mass. The {\em mass-weighted} median stellar mass $M_\mathrm{50}$ becomes insensitive to resolution once turbulent fragmentation is well-resolved. Without imposing Larson-like scaling laws, our simulations find $M_\mathrm{50} \propto M_\mathrm{0} \mathcal{M}^{-3} \alpha_\mathrm{turb} \mathrm{SFE}^{1/3}$ for GMC mass $M_\mathrm{0}$, sonic Mach number $\mathcal{M}$, virial parameter $\alpha_\mathrm{turb}$, and star formation efficiency $\mathrm{SFE}=M_\mathrm{\star}/M_\mathrm{0}$. This fit agrees well with previous IMF results from the RAMSES, ORION2, and SphNG codes. Although $M_\mathrm{50}$ has no significant dependence on the magnetic field strength at the cloud scale, MHD is necessary to prevent a fragmentation cascade that results in non-convergent stellar masses. For initial conditions and SFE similar to star-forming GMCs in our Galaxy, we predict $M_\mathrm{50}$ to be $>20 M_{\odot}$, an order of magnitude larger than observed ($\sim 2 M_\odot$), together with an excess of brown dwarfs. Moreover, $M_\mathrm{50}$ is sensitive to initial cloud properties and evolves strongly in time within a given cloud, predicting much larger IMF variations than are observationally allowed. We conclude that physics beyond MHD turbulence and gravity are necessary ingredients for the IMF.

Giant molecular clouds (GMCs) are well-studied in the local Universe, however, exactly how their properties vary during galaxy evolution is poorly understood due to challenging resolution requirements, both observational and computational. We present the first time-dependent analysis of giant molecular clouds in a Milky Way-like galaxy and an LMC-like dwarf galaxy of the FIRE-2 (Feedback In Realistic Environments) simulation suite, which have sufficient resolution to predict the bulk properties of GMCs in cosmological galaxy formation self-consistently. We show explicitly that the majority of star formation outside the galactic center occurs within self-gravitating gas structures that have properties consistent with observed bound GMCs. We find that the typical cloud bulk properties such as mass and surface density do not vary more than a factor of 2 in any systematic way after the first Gyr of cosmic evolution within a given galaxy from its progenitor. While the median properties are constant, the tails of the distributions can briefly undergo drastic changes, which can produce very massive and dense self-gravitating gas clouds. Once the galaxy forms, we identify only two systematic trends in bulk properties over cosmic time: a steady increase in metallicity produced by previous stellar populations and a weak decrease in bulk cloud temperatures. With the exception of metallicity we find no significant differences in cloud properties between the Milky Way-like and dwarf galaxies. These results have important implications for cosmological star and star cluster formation and put especially strong constraints on theories relating the stellar initial mass function to cloud properties.

We present the first measurement of the lifetimes of Giant Molecular Clouds (GMCs) in cosmological simulations at $z = 0$, using the Latte suite of FIRE-2 simulations of Milky Way-mass galaxies. We track GMCs with total gas mass $\gtrsim 10^5$ M$_\odot$ at high spatial ($\sim1$ pc), mass ($7100$ M$_{\odot}$), and temporal (1 Myr) resolution. Our simulated GMCs are consistent with the distribution of masses for massive GMCs in the Milky Way and nearby galaxies. We find GMC lifetimes of $5-7$ Myr, or 1-2 freefall times, on average, with less than 1\% of clouds living longer than 20 Myr. We find increasing GMC lifetimes with galactocentric radius, implying that environment affects the evolutionary cycle of GMCs. However, our GMC lifetimes show no systematic dependence on GMC mass or amount of star formation. These results are broadly consistent with inferences from the literature and provide an initial investigation into ultimately understanding the physical processes that govern GMC lifetimes in a cosmological setting.

One of the most robust observations of the stellar initial mass function (IMF) is its near-universality in the Milky Way and neighboring galaxies. But recent observations of early-type galaxies can be interpreted to imply a bottom-heavy IMF, while others of ultra-faint dwarfs could imply a top-heavy IMF. This would impose powerful constraints on star formation models. We explore what sort of cloud-scale IMF models could possibly satisfy these constraints. We utilize simulated galaxies which reproduce (broadly) the observed galaxy properties, while they also provide the detailed star formation history and properties of each progenitor star-forming cloud. We then consider generic models where the characteristic mass of the IMF is some arbitrary power-law function of progenitor cloud properties, along with well-known literature IMF models, which scale with Jeans mass, turbulent Bonnor-Ebert mass, temperature, the opacity limit, metallicity, or the protostellar heating mass. We show that no IMF models currently in the literature - nor any model where the turnover mass is an arbitrary power-law function of a combination of cloud temperature/density/size/metallicity/velocity dispersion/magnetic field - can reproduce the claimed IMF variation in ellipticals or dwarfs without severely violating observational constraints in the Milky Way. Specifically, they predict too much variation in the extreme environments of the Galaxy, compared to that observed. Either the IMF varies in a more complicated manner, or alternative interpretations of the extragalactic observations must be explored.

Young massive star clusters spanning $\sim 10^4 - 10^8 M_\odot$ in mass have been observed to have similar surface brightness profiles. Recent hydrodynamical simulations of star cluster formation have also produced star clusters with this structure. We argue analytically that this type of mass distribution arises naturally in the relaxation from a hierarchically-clustered distribution of stars into a monolithic star cluster through hierarchical merging. We show that arbitrary initial profiles will tend to converge to a universal profile under hierarchical merging, owing to phase-space mixing obeying certain conservation constraints. We perform $N$-body simulations of a pairwise merger of model star clusters and find that mergers readily produce the shallow surface brightness profiles observed in young massive clusters. Finally, we simulate the relaxation of a hierarchically-clustered mass distribution constructed from an idealized fragmentation model. Assuming only power-law spatial and kinematic scaling relations, these numerical experiments are able to reproduce the surface density profiles of observed young massive star clusters. Thus we provide physical motivation for the structure of young massive clusters within the paradigm of hierarchical star formation. This has important implications for the structure of nascent globular clusters.

The evolution of self-gravitating clouds of isothermal gas forms the basis of many star formation theories. Therefore it is important to know under what conditions such a cloud will undergo homologous collapse into a single, massive object, or will fragment into a spectrum of smaller ones. And if it fragments, do initial conditions (e.g. Jeans mass, sonic mass) influence the mass function of the fragments, as predicted by many theories of star formation? In this paper we show that the relevant parameter separating homologous collapse from fragmentation is not the Mach number of the initial turbulence (as suspected by many), but the infall Mach number $\mathcal{M}_{\rm infall}\sim\sqrt{G M/(R c_s^2)}$, equivalent to the number of Jeans masses in the initial cloud $N_J$. We also show that fragmenting clouds produce a power-law mass function with slopes close to the expected -2 (i.e. equal mass in all logarithmic mass intervals). However, the low-mass cut-off of this mass function is entirely numerical; the initial properties of the cloud have no effect on it. In other words, if $\mathcal{M}_{\rm infall}\gg 1$, fragmentation proceeds without limit to masses much smaller than the initial Jeans mass.

A large number of astronomical phenomena exhibit remarkably similar scaling relations. The most well-known of these is the mass distribution $\mathrm{d} N/\mathrm{d} M\propto M^{-2}$ which (to first order) describes stars, protostellar cores, clumps, giant molecular clouds, star clusters and even dark matter halos. In this paper we propose that this ubiquity is not a coincidence and that it is the generic result of scale-free structure formation where the different scales are uncorrelated. We show that all such systems produce a mass function proportional to $M^{-2}$ and a column density distribution with a power law tail of $\mathrm{d} A/\mathrm{d} \ln\Sigma\propto\Sigma^{-1}$. In the case where structure formation is controlled by gravity the two-point correlation becomes $\xi_{2D}\propto R^{-1}$. Furthermore, structures formed by such processes (e.g. young star clusters, DM halos) tend to a $\rho\propto R^{-3}$ density profile. We compare these predictions with observations, analytical fragmentation cascade models, semi-analytical models of gravito-turbulent fragmentation and detailed "full physics" hydrodynamical simulations. We find that these power-laws are good first order descriptions in all cases.

One of the key observations regarding the stellar initial mass function (IMF) is its near-universality in the Milky Way (MW), which provides a powerful way to constrain different star formation models that predict the IMF. However, those models are almost universally "cloud-scale" or smaller -- they take as input or simulate single molecular clouds (GMCs), clumps, or cores, and predict the resulting IMF as a function of the cloud properties. Without a model for the progenitor properties of all clouds which formed the stars at different locations in the MW (including ancient stellar populations formed in high-redshift, likely gas-rich dwarf progenitor galaxies that looked little like the Galaxy today), the predictions cannot be fully explored, nor safely applied to "live" cosmological calculations of the IMF in different galaxies at different cosmological times. We therefore combine a suite of high-resolution cosmological simulations (from the Feedback In Realistic Environments project), which form MW-like galaxies with reasonable star formation properties and explicitly resolve massive GMCs, with various proposed cloud-scale IMF models. We apply the models independently to every star particle formed in the simulations to synthesize the predicted IMF in the present-day galaxy. We explore models where the IMF depends on Jeans mass, sonic or turbulent Bonner-Ebert mass, fragmentation with a polytropic equation-of-state, or where it is self-regulated by protostellar feedback. We show that all of these models, except the feedback-regulated ones, predict far more variation ($\sim 0.6-1$ dex $1\,\sigma$ scatter in the IMF turnover mass) in the simulations than is observed in the MW.

Stars are strongly clustered on both large (~pc) and small (~binary) scales, but there are few analytic or even semi-analytic theories for the correlation function and multiplicity of stars. In this paper we present such a theory, based on our recently-developed semi-analytic framework called MISFIT, which models gravito-turbulent fragmentation, including the suppression of fragmentation by protostellar radiation feedback. We compare the results including feedback to a control model in which it is omitted. We show that both classes of models robustly reproduce the stellar correlation function at >0.01 pc scales, which is well approximated by a power-law that follows generally from scale-free physics (turbulence plus gravity) on large scales. On smaller scales protostellar disk fragmentation becomes dominant over common core fragmentation, leading to a steepening of the correlation function. Multiplicity is more sensitive to feedback: we found that a model with the protostellar heating reproduces the observed multiplicity fractions and mass ratio distributions for both Solar and sub-Solar mass stars (in particular the brown dwarf desert), while a model without feedback fails to do so. The model with feedback also produces an at-formation period distribution consistent with the one inferred from observations. However, it is unable to produce short-range binaries below the length scale of protostellar disks. We suggest that such close binaries are produced primarily by disk fragmentation and further decrease their separation through orbital decay.

Turbulence is thought to be a primary driving force behind the early stages of star formation. In this framework large, self gravitating, turbulent clouds fragment into smaller clouds which in turn fragment into even smaller ones. At the end of this cascade we find the clouds which collapse into protostars. Following this process is extremely challenging numerically due to the large dynamical range so in this paper we propose a semi analytic framework which is able to follow star formation from the largest, giant molecular cloud (GMC) scale, to the final protostellar size scale. Due to the simplicity of the framework it is ideal for theoretical experimentation to explore the principal processes behind different aspects of star formation, at the cost of strong assumptions. The basic version of the model discussed in this paper only contains turbulence, gravity and crude assumptions about feedback, nevertheless it can reproduce the observed core mass function (CMF) and provide the protostellar system mass function (PSMF), which shows a striking resemblance to the observed IMF. Furthermore we find that to produce a universal IMF protostellar feedback must be taken into account otherwise the PSMF peak shows a strong dependence on the background temperature.

A popular theory of star formation is gravito-turbulent fragmentation, in which self-gravitating structures are created by turbulence-driven density fluctuations. Simple theories of isothermal fragmentation successfully reproduce the core mass function (CMF) which has a very similar shape to the initial mass function (IMF) of stars. However, numerical simulations of isothermal turbulent fragmentation thus far have not succeeded in identifying a fragment mass scale that is independent of the simulation resolution. Moreover, the fluid equations for magnetized, self-gravitating, isothermal turbulence are scale-free, and do not predict any characteristic mass. In this paper we show that, although an isothermal self-gravitating flow does produce a CMF with a mass scale imposed by the initial conditions, this scale changes as the parent cloud evolves. In addition, the cores that form undergo further fragmentation and after sufficient time forget about their initial conditions, yielding a scale-free pure power-law distribution $\mathrm{d} N/\mathrm{d} M\propto M^{-2}$ for the stellar IMF. We show that this problem can be alleviated by introducing additional physics that provides a termination scale for the cascade. Our candidate for such physics is a simple model for stellar radiation feedback. Radiative heating, powered by accretion onto forming stars, arrests the fragmentation cascade and imposes a characteristic mass scale that is nearly independent of the time-evolution or initial conditions in the star-forming cloud, and that agrees well with the peak of the observed IMF. In contrast, models that introduce a stiff equation of state for denser clouds but that do not explicitly include the effects of feedback do not yield an invariant IMF.

It has been shown that fragmentation within self-gravitating, turbulent molecular clouds ("turbulent fragmentation") can naturally explain the observed properties of protostellar cores, including the core mass function (CMF). Here, we extend recently-developed analytic models for turbulent fragmentation to follow the time-dependent hierarchical fragmentation of self-gravitating cores, until they reach effectively infinite density (and form stars). We show that turbulent fragmentation robustly predicts two key features of the IMF. First, a high-mass power-law scaling very close to the Salpeter slope, which is a generic consequence of the scale-free nature of turbulence and self-gravity. We predict the IMF slope (-2.3) is slightly steeper then the CMF slope (-2.1), owing to the slower collapse and easier fragmentation of large cores. Second, a turnover mass, which is set by a combination of the CMF turnover mass (a couple solar masses, determined by the `sonic scale' of galactic turbulence, and so weakly dependent on galaxy properties), and the equation of state (EOS). A "soft" EOS with polytropic index $\gamma<1.0$ predicts that the IMF slope becomes "shallow" below the sonic scale, but fails to produce the full turnover observed. An EOS which becomes "stiff" at sufficiently low surface densities $\Sigma_{\rm gas} \sim 5000\,M_{\odot}\,{\rm pc^{-2}}$, and/or models where each collapsing core is able to heat and effectively stiffen the EOS of a modest mass ($\sim 0.02\,M_{\odot}$) of surrounding gas, are able to reproduce the observed turnover. Such features are likely a consequence of more detailed chemistry and radiative feedback.