Abstract
Observations of ionized carbon at 158 μm (C ii) from luminous star-forming galaxies at z ∼ 0 show that their ratios of C ii to far-infrared (FIR) luminosity are systematically lower than ...those of more modestly star-forming galaxies. In this paper, we provide a theory for the origin of this so-called C ii deficit in galaxies. Our model treats the interstellar medium as a collection of clouds with radially stratified chemical and thermal properties, which are dictated by the clouds’ volume and surface densities, as well as the interstellar radiation and cosmic ray fields to which they are exposed. C ii emission arises from the outer, H i-dominated layers of clouds, and from regions where the hydrogen is H2 but the carbon is predominantly C+. In contrast, the most shielded regions of clouds are dominated by CO, and produce little C ii emission. This provides a natural mechanism to explain the observed C ii–star formation relation: galaxies’ star formation rates are largely driven by the surface densities of their clouds. As this rises, so does the fraction of gas in the CO-dominated phase that produces little C ii emission. Our model further suggests that the apparent offset in the C ii–FIR relation for high-z sources compared to those at present epoch may arise from systematically larger gas masses at early times: a galaxy with a large gas mass can sustain a high star formation rate even with a relatively modest surface density, allowing copious C ii emission to coexist with rapid star formation.
How do bound star clusters form? Krumholz, Mark R; McKee, Christopher F
Monthly Notices of the Royal Astronomical Society,
05/2020, Volume:
494, Issue:
1
Journal Article
Peer reviewed
Open access
ABSTRACT
Gravitationally bound clusters that survive gas removal represent an unusual mode of star formation in the Milky Way and similar spiral galaxies. While forming, they can be distinguished ...observationally from unbound star formation by their high densities, virialized velocity structures, and star formation histories that accelerate towards the present, but extend multiple free-fall times into the past. In this paper, we examine several proposed scenarios for how such structures might form and evolve, and carry out a Bayesian analysis to test these models against observed distributions of protostellar age, counts of young stellar objects relative to gas, and the overall star formation rate of the Milky Way. We show that models in which the acceleration of star formation is due either to a large-scale collapse or a time-dependent increase in star formation efficiency are unable to satisfy the combined set of observational constraints. In contrast, models in which clusters form in a ‘conveyor belt’ mode where gas accretion and star formation occur simultaneously, but the star formation rate per free-fall time is low, can match the observations.
The star formation rates (SFRs) of low-metallicity galaxies depend sensitively on the gas metallicity, because metals are crucial to mediating the transition from intermediate-temperature atomic gas ...to cold molecular gas, a necessary precursor to star formation. We study the impact of this effect on the star formation history of galaxies. We incorporate metallicity-dependent star formation and metal enrichment in a simple model that follows the evolution of a halo main progenitor. We specify the expected dependence of sSFR and metallicity on stellar mass and redshift. At a given z, and below a critical mass, these relations are predicted to be flat and rising, respectively. Our model predictions qualitatively match some of the puzzling features in the observed star formation history.
Young stars typically form in star clusters, so the supernovae (SNe) they produce are clustered in space and time. This clustering of SNe may alter the momentum per SN deposited in the interstellar ...medium (ISM) by affecting the local ISM density, which in turn affects the cooling rate. We study the effect of multiple SNe using idealized 1D hydrodynamic simulations which explore a large parameter space of the number of SNe, and the background gas density and metallicity. The results are provided as a table and an analytic fitting formula. We find that for clusters with up to ~100 SNe, the asymptotic momentum scales superlinearly with the number of SNe, resulting in a momentum per SN which can be an order of magnitude larger than for a single SN, with a maximum efficiency for clusters with 10-100 SNe. We argue that additional physical processes not included in our simulations -- self-gravity, breakout from a galactic disc, and galactic shear -- can slightly reduce the momentum enhancement from clustering, but the average momentum per SN still remains a factor of 4 larger than the isolated SN value when averaged over a realistic cluster mass function for a star-forming galaxy. We conclude with a discussion of the possible role of mixing between hot and cold gas, induced by multidimensional instabilities or pre-existing density variations, as a limiting factor in the build-up of momentum by clustered SNe, and suggest future numerical experiments to explore these effects.
While the evolution of superbubbles driven by clustered supernovae has been studied by numerous authors, the resulting radial momentum yield is uncertain by as much as an order of magnitude depending ...on the computational methods and assumed properties of the surrounding interstellar medium (ISM). In this work, we study the origin of these discrepancies, and seek to determine the correct momentum budget for a homogeneous ISM. We carry out 3D hydrodynamic and magnetohydrodynamic (MHD) simulations of clustered supernova explosions, using a Lagrangian method and checking for convergence with respect to resolution. We find that the terminal momentum of a shell driven by clustered supernovae is dictated primarily by the mixing rate across the contact discontinuity between the hot and cold phases, and that this energy mixing rate is dominated by numerical diffusion even at the highest resolution we can complete, 0.03 M⊙. Magnetic fields also reduce the mixing rate, so that MHD simulations produce higher momentum yields than HD ones at equal resolution. As a result, we obtain only a lower limit on the momentum yield from clustered supernovae. Combining this with our previous 1D results, which provide an upper limit because they allow almost no mixing across the contact discontinuity, we conclude that the momentum yield per supernova from clustered supernovae in a homogeneous ISM is bounded between 2 × 105 and 3 × 106 M⊙ km s-1. A converged value for the simple homogeneous ISM remains elusive.
Observations of molecular gas in high-z star-forming galaxies
typically rely on emission from CO lines arising from states with rotational quantum numbers J > 1. Converting these observations to an ...estimate of the CO J = 1-0 intensity, and thus inferring H2 gas masses, requires knowledge of the CO excitation ladder or spectral line energy distribution (SLED). The few available multi-J CO observations of galaxies show a very broad range of SLEDs, even at fixed galaxy mass and star formation rate (SFR), making the conversion to J = 1-0 emission and hence molecular gas mass highly uncertain. Here, we combine numerical simulations of disc galaxies and galaxy mergers with molecular line radiative transfer calculations to develop a model for the physical parameters that drive variations in CO SLEDs in galaxies. An essential feature of our model is a fully self-consistent computation of the molecular gas temperature and excitation structure. We find that, while the shape of the SLED is ultimately determined by difficult-to-observe quantities such as the gas density, temperature and optical depth distributions, all of these quantities are well correlated with the galaxy's mean star formation rate surface density (ΣSFR), which is observable. We use this result to develop a model for the CO SLED in terms of ΣSFR, and show that this model quantitatively reproduces the SLEDs of galaxies over a dynamic range of ∼200 in SFR surface density, at redshifts from z = 0 to 6. This model should make it possible to significantly reduce the uncertainty in deducing molecular gas masses from observations of high-J CO emission.
The pressure exerted by the radiation of young stars may be an important feedback mechanism that drives turbulence and winds in forming star clusters and the disks of starburst galaxies. In this ...paper, we report a series of two-dimensional flux-limited diffusion radiation-hydrodynamics calculations performed with the code orion in which we drive strong radiation fluxes through columns of dusty matter confined by gravity in order to answer these questions. We consider both systems where the radiation flux is sub-Eddington throughout the gas column, and those where it is super-Eddington at the midplane but sub-Eddington in the atmosphere. However, the instability also produces a channel structure in which the radiation-matter interaction is reduced compared to time-steady analytic models because the radiation field is not fully trapped. We provide an approximation formula, appropriate for implementation in analytic models and non-radiative simulations, for the force exerted by the infrared radiation field in this regime.
On the Origin of Stellar Masses Krumholz, Mark R
Astrophysical journal/The Astrophysical journal,
12/2011, Volume:
743, Issue:
2
Journal Article
Peer reviewed
Open access
It has been a longstanding problem to determine, as far as possible, the characteristic masses of stars in terms of fundamental constants; the almost complete invariance of this mass as a function of ...the star-forming environment suggests that this should be possible. Here I provide such a calculation. The typical stellar mass is set by the characteristic fragment mass in a star-forming cloud, which depends on the cloud's density and temperature structure. Except in the very early universe, the latter is determined mainly by the radiation released as matter falls onto seed protostars. The energy yield from this process is ultimately set by the properties of deuterium burning in protostellar cores, which determines the stars' radii. I show that it is possible to combine these considerations to compute a characteristic stellar mass almost entirely in terms of fundamental constants, with an extremely weak residual dependence on the interstellar pressure and metallicity. This result not only explains the invariance of stellar masses, it resolves a second mystery: why fragmentation of a cold, low-density interstellar cloud, a process with no obvious dependence on the properties of nuclear reactions, happens to select a stellar mass scale such that stellar cores can ignite hydrogen. Finally, the weak residual dependence on the interstellar pressure and metallicity may explain recent observational hints of a smaller characteristic mass in the high-pressure, high-metallicity cores of giant elliptical galaxies.