26Al is a short-lived radioactive isotope thought to be injected into the interstellar medium (ISM) by massive stellar winds and supernovae. However, all-sky maps of 26Al emission show a distribution with a much larger scale height and faster rotation speed than either massive stars or the cold ISM. We investigate the origin of this discrepancy using an N-body+hydrodynamics simulation of a Milky-Way-like galaxy, self-consistently including self-gravity, star formation, stellar feedback, and 26Al production. We find no evidence that the Milky Way's spiral structure explains the 26Al anomaly. Stars and the 26Al bubbles they produce form along spiral arms, but, because our simulation produces material arms that arise spontaneously rather than propagating arms forced by an external potential, star formation occurs at arm centres rather than leading edges. As a result, we find a scale height and rotation speed for 26Al similar to that of the cold ISM. However, we also show that a synthetic 26Al emission map produced for a possible Solar position at the edge of a large 26Al bubble recovers many of the major qualitative features of the observed 26Al sky. This suggests that the observed anomalous 26Al distribution is the product of foreground emission from the 26Al produced by a nearby, recent supernova.

Magnetic fields play an important role in the dynamics of present-day molecular clouds. Recent work has shown that magnetic fields are equally important for primordial clouds, which form the first stars in the Universe. While the primordial magnetic field strength on cosmic scales is largely unconstrained, theoretical models strongly suggest that a weak seed field existed in the early Universe. We study how the amplification of such a weak field can influence the evolution of accretion discs around the first stars, and thus affect the primordial initial mass function (IMF). We perform a suite of 3D magneto-hydrodynamic (MHD) simulations with different initial field strengths and numerical resolutions. We find that, in simulations with sufficient spatial resolution to resolve the Jeans scale during the collapse, even initially weak magnetic fields grow exponentially to become dynamically important due to both the so-called \textit{small-scale turbulent dynamo} and the \textit{large-scale mean-field dynamo}. Capturing the small-scale dynamo action depends primarily on how well we resolve the Jeans length, while capturing the large-scale dynamo depends on the Jeans resolution as well as the maximum absolute resolution. Provided enough resolution, we find that fragmentation does not depend strongly on the initial field strength, because even weak fields grow to become strong. However, fragmentation in runs with magnetic fields differs significantly from those without magnetic fields. We conclude that the development of dynamically strong magnetic fields during the formation of the first stars is likely inevitable, and that these fields had a significant impact on the primordial IMF.

Cosmic rays (CRs) are a plausible mechanism for launching winds of cool material from the discs of star-forming galaxies. However, there is no consensus on what types of galaxies likely host CR-driven winds, or what role these winds might play in regulating galaxies' star formation rates. Using a detailed treatment of the transport and losses of hadronic CRs developed in the previous paper in this series, here we develop a semi-analytic model that allows us to assess the viability of using CRs to launch cool winds from galactic discs. In particular, we determine the critical CR fluxes -- and corresponding star formation rate surface densities -- above which hydrostatic equilibrium within a given galaxy is precluded because CRs drive the gas off in a wind or otherwise render it unstable. We show that, for star-forming galaxies with lower gas surface densities typical of the Galaxy and local dwarfs, the locus of this CR stability curve patrols the high side of the observed distribution of galaxies in the Kennicutt-Schmidt parameter space of star formation rate versus gas surface density. However, hadronic losses render CRs unable to drive winds in galaxies with higher surface densities. Our results show that quiescent, low surface density galaxies like the Milky Way are poised on the cusp of instability, such that small changes to ISM parameters can lead to the launching of CR-driven outflows, and we suggest that, as a result, CR feedback sets an ultimate limit to the star formation efficiency of most modern galaxies.

Magnetic fields play an important role for the formation of stars in both local and high-redshift galaxies. Recent studies of dynamo amplification in the first dark matter haloes suggest that significant magnetic fields were likely present during the formation of the first stars in the Universe at redshifts of 15 and above. In this work, we study how these magnetic fields potentially impact the initial mass function (IMF) of the first stars. We perform 200 high-resolution, three-dimensional (3D), magneto-hydrodynamic (MHD) simulations of the collapse of primordial clouds with different initial turbulent magnetic field strengths as predicted from turbulent dynamo theory in the early Universe, forming more than 1100 first stars in total. We detect a strong statistical signature of suppressed fragmentation in the presence of strong magnetic fields, leading to a dramatic reduction in the number of first stars with masses low enough that they might be expected to survive to the present day. Additionally, strong fields shift the transition point where stars go from being mostly single to mostly multiple to higher masses. However, irrespective of the field strength, individual simulations are highly chaotic, show different levels of fragmentation and clustering, and the outcome depends on the exact realisation of the turbulence in the primordial clouds. While these are still idealised simulations that do not start from cosmological initial conditions, our work shows that magnetic fields play a key role for the primordial IMF, potentially even more so than for the present-day IMF.

In the Milky Way, cosmic rays (CRs) are dynamically important in the interstellar medium, contribute to hydrostatic balance, and may help regulate star formation. However, we know far less about the importance of CRs in galaxies whose gas content or star formation rate differ significantly from those of the Milky Way. Here we construct self-consistent models for hadronic CR transport, losses, and contribution to pressure balance as a function of galaxy properties, covering a broad range of parameters from dwarfs to extreme starbursts. While the CR energy density increases from $\sim 1$ eV cm$^{-3}$ to $\sim 1$ keV cm$^{-3}$ over the range from sub-Milky Way dwarfs to bright starbursts, strong hadronic losses render CRs increasingly unimportant dynamically as the star formation rate surface density increases. In Milky Way-like systems, CR pressure is typically comparable to turbulent gas and magnetic pressure at the galactic midplane, but the ratio of CR pressure to gas pressure drops to $\sim 10^{-3}$ in dense starbursts. Galaxies also become increasingly CR calorimetric and gamma-ray bright in this limit. The degree of calorimetry at fixed galaxy properties is sensitive to the assumed model for CR transport, and in particular to the time CRs spend interacting with neutral ISM, where they undergo strong streaming losses. We also find that in some regimes of parameter space hydrostatic equilibrium discs cannot exist, and in Paper II of this series we use this result to derive a critical surface in the plane of star formation surface density and gas surface density beyond which CRs may drive large-scale galactic winds.

Stellar feedback in the form of radiation pressure and magnetically-driven collimated outflows may limit the maximum mass that a star can achieve and affect the star-formation efficiency of massive pre-stellar cores. Here we present a series of 3D adaptive mesh refinement radiation-magnetohydrodynamic simulations of the collapse of initially turbulent, massive pre-stellar cores. Our simulations include radiative feedback from both the direct stellar and dust-reprocessed radiation fields, and collimated outflow feedback from the accreting stars. We find that protostellar outflows punches holes in the dusty circumstellar gas along the star's polar directions, thereby increasing the size of optically thin regions through which radiation can escape. Precession of the outflows as the star's spin axis changes due to the turbulent accretion flow further broadens the outflow, and causes more material to be entrained. Additionally, the presence of magnetic fields in the entrained material leads to broader entrained outflows that escape the core. We compare the injected and entrained outflow properties and find that the entrained outflow mass is a factor of $\sim$3 larger than the injected mass and the momentum and energy contained in the entrained material are $\sim$25% and $\sim$5% of the injected momentum and energy, respectively. As a result, we find that, when one includes both outflows and radiation pressure, the former are a much more effective and important feedback mechanism, even for massive stars with significant radiative outputs.

We infer the intrinsic ionised gas kinematics for 383 star-forming galaxies across a range of integrated star-formation rates (SFR $\in [10^{-3}, 10^2]$ M$_\odot$ yr$^{-1}$) at $z \lesssim 0.1$ using a consistent 3D forward-modelling technique. The total sample is a combination of galaxies from the SAMI Galaxy Survey and DYNAMO survey. For typical low-$z$ galaxies taken from the SAMI Galaxy Survey, we find the vertical velocity dispersion ($\sigma_{v, z}$) to be positively correlated with measures of star-formation rate, stellar mass, HI gas mass, and rotational velocity. The greatest correlation is with star-formation rate surface density ($\Sigma_\text{SFR}$). Using the total sample, we find $\sigma_{v, z}$ increases slowly as a function of integrated star-formation rate in the range SFR $\in$ [$10^{-3}$, 1] M$_\odot$ yr$^{-1}$ from $17\pm3$ km s$^{-1}$ to $24\pm5$ km s$^{-1}$ followed by a steeper increase up to $\sigma_{v, z}$ $\sim 80$ km s$^{-1}$ for SFR $\gtrsim 1$ M$_\odot$ yr$^{-1}$. This is consistent with recent theoretical models that suggest a $\sigma_{v, z}$ floor driven by star-formation feedback processes with an upturn in $\sigma_{v, z}$ at higher SFR driven by gravitational transport of gas through the disc.

Understanding the evolution of the N/O ratio in the interstellar medium (ISM) of galaxies is essential if we are to complete our picture of the chemical evolution of galaxies at high redshift, since most observational calibrations of O/H implicitly depend upon the intrinsic N/O ratio. The observed N/O ratio, however, shows large scatter at low O/H, and is strongly dependent on galactic environment. We show that several heretofore unexplained features of the N/O distribution at low O/H can be explained by the N seen in metal-poor galaxies being mostly primary nitrogen that is returned to the ISM via pre-supernova winds from rapidly rotating massive stars ($M \gtrsim 10$ M$_\odot$, $v/v_{\rm crit} \gtrsim 0.4$). This mechanism naturally produces the observed N/O plateau at low O/H. We show that the large scatter in N/O at low O/H also arises naturally from variations in star-formation efficiency. By contrast, models in which the N and O come primarily from supernovae provide a very poor fit to the observed abundance distribution. We propose that the peculiar abundance patterns we observe at low O/H are a signature that dwarf galaxies retain little of their SN ejecta, leaving them with abundance patterns typical of winds.

The evolutionary paths taken by massive stars with $M \gtrsim 60 \, \mathrm{M}_\odot$ remain substantially uncertain. They begin their lives as main sequence (MS) O-stars. Depending on their masses, rotation rates, and metallicities, they can then encounter a wide range of evolutionary states with an equally broad set of possible surface compositions and spectral classifications. We present a new grid of calculations for the evolution of such stars that covers a broad range in mass, M/M$_\odot = 60$ to $150$, rotation rate, $v \, / \, v_{\rm crit} = 0$ to $0.6$, metallicity, $[\mathrm{Fe}/\mathrm{H}] = -4$ to $0$, and $\alpha$-element enhancement, $[\alpha/\mathrm{Fe}] = 0$ to $0.4$. We show that rotating stars undergo rotationally-induced dredge-up of nucleosynthetic products, mostly He and N, to their surfaces while still on the MS. Non-rotating metal-rich stars also reveal the products of nucleosynthesis on their surfaces because even modest amounts of mass loss expose their "fossil" convective cores: regions that are no longer convective, but which were part of the convective core at an early stage in the star's evolution. Thus surface enhancement of He and N is expected for rotating stars at all metallicities, and for non-rotating stars if they are relatively metal-rich. We calculate a stellar atmosphere for a representative model from our grid, properly accounting for He- and N-enhancement, and show that the resulting spectrum provides a good match to observed WNL stars, strongly suggesting that the physical mechanisms we have identified are the ultimate cause of the WNL phase.

Metallicity gradients are important diagnostics of galaxy evolution, because they record the history of events such as mergers, gas inflow and star-formation. However, the accuracy with which gradients can be measured is limited by spatial resolution and noise, and hence measurements need to be corrected for such effects. We use high resolution (~20 pc) simulation of a face-on Milky Way mass galaxy, coupled with photoionisation models, to produce a suite of synthetic high resolution integral field spectroscopy (IFS) datacubes. We then degrade the datacubes, with a range of realistic models for spatial resolution (2 to 16 beams per galaxy scale length) and noise, to investigate and quantify how well the input metallicity gradient can be recovered as a function of resolution and signal-to-noise ratio (SNR) with the intention to compare with modern IFS surveys like MaNGA and SAMI. Given appropriate propagation of uncertainties and pruning of low SNR pixels, we show that a resolution of 3-4 telescope beams per galaxy scale length is sufficient to recover the gradient to ~10-20% uncertainty. The uncertainty escalates to ~60% for lower resolution. Inclusion of the low SNR pixels causes the uncertainty in the inferred gradient to deteriorate. Our results can potentially inform future IFS surveys regarding the resolution and SNR required to achieve a desired accuracy in metallicity gradient measurements.

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, virialised velocity structures, and star formation histories that accelerate toward 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.

We obtained CO(2-1) observations of seven GRB hosts with the APEX and IRAM 30m telescopes. We analysed these data together with all other hosts with previous CO observations. We obtained detections for 3 GRB hosts (980425, 080207, and 111005A) and upper limits for the remaining 4 (031203, 060505, 060814, and 100316D). In our entire sample of 12 CO-observed GRB hosts, 3 are clearly deficient in molecular gas, even taking into account their metallicity (980425, 060814, and 080517). Four others are close to the best-fit line for other star-forming galaxies on the SFR-MH2 plot (051022, 060505, 080207, and 100316D). One host is clearly molecule rich (111005A). Finally, the data for 4 GRB hosts are not deep enough to judge whether they are molecule deficient (000418, 030329, 031203, and 090423). The median value of the molecular gas depletion time, MH2/SFR, of GRB hosts is ~0.3 dex below that of other star-forming galaxies, but this result has low statistical significance. A Kolmogorov-Smirnov test performed on MH2/SFR shows an only ~2sigma difference between GRB hosts and other galaxies. This difference can partly be explained by metallicity effects, since the significance decreases to ~1sigma for MH2/SFR versus~metallicity. We found that any molecular gas deficiency of GRB hosts has low statistical significance and that it can be attributed to their lower metallicities; and thus the sample of GRB hosts has molecular properties that are consistent with those of other galaxies, and they can be treated as representative star-forming galaxies. Given the concentration of atomic gas recently found close to GRB and supernova sites, indicating recent gas inflow, our results about the weak molecular deficiency imply that such an inflow does not enhance the SFRs significantly, or that atomic gas converts efficiently into the molecular phase, which fuels star formation.

Starburst galaxies are efficient $\gamma$-ray producers, because their high supernova rates generate copious cosmic ray (CR) protons, and their high gas densities act as thick targets off which these protons can produce neutral pions and thence $\gamma$-rays. In this paper we present a first-principles calculation of the mechanisms by which CRs propagate through such environments, combining astrochemical models with analysis of turbulence in weakly ionised plasma. We show that CRs cannot scatter off the strong large-scale turbulence found in starbursts, because efficient ion-neutral damping prevents such turbulence from cascading down to the scales of CR gyroradii. Instead, CRs stream along field lines at a rate determined by the competition between streaming instability and ion-neutral damping, leading to transport via a process of field line random walk. This results in an effective diffusion coefficient that is nearly energy-independent up to CR energies of $\sim 1$ TeV. We apply our computed diffusion coefficient to a simple model of CR escape and loss, and show that the resulting $\gamma$-ray spectra are in good agreement with the observed spectra of the starbursts NGC 253, M82, and Arp 220. In particular, our model reproduces these galaxies' relatively hard GeV $\gamma$-ray spectra and softer TeV spectra without the need for any fine-tuning of advective escape times or the shape of the CR injection spectrum.

We use the hydrodynamical simulation of our inner Galaxy presented in Armillotta et al. (2019) to study the gas distribution and kinematics within the CMZ. We use a resolution high enough to capture the gas emitting in dense molecular tracers such as NH3 and HCN, and simulate a time window of 50 Myr, long enough to capture phases during which the CMZ experiences both quiescent and intense star formation. We then post-process the simulated CMZ to calculate its spatially-dependent chemical and thermal state, producing synthetic emission data cubes and maps of both HI and the molecular gas tracers CO, NH3 and HCN. We show that, as viewed from Earth, gas in the CMZ is distributed mainly in two parallel and elongated features extending from positive longitudes and velocities to negative longitudes and velocities. The molecular gas emission within these two streams is not uniform, and it is mostly associated to the region where gas flowing towards the Galactic Center through the dust lanes collides with gas orbiting within the ring. Our simulated data cubes reproduce a number of features found in the observed CMZ. However, some discrepancies emerge when we use our results to interpret the position of individual molecular clouds. Finally, we show that, when the CMZ is near a period of intense star formation, the ring is mostly fragmented as a consequence of supernova feedback, and the bulk of the emission comes from star-forming molecular clouds. This correlation between morphology and star formation rate should be detectable in observations of extragalactic CMZs.

Transitional discs have central regions characterised by significant depletion of both dust and gas compared to younger, optically-thick discs. However, gas and dust are not depleted by equal amounts: gas surface densities are typically reduced by factors of $\sim 100$, but small dust grains are sometimes depleted by far larger factors, to the point of being undetectable. While this extreme dust depletion is often attributed to planet formation, in this paper we show that another physical mechanism is possible: expulsion of grains from the disc by radiation pressure. We explore this mechanism using 2D simulations of dust dynamics, simultaneously solving the equation of radiative transfer with the evolution equations for dust diffusion and advection under the combined effects of stellar radiation and hydrodynamic interaction with a turbulent, accreting background gas disc. We show that, in transition discs that are depleted in both gas and dust fraction by factors of $\sim 100-1000$ compared to minimum mass Solar nebular values, radiative clearing of any remaining $\sim 0.5$ $\mu$m and larger grains is both rapid and inevitable. The process is size-dependent, with smaller grains removed fastest and larger ones persisting for longer times. Our proposed mechanism thus naturally explains the extreme depletion of small grains commonly-found in transition discs. We further suggest that the dependence of this mechanism on grain size and optical properties may explain some of the unusual grain properties recently discovered in a number of transition discs. The simulation code we develop is freely available.

Using a 1D Lagrangian code specifically designed to assess the impact of multiple, time-resolved supernovae (SNe) from a single star cluster on the surrounding medium, we test three commonly used feedback recipes: delayed cooling (e.g., used in the GASOLINE-2 code), momentum-energy injection (a resolution-dependent transition between momentum-dominated feedback and energy-dominated feedback used, e.g., in the FIRE-2 code), and simultaneous energy injection (e.g., used in the EAGLE simulations). Our work provides an intermediary test for these recipes: we analyse a setting that is more complex than the simplified scenarios for which many were designed, but one more controlled than a full galactic simulation. In particular, we test how well these models reproduce the enhanced momentum efficiency seen for an 11 SN cluster simulated at high resolution (0.6 pc; a factor of 12 enhancement relative to the isolated SN case) when these subgrid recipes are implemented in low resolution (20 pc) runs. We find that: 1) the delayed cooling model performs well -- resulting in 9 times the momentum efficiency of the fiducial isolated SN value -- when SNe are clustered and $10^{51}$ erg are injected per SN, while clearly over-predicting the momentum efficiency in the single SN test case; 2) the momentum-energy model always achieves good results, with a factor of 5 boost in momentum efficiency; and 3) injecting the energy from all SNe simultaneously does little to prevent over-cooling and greatly under-produces the momentum deposited by clustered SNe, resulting in a factor of 3 decrease in momentum efficiency on the average.

Context. The diffuse gamma-ray emission of $^{26}{\rm Al}$ at 1.8 MeV reflects ongoing nucleosynthesis in the Milky Way, and traces massive-star feedback in the interstellar medium due to its 1 Myr radioactive lifetime. Interstellar-medium morphology and dynamics are investigated in astrophysics through 3D hydrodynamic simulations in fine detail, as only few suitable astronomical probes are available. Aims. We compare a galactic-scale hydrodynamic simulation of the Galaxy's interstellar medium, including feedback and nucleosynthesis, with gamma-ray data on $^{26}{\rm Al}$ emission in the Milky Way extracting constraints that are only weakly dependent on the particular realisation of the simulation or Galaxy structure. Methods. Due to constraints and biases in both the simulations and the gamma-ray observations, such comparisons are not straightforward. For a direct comparison, we perform maximum likelihood fits of simulated sky maps as well as observation-based maximum entropy maps to measurements with INTEGRAL/SPI. To study general morphological properties, we compare the scale heights of $^{26}{\rm Al}$ emission produced by the simulation to INTEGRAL/SPI measurements.} Results. The direct comparison shows that the simulation describes the observed inner Galaxy well, but differs significantly from the observed full-sky emission morphology. Comparing the scale height distribution, we see similarities for small scale height features and a mismatch at larger scale heights. We attribute this to the prominent foreground emission sites that are not captured by the simulation.

We present a study of the gas cycle and star formation history in the central 500 pc of the Milky Way, known as Central Molecular Zone (CMZ). Through hydrodynamical simulations of the inner 4.5 kpc of our Galaxy, we follow the gas cycle in a completely self-consistent way, starting from gas radial inflow due to the Galactic bar, the channelling of this gas into a dense, star-forming ring/stream at ~ 200 - 300 pc from the Galactic centre, and the launching of galactic outflows powered by stellar feedback. We find that star formation activity in the CMZ goes through oscillatory burst/quench cycles, with a period of tens to hundreds of Myr, characterised by roughly constant gas mass but order-of-magnitude level variations in the star formation rate. Comparison with the observed present-day star formation rate of the CMZ suggests that we are currently near a minimum of this cycle. Stellar feedback drives a mainly two-phase wind off the Galactic disc. The warm phase dominates the mass flux, and carries 100 - 200 % of the gas mass converted into stars. However, most of this gas goes into a fountain and falls back onto the disc rather than escaping the Galaxy. The hot phase carries most of the energy, with a time-averaged energy outflow rate of 10 - 20 % of the supernova energy budget.

The adiabatic index of H$_2\,$ ($\gamma_{\mathrm{H_2}}$) is non-constant at temperatures between $100-10^4\,\mathrm{K}$ due to the large energy spacing between its rotational and vibrational modes. For the formation of the first stars at redshifts 20 and above, this variation can be significant because primordial molecular clouds are in this temperature range due to the absence of efficient cooling by dust and metals. We study the possible importance of variations in $\gamma_{\mathrm{H_2}}$ for the primordial initial mass function by carrying out 80 3D gravito-hydrodynamic simulations of collapsing clouds with different random turbulent velocity fields, half using fixed $\gamma_{\rm H_2} = 7/5$ in the limit of classical diatomic gas (used in earlier works) and half using an accurate quantum mechanical treatment of $\gamma_{\mathrm{H_2}}$. We use the adaptive mesh refinement code FLASH with the primordial chemistry network from KROME for this study. The simulation suite produces almost 400 stars, with masses from $0.02 - 50$ M$_\odot$ (mean mass $\sim 10.5\,\mathrm{M_{\odot}}$ and mean multiplicity fraction $\sim 0.4$). While the results of individual simulations do differ when we change our treatment of $\gamma_{\mathrm{H_2}}$, we find no statistically significant differences in the overall mass or multiplicity distributions of the stars formed in the two sets of runs. We conclude that, at least prior to the onset of radiation feedback, approximating H$_2$ as a classical diatomic gas with $\gamma_{\rm H_2} = 7/5$ does not induce significant errors in simulations of the fragmentation of primordial gas. Nonetheless, we recommend using the accurate formulation of the H$_2\,$ adiabatic index in primordial star formation studies since it is not computationally more expensive and provides a better treatment of the thermodynamics.

Gaia DR2 provides an unprecedented sample of stars with full 6D phase-space measurements, creating the need for a self-consistent means of discovering and characterising the phase-space overdensities known as moving groups or associations. Here we present Chronostar, a new Bayesian analysis tool that meets this need. Chronostar uses the Expectation-Maximisation algorithm to remove the circular dependency between association membership lists and fits to their phase-space distributions, making it possible to discover unknown associations within a kinematic data set. It uses forward-modelling of orbits through the Galactic potential to overcome the problem of tracing backward stars whose kinematics have significant observational errors, thereby providing reliable ages. In tests using synthetic data sets with realistic measurement errors and complex initial distributions, Chronostar successfully recovers membership assignments and kinematic ages up to ~ 100 Myr. In tests on real stellar kinematic data in the phase-space vicinity of the ${\beta}$ Pictoris Moving Group, Chronostar successfully rediscovers the association without any human intervention, identifies 15 new likely members, corroborates 43 candidate members, and returns a kinematic age of 17.8+-1.2 Myr. In the process we also rediscover the Tucana-Horologium Moving Group, for which we obtain a kinematic age of 36.3+1.3-1.4 Myr.

Most gas in giant molecular clouds is relatively low-density and forms star inefficiently, converting only a small fraction of its mass to stars per dynamical time. However, star formation models generally predict the existence of a threshold density above which the process is efficient and most mass collapses to stars on a dynamical timescale. A number of authors have proposed observational techniques to search for a threshold density above which star formation is efficient, but it is unclear which of these techniques, if any, are reliable. In this paper we use detailed simulations of turbulent, magnetised star-forming clouds, including stellar radiation and outflow feedback, to investigate whether it is possible to recover star formation thresholds using current observational techniques. Using mock observations of the simulations at realistic resolutions, we show that plots of projected star formation efficiency per free-fall time $\epsilon_{\rm ff}$ can detect the presence of a threshold, but that the resolutions typical of current dust emission or absorption surveys are insufficient to determine its value. In contrast, proposed alternative diagnostics based on a change in the slope of the gas surface density versus star formation rate surface density (Kennicutt-Schmidt relation) or on the correlation between young stellar object counts and gas mass as a function of density are ineffective at detecting thresholds even when they are present. The signatures in these diagnostics sometimes taken as indicative of a threshold in observations, which we generally reproduce in our mock observations, do not prove to correspond to real physical features in the 3D gas distribution.

Direct comparisons between galaxy simulations and observations that both reach scales < 100 pc are strong tools to investigate the cloud-scale physics of star formation and feedback in nearby galaxies. Here we carry out such a comparison for hydrodynamical simulations of a Milky Way-like galaxy, including stochastic star formation, HII region and supernova feedback, and chemical post-processing at 8 pc resolution. Our simulation shows excellent agreement with almost all kpc-scale and larger observables, including total star formation rates, radial profiles of CO, HI, and star formation through the galactic disc, mass ratios of the ISM components, both whole-galaxy and resolved Kennicutt-Schmidt relations, and giant molecular cloud properties. However, we find that our simulation does not reproduce the observed de-correlation between tracers of gas and star formation on < 100 pc scales, known as the star formation 'uncertainty principle', which indicates that observed clouds undergo rapid evolutionary lifecycles. We conclude that the discrepancy is driven by insufficiently-strong pre-supernova feedback in our simulation, which does not disperse the surrounding gas completely, leaving star formation tracer emission too strongly associated with molecular gas tracer emission, inconsistent with observations. This result implies that the cloud-scale de-correlation of gas and star formation is a fundamental test for feedback prescriptions in galaxy simulations, one that can fail even in simulations that reproduce all other macroscopic properties of star-forming galaxies.

The distribution of gas-phase abundances in galaxy disks encodes the history of nucleosynthesis and transport through the interstellar medium (ISM) over cosmic time. Multi-object and high resolution integral-field spectroscopy have started to measure these distributions across hundreds of HII regions individually resolved at $\lesssim 100$ pc scales in a handful of objects, but in the coming decade these studies will expand to larger samples of galaxies. This will allow us to understand the role of feedback and turbulence in driving the mixing and diffusion of metals in the ISM, and statistically assess the role of galaxy environment and disk dynamics in modifying how mixing proceeds. Detailed searches for over- and under-enriched regions can address to what extent star formation is triggered by previous generations of star formation and by pristine and recycled gas flows. Local galaxies, for which these detailed measurements will be possible within the next decade, will inform the interpretation of integrated measurements at high-z, where very different dynamical gas-rich environments are found in early disk galaxies. Currently, progress in the field is severely hampered by the 0.2-0.3 dex level systematic uncertainties plaguing nebular abundance diagnostics. Improving our detailed understanding of ionized nebulae at $<$20 pc scales will help us find a solution to this problem, which will prove key to the study of metal enrichment and mixing across the galaxy population in the next decade.

We address a simple model where the Kennicutt-Schmidt (KS) relation between the macroscopic densities of star-formation rate (SFR, $\rho_{\rm sfr}$) and gas ($n$) in galactic discs emerges from self-regulation of the SFR via supernova feedback. It arises from the physics of supernova bubbles, insensitive to the microscopic SFR recipe and not explicitly dependent on gravity. The key is that the filling factor of SFR-suppressed supernova bubbles self-regulates to a constant, $f\sim 0.5$. Expressing the bubble fading radius and time in terms of $n$, the filling factor is $f \propto S\,n^{-s}$ with $s\sim 1.5$, where $S$ is the supernova rate density. A constant $f$ thus refers to $\rho_{\rm sfr} \propto n^{1.5}$, with a density-independent SFR efficiency per free-fall time $\sim 0.01$. The self-regulation to $f \sim 0.5$ and the convergence to a KS relation independent of the local SFR recipe are demonstrated in cosmological and isolated-galaxy simulations using different codes and recipes. In parallel, the spherical analysis of bubble evolution is generalized to clustered supernovae, analytically and via simulations, yielding $s \simeq 1.5 \pm 0.5$. An analysis of photo-ionized bubbles about pre-supernova stars yields a range of KS slopes but the KS relation is dominated by the supernova bubbles. Superbubble blowouts may lead to an alternative self-regulation by outflows and recycling. While the model is over-simplified, its simplicity and validity in the simulations may argue that it captures the origin of the KS relation.

Star-forming gas clouds are strongly magnetized, and their ionization fractions are high enough to place them close to the regime of ideal magnetohydrodyamics on all but the smallest size scales. In this review we discuss the effects of magnetic fields on the star formation rate (SFR) in these clouds, and on the mass spectrum of the fragments that are the outcome of the star formation process, the stellar initial mass function (IMF). Current numerical results suggest that magnetic fields by themselves are minor players in setting either the SFR or the IMF, changing star formation rates and median stellar masses only by factors of $\sim 2-3$ compared to non-magnetized flows. However, the indirect effects of magnetic fields, via their interaction with star formation feedback in the form of jets, photoionization, radiative heating, and supernovae, could have significantly larger effects. We explore evidence for this possibility in current simulations, and suggest avenues for future exploration, both in simulations and observations.

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 (HD) 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_\odot$. 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\times 10^5$ and $3\times 10^6$ $M_\odot$ km s$^{-1}$. A converged value for the simple homogeneous ISM remains elusive.

We present the star cluster catalogs for 17 dwarf and irregular galaxies in the $HST$ Treasury Program "Legacy ExtraGalactic UV Survey" (LEGUS). Cluster identification and photometry in this subsample are similar to that of the entire LEGUS sample, but special methods were developed to provide robust catalogs with accurate fluxes due to low cluster statistics. The colors and ages are largely consistent for two widely used aperture corrections, but a significant fraction of the clusters are more compact than the average training cluster. However, the ensemble luminosity, mass, and age distributions are consistent suggesting that the systematics between the two methods are less than the random errors. When compared with the clusters from previous dwarf galaxy samples, we find that the LEGUS catalogs are more complete and provide more accurate total fluxes. Combining all clusters into a composite dwarf galaxy, we find that the luminosity and mass functions can be described by a power law with the canonical index of $-2$ independent of age and global SFR binning. The age distribution declines as a power law, with an index of $\approx-0.80\pm0.15$, independent of cluster mass and global SFR binning. This decline of clusters is dominated by cluster disruption since the combined star formation histories and integrated-light SFRs are both approximately constant over the last few hundred Myr. Finally, we find little evidence for an upper-mass cutoff ($<2\sigma$) in the composite cluster mass function, and can rule out a truncation mass below $\approx10^{4.5}$M$_{\odot}$ but cannot rule out the existence of a truncation at higher masses.

Star clusters stand at the intersection of much of modern astrophysics: the interstellar medium, gravitational dynamics, stellar evolution, and cosmology. Here we review observations and theoretical models for the formation, evolution, and eventual disruption of star clusters. Current literature suggests a picture of this life cycle with several phases: (1) Clusters form in hierarchically-structured, accreting molecular clouds that convert gas into stars at a low rate per dynamical time until feedback disperses the gas. (2) The densest parts of the hierarchy resist gas removal long enough to reach high star formation efficiency, becoming dynamically-relaxed and well-mixed. These remain bound after gas removal. (3) In the first $\sim 100$ Myr after gas removal, clusters disperse moderately fast, through a combination of mass loss and tidal shocks by dense molecular structures in the star-forming environment. (4) After $\sim 100$ Myr, clusters lose mass via two-body relaxation and shocks by giant molecular clouds, processes that preferentially affect low-mass clusters and cause a turnover in the cluster mass function to appear on $\sim 1-10$ Gyr timescales. (5) Even after dispersal, some clusters remain coherent and thus detectable in chemical or action space for multiple galactic orbits. In the next decade a new generation of space- and AO-assisted ground-based telescopes will enable us to test and refine this picture.

Bathed in intense ionizing radiation, close-in gaseous planets undergo hydrodynamic atmospheric escape, which ejects the upper extent of their atmospheres into the interplanetary medium. Ultraviolet detections of escaping gas around transiting planets corroborate such a framework. Exposed to the stellar environment, the outflow is shaped by its interaction with the stellar wind and by the planet's orbit. We model these effects using Athena to perform 3-D radiative-hydrodynamic simulations of tidally-locked hydrogen atmospheres receiving large amounts of ionizing extreme-ultraviolet flux in various stellar environments for the low-magnetic-field case. Through a step-by-step exploration of orbital and stellar wind effects on the planetary outflow, we find three structurally distinct stellar wind regimes: weak, intermediate, and strong. We perform synthetic Lyman-$\alpha$ observations and find unique observational signatures for each regime. A weak stellar wind$\textrm{---}$which cannot confine the planetary outflow, leading to a torus of material around the star$\textrm{---}$has a pre-transit, red-shifted dayside arm and a slightly redward-skewed spectrum during transit. The intermediate regime truncates the dayside outflow at large distances from the planet and causes periodic disruptions of the outflow, producing observational signatures that mimic a double transit. The first of these dips is blue-shifted and precedes the optical transit. Finally, strong stellar winds completely confine the outflow into a cometary tail and accelerate the outflow outwards, producing large blue-shifted signals post-transit. Across all three regimes, large signals occur far outside of transit, offering motivation to continue ultraviolet observations outside of direct transit.

We use ALMA and IRAM 30-m telescope data to investigate the relationship between the spectroscopically-traced dense gas fraction and the cloud-scale (120 pc) molecular gas surface density in five nearby, star-forming galaxies. We estimate the dense gas mass fraction at 650 pc and 2800 pc scales using the ratio of HCN (1-0) to CO (1-0) emission. We then use high resolution (120 pc) CO (2-1) maps to calculate the mass-weighted average molecular gas surface density within 650 pc or 2770 pc beam where the dense gas fraction is estimated. On average, the dense gas fraction correlates with the mass-weighted average molecular gas surface density. Thus, parts of a galaxy with higher mean cloud-scale gas surface density also appear to have a larger fraction of dense gas. The normalization and slope of the correlation do vary from galaxy to galaxy and with the size of the regions studied. This correlation is consistent with a scenario where the large-scale environment sets the gas volume density distribution, and this distribution manifests in both the cloud-scale surface density and the dense gas mass fraction.

The advent of new sub-millimeter observational facilities has stimulated the desire to model the sub-mm line emission of galaxies within cosmological galaxy formation models. This is typically done in post-processing by applying sub-resolution recipes to describe the properties of the unresolved interstellar medium. At the same time, while there is freedom in how one implements these sorts of recipes, the impact of various choices has yet to be systematically explored in simulations. In this paper, we do just that. We combine a semi-analytic model of galaxy formation with chemical equilibrium networks and numerical radiative transfer models and explore how different choices for the sub-grid modeling affect the predicted CO, [CI], and [CII] emission of galaxies. We find that a key component for a successful model includes a molecular cloud mass-size relation and scaling for the ultraviolet and cosmic ray radiation field that varies based on local ISM properties. Our most successful model adopts a Plummer profile for the radial density distribution of gas within molecular clouds. On the other hand, different assumptions for the clumping of gas within molecular clouds and changes in the slope of the molecular cloud mass distribution function hardly affect the CO, [CI], and [CII] luminosities of galaxies. At fixed star-formation rate the [CII]-SFR ratio of galaxies scales inversely with the pressure acting on molecular clouds, increasing the molecular clouds density and hence decreasing the importance of [CII] line cooling. Overall we find that it is essential that a wide range of sub-mm emission lines arising in vastly different phases of the ISM are used as model constraints in order to limit the freedom in sub-grid choices and we present the details of the most successful model variant.

We present a flexible, detailed model for the evolution of galactic discs in a cosmological context since $z\approx 4$, including a physically-motivated model for radial transport of gas and stars within galactic discs. This expansion beyond traditional semi-analytic models that do not include radial structure, or include only a prescribed radial structure, enables us to study the internal structure of disc galaxies and the processes that drive it. In order to efficiently explore the large parameter space allowed by this model, we construct a neural network-based emulator that can quickly return a reasonable approximation for many observables we can extract from the model, e.g. the star formation rate or the half mass stellar radius, at different redshifts. We employ the emulator to constrain the model parameters with Bayesian inference by comparing its predictions to 11 observed galaxy scaling relations at a variety of redshifts. The constrained models agree well with observations, both those used to fit the data and those not included in the fitting procedure. These models will be useful theoretical tools for understanding the increasingly detailed observational datasets from IFUs.

The $\rm H_2$ mass of molecular clouds has traditionally been traced by the CO(J=1-0) rotational transition line. This said, CO is relatively easily photodissociated, and can also be destroyed by cosmic rays, thus rendering some fraction of molecular gas to be "CO-dark". We investigate the amount and physical properties of CO-dark gas in two $z \sim 0$ disc galaxies, and develop predictions for the expected intensities of promising alternative tracers ([CI 609 $\mu$m and [CII] 158 $\mu$m emission). We do this by combining cosmological zoom simulations of disc galaxies with thermal-radiative-chemical equilibrium interstellar medium (ISM) calculations to model the predicted H~\textsc{i} and $\rm H_2$ abundances and CO(J=1-0), [CI] 609 $\mu$m and [CII] 158 $\mu$m emission properties. Our model treats the ISM as a collection of radially stratified clouds whose properties are dictated by their volume and column densities, the gas-phase metallicity, and the interstellar radiation field and cosmic ray ionization rates. Our main results follow. Adopting an observationally motivated definition of CO-dark gas, i.e. $\rm H_2$ gas with $W_{\rm CO} < 0.1 $ K-km/s, we find that a significant amount ($> 50\%$) of the total $\rm H_2$ mass lies in CO-dark gas, most of which is diffuse gas, poorly shielded due to low dust column density. The CO-dark molecular gas tends to be dominated by [CII], though [CI] also serves as a bright tracer of the dark gas in many instances. At the same time, [CII] also tends to trace neutral atomic gas. As a result, when we quantify the conversion factors for the three carbon-based tracers of molecular gas, we find that [CI] suffers the least contamination from diffuse atomic gas, and is relatively insensitive to secondary parameters.

We recently presented a new statistical method to constrain the physics of star formation and feedback on the cloud scale. Fundamentally, this new method is only able to recover the relative durations of different evolutionary phases. It therefore requires a "reference time-scale" to convert the relative time-scales into absolute values. The phase over which the star formation rate (SFR) tracer is visible is the logical choice to assign this reference time-scale to, since it can be characterised using stellar population synthesis models. In this paper, we use hydrodynamical disc galaxy simulations to produce synthetic emission maps of several SFR tracers and apply our statistical method to measure the associated characteristic time-scale of each tracer. These cover 12 ultraviolet (UV) filters (from GALEX, Swift, and HST) covering a wavelength range 150-350 nm, as well as H$\alpha$. For solar-metallicity environments, we find the characteristic time-scales for H$\alpha$ with (without) continuum subtraction to be ~4.3 (6-16) Myr and for the UV filters to be in the range 17-33 Myr, monotonically increasing with wavelength. We find that the characteristic time-scale decreases towards higher metallicities, as well as to lower star formation rate surface densities, if the stellar initial mass function is not well-sampled. We provide fitting functions to the resulting reference time-scales to facilitate observational applications of our statistical method across a wide range of galactic environments. More generally, our results predict the time-scales over which photoionisation and UV heating take place around star-forming regions.

We describe a novel method for determining the demographics of a population of star clusters, for example distributions of cluster mass and age, from unresolved photometry. This method has a number of desirable properties: it fully exploits all the information available in a data set without any binning, correctly accounts for both measurement error and sample incompleteness, naturally handles heterogenous data (for example fields that have been imaged with different sets of filters or to different depths), marginalises over uncertain extinctions, and returns the full posterior distributions of the parameters describing star cluster demographics. We demonstrate the method using mock star cluster catalogs and show that our method is robust and accurate, and that it can recover the demographics of star cluster populations significantly better than traditional fitting methods. For realistic sample sizes, our method is sufficiently powerful that its accuracy is ultimately limited by the accuracy of the underlying physical models for stellar evolution and interstellar dust, rather than by statistical uncertainties. Our method is implemented as part of the Stochastically Lighting Up Galaxies (SLUG) stellar populations code, and is freely available.

The large columns of dusty gas enshrouding and fuelling star-formation in young, massive stellar clusters may render such systems optically thick to radiation well into the infrared. This raises the prospect that both "direct" radiation pressure produced by absorption of photons leaving stellar surfaces and "indirect" radiation pressure from photons absorbed and then re-emitted by dust grains may be important sources of feedback in such systems. Here we evaluate this possibility by deriving the conditions under which a spheroidal, self-gravitating, mixed gas-star cloud can avoid catastrophic disruption by the combined effects of direct and indirect radiation pressure. We show that radiation pressure sets a maximum star cluster formation efficiency of $\epsilon_{\rm max} \sim 0.9$ at a (very large) gas surface density of $\sim 10^5 M_\odot$ pc$^{-2} (Z_\odot/Z) \simeq 20$ g cm$^{-2} (Z_\odot/Z)$, but that gas clouds above this limit undergo significant radiation-driven expansion during star formation, leading to a maximum stellar surface density very near this value for all star clusters. Data on the central surface mass density of compact stellar systems, while sparse and partly confused by dynamical effects, are broadly consistent with the existence of a metallicity-dependent upper-limit comparable to this value. Our results imply that this limit may preclude the formation of the progenitors of intermediate-mass black holes for systems with $Z \gtrsim 0.2 Z_\odot$.

HCN(1$-$0) emission traces dense gas and correlates very strongly with star formation rates (SFRs) on scales from small Milky Way clouds to whole galaxies. The observed correlation offers strong constraints on the efficiency of star formation in dense gas, but quantitative interpretation of this constraint requires a mapping from HCN emission to gas mass and density. In this paper we provide the required calibration by postprocessing high-resolution simulations of dense, star-forming clouds to calculate their HCN emission ($L_{\rm HCN}$) and to determine how that emission is related to the underlying gas density distribution and star formation efficiency. We find that HCN emission traces gas with a luminosity-weighted mean number density of $0.8-1.7 \times 10^4\,{\rm cm}^{-3}$ and that HCN luminosity is related to mass of dense gas of $\gtrsim 10^4\,{\rm cm}^{-3}$ with a conversion factor of ${\alpha}_{\rm HCN} \approx 14\,\rm M_{\odot}/(K\,km\,s^{-1}\,{pc}^2)$. We also measure a new empirical relationship between the SFR per global mean freefall time (${\epsilon}_{\rm ff}$) and the SFR$-$HCN relationship, ${\rm SFR}/L_{\rm HCN} = 2.0 \times 10^{-7}\,({\epsilon}_{\rm ff}/0.01)^{1.1}\,\rm M_{\odot}\,{yr}^{-1}/(K\,km\,s^{-1}\,{pc}^2)$. The observed SFR$-$HCN correlation strongly constrains ${\epsilon}_{\rm ff} \approx 1\%$ with a factor of $\sim 3$ systematic uncertainty. The scatter in ${\epsilon}_{\rm ff}$ from cloud to cloud within the Milky Way is a factor of a few. We conclude that $L_{\rm HCN}$ is an effective tracer of dense gas and that the IR$-$HCN correlation is a significant diagnostic of the microphysics of star formation in dense gas.

We present an investigation of the relationship between giant molecular cloud (GMC) properties and the associated stellar clusters in the nearby flocculent galaxy NGC 7793. We combine the star cluster catalog from the HST LEGUS (Legacy ExtraGalactic UV Survey) program with the 15 parsec resolution ALMA CO(2-1) observations. We find a strong spatial correlation between young star clusters and GMCs such that all clusters still associated with a GMC are younger than 11 Myr and display a median age of 2 Myr. The age distribution increases gradually as the cluster-GMC distance increases, with star clusters that are spatially unassociated with molecular gas exhibiting a median age of 7 Myr. Thus, star clusters are able to emerge from their natal clouds long before the timescale required for clouds to disperse. To investigate if the hierarchy observed in the stellar components is inherited from the GMCs, we quantify the amount of clustering in the spatial distributions of the components and find that the star clusters have a fractal dimension slope of $-0.35 \pm 0.03$, significantly more clustered than the molecular cloud hierarchy with slope of $-0.18 \pm 0.04$ over the range 40-800 pc. We find, however, that the spatial clustering becomes comparable in strength for GMCs and star clusters with slopes of $-0.44\pm0.03$ and $-0.45\pm0.06$ respectively, when we compare massive ($>$10$^5$ M$_{\odot}$) GMCs to massive and young star clusters. This shows that massive star clusters trace the same hierarchy as their parent GMCs, under the assumption that the star formation efficiency is a few percent.

Meteoritic evidence shows that the Solar system at birth contained significant quantities of short-lived radioisotopes (SLRs) such as 60Fe and 26Al (with half-lives of 2.6 and 0.7 Myr respectively) produced in supernova explosions and in the Wolf-Rayet winds that precede them. Proposed explanations for the high SLR abundance include formation of the Sun in a supernova-triggered collapse or in a giant molecular cloud (GMC) that was massive enough to survive multiple supernovae (SNe) and confine their ejecta. However, the former scenario is possible only if the Sun is a rare outlier among massive stars, while the latter appears to be inconsistent with the observation that 26Al is distributed with a scale height significantly larger than GMCs. In this paper, we present a high-resolution chemo-hydrodynamical simulation of the entire Milky-Way Galaxy, including stochastic star formation, HII regions, SNe, and element injection, that allows us to measure for the distribution of 60Fe/56Fe and 26Al/27Al ratios over all stars in the Galaxy. We show that the Solar System's abundance ratios are well within the normal range, but that SLRs originate neither from triggering nor from confinement in long-lived clouds as previously conjectured. Instead, we find that SLRs are abundant in newborn stars because star formation is correlated on galactic scales, so that ejecta preferentially enrich atomic gas that will subsequently be accreted onto existing GMCs or will form new ones. Thus new generations of stars preferentially form in patches of the Galaxy contaminated by previous generations of stellar winds and supernovae.

We explore the possibility of using UV spectroscopy in combination with broad-band photometry as diagnostic tools for understanding the shape of the Initial Mass Function (IMF) in unresolved stellar populations. Building on our previous work, we extend the Stochastically Lighting Up Galaxies code (slug) to include a high-resolution UV spectral synthesiser and equivalent width calculation capabilities. We first gain a qualitative understanding of how UV spectral features behave as the parameters that define a star cluster in slug (mass, age, extinction, and IMF slope alpha) are changed. We then exploit Bayesian inference techniques to recover the alpha values for clusters simulated with slug, using mock observations of these clusters comprised of broad-band photometry and equivalent width measurements of a selection of UV spectral features. We find some improvement when compared to attempts using broad-band photometry alone (with the interquartile range of the alpha posterior PDF shrinking by 32%), although we still do not yet fully break the known degeneracy between the cluster mass and alpha. Finally, we make predictions about how effective real observations will be by quantifying our ability to constrain alpha as a function of limiting equivalent width. We find that observations sensitive to a modest equivalent width of 9 A are sufficient to improve the recovery of the IMF slope parameter by 32% (interquartile range of posterior PDF median residuals), moving to 39% when we include all the significant spectral features in the wavelength range 900-3000 A.

In recent years a number of authors have introduced methods to model the effects of radiation pressure feedback on flows of interstellar and intergalactic gas, and have posited that the forces exerted by stars' radiation output represents an important feedback mechanism capable of halting accretion and thereby regulating star formation. However, numerical simulations have reached widely varying conclusions about the effectiveness of this feedback. In this paper I show that much of the divergence in the literature is a result of failure to obey an important resolution criterion: whether radiation feedback is able to reverse an accretion flow is determined on scales comparable to the dust destruction radius, which is $\lesssim 1000$ AU even for the most luminous stellar sources. Simulations that fail to resolve this scale can produce unphysical results, in many cases leading to a dramatic overestimate of the effectiveness of radiation feedback. Most published simulations of radiation feedback on molecular cloud and galactic scales fail to satisfy this condition. I show how the problem can be circumvented by introducing a new subgrid model that explicitly accounts for momentum balance on unresolved scales, making it possible to simulate dusty accretion flows safely even at low resolution.

Ongoing surveys are in the process of measuring the chemical abundances in large numbers of stars, with the ultimate goal of reconstructing the formation history of the Milky Way using abundances as tracers. However, interpretation of these data requires that we understand the relationship between stellar distributions in chemical and physical space, i.e., how similar in chemical abundance do we expect a pair of stars to be as a function of the distance between their formation sites. We investigate this question by simulating the gravitational collapse of a turbulent molecular cloud extracted from a galaxy-scale simulation, seeded with chemical inhomogeneities with different initial spatial scales. We follow the collapse from galactic scales down to resolutions scales of $\approx 10^{-3}$ pc, and find that, during this process, turbulence mixes the metal patterns, reducing the abundance scatter initially present in the gas by an amount that depends on the initial scale of inhomogeneity of each metal field. However, we find that regardless of the initial spatial structure of the metals at the onset of collapse, the final stellar abundances are highly correlated on distances below a few pc, and nearly uncorrelated on larger distances. Consequently, the star formation process defines a natural size scale of $\sim 1$ pc for chemically-homogenous star clusters, suggesting that any clusters identified as homogenous in chemical space must have formed within $\sim 1$ pc of one another. However, in order to distinguish different star clusters in chemical space, observations across multiple elements will be required, and the elements that are likely to be most efficient at separating distinct clusters in chemical space are those whose correlation length in the ISM is of order tens of pc, comparable to the sizes of individual molecular clouds.

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.

We apply an analytic theory for environmentally-dependent molecular cloud lifetimes to the Central Molecular Zone of the Milky Way. Within this theory, the cloud lifetime in the Galactic centre is obtained by combining the time-scales for gravitational instability, galactic shear, epicyclic perturbations and cloud-cloud collisions. We find that at galactocentric radii $\sim 45$-$120$ pc, corresponding to the location of the '100-pc stream', cloud evolution is primarily dominated by gravitational collapse, with median cloud lifetimes between 1.4 and 3.9 Myr. At all other galactocentric radii, galactic shear dominates the cloud lifecycle, and we predict that molecular clouds are dispersed on time-scales between 3 and 9 Myr, without a significant degree of star formation. Along the outer edge of the 100-pc stream, between radii of 100 and 120 pc, the time-scales for epicyclic perturbations and gravitational free-fall are similar. This similarity of time-scales lends support to the hypothesis that, depending on the orbital geometry and timing of the orbital phase, cloud collapse and star formation in the 100-pc stream may be triggered by a tidal compression at pericentre. Based on the derived time-scales, this should happen in approximately 20 per cent of all accretion events onto the 100-pc stream.

The importance of radiation pressure feedback in galaxy formation has been extensively debated over the last decade. The regime of greatest uncertainty is in the most actively star-forming galaxies, where large dust columns can potentially produce a dust-reprocessed infrared radiation field with enough pressure to drive turbulence or eject material. Here we derive the conditions under which a self-gravitating, mixed gas-star disc can remain hydrostatic despite trapped radiation pressure. Consistently taking into account the self-gravity of the medium, the star- and dust-to-gas ratios, and the effects of turbulent motions not driven by radiation, we show that galaxies can achieve a maximum Eddington-limited star formation rate per unit area $\dot{\Sigma}_{\rm *,crit} \sim 10^3 M_{\odot}$ pc$^{-2}$ Myr$^{-1}$, corresponding to a critical flux of $F_{\rm *,crit} \sim 10^{13} L_{\odot}$ kpc$^{-2}$ similar to previous estimates; higher fluxes eject mass in bulk, halting further star formation. Conversely, we show that in galaxies below this limit, our one-dimensional models imply simple vertical hydrostatic equilibrium and that radiation pressure is ineffective at driving turbulence or ejecting matter. Because the vast majority of star-forming galaxies lie below the maximum limit for typical dust-to-gas ratios, we conclude that infrared radiation pressure is likely unimportant for all but the most extreme systems on galaxy-wide scales. Thus, while radiation pressure does not explain the Kennicutt-Schmidt relation, it does impose an upper truncation on it. Our predicted truncation is in good agreement with the highest observed gas and star formation rate surface densities found both locally and at high redshift.

We introduce a new model for the structure and evolution of the gas in galactic discs. In the model the gas is in vertical pressure and energy balance. Star formation feedback injects energy and momentum, and non-axisymmetric torques prevent the gas from becoming more than marginally gravitationally unstable. From these assumptions we derive the relationship between galaxies' bulk properties (gas surface density, stellar content, and rotation curve) and their star formation rates, gas velocity dispersions, and rates of radial inflow. We show that the turbulence in discs can be powered primarily by star formation feedback, radial transport, or a combination of the two. In contrast to models that omit either radial transport or star formation feedback, the predictions of this model yield excellent agreement with a wide range of observations, including the star formation law measured in both spatially resolved and unresolved data, the correlation between galaxies' star formation rates and velocity dispersions, and observed rates of radial inflow. The agreement holds across a wide range of galaxy mass and type, from local dwarfs to extreme starbursts to high-redshifts discs. We apply the model to galaxies on the star-forming main sequence, and show that it predicts a transition from mostly gravity-driven turbulence at high redshift to star formation-driven turbulence at low redshift. This transition, and the changes in mass transport rates that it produces, naturally explain why galaxy bulges tend to form at high redshift and discs at lower redshift, and why galaxies tend to quench inside-out.

We use new ALMA observations to investigate the connection between dense gas fraction, star formation rate, and local environment across the inner region of four local galaxies showing a wide range of molecular gas depletion times. We map HCN (1-0), HCO$^+$ (1-0), CS (2-1), $^{13}$CO (1-0), and C$^{18}$O (1-0) across the inner few kpc of each target. We combine these data with short spacing information from the IRAM large program EMPIRE, archival CO maps, tracers of stellar structure and recent star formation, and recent HCN surveys by Bigiel et al. and Usero et al. We test the degree to which changes in the dense gas fraction drive changes in the SFR. $I_{HCN}/I_{CO}$ (tracing the dense gas fraction) correlates strongly with $I_{CO}$ (tracing molecular gas surface density), stellar surface density, and dynamical equilibrium pressure, $P_{DE}$. Therefore, $I_{HCN}/I_{CO}$ becomes very low and HCN becomes very faint at large galactocentric radii, where ratios as low as $I_{HCN}/I_{CO} \sim 0.01$ become common. The apparent ability of dense gas to form stars, $\Sigma_{SFR}/\Sigma_{dense}$ (where $\Sigma_{dense}$ is traced by the HCN intensity and the star formation rate is traced by a combination of H$\alpha$ and 24$\mu$m emission), also depends on environment. $\Sigma_{SFR}/\Sigma_{dense}$ decreases in regions of high gas surface density, high stellar surface density, and high $P_{DE}$. Statistically, these correlations between environment and both $\Sigma_{SFR}/\Sigma_{dense}$ and $I_{HCN}/I_{CO}$ are stronger than that between apparent dense gas fraction ($I_{HCN}/I_{CO}$) and the apparent molecular gas star formation efficiency $\Sigma_{SFR}/\Sigma_{mol}$. We show that these results are not specific to HCN.

As part of our long-term campaign to understand how cold streams feed massive galaxies at high redshift, we study the Kelvin-Helmholtz instability (KHI) of a supersonic, cold, dense gas stream as it penetrates through a hot, dilute circumgalactic medium (CGM). A linear analysis (Paper I) showed that, for realistic conditions, KHI may produce nonlinear perturbations to the stream during infall. Therefore, we proceed here to study the nonlinear stage of KHI, still limited to a two-dimensional slab with no radiative cooling or gravity. Using analytic models and numerical simulations, we examine stream breakup, deceleration and heating via surface modes and body modes. The relevant parameters are the density contrast between stream and CGM ($\delta$), the Mach number of the stream velocity with respect to the CGM ($M_{\rm b}$) and the stream radius relative to the halo virial radius ($R_{\rm s}/R_{\rm v}$). We find that sufficiently thin streams disintegrate prior to reaching the central galaxy. The condition for breakup ranges from $R_{\rm s} < 0.03 R_{\rm v}$ for $(M_{\rm b} \sim 0.75, \delta \sim 10)$ to $R_{\rm s} < 0.003 R_{\rm v}$ for $(M_{\rm b} \sim 2.25, \delta \sim 100)$. However, due to the large stream inertia, KHI has only a small effect on the stream inflow rate and a small contribution to heating and subsequent Lyman-$\alpha$ cooling emission.

We present a study of the dust-to-gas ratios in five nearby galaxies NGC 628 (M74), NGC 6503, NGC 7793, UGC 5139 (Holmberg I), and UGC 4305 (Holmberg II). Using Hubble Space Telescope broad band WFC3/UVIS UV and optical images from the Treasury program LEGUS (Legacy ExtraGalactic UV Survey) combined with archival HST/ACS data, we correct thousands of individual stars for extinction across these five galaxies using an isochrone-matching (reddening-free Q) method. We generate extinction maps for each galaxy from the individual stellar extinctions using both adaptive and fixed resolution techniques, and correlate these maps with neutral HI and CO gas maps from literature, including The HI Nearby Galaxy Survey (THINGS) and the HERA CO-Line Extragalactic Survey (HERACLES). We calculate dust-to-gas ratios and investigate variations in the dust-to-gas ratio with galaxy metallicity. We find a power law relationship between dust-to-gas ratio and metallicity, consistent with other studies of dust-to-gas ratio compared to metallicity. We find a change in the relation when H$_2$ is not included. This implies that underestimation of $N_{H_2}$ in low-metallicity dwarfs from a too-low CO-to-H$_2$ conversion factor $X_{CO}$ could have produced too low a slope in the derived relationship between dust-to-gas ratio and metallicity. We also compare our extinctions to those derived from fitting the spectral energy distribution (SED) using the Bayesian Extinction and Stellar Tool (BEAST) for NGC 7793 and find systematically lower extinctions from SED-fitting as compared to isochrone matching.

We investigate the spatially-resolved star formation relation using a galactic disk formed in a comprehensive high-resolution (3.8 pc) simulation. Our new implementation of stellar feedback includes ionizing radiation as well as supernova explosions, and we handle ionizing radiation by solving the radiative transfer equation rather than by a subgrid model. Photoheating by stellar radiation stabilizes gas against Jeans fragmentation, reducing the star formation rate. Because we have self-consistently calculated the location of ionized gas, we are able to make spatially-resolved mock observations of star formation tracers, such as H-alpha emission. We can also observe how stellar feedback manifests itself in the correlation between ionized and molecular gas. Applying our techniques to the disk in a galactic halo of 2.3e11 Msun, we find that the correlation between star formation rate density (estimated from mock H-alpha emission) and molecular hydrogen density shows large scatter, especially at high resolutions of <~ 75 pc that are comparable to the size of giant molecular clouds (GMCs). This is because an aperture of GMC size captures only particular stages of GMC evolution, and because H-alpha traces hot gas around star-forming regions and is displaced from the molecular hydrogen peaks themselves. By examining the evolving environment around star clusters, we speculate that the breakdown of the traditional star formation laws of the Kennicutt-Schmidt type at small scales is further aided by a combination of stars drifting from their birthplaces, and molecular clouds being dispersed via stellar feedback.

We describe a new method for simulating ionizing radiation and supernova feedback in the analogues of low-redshift galactic disks. In this method, which we call star-forming molecular cloud (SFMC) particles, we use a ray-tracing technique to solve the radiative transfer equation for ultraviolet photons emitted by thousands of distinct particles on the fly. Joined with high numerical resolution of 3.8 pc, the realistic description of stellar feedback helps to self-regulate star formation. This new feedback scheme also enables us to study the escape of ionizing photons from star-forming clumps and from a galaxy, and to examine the evolving environment of star-forming gas clumps. By simulating a galactic disk in a halo of 2.3e11 Msun, we find that the average escape fraction from all radiating sources on the spiral arms (excluding the central 2.5 kpc) fluctuates between 0.08% and 5.9% during a ~20 Myr period with a mean value of 1.1%. The flux of escaped photons from these sources is not strongly beamed, but manifests a large opening angle of more than 60 degree from the galactic pole. Further, we investigate the escape fraction per SFMC particle, f_esc(i), and how it evolves as the particle ages. We discover that the average escape fraction f_esc is dominated by a small number of SFMC particles with high f_esc(i). On average, the escape fraction from a SFMC particle rises from 0.27% at its birth to 2.1% at the end of a particle lifetime, 6 Myrs. This is because SFMC particles drift away from the dense gas clumps in which they were born, and because the gas around the star-forming clumps is dispersed by ionizing radiation and supernova feedback. The framework established in this study brings deeper insight into the physics of photon escape fraction from an individual star-forming clump, and from a galactic disk.

Interstellar gas clouds are often both highly magnetized and supersonically turbulent, with velocity dispersions set by a competition between driving and dissipation. This balance has been studied extensively in the context of gases with constant mean density. However, many astrophysical systems are contracting under the influence of external pressure or gravity, and the balance between driving and dissipation in a contracting, magnetized medium has yet to be studied. In this paper we present three-dimensional (3D) magnetohydrodynamic (MHD) simulations of compression in a turbulent, magnetized medium that resembles the physical conditions inside molecular clouds. We find that in some circumstances the combination of compression and magnetic fields leads to a rate of turbulent dissipation far less than that observed in non-magnetized gas, or in non-compressing magnetized gas. As a result, a compressing, magnetized gas reaches an equilibrium velocity dispersion much greater than would be expected for either the hydrodynamic or the non-compressing case. We use the simulation results to construct an analytic model that gives an effective equation of state for a coarse-grained parcel of the gas, in the form of an ideal equation of state with a polytropic index that depends on the dissipation and energy transfer rates between the magnetic and turbulent components. We argue that the reduced dissipation rate and larger equilibrium velocity dispersion produced by compressing, magnetized turbulence has important implications for the driving and maintenance of turbulence in molecular clouds, and for the rates of chemical and radiative processes that are sensitive to shocks and dissipation.

The distributions of a galaxy's gas and stars in chemical space encodes a tremendous amount of information about that galaxy's physical properties and assembly history. However, present methods for extracting information from chemical distributions are based either on coarse averages measured over galactic scales (e.g., metallicity gradients) or on searching for clusters in chemical space that can be identified with individual star clusters or gas clouds on $\sim 1$ pc scales. These approaches discard most of the information, because in galaxies gas and young stars are observed to be distributed fractally, with correlations on all scales, and the same is likely to be true of metals. In this paper we introduce a first theoretical model, based on stochastically-forced diffusion, capable of predicting the multi-scale statistics of metal fields. We derive the variance, correlation function, and power spectrum of the metal distribution from first principles, and determine how these quantities depend on elements' astrophysical origin sites and on the large-scale properties of galaxies. Among other results, we explain for the first time why the typical abundance scatter observed in the interstellar media of nearby galaxies is $\approx 0.1$ dex, and we predict that this scatter will be correlated on spatial scales of $\sim 0.5-1$ kpc, and over time scales of $\sim 100-300$ Myr. We discuss the implications of our results for future chemical tagging studies.

We present a large suite of simulations of the formation of low-mass star clusters. Our simulations include an extensive set of physical processes -- magnetohydrodynamics, radiative transfer, and protostellar outflows -- and span a wide range of virial parameters and magnetic field strengths. Comparing the outcomes of our simulations to observations, we find that simulations remaining close to virial balance throughout their history produce star formation efficiencies and initial mass function (IMF) peaks that are stable in time and in reasonable agreement with observations. Our results indicate that small-scale dissipation effects near the protostellar surface provide a feedback loop for stabilizing the star formation efficiency. This is true regardless of whether the balance is maintained by input of energy from large scale forcing or by strong magnetic fields that inhibit collapse. In contrast, simulations that leave virial balance and undergo runaway collapse form stars too efficiently and produce an IMF that becomes increasingly top-heavy with time. In all cases we find that the competition between magnetic flux advection toward the protostar and outward advection due to magnetic interchange instabilities, and the competition between turbulent amplification and reconnection close to newly-formed protostars renders the local magnetic field structure insensitive to the strength of the large-scale field, ensuring that radiation is always more important than magnetic support in setting the fragmentation scale and thus the IMF peak mass. The statistics of multiple stellar systems are similarly insensitive to variations in the initial conditions and generally agree with observations within the range of statistical uncertainty.

The Legacy ExtraGalactic UV Survey (LEGUS) is a multiwavelength Cycle 21 Treasury program on the Hubble Space Telescope. It studied 50 nearby star-forming galaxies in five bands from the near UV to the I-band, combining new Wide Field Camera 3 observations with archival Advanced Camera for Surveys data. LEGUS was designed to investigate how star formation occurs and develops on both small and large scales, and how it relates to the galactic environments. In this paper we present the photometric catalogs for all the apparently single stars identified in the 50 LEGUS galaxies. Photometric catalogs and mosaicked images for all filters are available for download. We present optical and near UV color-magnitude diagrams for all the galaxies. For each galaxy we derived the distance from the tip of the red giant branch. We then used the NUV color-magnitude diagrams to identify stars more massive than 14 Mo, and compared their number with the number of massive stars expected from the GALEX FUV luminosity. Our analysis shows that the fraction of massive stars forming in star clusters and stellar associations is about constant with the star formation rate. This lack of a relation suggests that the time scale for evaporation of unbound structures is comparable or longer than 10 Myr. At low star formation rates this translates to an excess of mass in clustered environments as compared to model predictions of cluster evolution, suggesting that a significant fraction of stars form in unbound systems.

Carbon monoxide (CO) provides crucial information about the molecular gas properties of galaxies. While $^{12}$CO has been targeted extensively, isotopologues such as $^{13}$CO have the advantage of being less optically thick and observations have recently become accessible across full galaxy discs. We present a comprehensive new dataset of $^{13}$CO(1-0) observations with the IRAM 30-m telescope of the full discs of 9 nearby spiral galaxies from the EMPIRE survey at a spatial resolution of $\sim$1.5kpc. $^{13}$CO(1-0) is mapped out to $0.7-1r_{25}$ and detected at high signal-to-noise throughout our maps. We analyse the $^{12}$CO(1-0)-to-$^{13}$CO(1-0) ratio ($\Re$) as a function of galactocentric radius and other parameters such as the $^{12}$CO(2-1)-to-$^{12}$CO(1-0) intensity ratio, the 70-to-160$\mu$m flux density ratio, the star-formation rate surface density, the star-formation efficiency, and the CO-to-H$_2$ conversion factor. We find that $\Re$ varies by a factor of 2 at most within and amongst galaxies, with a median value of 11 and larger variations in the galaxy centres than in the discs. We argue that optical depth effects, most likely due to changes in the mixture of diffuse/dense gas, are favored explanations for the observed $\Re$ variations, while abundance changes may also be at play. We calculate a spatially-resolved $^{13}$CO(1-0)-to-H$_2$ conversion factor and find an average value of $1.0\times10^{21}$ cm$^{-2}$ (K.km/s)$^{-1}$ over our sample with a standard deviation of a factor of 2. We find that $^{13}$CO(1-0) does not appear to be a good predictor of the bulk molecular gas mass in normal galaxy discs due to the presence of a large diffuse phase, but it may be a better tracer of the mass than $^{12}$CO(1-0) in the galaxy centres where the fraction of dense gas is larger.

The radiation force on dust grains may be dynamically important in driving turbulence and outflows in rapidly star-forming galaxies. Recent studies focus on the highly optically-thick limit relevant to the densest ultra-luminous galaxies and super star clusters, where reprocessed infrared photons provide the dominant source of electromagnetic momentum. However, even among starburst galaxies, the great majority instead lie in the so-called "single-scattering" limit, where the system is optically-thick to the incident starlight, but optically-thin to the re-radiated infrared. In this paper we present a stability analysis and multidimensional radiation-hydrodynamic simulations exploring the stability and dynamics of isothermal dusty gas columns in this regime. We describe our algorithm for full angle-dependent radiation transport based on the discontinuous Galerkin finite element method. For a range of near-Eddington fluxes, we show that the medium is unstable, producing convective-like motions in a turbulent atmosphere with a scale height significantly inflated compared to the gas pressure scale height and mass-weighted turbulent energy densities of $\sim 0.01-0.1$ of the midplane radiation energy density, corresponding to mass-weighted velocity dispersions of Mach number $\sim 0.5-2$. Extrapolation of our results to optical depths of $10^3$ implies maximum turbulent Mach numbers of $\sim20$. Comparing our results to galaxy-averaged observations, and subject to the approximations of our calculations, we find that radiation pressure does not contribute significantly to the effective supersonic pressure support in star-forming disks, which in general are substantially sub-Eddington. We further examine the time-averaged vertical density profiles in dynamical equilibrium and comment on implications for radiation-pressure-driven galactic winds.

Winds arising from galaxies, star clusters, and active galactic nuclei are crucial players in star and galaxy formation, but it has proven remarkably difficult to use observations of them to determine physical properties of interest, particularly mass fluxes. Much of the difficulty stems from a lack of a theory that links a physically-realistic model for winds' density, velocity, and covering factors to calculations of light emission and absorption. In this paper we provide such a model. We consider a wind launched from a turbulent region with a range of column densities, derive the differential acceleration of gas as a function of column density, and use this result to compute winds' absorption profiles, emission profiles, and emission intensity maps in both optically thin and optically thick species. The model is sufficiently simple that all required computations can be done analytically up to straightforward numerical integrals, rendering it suitable for the problem of deriving physical parameters by fitting models to observed data. We show that our model produces realistic absorption and emission profiles for some example cases, and argue that the most promising methods of deducing mass fluxes are based on combinations of absorption lines of different optical depths, or on combining absorption with measurements of molecular line emission. In the second paper in this series, we expand on these ideas by introducing a set of observational diagnostics that are significantly more robust that those commonly in use, and that can be used to obtain improved estimates of wind properties.

We present the first images of the 28SiO v=1, J=2-1 maser emission around the closest known massive young stellar object Orion Source I observed at 86 GHz (3mm) with the VLBA. These images have high spatial (~0.3 mas) and spectral (~0.054 km/s) resolutions. We find that the 3mm masers lie in an X-shaped locus consisting of four arms, with blue-shifted emission in the south and east arms and red-shifted emission in the north and west arms. Comparisons with previous images of the 28SiO v=1,2, J=1-0 transitions at 7mm (observed in 2001-2002) show that the bulk of the J=2-1 transition emission follows the streamlines of the J=1-0 emission and exhibits an overall velocity gradient consistent with the gradient at 7mm. While there is spatial overlap between the 3mm and 7mm transitions, the 3mm emission, on average, lies at larger projected distances from Source I (~44 AU compared with ~35 AU for 7mm). The spatial overlap between the v=1, J=1-0 and J=2-1 transitions is suggestive of a range of temperatures and densities where physical conditions are favorable for both transitions of a same vibrational state. However, the observed spatial offset between the bulk of emission at 3mm and 7mm possibly indicates different ranges of temperatures and densities for optimal excitation of the masers. We discuss different maser pumping models that may explain the observed offset. We interpret the 3mm and 7mm masers as being part of a single wide-angle outflow arising from the surface of an edge-on disk rotating about a northeast-southwest axis, with a continuous velocity gradient indicative of differential rotation consistent with a Keplerian profile in a high-mass proto-binary.

We present an analysis of the positions and ages of young star clusters in eight local galaxies to investigate the connection between the age difference and separation of cluster pairs. We find that star clusters do not form uniformly but instead are distributed such that the age difference increases with the cluster pair separation to the 0.25-0.6 power, and that the maximum size over which star formation is physically correlated ranges from ~200 pc to ~1 kpc. The observed trends between age difference and separation suggest that cluster formation is hierarchical both in space and time: clusters that are close to each other are more similar in age than clusters born further apart. The temporal correlations between stellar aggregates have slopes that are consistent with turbulence acting as the primary driver of star formation. The velocity associated with the maximum size is proportional to the galaxy's shear, suggesting that the galactic environment influences the maximum size of the star-forming structures.

We present the implementation of a Bayesian formalism within the Stochastically Lighting Up Galaxies (SLUG) stellar population synthesis code, which is designed to investigate variations in the initial mass function (IMF) of star clusters. By comparing observed cluster photometry to large libraries of clusters simulated with a continuously varying IMF, our formalism yields the posterior probability distribution function (PDF) of the cluster mass, age, and extinction, jointly with the parameters describing the IMF. We apply this formalism to a sample of star clusters from the nearby galaxy NGC 628, for which broad-band photometry in five filters is available as part of the Legacy ExtraGalactic UV Survey (LEGUS). After allowing the upper-end slope of the IMF ($\alpha_3$) to vary, we recover PDFs for the mass, age, and extinction that are broadly consistent with what is found when assuming an invariant Kroupa IMF. However, the posterior PDF for $\alpha_3$ is very broad due to a strong degeneracy with the cluster mass, and it is found to be sensitive to the choice of priors, particularly on the cluster mass. We find only a modest improvement in the constraining power of $\alpha_3$ when adding H$\alpha$ photometry from the companion H$\alpha$-LEGUS survey. Conversely, H$\alpha$ photometry significantly improves the age determination, reducing the frequency of multi-modal PDFs. With the aid of mock clusters we quantify the degeneracy between physical parameters, showing how constraints on the cluster mass that are independent of photometry can be used to pin down the IMF properties of star clusters.

We report the large effort which is producing comprehensive high-level young star cluster (YSC) catalogues for a significant fraction of galaxies observed with the Legacy ExtraGalactic UV Survey (LEGUS) Hubble treasury program. We present the methodology developed to extract cluster positions, verify their genuine nature, produce multiband photometry (from NUV to NIR), and derive their physical properties via spectral energy distribution fitting analyses. We use the nearby spiral galaxy NGC628 as a test case for demonstrating the impact that LEGUS will have on our understanding of the formation and evolution of YSCs and compact stellar associations within their host galaxy. Our analysis of the cluster luminosity function from the UV to the NIR finds a steepening at the bright end and at all wavelengths suggesting a dearth of luminous clusters. The cluster mass function of NGC628 is consistent with a power-law distribution of slopes $\sim -2$ and a truncation of a few times $10^5$ M$_\odot$. After their formation YSCs and compact associations follow different evolutionary paths. YSCs survive for a longer timeframe, confirming their being potentially bound systems. Associations disappear on time scales comparable to hierarchically organized star-forming regions, suggesting that they are expanding systems. We find mass-independent cluster disruption in the inner region of NGC628, while in the outer part of the galaxy there is little or no disruption. We observe faster disruption rates for low mass ($\leq$ $10^4$ M$_\odot$) clusters suggesting that a mass-dependent component is necessary to fully describe the YSC disruption process in NGC628.

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.

Observations of ionised carbon at 158 micron ([CII]) from luminous star-forming galaxies at z~0 show that their ratios of [CII] 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 "[CII] 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. [CII] emission arises from the outer, HI 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 [CII] emission. This provides a natural mechanism to explain the observed [CII]-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 [CII] emission. Our model further suggests that the apparent offset in the [CII]-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 relatively modest surface density, allowing copious [CII] emission to coexist with rapid star formation.

We use the IRAM Large Program EMPIRE and new high-resolution ALMA data to measure 13CO(1-0)/C18O(1-0) intensity ratios across nine nearby spiral galaxies. These isotopologues of CO are typically optically thin across most of the area in galaxy disks, and this ratio allows us to gauge their relative abundance due to chemistry or stellar nucleosynthesis effects. Resolved 13CO/C18O gradients across normal galaxies have been rare due to the faintness of these lines. We find a mean 13CO/C18O ratio of 6.0$\pm$0.9 for the central regions of our galaxies. This agrees well with results in the Milky Way, but differs from results for starburst galaxies (3.4$\pm$0.9) and ultraluminous infrared galaxies (1.1$\pm$0.4). In our sample, the 13CO/C18O ratio consistently increases with increasing galactocentric radius and decreases with increasing star formation rate surface density. These trends qualitatively agree with expectations for carbon and oxygen isotopic abundance variations due to stellar nucleosynthesis, with a possible effect of fractionation.

High critical density molecular lines like HCN(1-0) or HCO+(1-0) represent our best tool to study currently star-forming, dense molecular gas at extragalactic distances. The optical depth of these lines is a key ingredient to estimate the effective density required to excite emission. However, constraints on this quantity are even scarcer in the literature than measurements of the high density tracers themselves. Here, we combine new observations of HCN, HCO+ and HNC(1-0) and their optically thin isotopologues H13CN, H13CO+ and HN13C(1-0) to measure isotopologue line ratios. We use IRAM 30-m observations from the large program EMPIRE and new ALMA observations, which together target 6 nearby star-forming galaxies. Using spectral stacking techniques, we calculate or place strong upper limits on the HCN/H13CN, HCO+/H13CO+ and HNC/HN13C line ratios in the inner parts of these galaxies. Under simple assumptions, we use these to estimate the optical depths of HCN(1-0) and HCO+(1-0) to be \tau ~2-11 in the active, inner regions of our targets. The critical densities are consequently lowered to values between 5-20$\times 10^5$, 1-3$\times 10^5$ and 9$\times 10^4$ cm-3 for HCN, HCO+ and HNC, respectively. We study the impact of having different beam-filling factors, $\eta$, on these estimates and find that the effective critical densities decrease by a factor of $\frac{\eta_{12}}{\eta_{13}}\,\tau_{12}$. A comparison to existing work in NGC 5194 and NGC 253 shows HCN/H13CN and HCO+/H13CO+ ratios in agreement with our measurements within the uncertainties. The same is true for studies in other environments such as the Galactic Centre or nuclear regions of AGN-dominated nearby galaxies.

This book provides an introduction to the field of star formation at a level suitable for graduate students or advanced undergraduates in astronomy or physics. The structure of the book is as follows. The first two chapters begin with a discussion of observational techniques, and the basic phenomenology they reveal. The goal is to familiarize students with the basic techniques that will be used throughout, and to provide a common vocabulary for the rest of the book. The next five chapters provide a similar review of the basic physical processes that are important for star formation. Again, the goal is to provide a basis for what follows. The remaining chapters discuss star formation over a variety of scales, starting with the galactic scale and working down to the scales of individual stars and their disks. The book concludes with a brief discussion of the clearing of disks and the transition to planet formation. The book includes five problem sets, complete with solutions.

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 super-linearly with the number of SNe, resulting in a momentum per SN that 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 disk, 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 multi-dimensional instabilities or preexisting density variations, as a limiting factor in the buildup of momentum by clustered SNe, and suggest future numerical experiments to explore these effects.

We investigate the burstiness of star formation histories (SFHs) of galaxies at $0.4<z<1$ by using the ratio of star formation rates (SFRs) measured from H$\beta$ and FUV (1500 \AA) (H$\beta$--to--FUV ratio). Our sample contains 164 galaxies down to stellar mass (M*) of $10^{8.5} M_\odot$ in the CANDELS GOODS-N region, where Team Keck Redshift Survey DEIMOS spectroscopy and HST/WFC3 F275W images from CANDELS and Hubble Deep UV Legacy Survey are available. When the {\it ratio} of H$\beta$- and FUV-derived SFRs is measured, dust extinction correction is negligible (except for very dusty galaxies) with the Calzetti attenuation curve. The H$\beta$--to--FUV ratio of our sample increases with M* and SFR. The median ratio is $\sim$0.7 at M*$\sim10^{8.5} M_\odot$ (or SFR$\sim 0.5 M_\odot/yr$) and increases to $\sim$1 at M*$\sim10^{10} M_\odot$ (or SFR $\sim 10 M_\odot/yr$). At M*$<10^{9.5} M_\odot$, our median H$\beta$--to--FUV ratio is lower than that of local galaxies at the same M*, implying a redshift evolution. Bursty SFH on a timescale of a few tens of megayears on galactic scales provides a plausible explanation of our results, and the importance of the burstiness increases as M* decreases. Due to sample selection effects, our H$\beta$--to--FUV ratio may be an upper limit of the true value of a complete sample, which strengthens our conclusions. Other models, e.g., non-universal initial mass function or stochastic star formation on star cluster scales, are unable to plausibly explain our results.

Massive galaxies at high redshift are predicted to be fed from the cosmic web by narrow, dense, cold streams. These streams penetrate supersonically through the hot medium encompassed by a stable shock near the virial radius of the dark-matter halo. Our long-term goal is to explore the heating and dissipation rate of the streams and their fragmentation and possible breakup, in order to understand how galaxies are fed, and how this affects their star-formation rate and morphology. We present here the first step, where we analyze the linear Kelvin-Helmholtz instability (KHI) of a cold, dense slab or cylinder flowing through a hot, dilute medium in the transonic regime. The current analysis is limited to the adiabatic case with no gravity and assuming equal pressure in the stream and the medium. By analytically solving the linear dispersion relation, we find a transition from a dominance of the familiar rapidly growing surface modes in the subsonic regime to more slowly growing body modes in the supersonic regime. The system is parameterized by three parameters: the density contrast between the stream and the medium, the Mach number of stream velocity with respect to the medium, and the stream width with respect to the halo virial radius. We find that a realistic choice for these parameters places the streams near the mode transition, with the KHI exponential-growth time in the range 0.01-10 virial crossing times for a perturbation wavelength comparable to the stream width. We confirm our analytic predictions with idealized hydrodynamical simulations. Our linear-KHI estimates thus indicate that KHI may in principle be effective in the evolution of streams by the time they reach the galaxy. More definite conclusions await the extension of the analysis to the nonlinear regime and the inclusion of cooling, thermal conduction, the halo potential well, self-gravity and magnetic fields.

The pressure exerted by massive stars' radiation fields is an important mechanism regulating their formation. Detailed simulation of massive star formation therefore requires an accurate treatment of radiation. However, all published simulations have either used a diffusion approximation of limited validity; have only been able to simulate a single star fixed in space, thereby suppressing potentially-important instabilities; or did not provide adequate resolution at locations where instabilities may develop. To remedy this we have developed a new, highly accurate radiation algorithm that properly treats the absorption of the direct radiation field from stars and the re-emission and processing by interstellar dust. We use our new tool to perform three-dimensional radiation-hydrodynamic simulations of the collapse of massive pre-stellar cores with laminar and turbulent initial conditions and properly resolve regions where we expect instabilities to grow. We find that mass is channeled to the stellar system via gravitational and Rayleigh-Taylor (RT) instabilities, in agreement with previous results using stars capable of moving, but in disagreement with methods where the star is held fixed or with simulations that do not adequately resolve the development of RT instabilities. For laminar initial conditions, proper treatment of the direct radiation field produces later onset of instability, but does not suppress it entirely provided the edges of radiation-dominated bubbles are adequately resolved. Instabilities arise immediately for turbulent pre-stellar cores because the initial turbulence seeds the instabilities. Our results suggest that RT features are significant and should be present around accreting massive stars throughout their formation.

We present a highly-parallel multi-frequency hybrid radiation hydrodynamics algorithm that combines a spatially-adaptive long characteristics method for the radiation field from point sources with a moment method that handles the diffuse radiation field produced by a volume-filling fluid. Our Hybrid Adaptive Ray-Moment Method (HARM$^2$) operates on patch-based adaptive grids, is compatible with asynchronous time stepping, and works with any moment method. In comparison to previous long characteristics methods, we have greatly improved the parallel performance of the adaptive long-characteristics method by developing a new completely asynchronous and non-blocking communication algorithm. As a result of this improvement, our implementation achieves near-perfect scaling up to $\mathcal{O}(10^3)$ processors on distributed memory machines. We present a series of tests to demonstrate the accuracy and performance of the method.

To understand the conditions under which dense, molecular gas is able to form within a galaxy, we post-process a series of three-dimensional galactic-disk-scale simulations with ray-tracing based radiative transfer and chemical network integration to compute the equilibrium chemical and thermal state of the gas. In performing these simulations we vary a number of parameters, such as the ISRF strength, vertical scale height of stellar sources, cosmic ray flux, to gauge the sensitivity of our results to these variations. Self-shielding permits significant molecular hydrogen (H2) abundances in dense filaments around the disk midplane, accounting for approximately ~10-15% of the total gas mass. Significant CO fractions only form in the densest, n>~10^3 cm^-3, gas where a combination of dust, H2, and self-shielding attenuate the FUV background. We additionally compare these ray-tracing based solutions to photochemistry with complementary models where photo-shielding is accounted for with locally computed prescriptions. With some exceptions, these local models for the radiative shielding length perform reasonably well at reproducing the distribution and amount of molecular gas as compared with a detailed, global ray tracing calculation. Specifically, an approach based on the Jeans Length with a T=40K temperature cap performs the best in regards to a number of different quantitative measures based on the H2 and CO abundances.

As star-forming clouds collapse, the gas within them fragments to ever-smaller masses. Naively one might expect this process to continue down to the smallest mass that is able to radiate away its binding energy on a dynamical timescale, the opacity limit for fragmentation, at $\sim 0.01$ $M_\odot$. However, the observed peak of the initial mass function (IMF) lies a factor of $20-30$ higher in mass, suggesting that some other mechanism halts fragmentation before the opacity limit is reached. In this paper we analyse radiation-magnetohydrodynamic simulations of star cluster formation in typical Milky Way environments in order to determine what physical process limits fragmentation in them. We examine the regions in the vicinity of stars that form in the simulations to determine the amounts of mass that are prevented from fragmenting by thermal and magnetic pressure. We show that, on small scales, thermal pressure enhanced by stellar radiation heating is the dominant mechanism limiting the ability of the gas to further fragment. In the brown dwarf mass regime, $\sim 0.01$ $M_\odot$, the typical object that forms in the simulations is surrounded by gas whose mass is several times its own that is unable to escape or fragment, and instead is likely to accrete. This mechanism explains why $\sim 0.01$ $M_\odot$ objects are rare: unless an outside agent intervenes (e.g., a shock strips away the gas around them), they will grow by accreting the warmed gas around them. In contrast, by the time stars grow to masses of $\sim 0.2$ $M_\odot$, the mass of heated gas is only tens of percent of the central star mass, too small to alter its final mass by a large factor. This naturally explains why the IMF peak is at $\sim 0.2$ $M_\odot$.

Self-gravity and stellar feedback are capable of driving turbulence and transporting mass and angular momentum in disk galaxies, but the balance between them is not well understood. In the previous paper in this series, we showed that gravity alone can drive turbulence in galactic disks, regulate their Toomre $Q$ parameters to $\sim$ 1, and transport mass inwards at a rate sufficient to fuel star formation in the centers of present-day galaxies. In this paper we extend our models to include the effects of star formation feedback. We show that feedback suppresses galaxies' star formation rates by a factor of $\sim$ 5 and leads to the formation of a multi-phase atomic and molecular ISM. Both the star formation rate and the phase balance produced in our simulations agree well with observations of nearby spirals. After our galaxies reach steady state, we find that the inclusion of feedback actually lowers the gas velocity dispersion slightly compared to the case of pure self-gravity, and also slightly reduces the rate of inward mass transport. Nevertheless, we find that, even with feedback included, our galactic disks self-regulate to $Q$ $\sim$ 1, and transport mass inwards at a rate sufficient to supply a substantial fraction of the inner disk star formation. We argue that gravitational instability is therefore likely to be the dominant source of turbulence and transport in galactic disks, and that it is responsible for fueling star formation in the inner parts of galactic disks over cosmological times.

Photoelectric heating has long been recognized as the primary source of heating for the neutral interstellar medium. Simulations of spiral galaxies found some indication that photoelectric heating could suppress star formation. However, simulations that include photoelectric heating have typically found that it has little effect on the rate of star formation in either spiral galaxies or dwarfs suggesting that supernovae and not photoelectric heating are responsible for setting the star formation law in galaxies. This result is in tension with recent work indicating that a star formation law that depends on galaxy metallicity, as expected for photoelectric heating but not for supernovae, reproduces the present-day galaxy population better than a metallicity-independent one. Here we report a series of simulations of dwarf galaxies, where the effects of both photoelectric heating and supernovae are expected to be strongest. We simultaneously include space- and time-dependent photoelectric heating, and we resolve the Sedov phase of every supernova blast wave, allowing us to make a direct measurement of the relative importance of momentum injection by supernovae and dust heating by far ultraviolet (FUV) photons in suppressing star formation. We find that supernovae are unable to account for the long observed gas depletion times in dwarf galaxies. Instead, ordinary photoelectric heating is the dominant means by which dwarf galaxies regulate their star formation rate at any given time, suppressing the star formation rate by more than an order of magnitude relative to simulations with only supernovae.

Galaxies' interstellar media (ISM) are observed to be supersonically-turbulent, but the ultimate power source that drives turbulent motion remains uncertain. The two dominant models are that the turbulence is driven by star formation feedback and/or that it is produced by gravitational instability in the gas. Here we show that, while both models predict that the galaxies' ISM velocity dispersions will be positively correlated with their star formation rates, the forms of the correlation predicted by these two models are subtly but measurably different. A feedback-driven origin for the turbulence predicts a velocity dispersion that rises more sharply with star formation rate, and that does not depend on the gas fraction (i.e. $\dot{M}_* \propto \sigma^2$), while a gravity-driven model yields a shallower rise and a strong dependence on gas fraction (i.e. $\dot{M}_* \propto f_g^2 \sigma$). We compare the models to a collection of data on local and high-redshift galaxies culled from the literature, and show that the correlation expected for gravity-driven turbulence is a better match to the observations than a feedback-driven model. This suggests that gravity is the ultimate source of ISM turbulence, at least in the rapidly-star-forming, high velocity dispersion galaxies for which our test is most effective. We conclude by discussing the limitations of the present data set, and the prospects for future measurements to enable a more definitive test of the two models.

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.

We show that the turbulent gas in the star-forming regions of galaxies is unstable to wind formation via momentum deposition by radiation pressure or other momentum sources like supernova explosions, even if the system is below the average Eddington limit. This conclusion follows from the fact that the critical momentum injection rate per unit mass for unbinding gas from a self-gravitating system is proportional to the gas surface density and that a turbulent medium presents a broad distribution of column densities to the sources. For an average Eddington ratio of <Gamma>~0.1 and for turbulent Mach numbers >30, we find that ~1% of the gas is ejected per dynamical timescale at velocities larger than the local escape velocity. Because of the lognormal shape of the surface density distribution, the mass loss rate is highly sensitive to the average Eddington ratio, reaching 20-40% of the gas mass per dynamical time for <Gamma>=1. Using this model we find a large scatter in the mass-loading factor for star-forming galaxies, ranging from 0.001-10, but with significant uncertainties. Implications for the efficiency of star formation in giant molecular clouds are highlighted. For radiation pressure feedback alone, we find an increasing star formation efficiency as a function of initial gas surface density. Uncertainties are discussed.

The role of gravitational instability-driven turbulence in determining the structure and evolution of disk galaxies, and the extent to which gravity rather than feedback can explain galaxy properties, remains an open question. To address it, we present high resolution adaptive mesh refinement simulations of Milky Way-like isolated disk galaxies, including realistic heating and cooling rates and a physically motivated prescription for star formation, but no form of star formation feedback. After an initial transient, our galaxies reach a state of fully-nonlinear gravitational instability. In this state, gravity drives turbulence and radial inflow. Despite the lack of feedback, the gas in our galaxy models shows substantial turbulent velocity dispersions, indicating that gravitational instability alone may be able to power the velocity dispersions observed in nearby disk galaxies on 100 pc scales. Moreover, the rate of mass transport produced by this turbulence approaches $\sim 1$ $M_\odot$ yr$^{-1}$ for Milky Way-like conditions, sufficient to fully fuel star formation in the inner disks of galaxies. In a companion paper we add feedback to our models, and use the comparison between the two cases to understand what galaxy properties depend sensitively on feedback, and which can be understood as the product of gravity alone. All of the code, initial conditions, and simulation data for our model are publicly available.

We investigate a novel Bayesian analysis method, based on the Stochastically Lighting Up Galaxies (slug) code, to derive the masses, ages, and extinctions of star clusters from integrated light photometry. Unlike many analysis methods, slug correctly accounts for incomplete IMF sampling, and returns full posterior probability distributions rather than simply probability maxima. We apply our technique to 621 visually-confirmed clusters in two nearby galaxies, NGC 628 and NGC 7793, that are part of the Legacy Extragalactic UV Survey (LEGUS). LEGUS provides Hubble Space Telescope photometry in the NUV, U, B, V, and I bands. We analyze the sensitivity of the derived cluster properties to choices of prior probability distribution, evolutionary tracks, IMF, metallicity, treatment of nebular emission, and extinction curve. We find that slug's results for individual clusters are insensitive to most of these choices, but that the posterior probability distributions we derive are often quite broad, and sometimes multi-peaked and quite sensitive to the choice of priors. In contrast, the properties of the cluster population as a whole are relatively robust against all of these choices. We also compare our results from slug to those derived with a conventional non-stochastic fitting code, Yggdrasil. We show that slug's stochastic models are generally a better fit to the observations than the deterministic ones used by Yggdrasil. However, the overall properties of the cluster populations recovered by both codes are qualitatively similar.

The nearby dwarf starburst galaxy NGC5253 hosts a number of young, massive star clusters, the two youngest of which are centrally concentrated and surrounded by thermal radio emission (the `radio nebula'). To investigate the role of these clusters in the starburst energetics, we combine new and archival Hubble Space Telescope images of NGC5253 with wavelength coverage from 1500 Ang to 1.9 micron in 13 filters. These include H-alpha, P-beta, and P-alpha, and the imaging from the Hubble Treasury Program LEGUS (Legacy Extragalactic UV Survey). The extraordinarily well-sampled spectral energy distributions enable modeling with unprecedented accuracy the ages, masses, and extinctions of the 9 optically brightest clusters (M_V < -8.8) and the two young radio nebula clusters. The clusters have ages ~1-15 Myr and masses ~1x10^4 - 2.5x10^5 M_sun. The clusters' spatial location and ages indicate that star formation has become more concentrated towards the radio nebula over the last ~15 Myr. The most massive cluster is in the radio nebula; with a mass 2.5x10^5 M_sun and an age ~1 Myr, it is 2-4 times less massive and younger than previously estimated. It is within a dust cloud with A_V~50 mag, and shows a clear nearIR excess, likely from hot dust. The second radio nebula cluster is also ~1 Myr old, confirming the extreme youth of the starburst region. These two clusters account for about half of the ionizing photon rate in the radio nebula, and will eventually supply about 2/3 of the mechanical energy in present-day shocks. Additional sources are required to supply the remaining ionizing radiation, and may include very massive stars.

Long gamma-ray bursts (GRBs), among the most energetic events in the Universe, are explosions of massive and short-lived stars, so they pinpoint locations of recent star formation. However, several GRB host galaxies have recently been found to be deficient in molecular gas (H2), believed to be the fuel of star formation. Moreover, optical spectroscopy of GRB afterglows implies that the molecular phase constitutes only a small fraction of the gas along the GRB line-of-sight. Here we report the first ever 21 cm line observations of GRB host galaxies, using the Australia Telescope Compact Array, implying high levels of atomic hydrogen (HI), which suggests that the connection between atomic gas and star formation is stronger than previously thought, with star formation being potentially directly fuelled by atomic gas (or with very efficient HI-to-H2 conversion and rapid exhaustion of molecular gas), as has been theoretically shown to be possible. This can happen in low metallicity gas near the onset of star formation, because cooling of gas (necessary for star formation) is faster than the HI-to-H2 conversion. Indeed, large atomic gas reservoirs, together with low molecular gas masses, stellar and dust masses are consistent with GRB hosts being preferentially galaxies which have very recently started a star formation episode after accreting metal-poor gas from the intergalactic medium. This provides a natural route for forming GRBs in low-metallicity environments. The gas inflow scenario is also consistent with the existence of the companion HI object with no optical counterpart ~19 kpc from the GRB 060505 host, and with the fact that the HI centroids of the GRB 980425 and 060505 hosts do not coincide with optical centres of these galaxies, but are located close to the GRB positions.

Ionizing stellar photons heat the upper regions of planetary atmospheres, driving atmospheric mass loss. Gas escaping from several hot, hydrogen-rich planets has been detected using UV and X-ray transmission spectroscopy. Because these planets are tidally locked, and thus asymmetrically irradiated, escaping gas is unlikely to be spherically symmetric. In this paper, we focus on the effects of asymmetric heating on local outflow structure. We use the Athena code for hydrodynamics to produce 3D simulations of hot Jupiter mass loss that jointly model wind launching and stellar heating via photoionization. Our fiducial planet is an inflated, hot Jupiter with radius $R_p=2.14 R_{\rm Jup}$ and mass $M_p = 0.53 M_{\rm Jup}$. We irradiate the initially neutral, atomic hydrogen atmosphere with 13.6 eV photons and compute the outflow's ionization structure. There are clear asymmetries in the atmospheric outflow, including a neutral shadow on the planet's nightside. Given an incident ionizing UV flux comparable to that of the Sun, we find a steady-state mass loss rate of ~$2\times10^{10}$ g s$^{-1}$. The total mass loss rate and the outflow substructure along the substellar ray show good agreement with earlier 1D models, for two different fluxes. Our 3D data cube can be used to generate the outflow's extinction spectrum during transit. As a proof of concept, we find absorption of stellar Lyman-alpha at Doppler-shifted velocities of up to $\pm 50$ km s$^{-1}$. Our work provides a starting point for further 3D models that can be used to predict observable signatures of hot Jupiter mass loss.

We develop a simple dynamical model for the evolution of gas in the centres of barred spiral galaxies, using the Milky Way's Central Molecular Zone (CMZ, i.e., the central few hundred pc) as a case study. We show that, in the presence of a galactic bar, gas in a disc in the central regions of a galaxy will be driven inwards by angular momentum transport induced by acoustic instabilities within the bar's inner Lindblad resonance. This transport process drives turbulence within the gas that temporarily keeps it strongly gravitationally stable and prevents the onset of rapid star formation. However, at some point the rotation curve must transition from approximately flat to approximately solid body, and the resulting reduction in shear reduces the transport rates and causes gas to build up, eventually producing a gravitationally-unstable region that is subject to rapid and violent star formation. For the observed rotation curve of the Milky Way, the accumulation happens $\sim 100$ pc from the centre of the Galaxy, in good agreement with the observed location of gas clouds and young star clusters in the CMZ. The characteristic timescale for gas accumulation and star formation is of order $10-20$ Myr. We argue that similar phenomena should be ubiquitous in other barred spiral galaxies.

Stellar population synthesis techniques for predicting the observable light emitted by a stellar population have extensive applications in numerous areas of astronomy. However, accurate predictions for small populations of young stars, such as those found in individual star clusters, star-forming dwarf galaxies, and small segments of spiral galaxies, require that the population be treated stochastically. Conversely, accurate deductions of the properties of such objects also requires consideration of stochasticity. Here we describe a comprehensive suite of modular, open-source software tools for tackling these related problems. These include: a greatly-enhanced version of the slug code introduced by da Silva et al. (2012), which computes spectra and photometry for stochastically- or deterministically-sampled stellar populations with nearly-arbitrary star formation histories, clustering properties, and initial mass functions; cloudy_slug, a tool that automatically couples slug-computed spectra with the cloudy radiative transfer code in order to predict stochastic nebular emission; bayesphot, a general-purpose tool for performing Bayesian inference on the physical properties of stellar systems based on unresolved photometry; and cluster_slug and sfr_slug, a pair of tools that use bayesphot on a library of slug models to compute the mass, age, and extinction of mono-age star clusters, and the star formation rate of galaxies, respectively. The latter two tools make use of an extensive library of pre-computed stellar population models, which are included the software. The complete package is available at http://www.slugsps.com.

In this paper we present the MOSFIRE Deep Evolution Field (MOSDEF) survey. The MOSDEF survey aims to obtain moderate-resolution (R=3000-3650) rest-frame optical spectra (~3700-7000 Angstrom) for ~1500 galaxies at 1.37<z<3.80 in three well-studied CANDELS fields: AEGIS, COSMOS, and GOODS-N. Targets are selected in three redshift intervals: 1.37<z<1.70, 2.09<z<2.61, and 2.95<z<3.80, down to fixed H_AB (F160W) magnitudes of 24.0, 24.5 and 25.0, respectively, using the photometric and spectroscopic catalogs from the 3D-HST survey. We target both strong nebular emission lines (e.g., [OII], Hbeta, [OIII], 5008, Halpha, [NII], and [SII]) and stellar continuum and absorption features (e.g., Balmer lines, Ca-II H and K, Mgb, 4000 Angstrom break). Here we present an overview of our survey, the observational strategy, the data reduction and analysis, and the sample characteristics based on spectra obtained during the first 24 nights. To date, we have completed 21 masks, obtaining spectra for 591 galaxies. For ~80% of the targets we derive a robust redshift from either emission or absorption lines. In addition, we confirm 55 additional galaxies, which were serendipitously detected. The MOSDEF galaxy sample includes unobscured star-forming, dusty star-forming, and quiescent galaxies and spans a wide range in stellar mass (~10^9-10^11.5 Msol) and star formation rate (~10^0-10^3 Msol/yr). The spectroscopically confirmed sample is roughly representative of an H-band limited galaxy sample at these redshifts. With its large sample size, broad diversity in galaxy properties, and wealth of available ancillary data, MOSDEF will transform our understanding of the stellar, gaseous, metal, dust, and black hole content of galaxies during the time when the universe was most active.

Metal production in galaxies traces star formation, and is highly concentrated toward the centers of galactic discs. This suggests that galaxies should have inhomogeneous metal distributions with strong radial gradients, but observations of present-day galaxies show only shallow gradients with little azimuthal variation, implying the existence of a redistribution mechanism. We study the role of gravitational instability-driven turbulence as a mixing mechanism by simulating an isolated galactic disc at high resolution, including metal fields treated as passive scalars. Since any cylindrical field can be decomposed into a sum of Fourier-Bessel basis functions, we set up initial metal fields characterized by these functions and study how different modes mix. We find both shear and turbulence contribute to mixing, but the mixing strongly depends on the symmetries of the mode. Non-axisymmetric modes have decay times smaller than the galactic orbital period because shear winds them up to small spatial scales, where they are erased by turbulence. The decay timescales for axisymmetric modes are much greater, though for all but the largest-scale inhomogeneities the mixing timescale is still short enough to erase chemical inhomogeneities over cosmological times. These different timescales provide an explanation for why galaxies retain metallicity gradients while there is almost no variation at a fixed radius. Moreover, the comparatively long timescales required for mixing axisymmetric modes may explain the greater diversity of metallicity gradients observed in high redshift galaxies as compared to local ones: these systems have not yet reached equilibrium between metal production and diffusion.

The evolution of thin axisymmetric viscous accretion disks is a classic problem in astrophysics. While models based on this simplified geometry provide only approximations to the true processes of instability-driven mass and angular momentum transport, their simplicity makes them invaluable tools for both semi-analytic modeling and simulations of long-term evolution where two- or three-dimensional calculations are too computationally costly. Despite the utility of these models, the only publicly-available frameworks for simulating them are rather specialized and non-general. Here we describe a highly flexible, general numerical method for simulating viscous thin disks with arbitrary rotation curves, viscosities, boundary conditions, grid spacings, equations of state, and rates of gain or loss of mass (e.g., through winds) and energy (e.g., through radiation). Our method is based on a conservative, finite-volume, second-order accurate discretization of the equations, which we solve using an unconditionally-stable implicit scheme. We implement Anderson acceleration to speed convergence of the scheme, and show that this leads to factor of $\sim 5$ speed gains over non-accelerated methods in realistic problems, though the amount of speedup is highly problem-dependent. We have implemented our method in the new code Viscous Accretion Disk Evolution Resource (VADER), which is freely available for download from https://bitbucket.org/krumholz/vader/ under the terms of the GNU General Public License.

The Legacy ExtraGalactic UV Survey (LEGUS) is a Cycle 21 Treasury program on the Hubble Space Telescope, aimed at the investigation of star formation and its relation with galactic environment in nearby galaxies, from the scales of individual stars to those of ~kpc-size clustered structures. Five-band imaging, from the near-ultraviolet to the I-band, with the Wide Field Camera 3, plus parallel optical imaging with the Advanced Camera for Surveys, is being collected for selected pointings of 50 galaxies within the local 12 Mpc. The filters used for the observations with the Wide Field Camera 3 are: F275W(2,704 A), F336W(3,355 A), F438W(4,325 A), F555W(5,308 A), and F814W(8,024 A); the parallel observations with the Advanced Camera for Surveys use the filters: F435W(4,328 A), F606W(5,921 A), and F814W(8,057 A). The multi-band images are yielding accurate recent (<~50 Myr) star formation histories from resolved massive stars and the extinction-corrected ages and masses of star clusters and associations. The extensive inventories of massive stars and clustered systems will be used to investigate the spatial and temporal evolution of star formation within galaxies. This will, in turn, inform theories of galaxy evolution and improve the understanding of the physical underpinning of the gas-star formation relation and the nature of star formation at high redshift. This paper describes the survey, its goals and observational strategy, and the initial science results. Because LEGUS will provide a reference survey and a foundation for future observations with JWST and with ALMA, a large number of data products are planned for delivery to the community.

The enormous radiative and mechanical luminosities of massive stars impact a vast range of scales and processes, from the reionization of the universe, to the evolution of galaxies, to the regulation of the interstellar medium, to the formation of star clusters, and even to the formation of planets around stars in such clusters. Two main classes of massive star formation theory are under active study, Core Accretion and Competitive Accretion. In Core Accretion, the initial conditions are self-gravitating, centrally concentrated cores that condense with a range of masses from the surrounding, fragmenting clump environment. They then undergo relatively ordered collapse via a central disk to form a single star or a small-N multiple. In this case, the pre-stellar core mass function has a similar form to the stellar initial mass function. In Competitive Accretion, the material that forms a massive star is drawn more chaotically from a wider region of the clump without passing through a phase of being in a massive, coherent core. In this case, massive star formation must proceed hand in hand with star cluster formation. If stellar densities become very high near the cluster center, then collisions between stars may also help to form the most massive stars. We review recent theoretical and observational progress towards understanding massive star formation, considering physical and chemical processes, comparisons with low and intermediate-mass stars, and connections to star cluster formation.

Stellar feedback is often cited as the biggest uncertainty in galaxy formation models today. This uncertainty stems from a dearth of observational constraints as well as the great dynamic range between the small scales (<1 pc) where the feedback occurs and the large scales of galaxies (>1 kpc) that are shaped by this feedback. To bridge this divide, in this paper we aim to assess observationally the role of stellar feedback at the intermediate scales of HII regions. In particular, we employ multiwavelength data to examine several stellar feedback mechanisms in a sample of 32 HII regions in the Large and Small Magellanic Clouds (LMC and SMC, respectively). Using optical, infrared, radio, and X-ray images, we measure the pressures exerted on the shells from the direct stellar radiation, the dust-processed radiation, the warm ionized gas, and the hot X-ray emitting gas. We find that the warm ionized gas dominates over the other terms in all of the sources, although two have comparable dust-processed radiation pressures to their warm gas pressures. The hot gas pressures are comparatively weak, while the direct radiation pressures are 1-2 orders of magnitude below the other terms. We discuss the implications of these results, particularly highlighting evidence for hot gas leakage from the HII shells and regarding the momentum deposition from the dust-processed radiation to the warm gas. Furthermore, we emphasize that similar observational work should be done on very young HII regions to test whether direct radiation pressure and hot gas can drive the dynamics at early times.

The abundances of elements in stars are a critical clue to their origins. Observed star-to-star variations in logarithmic abundance within an open cluster are typically only $\sim 0.01-0.05$ over many elements, significantly smaller than the variation of $\sim 0.06-0.3$ seen in the interstellar medium from which the stars form. It is unknown why clusters are so homogenous, and whether homogeneity should also prevail in regions of lower star formation efficiency that do not produce bound clusters. Here we report adaptive mesh simulations using passively-advected scalars in order to trace the mixing of chemical elements as star-forming clouds form and collapse. We show that turbulent mixing during cloud assembly naturally produces a stellar abundance scatter at least ~5 times smaller than that in the gas, sufficient to fully explain the observed chemical homogeneity of stars. Moreover, mixing occurs very early, so that regions with efficiencies $\varepsilon \sim 10\%$ are nearly as well-mixed as those with $\varepsilon\sim 50\%$. This implies that even regions that do not form bound clusters are likely to be well-mixed, and enhances the prospects for using chemical tagging to reconstruct dissolved star clusters via their unique chemical signatures.

The integrated light of a stellar population, measured through photometric filters that are sensitive to the presence of young stars, is often used to infer the star formation rate (SFR) for that population. However, these techniques rely on an assumption that star formation is a continuous process, whereas in reality stars form in discrete spatially- and temporally-correlated structures. This discreteness causes the light output to undergo significant time-dependent fluctuations, which, if not accounted for, introduce systematic errors in the inferred SFRs due to the intrinsic distribution of luminosities at any fix SFR. We use SLUG, a code that Stochastically Lights Up Galaxies, to simulate galaxies undergoing stochastic star formation. We then use these simulations to present a quantitative analysis of these effects and provide tools for calculating probability distribution functions of SFRs given a set of observations. We show that, depending on the SFR tracer used, stochastic fluctuations can produce non-trivial errors at SFRs as high as 1 M/yr, and biases >0.5 dex at the lowest SFRs. We emphasize that due to the stochastic behavior of blue SFR tracers, one cannot assign a deterministic single value to the SFR of an individual galaxy, but we suggest methods by which future analyses that rely on integrated-light indicators can properly account for these stochastic effects.

We have used multi-wavelength Hubble Space Telescope WFC3 data of the starbursting spiral galaxy M83 in order to measure variations in the upper end of the stellar initial mass function (uIMF) using the production rate of ionizing photons in unresolved clusters with ages $\leq$ 8 Myr. As in earlier papers on M51 and NGC 4214, the upper end of the stellar IMF in M83 is consistent with an universal IMF, and stochastic sampling of the stellar populations in the $\lessapprox$ 10$^{3}$ Msun clusters are responsible for any deviations in this universality. The ensemble cluster population, as well as individual clusters, also imply that the most massive star in a cluster does not depend on the cluster mass. In fact, we have found that these small clusters seem to have an over-abundance of ionizing photons when compared to an expected universal or truncated IMF. This also suggests that the presence of massive stars in these clusters does not affect the star formation in a destructive way.

Star formation lies at the center of a web of processes that drive cosmic evolution: generation of radiant energy, synthesis of elements, formation of planets, and development of life. Decades of observations have yielded a variety of empirical rules about how it operates, but at present we have no comprehensive, quantitative theory. In this review I discuss the current state of the field of star formation, focusing on three central questions: what controls the rate at which gas in a galaxy converts to stars? What determines how those stars are clustered, and what fraction of the stellar population ends up in gravitationally-bound structures? What determines the stellar initial mass function, and does it vary with star-forming environment? I use these three question as a lens to introduce the basics of star formation, beginning with a review of the observational phenomenology and the basic physical processes. I then review the status of current theories that attempt to solve each of the three problems, pointing out links between them and opportunities for theoretical and numerical work that crosses the scale between them. I conclude with a discussion of prospects for theoretical progress in the coming years.

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, making the conversion to J=1-0 emission and hence molecular gas mass highly uncertain. Here, we combine numerical simulations of disk 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 (Sigma_SFR), which is observable. We use this result to develop a model for the CO SLED in terms of Sigma_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-6. This model should make it possible to significantly reduce the uncertainty in deducing molecular gas masses from observations of high-J CO emission.

Star clusters larger than $\sim 10^{3}$ $M_\odot$ contain multiple hot stars that launch fast stellar winds. The integrated kinetic energy carried by these winds is comparable to that delivered by supernova explosions, suggesting that at early times winds could be an important form of feedback on the surrounding cold material from which the star cluster formed. However, the interaction of these winds with the surrounding clumpy, turbulent, cold gas is complex and poorly understood. Here we investigate this problem via an accounting exercise: we use empirically determined properties of four well-studied massive star clusters to determine where the energy injected by stellar winds ultimately ends up. We consider a range of kinetic energy loss channels, including radiative cooling, mechanical work on the cold interstellar medium, thermal conduction, heating of dust via collisions by the hot gas, and bulk advection of thermal energy by the hot gas. We show that, for at least some of the clusters, none of these channels can account for more than a small fraction of the injected energy. We suggest that turbulent mixing at the hot-cold interface or physical leakage of the hot gas from the HII region can efficiently remove the kinetic energy injected by the massive stars in young star clusters. Even for the clusters where we are able to account for all the injected kinetic energy, we show that our accounting sets strong constraints on the importance of stellar winds as a mechanism for feedback on the cold interstellar medium.

In this chapter I review theoretical models for the formation of very massive stars. After a brief overview of some relevant observations, I spend the bulk of the chapter describing two possible routes to the formation of very massive stars: formation via gas accretion, and formation via collisions between smaller stars. For direct accretion, I discuss the problems of how interstellar gas may be prevented from fragmenting so that it is available for incorporation into a single very massive star, and I discuss the problems presented for massive star formation by feedback in the form of radiation pressure, photoionization, and stellar winds. For collision, I discuss several mechanisms by which stars might be induced to collide, and I discuss what sorts of environments are required to enable each of these mechanisms to function. I then compare the direct accretion and collision scenarios, and discuss possible observational signatures that could be used to distinguish between them. Finally, I come to the question of whether the process of star formation sets any upper limits on the masses of stars that can form.

We present three Orion simulations of star cluster formation in a 1000 Msun, turbulent molecular cloud clump, including the effects of radiative transfer, protostellar outflows, and magnetic fields. Our simulations all use self-consistent turbulent initial conditions and vary the mean mass-to-flux ratio relative to the critical value over 2, 10, and infinity to gauge the influence of magnetic fields on star cluster formation. We find, in good agreement with previous studies, that magnetic fields of typically observed strengths lower the star formation rate by a factor of 2.4 and reduce the amount of fragmentation by a factor of 2 relative to the zero-field case. We also find that the field increases the characteristic sink particle mass, again by a factor of 2.4. The magnetic field also increases the degree of clustering in our simulations, such that the maximum stellar densities in the strong field case are higher than the others by again a factor of 2. This clustering tends to encourage the formation of multiple systems, which are more common in the rad-MHD runs than the rad-hydro run. The companion frequency in our simulations is consistent with observations of multiplicity in Class I sources, particularly for the strong field case. Finally, we find evidence of primordial mass segregation in our simulations reminiscent of that observed in star clusters like the Orion Nebula Cluster.

Stars do not generally form in isolation. Instead, they form in clusters, and in these clustered environments newborn stars can have profound effects on one another and on their parent gas clouds. Feedback from clustered stars is almost certainly responsible for a number of otherwise puzzling facts about star formation: that it is an inefficient process that proceeds slowly when averaged over galactic scales; that most stars disperse from their birth sites and dissolve into the galactic field over timescales $\ll 1$ Gyr; and that newborn stars follow an initial mass function (IMF) with a distinct peak in the range $0.1 - 1$ $M_\odot$, rather than an IMF dominated by brown dwarfs. In this review we summarize current observational constraints and theoretical models for the complex interplay between clustered star formation and feedback.

We introduce the AGORA project, a comprehensive numerical study of well-resolved galaxies within the LCDM cosmology. Cosmological hydrodynamic simulations with force resolutions of ~100 proper pc or better will be run with a variety of code platforms to follow the hierarchical growth, star formation history, morphological transformation, and the cycle of baryons in and out of 8 galaxies with halo masses M_vir ~= 1e10, 1e11, 1e12, and 1e13 Msun at z=0 and two different ("violent" and "quiescent") assembly histories. The numerical techniques and implementations used in this project include the smoothed particle hydrodynamics codes GADGET and GASOLINE, and the adaptive mesh refinement codes ART, ENZO, and RAMSES. The codes will share common initial conditions and common astrophysics packages including UV background, metal-dependent radiative cooling, metal and energy yields of supernovae, and stellar initial mass function. These are described in detail in the present paper. Subgrid star formation and feedback prescriptions will be tuned to provide a realistic interstellar and circumgalactic medium using a non-cosmological disk galaxy simulation. Cosmological runs will be systematically compared with each other using a common analysis toolkit, and validated against observations to verify that the solutions are robust - i.e., that the astrophysical assumptions are responsible for any success, rather than artifacts of particular implementations. The goals of the AGORA project are, broadly speaking, to raise the realism and predictive power of galaxy simulations and the understanding of the feedback processes that regulate galaxy "metabolism." The proof-of-concept dark matter-only test of the formation of a galactic halo with a z=0 mass of M_vir ~= 1.7e11 Msun by 9 different versions of the participating codes is also presented to validate the infrastructure of the project.

Giant molecular clouds (GMCs) are the primary reservoirs of cold, star-forming molecular gas in the Milky Way and similar galaxies, and thus any understanding of star formation must encompass a model for GMC formation, evolution, and destruction. These models are necessarily constrained by measurements of interstellar molecular and atomic gas, and the emergent, newborn stars. Both observations and theory have undergone great advances in recent years, the latter driven largely by improved numerical simulations, and the former by the advent of large-scale surveys with new telescopes and instruments. This chapter offers a thorough review of the current state of the field.

We present a simple toy model to understand what sets the scatter in star formation and metallicity of galaxies at fixed mass. The scatter ultimately arises from the intrinsic scatter in the accretion rate, but may be substantially reduced depending on the timescale on which the accretion varies compared to the timescale on which the galaxy loses gas mass. This model naturally produces an anti-correlation between star formation and metallicity at a fixed mass, the basis of the fundamental metallicity relation. We show that observational constraints on the scatter in galaxy scaling relations can be translated into constraints on the galaxy-to-galaxy variation in the mass loading factor, and the timescales and magnitude of stochastic accretion onto star-forming galaxies. We find a remarkably small scatter in the mass loading factor, < 0.1 dex, and that the scatter in accretion rates is smaller than expected from N-body simulations.

I describe DESPOTIC, a code to Derive the Energetics and SPectra of Optically Thick Interstellar Clouds. DESPOTIC represents such clouds using a one-zone model, and can calculate line luminosities, line cooling rates, and in restricted cases line profiles using an escape probability formalism. It also includes approximate treatments of the dominant heating, cooling, and chemical processes for the cold interstellar medium, including cosmic ray and X-ray heating, grain photoelectric heating, heating of the dust by infrared and ultraviolet radiation, thermal cooling of the dust, collisional energy exchange between dust and gas, and a simple network for carbon chemistry. Based on these heating, cooling, and chemical rates, DESPOTIC can calculate clouds' equilibrium gas and dust temperatures, equilibrium carbon chemical state, and time-dependent thermal and chemical evolution. The software is intended to allow rapid and interactive calculation of clouds' characteristic temperatures, identification of their dominant heating and cooling mechanisms, and prediction of their observable spectra across a wide range of interstellar environments. DESPOTIC is implemented as a Python package, and is released under the GNU General Public License.

In this paper, I investigate the processes that regulate the rate of star formation in regions of galaxies where the neutral interstellar medium is predominantly composed of non-star-forming HI. In such regions, found today predominantly in low-metallicity dwarf galaxies and in the outer parts of large spirals, the star formation rate per unit area and per unit mass is much smaller than in more molecule-rich regions. While in molecule-rich regions the ultraviolet radiation field produced by efficient star formation forces the density of the cold neutral medium to a value set by two-phase equilibrium, I show that the low rates of star formation found in molecule-poor regions preclude this condition. Instead, the density of the cold neutral gas is set by the requirements of hydrostatic balance. Using this result, I extend the Krumholz, McKee, & Tumlinson model for star formation and the atomic to molecular transition to the molecule-poor regime. This "KMT+" model matches a wide range of observations of the star formation rate and the balance between the atomic and molecular phases in dwarfs and in the outer parts of spirals, and is well-suited to implementation as a subgrid recipe for star formation in cosmological simulations and semi-analytic models. I discuss the implications of this model for star formation over cosmological times.

We present a cosmological hydrodynamic simulation of the formation of dwarf galaxies at redshifts z>~2.5 using a physically-motivated model for H2-regulated star formation. Our simulation, performed using the Enzo code and reaching a peak resolution of 109 proper parsecs at z=2.5, extends the results of Kuhlen et al. (2012) to significantly lower redshifts. We show that a star formation prescription regulated by the local H2 abundance leads to the suppression of star formation in dwarf galaxy halos with M_h <~ 10^10 Msun and to a large population of gas-rich "dark galaxies" at z=2.5 with low star formation efficiencies and gas depletion timescales >20 Gyr. The fraction of dark galaxies is 60% at M_h ~ 10^10 Msun and increases rapidly with decreasing halo mass. Dark galaxies form late and their gaseous disks never reach the surface densities, > ~5700 Msun / pc^2 (Z/10^-3 Zsun)^(-0.88), that are required to build a substantial molecular fraction. Despite this large population of dark galaxies, we show that our H2-regulated simulation is consistent with both the observed luminosity function of galaxies and the cosmological mass density of neutral gas at z>~2.5. Moreover, our results provide a theoretical explanation for the recent detection in fluorescent Ly-alpha emission of gaseous systems at high redshift with little or no associated star formation. We further propose that H2-regulation may offer a fresh solution to a number of outstanding "dwarf galaxy problems" in LambdaCDM. In particular, H2-regulation leads galaxy formation to become effectively stochastic on mass scales of M_h ~ 10^10 Msun, and thus these massive dwarfs are not "too big to fail".

[abridged] Radiation pressure on dust grains may be an important mechanism in driving winds in a wide variety of astrophysical systems. However, the efficiency of the coupling between the radiation field and the dusty gas is poorly understood in environments characterized by high optical depths. We present a series of idealized numerical experiments, performed with the radiation-hydrodynamic code ORION, in which we study the dynamics of such winds and quantify their properties. We find that, after wind acceleration begins, radiation Rayleigh-Taylor instability forces the gas into a configuration that reduces the rate of momentum transfer from the radiation field to the gas by a factor ~ 10 - 100 compared to an estimate based on the optical depth at the base of the atmosphere; instead, the rate of momentum transfer from a driving radiation field of luminosity L to the gas is roughly L/c multiplied by one plus half the optical depth evaluated using the photospheric temperature, which is far smaller than the optical depth one would obtain using the interior temperature. When we apply our results to conditions appropriate to ULIRGs and star clusters, we find that the asymptotic wind momentum flux from such objects should not significantly exceed that carried by the direct radiation field, L/c. This result constrains the expected mass loss rates from systems that exceed the Eddington limit to be of order the so-called "single-scattering" limit, and not significantly higher. We present an approximate fitting formula for the rate of momentum transfer from radiation to dusty gas through which it passes, which is suitable for implementation in sub-grid models of galaxy formation. Finally, we provide a first map of the column density distribution of gas in a radiatively-driven wind as a function of velocity, and velocity dispersion.

Over the past 10 Gyr, star-forming galaxies have changed dramatically, from clumpy and gas rich, to rather quiescent stellar-dominated disks with specific star formation rates lower by factors of a few tens. We present a general theoretical model for how this transition occurs, and what physical processes drive it, making use of 1D axisymmetric thin disk simulations with an improved version of the Gravitational Instability-Dominated Galaxy Evolution Tool (GIDGET) code. We show that at every radius galaxies tend to be in a slowly evolving equilibrium state wherein new accretion is balanced by star formation, galactic winds, and radial transport of gas through the disk by gravitational instability (GI) -driven torques. The gas surface density profile is determined by which of these terms are in balance at a given radius - direct accretion is balanced by star formation and galactic winds near galactic centers, and by transport at larger radii. We predict that galaxies undergo a smooth transition from a violent disk instability phase to secular evolution. This model provides a natural explanation for the high velocity dispersions and large clumps in z~2 galaxies, the growth and subsequent quenching of bulges, and features of the neutral gas profiles of local spiral galaxies.

Both simulations and observations indicate that stars form in filamentary, hierarchically clustered associations, most of which disperse into their galactic field once feedback destroys their parent clouds. However, during their early evolution in these substructured environments, stars can undergo close encounters with one another that might have significant impacts on their protoplanetary disks or young planetary systems. We perform N-body simulations of the early evolution of dissolving, substructured clusters with a wide range of properties, with the aim of quantifying the expected number and orbital element distributions of encounters as a function of cluster properties. We show that the presence of substructure both boosts the encounter rate and modifies the distribution of encounter velocities compared to what would be expected for a dynamically relaxed cluster. However, the boost only lasts for a dynamical time, and as a result the overall number of encounters expected remains low enough that gravitational stripping is unlikely to be a significant effect for the vast majority of star-forming environments in the Galaxy. We briefly discuss the implications of this result for models of the origin of the Solar System, and of free-floating planets. We also provide tabulated encounter rates and orbital element distributions suitable for inclusion in population synthesis models of planet formation in a clustered environment.

Late in their evolution, massive stars may undergo periods of violent instability and mass loss, but the mechanism responsible for these episodes has not been identified. We study one potential contributor: the development of local radiation-driven instabilities in the outer layers of main sequence (MS) and post-MS massive stars. We construct a sequence of massive stellar evolution models and investigate where they are subject to local radiative instabilities, both in the presence of magnetic fields and without them,and at a range of metallicities. We find that these types of instabilities do not occur in solar-metallicity MS stars up to 100\,M$_\odot$, but they set in immediately post-MS for stars heavier than $\sim 25$\,M$_\odot$. Once an instability appears, it involves a significant amount of mass in the star's upper layers (up to $\sim 1$ per cent of the initial stellar mass), suggesting that radiation-driven instabilities are a potentially viable mechanism for dynamic mass loss. We find that the presence of magnetic fields at strengths low enough not to disturb the hydrostatic balance of the star does not alter these results. Stars with sub-Solar metallicity also show instability, but their instabilities involve less mass and appear later in the star's evolution.

Optical and infrared emission lines from HII regions are an important diagnostic used to study galaxies, but interpretation of these lines requires significant modeling of both the internal structure and dynamical evolution of the emitting regions. Most of the models in common use today assume that HII region dynamics are dominated by the expansion of stellar wind bubbles, and have neglected the contribution of radiation pressure to the dynamics, and in some cases also to the internal structure. However, recent observations of nearby galaxies suggest that neither assumption is justified, motivating us to revisit the question of how HII region line emission depends on the physics of winds and radiation pressure. In a companion paper we construct models of single HII regions including and excluding radiation pressure and winds, and in this paper we describe a population synthesis code that uses these models to simulate galactic collections of HII regions with varying physical parameters. We show that the choice of physical parameters has significant effects on galactic emission line ratios, and that in some cases the line ratios can exceed previously claimed theoretical limits. Our results suggest that the recently-reported offset in line ratio values between high-redshift star-forming galaxies and those in the local universe may be partially explained by the presence of large numbers of radiation pressured-dominated HII regions within them.

The emission line ratios [OIII]5007/H-beta and [NII]6584/H-alpha have been adopted as an empirical way to distinguish between the fundamentally different mechanisms of ionization in emission-line galaxies. However, detailed interpretation of these diagnostics requires calculations of the internal structure of the emitting HII regions, and these calculations depend on the assumptions one makes about the relative importance of radiation pressure and stellar winds. In this paper we construct a grid of quasi-static HII region models to explore how choices about these parameters alter HII regions' emission line ratios. We find that, when radiation pressure is included in our models, HII regions reach a saturation point beyond which further increases in the luminosity of the driving stars does not produce any further increase in effective ionization parameter, and thus does not yield any further alteration in an HII region's line ratio. We also show that, if stellar winds are assumed to be strong, the maximum possible ionization parameter is quite low. As a result of this effect, it is inconsistent to simultaneously assume that HII regions are wind-blown bubbles and that they have high ionization parameters; some popular HII region models suffer from this inconsistency. Our work in this paper provides a foundation for a companion paper in which we embed the model grids we compute here within a population synthesis code that enables us to compute the integrated line emission from galactic populations of HII regions.

We predict the evolution of giant clumps undergoing star-driven outflows in high-z gravitationally unstable disk galaxies. We find that the mass loss is expected to occur through a steady wind over many tens of free-fall times (t_ff ~ 10 Myr) rather than by an explosive disruption in one or a few t_ff. Our analysis is based on the finding from simulations that radiation trapping is negligible because it destabilizes the wind (Krumholz & Thompson 2012, 2013). Each photon can therefore contribute to the wind momentum only once, so the radiative force is limited to L/c. When combining radiation, protostellar and main-sequence winds, and supernovae, we estimate the total direct injection rate of momentum into the outflow to be 2.5 L/c. The adiabatic phase of supernovae and main-sequence winds can double this rate. The resulting outflow mass-loading factor is of order unity, and if the clumps were to deplete their gas the timescale would have been a few disk orbital times, to end with half the original clump mass in stars. However, the clump migration time to the disk center is on the order of an orbital time, about 250 Myr, so the clumps are expected to complete their migration prior to depletion. Furthermore, the clumps are expected to double their mass in a disk orbital time by accretion from the disk and clump-clump mergers, so their mass actually grows in time and with decreasing radius. From the 6-7 giant clumps with observed outflows, 5 are consistent with these predictions, and one has a much higher mass-loading factor and momentum injection rate. The latter either indicates that the estimated outflow is an overestimate (within the 1-sigma error), that the SFR has dropped since the time when the outflow was launched, or that the driving mechanism is different, e.g. supernova feedback in a cavity generated by the other feedbacks.

Recent studies have claimed the existence of very massive stars (VMS) up to 300 solar masses in the local Universe. As this finding may represent a paradigm shift for the canonical stellar upper-mass limit of 150 Msun, it is timely to discuss the status of the data, as well as the far-reaching implications of such objects. We held a Joint Discussion at the General Assembly in Beijing to discuss (i) the determination of the current masses of the most massive stars, (ii) the formation of VMS, (iii) their mass loss, and (iv) their evolution and final fate. The prime aim was to reach broad consensus between observers and theorists on how to identify and quantify the dominant physical processes.

The production rate of ionizing photons in young < 8 Myr, unresolved stellar clusters in the nearby irregular galaxy NGC 4214 is probed using multi-wavelength Hubble Space Telescope WFC3 data. We normalize the ionizing photon rate by the cluster mass to investigate the upper end of the stellar initial mass function (IMF). We have found that within the uncertainties the upper end of the stellar IMF appears to be universal in this galaxy, and that deviations from a universal IMF can be attributed to stochastic sampling of stars in clusters with masses < 10^3 M_sun. Furthermore, we have found that there does not seem to be a dependence of the maximum stellar mass on the cluster mass. We have also found that for massive clusters, feedback may cause an underrepresentation in Ha luminosities, which needs to be taken into account when conducting this type of analysis.

We present a set of 3-dimensional, radiation-magnetohydrodynamic calculations of the gravitational collapse of massive (300 Msun), star-forming molecular cloud cores. We show that the combined effects of magnetic fields and radiative feedback strongly suppress core fragmentation, leading to the production of single star systems rather than small clusters. We find that the two processes are efficient at suppressing fragmentation in different regimes, with the feedback most effective in the dense, central region and the magnetic field most effective in more diffuse, outer regions. Thus, the combination of the two is much more effective at suppressing fragmentation than either one considered in isolation. Our work suggests that typical massive cores, which have mass-to-flux ratios of about 2 relative to critical, likely form a single star system, but that cores with weaker fields may form a small star cluster. This result helps us understand why the observed relationship between the core mass function and the stellar initial mass function holds even for ~100 Msun cores with many thermal Jeans masses of material. We also demonstrate that a ~40 AU Keplerian disk is able to form in our simulations, despite the braking effect caused by the strong magnetic field.

The gas from which stars form is magnetized, and strong magnetic fields can efficiently transport angular momentum. Most theoretical models of this phenomenon find that it should prevent formation of large (>100 AU), rotationally-supported disks around most protostars, even when non-ideal magnetohydrodynamic (MHD) effects that allow the field and gas to decouple are taken into account. Using recent observations of magnetic field strengths and orientations in protostellar cores, we show that this conclusion is incorrect. The distribution of magnetic field strengths is very broad, and alignments between fields and angular momentum vectors within protostellar cores are essentially random. By combining the field strength and misalignment data with MHD simulations showing that disk formation is expected for both weak and misaligned fields, we show that these observations imply that we should expect disk fractions of ~10 - 50% even when protostars are still deeply embedded in their parent cores, and even if the gas is governed by ideal MHD.

[abridged] The pressure exerted by the radiation of young stars may be an important feedback mechanism in forming star clusters and the disks of starburst galaxies. However, there is great uncertainty in how efficiently radiation couples to matter in these high optical depth environments. In particular, it is unclear what levels of turbulence the radiation can produce, and whether the infrared radiation trapped by the dust opacity can give rise to heavily mass-loaded winds. 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. We consider both systems where the radiation flux is sub-Eddington throughout the gas column, and where it is super-Eddington at the midplane but sub-Eddington in the atmosphere. In the latter, we find that the radiation-matter interaction gives rise to radiation-driven Rayleigh-Taylor instability, which drives supersonic turbulence at a level sufficient to fully explain the turbulence seen in Galactic protocluster gas clouds, and to make a non-trivial contribution to the turbulence observed in starburst galaxy disks. However, the instability also produces a channel structure in which the radiation-matter interaction is reduced because the radiation field is not fully trapped. For astrophysical parameters relevant to forming star clusters and starburst galaxies, we find that this effect reduces the net momentum deposition rate in the dusty gas by a factor of ~2-6 compared to simple analytic estimates, and that in steady state the Eddington ratio reaches unity and there are no strong winds. 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.

Observations of nearby galaxies have firmly established, over a broad range of galactic environments and metallicities, that star formation occurs exclusively in the molecular phase of the interstellar medium (ISM). Theoretical models show that this association results from the correlation between chemical phase, shielding, and temperature. Interstellar gas converts from atomic to molecular only in regions that are well shielded from interstellar ultraviolet (UV) photons, and since UV photons are also the dominant source of interstellar heating, only in these shielded regions does the gas become cold enough to be subject to Jeans instability. However, while the equilibrium temperature and chemical state of interstellar gas are well-correlated, the time scale required to reach chemical equilibrium is much longer than that required to reach thermal equilibrium, and both timescales are metallicity-dependent. Here I show that the difference in time scales implies that, at metallicities below a few percent of the Solar value, well-shielded gas will reach low temperatures and proceed to star formation before the bulk of it is able to convert from atomic to molecular. As a result, at extremely low metallicities, star formation will occur in a cold atomic phase of the ISM rather than a molecular phase. I calculate the observable consequences of this result for star formation in low metallicity galaxies, and I discuss how some current numerical models for H2-regulated star-formation may need to be modified.

Observations show that radial metallicity gradients in disk galaxies are relatively shallow, if not flat, especially at large galactocentric distances and for galaxies in the high-redshift universe. Given that star formation and metal production are centrally concentrated, this requires a mechanism to redistribute metals. However, the nature of this mechanism is poorly understood, let alone quantified. To address this problem, we conduct magnetohydrodynamical simulations of a local shearing sheet of a thin, thermally unstable, gaseous disk driven by a background stellar spiral potential, including metals modeled as passive scalar fields. Contrary to what a simple \alpha\ prescription for the gas disk would suggest, we find that turbulence driven by thermal instability is very efficient at mixing metals, regardless of the presence or absence of stellar spiral potentials or magnetic fields. The timescale for homogenizing randomly distributed metals is comparable to or less than the local orbital time in the disk. This implies that turbulent mixing of metals is a significant process in the history of chemical evolution of disk galaxies.

The Hubble Deep Field (HDF) is a region in the sky that provides one of the deepest multi-wavelength views of the distant universe and has led to the detection of thousands of galaxies seen throughout cosmic time. An early map of the HDF at a wavelength of 850 microns that is sensitive to dust emission powered by star formation revealed the brightest source in the field, dubbed HDF850.1. For more than a decade, this source remained elusive and, despite significant efforts, no counterpart at shorter wavelengths, and thus no redshift, size or mass, could be identified. Here we report, using a millimeter wave molecular line scan, an unambiguous redshift determination for HDF850.1 of z=5.183. This places HDF850.1 in a galaxy overdensity at z~5.2 in the HDF, corresponding to a cosmic age of only 1.1 Gyr after the Big Bang. This redshift is significantly higher than earlier estimates and higher than most of the >100 sub-millimeter bright galaxies identified to date. The source has a star formation rate of 850 M_sun/yr and is spatially resolved on scales of 5 kpc, with an implied dynamical mass of ~1.3x10^11 M_sun, a significant fraction of which is present in the form of molecular gas. Despite our accurate redshift and position, a counterpart arising from starlight remains elusive.

[abridged] We report a series of simulations of the formation of a star cluster similar to the Orion Nebula Cluster (ONC), including both radiative transfer and protostellar outflows, and starting from both smooth and self-consistently turbulent initial conditions. Each simulation forms >150 stars and brown dwarfs, yielding a stellar mass distribution from < 0.1 to > 10 Msun. We show that a simulation that begins with self-consistently turbulence embedded in a larger turbulent volume, and that includes protostellar outflows, produces an initial mass function (IMF) consistent both with that of the ONC and the Galactic field. This is the first simulation published to date that reproduces the observed IMF in a cluster large enough to contain massive stars, and where the result is determined by a fully self-consistent calculation of gas thermodynamics. This simulation also produces a star formation rate that, while still somewhat too high, is much closer to observed values than if we omit either the larger turbulent volume or the outflows. Moreover, we show that the combination of outflows, self-consistently turbulent initial conditions, and turbulence continually fed by motions on scales larger than that of the protocluster yields an IMF that is in agreement with observations and invariant with time, resolving the "overheating" problem in which simulations without these features have an IMF peak that shifts to progressively higher masses over time. The simulation that matches the observed IMF also reproduces the observed trend of stellar multiplicity strongly increasing with mass. We show that this simulation produces massive stars from distinct massive cores whose properties are consistent with those of observed massive cores. However, the stars formed in these cores also undergo dynamical interactions that naturally produce Trapezium-like hierarchical multiple systems.

A number of authors have argued that the Sun must have been born in a cluster of no more than about 1000 stars, on the basis that, in a larger cluster, close encounters between the Sun and other stars would have truncated the outer Solar System or excited the outer planets into eccentric orbits. However, this dynamical limit is in tension with meteoritic evidence that the Solar System was exposed to a nearby supernova during or shortly after its formation; a 1000-star cluster is much too small for supernova contamination to be likely. In this paper we revisit the dynamical limit in the light of improved observations of the properties of young clusters. We use a series of scattering simulations to measure the velocity-dependent cross-section for disruption of the outer Solar System by stellar encounters, and use this cross-section to compute the probability of a disruptive encounter as a function of birth cluster properties. We find that, contrary to prior work, the probability of disruption is small regardless of the cluster mass, and that it actually decreases rather than increases with cluster mass. Our results differ from prior work for three main reasons: (1) unlike in most previous work, we compute a velocity-dependent cross section and properly integrate over the cluster mass-dependent velocity distribution of incoming stars; (2) we adopt realistically-short cluster lifetimes of a few crossing times, rather than assuming lifetimes of 10 to 100 Myr; and (3) following recent observations, we adopt a mass-independent surface density for embedded clusters, rather than a mass-independent radius as assumed many earlier papers. Our results remove the tension between the dynamical limit and the meteoritic evidence, and suggest that the Sun was born in a massive cluster.

Observations of disk galaxies at z~2 have demonstrated that turbulence driven by gravitational instability can dominate the energetics of the disk. We present a 1D simulation code, which we have made publicly available, that economically evolves these galaxies from z~2 to z~0 on a single CPU in a matter of minutes, tracking column density, metallicity, and velocity dispersions of gaseous and multiple stellar components. We include an H$_2$ regulated star formation law and the effects of stellar heating by transient spiral structure. We use this code to demonstrate a possible explanation for the existence of a thin and thick disk stellar population and the age-velocity dispersion correlation of stars in the solar neighborhood: the high velocity dispersion of gas in disks at z~2 decreases along with the cosmological accretion rate, while at lower redshift, the dynamically colder gas forms the low velocity dispersion stars of the thin disk.

[abridged] The star formation rates (SFR) 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. Our model shows that including the effect of metallicity leads to suppression of star formation at redshift z>2 in dark halos with masses <~ 10^11 Msun, with the suppression becoming near total for halos below ~10^9.5-10 Msun. We find that at high redshift the SFR cannot catch up with the gas inflow rate (IR), because the SFR is limited by the free-fall time, and because it is suppressed further by a lack of metals. As a result, in each galaxy the SFR is growing in time faster than the IR, and the integrated cosmic SFR density is rising with time. The suppressed in situ SFR at high z makes the growth of stellar mass dominated by ex situ SFR which implies that the specific SFR (sSFR) remains constant with time. The intensely accreted gas at high z is accumulating as an atomic gas reservoir. This provides additional fuel for star formation in 10^10 - 10^12 Msun halos at z ~ 1-3, which allows the SFR to exceed the instantaneous IR, and may enable an even higher outflow rate. At z<1, following the natural decline in IR with time due to the universal expansion, the SFR and sSFR are expected to drop. 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.

We describe cosmological galaxy formation simulations with the adaptive mesh refinement code Enzo that incorporate a star formation prescription regulated by the local abundance of molecular hydrogen. We show that this H2-regulated prescription leads to a suppression of star formation in low mass halos (M_h < ~10^10 M_sun) at z>4, alleviating some of the dwarf galaxy problems faced by theoretical galaxy formation models. H2 regulation modifies the efficiency of star formation of cold gas directly, rather than indirectly reducing the cold gas content with "supernova feedback". We determine the local H2 abundance in our most refined grid cells (76 proper parsec in size at z=4) by applying the model of Krumholz, McKee, & Tumlinson, which is based on idealized 1D radiative transfer calculations of H2 formation-dissociation balance in ~100 pc atomic--molecular complexes. Our H2-regulated simulations are able to reproduce the empirical (albeit lower z) Kennicutt-Schmidt relation, including the low Sigma_gas cutoff due to the transition from atomic to molecular phase and the metallicity dependence thereof, without the use of an explicit density threshold in our star formation prescription. We compare the evolution of the luminosity function, stellar mass density, and star formation rate density from our simulations to recent observational determinations of the same at z=4-8 and find reasonable agreement between the two.

The physical mechanisms that set the initial rotation rates in massive stars are a crucial unknown in current star formation theory. Observations of young, massive stars provide evidence that they form in a similar fashion to their low-mass counterparts. The magnetic coupling between a star and its accretion disk may be sufficient to spin down low-mass pre-main sequence (PMS) stars to well below breakup at the end stage of their formation when the accretion rate is low. However, we show that these magnetic torques are insufficient to spin down massive PMS stars due to their short formation times and high accretion rates. We develop a model for the angular momentum evolution of stars over a wide range in mass, considering both magnetic and gravitational torques. We find that magnetic torques are unable to spin down either low or high mass stars during the main accretion phase, and that massive stars cannot be spun down significantly by magnetic torques during the end stage of their formation either. Spin-down occurs only if massive stars' disk lifetimes are substantially longer or their magnetic fields are much stronger than current observations suggest.

The most common means of converting an observed CO line intensity into a molecular gas mass requires the use of a conversion factor (Xco). While in the Milky Way this quantity does not appear to vary significantly, there is good reason to believe that Xco will depend on the larger-scale galactic environment. Utilising numerical models, we investigate how varying metallicities, gas temperatures and velocity dispersions in galaxies impact the way CO line emission traces the underlying H2 gas mass, and under what circumstances Xco may differ from the Galactic mean value. We find that, due to the combined effects of increased gas temperature and velocity dispersion, Xco is depressed below the Galactic mean in high surface density environments such as ULIRGs. In contrast, in low metallicity environments, Xco tends to be higher than in the Milky Way, due to photodissociation of CO in metal-poor clouds. At higher redshifts, gas-rich discs may have gravitationally unstable clumps which are warm (due to increased star formation) and have elevated velocity dispersions. These discs tend to have Xco values ranging between present-epoch gas-rich mergers and quiescent discs at low-z. This model shows that on average, mergers do have lower Xco values than disc galaxies, though there is significant overlap. Xco varies smoothly with the local conditions within a galaxy, and is not a function of global galaxy morphology. We combine our results to provide a general fitting formula for Xco as a function of CO line intensity and metallicity. We show that replacing the traditional approach of using one constant Xco for starbursts and another for discs with our best-fit function produces star formation laws that are continuous rather than bimodal, and that have significantly reduced scatter.

We have carried out a numerical study of the effect of large scale magnetic fields on the rate of accretion from a uniform, isothermal gas onto a resistive, stationary point mass. Only mass, not magnetic flux, accretes onto the point mass. The simulations for this study avoid complications arising from boundary conditions by keeping the boundaries far from the accreting object. Our simulations leverage adaptive refinement methodology to attain high spatial fidelity close to the accreting object. Our results are particularly relevant to the problem of star formation from a magnetized molecular cloud in which thermal energy is radiated away on time scales much shorter than the dynamical time scale. Contrary to the adiabatic case, our simulations show convergence toward a finite accretion rate in the limit in which the radius of the accreting object vanishes, regardless of magnetic field strength. For very weak magnetic fields, the accretion rate first approaches the Bondi value and then drops by a factor ~ 2 as magnetic flux builds up near the point mass. For strong magnetic fields, the steady-state accretion rate is reduced by a factor ~ 0.2 \beta^{1/2} compared to the Bondi value, where \beta is the ratio of the gas pressure to the magnetic pressure. We give a simple expression for the accretion rate as a function of the magnetic field strength. Approximate analytic results are given in the Appendixes for both time-dependent accretion in the limit of weak magnetic fields and steady-state accretion for the case of strong magnetic fields.

[abridged] While observations of Local Group galaxies show a very simple, local star formation law in which the star formation rate per unit area in each patch of a galaxy scales linearly with the molecular gas surface density, recent observations of both Milky Way molecular clouds and high redshift galaxies apparently show a more complicated relationship, in which regions of equal surface density can form stars at quite different rates. These data have been interpreted as implying either that different star formation laws apply in different circumstances, that the star formation law is sensitive to large-scale galaxy properties rather than local properties, or that there are high density thresholds for star formation. Here we collate resolved observations of Milky Way molecular clouds, kpc-scale observations of Local Group galaxies, and unresolved observations of both disk and starburst galaxies in the local universe and at high redshift. We show that all of these data are in fact consistent with a simple, local, volumetric star formation law. The apparent variations stem from the fact that the observed objects have a wide variety of 3D size scales and degrees of internal clumping, so even at fixed gas column density the regions being observed can have wildly varying volume densities. We provide a simple theoretical framework to remove this projection effect, and we show that all the data, from small Solar neighborhood clouds with masses ~10^3 Msun to sub-mm galaxies with masses ~10^11 Msun, fall on a single star formation law in which the SFR is simply ~1% of the molecular gas mass per local free-fall time. In contrast, proposed star formation laws in which the star formation timescale is set by the galactic rotation period or the SFR is linearly proportional to the gas mass above some density threshold fail to match at least some of the data.

We use the three-dimensional Athena ionizing radiation-magnetohydrodynamics (IRMHD) code to simulate blister-type HII regions driven by stars on the edge of magnetized gas clouds. We compare these to simulations of spherical HII regions where the star is embedded deep within a cloud, and to non-magnetized simulations of both types, in order to compare their ability to drive turbulence and influence star formation. We find that magnetized blister HII regions can be very efficient at injecting energy into clouds. This is partly a magnetic effect: the magnetic energy added to a cloud by an HII region is comparable to or larger than the kinetic energy, and magnetic fields can also help collimate the ejected gas, increasing its energy yield. As a result of these effects, a blister HII region expanding into a cloud with a magnetic field perpendicular to its edge injects twice as much energy by 5 Myr as a non-magnetized blister HII region driven by a star of the same luminosity. Blister HII regions are also more efficient at injecting kinetic energy than spherical HII regions, due to the recoil provided by escaping gas, but not by as much as predicted by some analytic approximations.

Recent observations indicate a lower Halpha to FUV ratio in dwarf galaxies than in brighter systems, a trend that could be explained by a truncated and/or steeper IMF in small galaxies. However, at low star formation rates (SFRs), the Halpha to FUV ratio can vary due to stochastic sampling even for a universal IMF, a hypothesis that has, prior to this work, received limited investigation. Using SLUG, a fully stochastic code for synthetic photometry in star clusters and galaxies, we compare the Halpha and FUV luminosity in a sample of ~450 nearby galaxies with models drawn from a universal Kroupa IMF and a modified IMF, the integrated galactic initial mass function (IGIMF). Once random sampling and time evolution are included, a Kroupa IMF convolved with the cluster mass function reproduces the observed Halpha distribution at all FUV luminosities, while a truncated IMF as implemented in current IGIMF models underpredicts the Halpha luminosity by more than an order of magnitude at the lowest SFRs. We conclude that the observed luminosity is the result of the joint probability distribution function of the SFR, cluster mass function, and a universal IMF, consistent with parts of the IGIMF theory, but that a truncation in the IMF in clusters is inconsistent with the observations. Future work will examine stochastic star formation and its time dependence in detail to study whether random sampling can explain other observations that suggest a varying IMF.

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.

We report the results of a series of AMR radiation-hydrodynamic simulations of the collapse of massive star forming clouds using the ORION code. These simulations are the first to include the feedback effects protostellar outflows, as well as protostellar radiative heating and radiation pressure exerted on the infalling, dusty gas. We find that that outflows evacuate polar cavities of reduced optical depth through the ambient core. These enhance the radiative flux in the poleward direction so that it is 1.7 to 15 times larger than that in the midplane. As a result the radiative heating and outward radiation force exerted on the protostellar disk and infalling cloud gas in the equatorial direction are greatly diminished. The simultaneously reduces the Eddington radiation pressure barrier to high-mass star formation and increases the minimum threshold surface density for radiative heating to suppress fragmentation compared to models that do not include outflows. The strength of both these effects depends on the initial core surface density. Lower surface density cores have longer free-fall times and thus massive stars formed within them undergo more Kelvin contraction as the core collapses, leading to more powerful outflows. Furthermore, in lower surface density clouds the ratio of the time required for the outflow to break out of the core to the core free-fall time is smaller, so that these clouds are consequently influenced by outflows at earlier stages of collapse. As a result, outflow effects are strongest in low surface density cores and weakest in high surface density one. We also find that radiation focusing in the direction of outflow cavities is sufficient to prevent the formation of radiation pressure-supported circumstellar gas bubbles, in contrast to models which neglect protostellar outflow feedback.

Relating the observed CO emission from giant molecular clouds (GMCs) to the underlying H2 column density is a long-standing problem in astrophysics. While the Galactic CO-H2 conversion factor (Xco) appears to be reasonably constant, observations indicate that Xco may be depressed in high-surface density starburst environments. Using a multi-scale approach, we investigate the dependence of Xco on the galactic environment in numerical simulations of disc galaxies and galaxy mergers. Xco is proportional to the GMC surface density divided by the integrated CO intensity, Wco, and Wco is related to the kinetic temperature and velocity dispersion in the cloud. In disc galaxies (except within the central ~ kpc), the galactic environment is largely unimportant in setting the physical properties of GMCs provided they are gravitationally bound. The temperatures are roughly constant at ~10 K due to the balance of CO cooling and cosmic ray heating, giving a nearly constant CO-H2 conversion factor in discs. In mergers, the velocity dispersion of the gas rises dramatically during coalescence. The gas temperature also rises as it couples well to the warm (~50 K) dust at high densities (n > 10^4 cm^-3). The rise in velocity dispersion and temperature combine to offset the rise in surface density in mergers, causing Xco to drop by a factor of ~2-10 compared to the disc simulation. This model predicts that high-resolution ALMA observations of nearby ULIRGs should show velocity dispersions of ~10-100 km/s, and brightness temperatures comparable to the dust temperatures.

One model for the origin of typical galactic star clusters such as the Orion Nebula Cluster (ONC) is that they form via the rapid, efficient collapse of a bound gas clump within a larger, gravitationally-unbound giant molecular cloud. However, simulations in support of this scenario have thus far have not included the radiation feedback produced by the stars; radiative simulations have been limited to significantly smaller or lower density regions. Here we use the ORION adaptive mesh refinement code to conduct the first ever radiation-hydrodynamic simulations of the global collapse scenario for the formation of an ONC-like cluster. We show that radiative feedback has a dramatic effect on the evolution: once the first ~10-20% of the gas mass is incorporated into stars, their radiative feedback raises the gas temperature high enough to suppress any further fragmentation. However, gas continues to accrete onto existing stars, and, as a result, the stellar mass distribution becomes increasingly top-heavy, eventually rendering it incompatible with the observed IMF. Systematic variation in the location of the IMF peak as star formation proceeds is incompatible with the observed invariance of the IMF between star clusters, unless some unknown mechanism synchronizes the IMFs in different clusters by ensuring that star formation is always truncated when the IMF peak reaches a particular value. We therefore conclude that the global collapse scenario, at least in its simplest form, is not compatible with the observed stellar IMF. We speculate that processes that slow down star formation, and thus reduce the accretion luminosity, may be able to resolve the problem.

We revisit the problem of low-mass pre-main-sequence (PMS) stellar evolution and its observational consequences for where stars fall on the Hertzsprung-Russell diagram (HRD). In contrast to most previous work, our models follow stars as they grow from small masses via accretion, and we perform a systematic study of how the stars' HRD evolution is influenced by their initial radius, by the radiative properties of the accretion flow, and by the accretion history, using both simple idealized accretion histories and histories taken from numerical simulations of star cluster formation. We compare our numerical results to both non-accreting isochrones and to the positions of observed stars in the HRD, with a goal of determining whether both the absolute ages and the age dispersions inferred from non-accreting isochrones are reliable. We show that non-accreting isochrones can sometimes overestimate stellar ages for more massive stars (those with effective temperatures above \sim 3500 K), thereby explaining why non-accreting isochrones often suggest a systematic age difference between more and less massive stars in the same cluster. However, we also find the only way to produce a similar overestimate for the ages of cooler stars is by selecting parameters that are strongly inconsistent with both theoretical models and observational constraints. We therefore conclude that inferred ages and age spreads for cool stars are reliable, at least to the extent that the observed bolometric luminosities and temperatures are accurate. Finally, we note that the time-dependence of the mass accretion rate has remarkably little effect on low-mass stars' evolution on the HRD, and that such time-dependence may be neglected for all stars except those with effective temperatures above \sim 4000 K.

The effects of stochasticity on the luminosities of stellar populations are an often neglected but crucial element for understanding populations in the low mass or low star formation rate regime. To address this issue, we present SLUG, a new code to "Stochastically Light Up Galaxies". SLUG synthesizes stellar populations using a Monte Carlo technique that treats stochastic sampling properly including the effects of clustering, the stellar initial mass function, star formation history, stellar evolution, and cluster disruption. This code produces many useful outputs, including i) catalogs of star clusters and their properties, such as their stellar initial mass distributions and their photometric properties in a variety of filters, ii) two dimensional histograms of color-magnitude diagrams of every star in the simulation, iii) and the photometric properties of field stars and the integrated photometry of the entire simulated galaxy. After presenting the SLUG algorithm in detail, we validate the code through comparisons with starburst99 in the well-sampled regime, and with observed photometry of Milky Way clusters. Finally, we demonstrate the SLUG's capabilities by presenting outputs in the stochastic regime. SLUG is publicly distributed through the website http://sites.google.com/site/runslug/.

We present virial models for the global evolution of giant molecular clouds. Focusing on the presence of an accretion flow, and accounting for the amount of mass, momentum, and energy supplied by accretion and star formation feedback, we are able to follow the growth, evolution, and dispersal of individual giant molecular clouds. Our model clouds reproduce the scaling relations observed in both galactic and extragalactic clouds. We find that accretion and star formation contribute contribute roughly equal amounts of turbulent kinetic energy over the lifetime of the cloud. Clouds attain virial equilibrium and grow in such a way as to maintain roughly constant surface densities, with typical surface densities of order 50 - 200 Msun pc^-2, in good agreement with observations of giant molecular clouds in the Milky Way and nearby external galaxies. We find that as clouds grow, their velocity dispersion and radius must also increase, implying that the linewidth-size relation constitutes an age sequence. Lastly, we compare our models to observations of giant molecular clouds and associated young star clusters in the LMC and find good agreement between our model clouds and the observed relationship between H ii regions, young star clusters, and giant molecular clouds.

Young protostars embedded in circumstellar discs accrete from an angular momentum-rich mass reservoir. Without some braking mechanism, all stars should be spinning at or near break-up velocity. In this paper, we perform simulations of the self-gravitational collapse of an isothermal cloud using the ORION adaptive mesh refinement code and investigate the role that gravitational torques might play in the spin-down of the dense central object. While magnetic effects likely dominate for low mass stars, high mass and Population III stars might be less well magnetised. We find that gravitational torques alone prevent the central object from spinning up to more than half of its breakup velocity, because higher rotation rates lead to bar-like deformations that enable efficient angular momentum transfer to the surrounding medium. We also find that the long-term spin evolution of the central object is dictated by the properties of the surrounding disc. In particular, spiral modes with azimuthal wavenumber $m=2$ couple more effectively to its spin than the lopsided $m=1$ mode, which was found to inhibit spin evolution. We suggest that even in the absence of magnetic fields, gravitational torques may provide an upper limit on stellar spin, and that moderately massive circumstellar discs can cause long-term spin down.

The stellar initial mass function (IMF), along with the star formation rate, is one of the fundamental properties that any theory of star formation must explain. An interesting feature of the IMF is that it appears to be remarkably universal across a wide range of environments. Particularly, there appears to be little variation in either the characteristic mass of the IMF or its high-mass tail between clusters with different metallicities. Previous attempts to understand this apparent independence of metallicity have not accounted for radiation feedback from high-mass protostars, which can dominate the energy balance of the gas in star-forming regions. We extend this work, showing that the fragmentation of molecular gas should depend only weakly on the amount of dust present, even when the primary heating source is radiation from massive protostars. First, we report a series of core collapse simulations using the ORION AMR code that systematically vary the dust opacity and show explicitly that this has little effect on the temperature or fragmentation of the gas. Then, we provide an analytic argument for why the IMF varies so little in observed star clusters, even as the metallicity varies by a factor of 100.

We study the relation between molecular gas and star formation in a volume-limited sample of 222 galaxies from the COLD GASS survey, with measurements of the CO(1-0) line from the IRAM 30m telescope. The galaxies are at redshifts 0.025<z<0.05 and have stellar masses in the range 10.0<log(M*/Msun)<11.5. The IRAM measurements are complemented by deep Arecibo HI observations and homogeneous SDSS and GALEX photometry. A reference sample that includes both UV and far-IR data is used to calibrate our estimates of star formation rates from the seven optical/UV bands. The mean molecular gas depletion timescale, tdep(H2), for all the galaxies in our sample is 1 Gyr, however tdep(H2) increases by a factor of 6 from a value of ~0.5 Gyr for galaxies with stellar masses of 10^10 Msun to ~3 Gyr for galaxies with masses of a few times 10^11 Msun. In contrast, the atomic gas depletion timescale remains contant at a value of around 3 Gyr. This implies that in high mass galaxies, molecular and atomic gas depletion timescales are comparable, but in low mass galaxies, molecular gas is being consumed much more quickly than atomic gas. The strongest dependences of tdep(H2) are on the stellar mass of the galaxy (parameterized as log tdep(H2)= (0.36+/-0.07)(log M* - 10.70)+(9.03+/-0.99)), and on the specific star formation rate. A single tdep(H2) versus sSFR relation is able to fit both "normal" star-forming galaxies in our COLD GASS sample, as well as more extreme starburst galaxies (LIRGs and ULIRGs), which have tdep(H2) < 10^8 yr. Normal galaxies at z=1-2 are displaced with respect to the local galaxy population in the tdep(H2) versus sSFR plane and have molecular gas depletion times that are a factor of 3-5 times longer at a given value of sSFR due to their significantly larger gas fractions.

We are conducting COLD GASS, a legacy survey for molecular gas in nearby galaxies. Using the IRAM 30m telescope, we measure the CO(1-0) line in a sample of ~350 nearby (D=100-200 Mpc), massive galaxies (log(M*/Msun)>10.0). The sample is selected purely according to stellar mass, and therefore provides an unbiased view of molecular gas in these systems. By combining the IRAM data with SDSS photometry and spectroscopy, GALEX imaging and high-quality Arecibo HI data, we investigate the partition of condensed baryons between stars, atomic gas and molecular gas in 0.1-10L* galaxies. In this paper, we present CO luminosities and molecular hydrogen masses for the first 222 galaxies. The overall CO detection rate is 54%, but our survey also uncovers the existence of sharp thresholds in galaxy structural parameters such as stellar mass surface density and concentration index, below which all galaxies have a measurable cold gas component but above which the detection rate of the CO line drops suddenly. The mean molecular gas fraction MH2/M* of the CO detections is 0.066+/-0.039, and this fraction does not depend on stellar mass, but is a strong function of NUV-r colour. Through stacking, we set a firm upper limit of MH2/M*=0.0016+/-0.0005 for red galaxies with NUV-r>5.0. The average molecular-to-atomic hydrogen ratio in present-day galaxies is 0.3, with significant scatter from one galaxy to the next. The existence of strong detection thresholds in both the HI and CO lines suggests that "quenching" processes have occurred in these systems. Intriguingly, atomic gas strongly dominates in the minority of galaxies with significant cold gas that lie above these thresholds. This suggests that some re-accretion of gas may still be possible following the quenching event.

Observations show that star formation is an inefficient and slow process. This result can be attributed to the injection of energy and momentum by stars that prevents free-fall collapse of molecular clouds. The mechanism of this stellar feedback is debated theoretically: possible sources of pressure include the classical warm HII gas, the hot gas generated by shock-heating from stellar winds and supernovae, direct radiation of stars, and the dust-processed radiation field trapped inside the HII shell. In this paper, we measure observationally the pressures associated with each component listed above across the giant HII region 30 Doradus in the Large Magellanic Cloud. We exploit high-resolution, multi-wavelengh images (radio, infrared, optical, and X-ray) to map these pressures as a function of position. We find that radiation pressure dominates within 75 pc of the central star cluster, R136, while the HII gas pressure dominates at larger radii. By contrast, the dust-processed radiation pressure and hot gas pressure are generally weak and not dynamically important, although the hot gas pressure may have played a more significant role at early times. Based on the low X-ray gas pressures, we demonstrate that the hot gas is only partially confined and must be leaking out the HII shell. Additionally, we consider the implications of a dominant radiation pressure on the early dynamics of 30 Doradus.

Stars and more particularly massive stars, have a drastic impact on galaxy evolution. Yet the conditions in which they form and collapse are still not fully understood. In particular, the influence of the magnetic field on the collapse of massive clumps is relatively unexplored, it is thus of great relevance in the context of the formation of massive stars to investigate its impact. We perform high resolution, MHD simulations of the collapse of hundred solar masses, turbulent and magnetized clouds, using the adaptive mesh refinement code RAMSES. We compute various quantities such as mass distribution, magnetic field and angular momentum within the collapsing core and study the episodic outflows and the fragmentation that occurs during the collapse. The magnetic field has a drastic impact on the cloud evolution. We find that magnetic braking is able to substantially reduce the angular momentum in the inner part of the collapsing cloud. Fast and episodic outflows are being launched with typical velocities of the order of 3-5 km s$^{-1}$ although the highest velocities can be as high as 30-40 km s$^{-1}$. The fragmentation in several objects, is reduced in substantially magnetized clouds with respect to hydrodynamical ones by a factor of the order of 1.5-2. We conclude that magnetic fields have a significant impact on the evolution of massive clumps. In combination with radiation, magnetic fields largely determine the outcome of massive core collapse. We stress that numerical convergence of MHD collapse is a challenging issue. In particular, numerical diffusion appears to be important at high density therefore possibly leading to an over-estimation of the number of fragments.

Nearby spiral galaxies show an extremely tight correlation between tracers of molecular hydrogen (H_2) in the interstellar medium (ISM) and tracers of recent star formation, but it is unclear whether this correlation is fundamental or accidental. In the galaxies that have been surveyed to date, H_2 resides predominantly in gravitationally bound clouds cooled by carbon monoxide (CO) molecules, but in galaxies of low metal content the correlations between bound clouds, CO, and H_2 break down, and it is unclear if the star formation rate will then correlate with H_2 or with some other quantity. Here we show that star formation will continue to follow H_2 independent of metallicity. This is not because H_2 is directly important for cooling, but instead because the transition from predominantly atomic hydrogen (HI) to H_2 occurs under the same conditions as a dramatic drop in gas temperature and Bonnor-Ebert mass that destabilizes clouds and initiates collapse. We use this model to compute how star formation rate will correlate with total gas mass, with mass of gas where the hydrogen is H_2, and with mass of gas where the carbon is CO in galaxies of varying metallicity, and show that preliminary observations match the trend we predict.

We perform analytic linear stability analyses of an interface separating two stratified media threaded by a radiation flux, a configuration relevant in several astrophysical contexts. We develop a general framework for analyzing such systems, and obtain exact stability conditions in several limiting cases. In the optically thin, isothermal regime, where the discontinuity is chemical in nature (e.g.\ at the boundary of a radiation pressure-driven H \textsc{ii} region), radiation acts as part of an effective gravitational field, and instability arises if the effective gravity per unit volume toward the interface overcomes that away from it. In the optically thick "adiabatic" regime where the total (gas plus radiation) specific entropy of a Lagrangian fluid element is conserved,for example at the edge of radiation pressure-driven bubble around a young massive star, we show that radiation acts like a modified equation of state, and we derive a generalized version of the classical Rayleigh-Taylor stability condition.

Star formation is one of the least understood processes in cosmic evolution. It is difficult to formulate a general theory for star formation in part because of the wide range of physical processes involved. The interstellar gas out of which stars form is a supersonically turbulent plasma governed by magnetohydrodynamics. This is hard enough by itself, since we do not understand even subsonic hydrodynamic turbulence very well, let alone supersonic non-ideal MHD turbulence. However, the behavior of star-forming clouds in the ISM is also obviously influenced by gravity, which adds complexity, and by both continuum and line radiative processes. Finally, the behavior of star-forming clouds is influenced by a wide variety of chemical processes, including formation and destruction of molecules and dust grains (which changes the thermodynamic behavior of the gas) and changes in ionization state (which alter how strongly the gas couples to magnetic fields). As a result of these complexities, there is nothing like a generally agreed-upon theory of star formation, as there is for stellar structure. Instead, we are forced to take a much more phenomenological approach. These notes provide an introduction to our current thinking about how star formation works.

Recent observations indicate that star formation occurs only in the molecular phase of a galaxy's interstellar medium. A realistic treatment of star formation in simulations and analytic models of galaxies therefore requires that one determine where the transition from the atomic to molecular gas occurs. In this paper we compare two methods for making this determination in cosmological simulations where the internal structures of molecular clouds are unresolved: a complex time-dependent chemistry network coupled to a radiative transfer calculation of the dissociating ultraviolet (UV) radiation field, and a simple time-independent analytic approximation. We show that these two methods produce excellent agreement at all metallicities >~10^-2 of the Milky Way value across a very wide range of UV fields. At lower metallicities the agreement is worse, likely because time-dependent effects become important; however, there are no observational calibrations of molecular gas content at such low metallicities, so it is unclear if either method is accurate. The comparison suggests that, in many but not all applications, the analytic approximation provides a viable and nearly cost-free alternative to full time-dependent chemistry and radiative transfer.

We characterize the infall rate onto protostellar systems forming in self-gravitating radiation-hydrodynamic simulations. Using two dimensionless parameters to determine disks' susceptability to gravitational fragmentation, we infer limits on protostellar system multiplicity and the mechanism of binary formation. We show that these parameters give robust predictions even in the case of marginally resolved protostellar disks. We find that protostellar systems with radiation feedback predominately form binaries via turbulent fragmentation, not disk instability, and we predict turbulent fragmentation is the dominant channel for binary formation for low-mass stars. We clearly demonstrate that systems forming in simulations including radiative feedback have fundamentally different parameters than those in purely hydrodynamic simulations.

We derive the evolution equations describing a thin axisymmetric disk of gas and stars with an arbitrary rotation curve that is kept in a state of marginal gravitational instability and energy equilibrium due to the balance between energy released by accretion and energy lost due to decay of turbulence. Rather than adopt a parameterized alpha prescription, we instead use the condition of marginal gravitational instability to self-consistently determine the position- and time-dependent transport rates. We show that there is a steady-state configuration for disks dominated by gravitational instability, and that this steady state persists even when star formation is taken into account if the accretion rate is sufficiently large. For disks in this state we analytically determine the velocity dispersion, surface density, and rates of mass and angular momentum transport as a function of the gas mass fraction, the rotation curve, and the rate of external accretion onto the disk edge. We show that disks that are initially out of steady state will evolve into it on the viscous timescale of the disk, which is comparable to the orbital period if the accretion rate is high. Finally, we discuss the implications of these results for the structure of disks in a broad range of environments, including high redshift galaxies, the outer gaseous disks of local galaxies, and accretion disks around protostars.

We present SLUG, a new code to "Stochastically Light Up Galaxies". SLUG populates star clusters by randomly drawing stars from an initial mass function (IMF) and then following their time evolution with stellar models and an observationally-motivated prescription for cluster disruption. For a choice of star formation history, metallicity, and IMF, SLUG outputs synthetic photometry for clusters and field stars with a proper treatment of stochastic star formation. SLUG generates realistic distributions of star clusters, demonstrating the range of properties that result from finite sampling of an IMF and a random distribution of ages. The simulated data sets provide a quantitative means to address open problems in studies of star formation in galaxies and clusters, such as a test for IMF variations that are suggested by the systematic deficiency in the H-alpha/UV ratio in outer disks or in dwarf galaxies. SLUG will be made publicly available through the website http://sites.google.com/site/runslug/.

I provide a pedagogic review of adaptive mesh refinement (AMR) radiation hydrodynamics (RHD) methods and codes used in simulations of star formation, at a level suitable for researchers who are not computational experts. I begin with a brief overview of the types of RHD processes that are most important to star formation, and then I formally introduce the equations of RHD and the approximations one uses to render them computationally tractable. I discuss strategies for solving these approximate equations on adaptive grids, with particular emphasis on identifying the main advantages and disadvantages of various approximations and numerical approaches. Finally, I conclude by discussing areas ripe for improvement.

Stars in galaxies form in giant molecular clouds that coalesce when the atomic hydrogen is converted into molecules. There are currently two dominant models for what property of the galactic disk determines its molecular fraction: either hydrostatic pressure driven by the gravity of gas and stars, or a combination of gas column density and metallicity. To assess the validity of these models, we compare theoretical predictions to the observed atomic gas content of low-metallicity dwarf galaxies with high stellar densities. The extreme conditions found in these systems are optimal to distinguish the two models, otherwise degenerate in nearby spirals. Locally, on scales <100 pc, we find that the state of the interstellar medium is mostly sensitive to the gas column density and metallicity rather than hydrostatic pressure. On larger scales where the average stellar density is considerably lower, both pressure and shielding models reproduce the observations, even at low metallicity. We conclude that models based on gas and dust shielding more closely describe the process of molecular formation, especially at the high resolution that can be achieved in modern galaxy simulations or with future radio/millimeter arrays.

The chemical abundance patterns of the oldest stars in the Galaxy are expected to contain residual signatures of the first stars in the early universe. Numerous studies attempt to explain the intrinsic abundance scatter observed in some metal-poor populations in terms of chemical inhomogeneities dispersed throughout the early Galactic medium due to discrete enrichment events. Just how the complex data and models are to be interpreted with respect to "progenitor yields" remains an open question. Here we show that stochastic chemical evolution models to date have overlooked a crucial fact. Essentially all stars today are born in highly homogeneous star clusters and it is likely that this was also true at early times. When this ingredient is included, the overall scatter in the abundance plane [Fe/H] vs. [X/Fe] (C-space), where X is a nucleosynthetic element, can be much less than derived from earlier models. Moreover, for moderately flat cluster mass functions (gamma < 2), and/or for mass functions with a high mass cut-off (M_max > 10^5 M_sun), stars exhibit a high degree of clumping in C-space that can be identified even in relatively small data samples. Since stellar abundances can be modified by mass transfer in close binaries, clustered signatures are essential for deriving the yields of the first supernovae. We present a statistical test to determine whether a given set of observations exhibit such behaviour. Our initial work focusses on two dimensions in C-space, but we show that the clustering signal can be greatly enhanced by additional abundance axes. The proposed experiment will be challenging on existing 8-10m telescopes, but relatively straightforward for a multi-object echelle spectrograph mounted on a 25-40m telescope.

We extend existing semi-analytic models of galaxy formation to track atomic and molecular gas in disk galaxies. Simple recipes for processes such as cooling, star formation, supernova feedback, and chemical enrichment of the stars and gas are grafted on to dark matter halo merger trees derived from the Millennium Simulation. Each galactic disk is represented by a series of concentric rings. We assume that surface density profile of infalling gas in a dark matter halo is exponential, with scale radius r_d that is proportional to the virial radius of the halo times its spin parameter $\lambda$. As the dark matter haloes grow through mergers and accretion, disk galaxies assemble from the inside out. We include two simple prescriptions for molecular gas formation processes in our models: one is based on the analytic calculations by Krumholz, McKee & Tumlinson (2008), and the other is a prescription where the H_2 fraction is determined by the kinematic pressure of the ISM. Motivated by the observational results of Leroy et al. (2008), we adopt a star formation law in which $\Sigma_{SFR}\propto\Sigma_{H_2}$ in the regime where the molecular gas dominates the total gas surface density, and $\Sigma_{SFR}\propto \Sigma_{gas}^2$ where atomic hydrogen dominates. We then fit these models to the radial surface density profiles of stars, HI and H_2 drawn from recent high resolution surveys of stars and gas in nearby galaxies. We explore how the ratios of atomic gas, molecular gas and stellar mass vary as a function of global galaxy scale parameters, including stellar mass, stellar surface density, and gas surface density. We elucidate how the trends can be understood in terms of three variables that determine the partition of baryons in disks: the mass of the dark matter halo, the spin parameter of the halo, and the amount of gas recently accreted from the external environment.

We investigate the effects of radiation pressure from stars on the survival of the star-forming giant clumps in high-redshift massive disc galaxies, during the most active phase of galaxy formation. The clumps, typically of mass ~10^8-10^9 Msun and radius ~0.5-1, are formed in the turbulent gas-rich discs by violent gravitational instability and then migrate into a central bulge in ~10 dynamical times. We show that the survival or disruption of these clumps under the influence of stellar feedback depends critically on the rate at which they form stars. If they convert a few percent of their gas mass to stars per free-fall time, as observed for all local star-forming systems and implied by the Kennicutt-Schmidt law, they cannot be disrupted. Only if clumps convert most of their mass to stars in a few free-fall times can feedback produce significant gas expulsion. We consider whether such rapid star formation is likely in high-redshift giant clumps.

The fragmentation of star-forming interstellar clouds, and the resulting stellar initial mass function (IMF), is strongly affected by the temperature structure of the collapsing gas. Since radiation feedback from embedded stars can modify this as collapse proceeds, feedback plays an important role in determining the IMF. However, the effects and importance of radiative heating are likely to depend strongly on the surface density of the collapsing clouds, which determines both their effectiveness at trapping radiation and the accretion luminosities of the stars forming within them. In this paper we report a suite of adaptive mesh refinement radiation-hydrodynamic simulations using the ORION code in which we isolate the effect of column density on fragmentation by following the collapse of clouds of varying column density while holding the mass, initial density and velocity structure, and initial virial ratio fixed. We find that radiation does not significantly modify the overall star formation rate or efficiency, but that it suppresses fragmentation more and more as cloud surface densities increase from those typical of low mass star-forming regions like Taurus, through the typical surface density of massive star-forming clouds in the Galaxy, up to conditions found only in super star clusters. In regions of low surface density, fragmentation during collapse leads to the formation of small clusters rather than individual massive star systems, greatly reducing the fraction of the stellar population with masses >~ 10 Msun. Our simulations have important implications for the formation of massive stars and the universality of the IMF.

The Galactic disk retains a vast amount of information about how it came to be, and how it evolved over cosmic time. However, we know very little about the secular processes associated with disk evolution. One major uncertainty is the extent to which stars migrate radially through the disk, thereby washing out signatures of their past (e.g. birth sites). Recent theoretical work finds that such "blurring" of the disk can be important if spiral arms are transient phenomena. Here we describe an experiment to determine the importance of diffusion from the Solar circle with cosmic time. Consider a star cluster that has been placed into a differentially rotating, stellar fluid. We show that all clusters up to ~10^4 solar masses, and a significant fraction of those up to ~10^5 solar masses, are expected to be chemically homogeneous, and that clusters of this size can be assigned a unique "chemical tag" by measuring the abundances of <~10 independent element groups, with better age and orbit determinations allowing fewer abundance measurements. The star cluster therefore acts like a "tracer dye", and the present-day distribution of its stars provides a strong constraint on the rate of radial diffusion or migration in the Galactic disk. Sellwood & Binney have argued for strong radial transport driven by transient spiral perturbations: in principle, we could measure the strength of this migration directly.

We investigate how the removal of interstellar material by stellar feedback limits the efficiency of star formation in molecular clouds and how this determines the shape of the mass function of young star clusters. In particular, we derive relations between the power-law exponents of the mass functions of the clouds and clusters in the limiting regimes in which the feedback is energy-driven and momentum-driven, corresponding to minimum and maximum radiative losses and likely to bracket all realistic cases. We find good agreement between the predicted and observed exponents, especially for momentum-driven feedback, provided the protoclusters have roughly constant mean surface density, as indicated by observations of the star-forming clumps within molecular clouds. We also consider a variety of specific feedback mechanisms, concluding that HII regions inflated by radiation pressure predominate in massive protoclusters, a momentum-limited process when photons can escape after only a few interactions with dust grains. We then present a first estimate of the star formation efficiency in this case, finding that it depends on the masses and sizes of the protoclusters only through their mean surface density, thus ensuring consistency between the observed power-law exponents of the mass functions of the clouds and clusters. The numerical value of this efficiency is also consistent with observations.

Understanding the molecular content of galaxies is a critical problem in star formation and galactic evolution. Here we present a new method, based on a Stromgren-type analysis, to calculate the amount of HI that surrounds a molecular cloud irradiated by an isotropic radiation field. We consider both planar and spherical clouds, and H_2 formation either in the gas phase or catalyzed by dust grains. Under the assumption that the transition from atomic to molecular gas is sharp, our method gives the solution without any reference to the photodissociation cross section. We test our results for the planar case against those of a PDR code, and find typical accuracies of about 10%. Our results are also consistent with the scaling relations found in Paper I of this series, but they apply to a wider range of physical conditions. We present simple, accurate analytic fits to our results that are suitable for comparison to observations and to implementation in numerical and semi-analytic models.

The kinematics of newly-formed star clusters are interesting both as a probe of the state of the gas clouds from which the stars form, and because they influence planet formation, stellar mass segregation, cluster disruption, and other processes controlled in part by dynamical interactions in young clusters. However, to date there have been no attempts to use simulations of star cluster formation to investigate how the kinematics of young stars change in response to variations in the properties of their parent molecular clouds. In this letter we report the results of turbulent self-gravitating simulations of cluster formation in which we consider both clouds in virial balance and those undergoing global collapse. We find that stars in these simulations generally have velocity dispersions smaller than that of the gas by a factor of ~ 5, independent of the dynamical state of the parent cloud, so that subvirial stellar velocity dispersions arise naturally even in virialized molecular clouds. The simulated clusters also show large-scale stellar velocity gradients of ~0.2-2 km s$^{-1}$ pc$^{-1}$ and strong correlations between the centroid velocities of stars and gas, both of which are observed in young clusters. We conclude that star clusters should display subvirial velocity dispersions, large-scale velocity gradients, and strong gas-star velocity correlations regardless of whether their parent clouds are in virial balance, and, conversely, that observations of these features cannot be used to infer the dynamical state of the parent gas clouds.

We evaluate the role of radiation pressure in the dynamics of HII regions. We first determine under what conditions radiation pressure is significant in comparison to gas pressure and show that, while radiation pressure is generally unimportant for HII regions driven by a handful of massive stars, it is dominant for the larger HII regions produced by the massive star clusters found near the Galactic center and in starburst environments. We then provide a solution for the problem of how HII regions expand when radiation pressure influences their behavior. Finally, we compare radiation-dominated HII regions to other sources of stellar feedback, and argue that HII regions are probably the primary mechanism for regulating the formation of massive star clusters.

Forming stars emit a substantial amount of radiation into their natal environment. We use ORION, an adaptive mesh refinement (AMR) three-dimensional gravito-radiation-hydrodynamics code, to simulate low-mass star formation in a turbulent molecular cloud. We compare the distribution of stellar masses, accretion rates, and temperatures in the cases with and without radiative transfer, and we demonstrate that radiative feedback has a profound effect on accretion, multiplicity, and mass by reducing the number of stars formed and the total rate at which gas turns into stars. We also show, that once star formation reaches a steady state, protostellar radiation is by far the dominant source of energy in the simulation, exceeding viscous dissipation and compressional heating by at least an order of magnitude. Calculations that omit radiative feedback from protstars significantly underestimate the gas temperature and the strength of this effect. Although heating from protostars is mainly confined to the protostellar cores, we find that it is sufficient to suppress disk fragmentation that would otherwise result in very low-mass companions or brown dwarfs. We demonstrate that the mean protostellar accretion rate increases with the final stellar mass so that the star formation time is only a weak function of mass.

We study rapidly accreting, gravitationally unstable disks with a series of global, three dimensional, numerical experiments using the code ORION. In this paper we conduct a numerical parameter study focused on protostellar disks, and show that one can predict disk behavior and the multiplicity of the accreting star system as a function of two dimensionless parameters which compare the disk's accretion rate to its sound speed and orbital period. Although gravitational instabilities become strong, we find that fragmentation into binary or multiple systems occurs only when material falls in several times more rapidly than the canonical isothermal limit. The disk-to-star accretion rate is proportional to the infall rate, and governed by gravitational torques generated by low-m spiral modes. We also confirm the existence of a maximum stable disk mass: disks that exceed ~50% of the total system mass are subject to fragmentation and the subsequent formation of binary companions.

We argue that the lack of observed damped Lyman alpha (DLA) systems that simultaneously have high HI columns densities and high metallicities results naturally from the formation of molecules in the cold phase of a two-phase atomic medium in pressure balance. Our result applies equally well in diffuse systems where the ultraviolet radiation field is dominated by the extragalactic background and in dense star-forming ones, where the local radiation field is likely to be orders of magnitude higher. We point out that such a radiation-insensitive model is required to explain the absence of high column - high metallicity systems among DLAs observed using gamma-ray burst afterglows, since these are likely subjected to strong radiation fields created by active star formation in the GRB host galaxy. Moreover, we show that the observed relationship between the maximum atomic gas column in DLAs sets a firm upper limit on the fraction of the mass in these systems that can be in the warm, diffuse phase. Finally, we argue that our result explains the observed lack of in situ star formation in DLA systems.

The formation of stars is a key process in astrophysics. Detailed knowledge of the physical mechanisms that govern stellar birth is a prerequisite for understanding the formation and evolution of our galactic home, the Milky Way. A theory of star formation is an essential part of any model for the origin of our solar system and of planets around other stars. Despite this pivotal importance, and despite many decades of research, our understanding of the processes that initiate and regulate star formation is still limited. Stars are born in cold interstellar clouds of molecular hydrogen gas. Star formation in these clouds is governed by the complex interplay between the gravitational attraction in the gas and agents such as turbulence, magnetic fields, radiation and thermal pressure that resist compression. The competition between these processes determines both the locations at which young stars form and how much mass they ultimately accrete. It plays out over many orders of magnitude in space and time, ranging from galactic to stellar scales. In addition, star formation is a highly stochastic process in which rare and hard-to-predict events, such as the formation of very massive stars and the resulting feedback, can play a dominant role in determining the evolution of a star-forming cloud. As a consequence of the wide range of scales and processes that control star formation, analytic models are usually restricted to highly idealized cases. These can yield insight, but the complexity of the problem means that they must be used in concert with large-scale numerical simulations. Here we summarize the state of modern star formation theory and review the recent advances in numerical simulation techniques.

We propose a simple theoretical model for star formation in which the local star formation rate in a galaxy is determined by three factors. First, the interplay between the interstellar radiation field and molecular self-shielding determines what fraction of the gas is in molecular form and thus eligible to form stars. Second, internal feedback determines the properties of the molecular clouds that form, which are nearly independent of galaxy properties until the galactic ISM pressure becomes comparable to the internal GMC pressure. Above this limit, galactic ISM pressure determines molecular gas properties. Third, the turbulence driven by feedback processes in GMCs makes star formation slow, allowing a small fraction of the gas to be converted to stars per free-fall time within the molecular clouds. We combine analytic estimates for each of these steps to formulate a single star formation law, and show that the predicted correlation between star formation rate, metallicity, and surface densities of atomic, molecular, and total gas agree well with observations.

We use a sample of 47 homogeneous and high sensitivity CO images taken from the Nobeyama and BIMA surveys to demonstrate that, contrary to common belief, a significant number (~40%) of HI-deficient nearby spiral galaxies are also depleted in molecular hydrogen. While HI-deficiency by itself is not a sufficient condition for molecular gas depletion, we find that H2 reduction is associated with the removal of HI inside the galaxy optical disk. Those HI-deficient galaxies with normal H2 content have lost HI mainly from outside their optical disks, where the H2 content is low in all galaxies. This finding is consistent with theoretical models in which the molecular fraction in a galaxy is determined primarily by its gas column density. Our result is supported by indirect evidence that molecular deficient galaxies form stars at a lower rate or have dimmer far infrared fluxes than gas rich galaxies, as expected if the star formation rate is determined by the molecular hydrogen content. Our result is consistent with a scenario in which, when the atomic gas column density is lowered inside the optical disk below the critical value required to form molecular hydrogen and stars, spirals become quiescent and passive evolving systems. We speculate that this process would act on the time-scale set by the gas depletion rate and might be a first step for the transition between the blue and red sequence observed in the color-magnitude diagram.

Great strides have been made in the last two decades in determining how galaxies evolve from their initial dark matter seeds to the complex structures we observe at z=0. The role of mergers has been documented through both observations and simulations, numerous satellites that may represent these initial dark matter seeds have been discovered in the Local Group, high redshift galaxies have been revealed with monstrous star formation rates, and the gaseous cosmic web has been mapped through absorption line experiments. Despite these efforts, the dark matter simulations that include baryons are still unable to accurately reproduce galaxies. One of the major problems is our incomplete understanding of how a galaxy accretes its baryons and subsequently forms stars. Galaxy formation simulations have been unable to accurately represent the required gas physics on cosmological timescales, and observations have only just begun to detect the star formation fuel over a range of redshifts and environments. How galaxies obtain gas and subsequently form stars is a major unsolved, yet tractable problem in contemporary extragalactic astrophysics. In this paper we outline how progress can be made in this area in the next decade.

Massive stars produce so much light that the radiation pressure they exert on the gas and dust around them is stronger than their gravitational attraction, a condition that has long been expected to prevent them from growing by accretion. We present three-dimensional radiation-hydrodynamic simulations of the collapse of a massive prestellar core and find that radiation pressure does not halt accretion. Instead, gravitational and Rayleigh-Taylor instabilities channel gas onto the star system through non-axisymmetric disks and filaments that self-shield against radiation, while allowing radiation to escape through optically-thin bubbles. Gravitational instabilities cause the disk to fragment and form a massive companion to the primary star. Radiation pressure does not limit stellar masses, but the instabilities that allow accretion to continue lead to small multiple systems.

In this study, we investigate the shapes of starless and protostellar cores using hydrodynamic, self-gravitating adaptive mesh refinement simulations of turbulent molecular clouds. We simulate observations of these cores in dust emission, including realistic noise and telescope resolution, and compare to the observed core shapes measured in Orion by Nutter & Ward-Thompson (2007). The simulations and the observations have generally high statistical similarity, with particularly good agreement between simulations and Orion B. Although protostellar cores tend to have semi-major axis to semi-minor axis ratios closer to one, the distribution of axis ratios for starless and protostellar cores are not significantly different for either the actual observations of Orion or the simulated observations. Because of the high level of agreement between the non-magnetic hydrodynamic simulations and observation, contrary to a number of previous authors, one cannot infer the presence of magnetic fields from core shape distributions.

Gas in galactic disks is collected by gravitational instabilities into giant atomic-molecular complexes, but only the inner, molecular parts of these structures are able to collapse to form stars. Determining what controls the ratio of atomic to molecular hydrogen in complexes is therefore a significant problem in star formation and galactic evolution. In this paper we use the model of H_2 formation, dissociation, and shielding developed in the previous paper in this series to make theoretical predictions for atomic to molecular ratios as a function of galactic properties. We find that the molecular fraction in a galaxy is determined primarily by its column density and secondarily by its metallicity, and is to good approximation independent of the strength of the interstellar radiation field. We show that the column of atomic hydrogen required to shield a molecular region against dissociation is ~10 Msun pc^-2 at solar metallicity. We compare our model to data from recent surveys of the Milky Way and of nearby galaxies, and show that the both the primary dependence of molecular fraction on column density and the secondary dependence on metallicity that we predict are in good agreement with observed galaxy properties.

In this series of papers we study the structure of the atomic to molecular transition in the giant atomic-molecular complexes that are the repositories of most molecular gas in galaxies, with the ultimate goal of attaining a better understanding of what determines galaxies' molecular content. Here we derive an approximate analytic solution for the structure of a photodissociation region (PDR) in a cloud of finite size that is bathed in an external dissociating radiation field. Our solution extends previous work, which with few exceptions has been restricted to a one-dimensional treatment of the radiation field. We show that our analytic results compare favorably to exact numerical calculations in the one-dimensional limit. However, our more general geometry provides a more realistic representation than a semi-infinite slab for atomic-molecular complexes exposed to the interstellar radiation field, particularly in environments such as low-metallicity dwarf galaxies where the curvature and finite size of the atomic envelope cannot be neglected. For clouds that are at least 20% molecular we obtain analytic expressions for the molecular fraction in terms of properties of the gas and radiation field that are accurate to tens of percent, while for clouds of lower molecular content we obtain upper limits. As a side benefit, our analysis helps clarify when self-shielding is the dominant process in H_2 formation, and under what circumstances shielding by dust makes a significant contribution.

In this study we investigate the formation and properties of prestellar and protostellar cores using hydrodynamic, self-gravitating Adaptive Mesh Refinement simulations, comparing the cases where turbulence is continually driven and where it is allowed to decay. We model observations of these cores in the C$^{18}$O$(2\to 1)$, NH$_3(1,1)$, and N$_2$H$^+(1\to 0)$ lines, and from the simulated observations we measure the linewidths of individual cores, the linewidths of the surrounding gas, and the motions of the cores relative to one another. Some of these distributions are significantly different in the driven and decaying runs, making them potential diagnostics for determining whether the turbulence in observed star-forming clouds is driven or decaying. Comparing our simulations with observed cores in the Perseus and $\rho$ Ophiuchus clouds shows reasonably good agreement between the observed and simulated core-to-core velocity dispersions for both the driven and decaying cases. However, we find that the linewidths through protostellar cores in both simulations are too large compared to the observations. The disagreement is noticably worse for the decaying simulation, in which cores show highly supersonic infall signatures in their centers that decrease toward their edges, a pattern not seen in the observed regions. This result gives some support to the use of driven turbulence for modeling regions of star formation, but reaching a firm conclusion on the relative merits of driven or decaying turbulence will require more complete data on a larger sample of clouds as well as simulations that include magnetic fields, outflows, and thermal feedback from the protostars.

Most analytic work to date on protostellar disks has focused on those in isolation from their environments. However, observations are now beginning to probe the earliest, most embedded phases of star formation, during which disks are rapidly accreting from their parent cores and cannot be modeled in isolation. We present a simple, one-zone model of protostellar accretion disks with high mass infall rates. Our model combines a self-consistent calculation of disk temperatures with an approximate treatment of angular momentum transport via two mechanisms. We use this model to survey the properties of protostellar disks across a wide range of stellar masses and evolutionary times, and make predictions for disks' masses, sizes, spiral structure, and fragmentation that will be directly testable by future large-scale surveys of deeply embedded disks. We define a dimensionless accretion-rotation parameter which, in conjunction with the disk's temperature, controls the disk evolution. We track the dominant mode of angular momentum transport, and demonstrate that for stars with final masses greater than roughly one solar mass, gravitational instabilities are the most important mechanism as most of the mass accumulates. We predict that binary formation through disk fission, fragmentation of the disk into small objects, and spiral arm strength all increase in importance to higher stellar masses.

Massive stars are very rare, but their extreme luminosities make them both the only type of young star we can observe in distant galaxies and the dominant energy sources in the universe today. They form rarely because efficient radiative cooling keeps most star-forming gas clouds close to isothermal as they collapse, and this favors fragmentation into stars <~1 Msun. Heating of a cloud by accreting low-mass stars within it can prevent fragmentation and allow formation of massive stars, but what properties a cloud must have to form massive stars, and thus where massive stars form in a galaxy, has not yet been determined. Here we show that only clouds with column densities >~ 1 g cm^-2 can avoid fragmentation and form massive stars. This threshold, and the environmental variation of the stellar initial mass function (IMF) that it implies, naturally explain the characteristic column densities of massive star clusters and the difference between the radial profiles of Halpha and UV emission in galactic disks. The existence of a threshold also implies that there should be detectable variations in the IMF with environment within the Galaxy and in the characteristic column densities of massive star clusters between galaxies, and that star formation rates in some galactic environments may have been systematically underestimated.

The formation of massive stars is currently an unsolved problems in astrophysics. Understanding the formation of massive stars is essential because they dominate the luminous, kinematic, and chemical output of stars. Furthermore, their feedback is likely to play a dominant role in the evolution of molecular clouds and any subsequent star formation therein. Although significant progress has been made observationally and theoretically, we still do not have a consensus as to how massive stars form. There are two contending models to explain the formation of massive stars, Core Accretion and Competitive Accretion. They differ primarily in how and when the mass that ultimately makes up the massive star is gathered. In the core accretion model, the mass is gathered in a prestellar stage due to the overlying pressure of a stellar cluster or a massive pre-cluster cloud clump. In contrast, competitive accretion envisions that the mass is gathered during the star formation process itself, being funneled to the centre of a stellar cluster by the gravitational potential of the stellar cluster. Although these differences may not appear overly significant, they involve significant differences in terms of the physical processes involved. Furthermore, the differences also have important implications in terms of the evolutionary phases of massive star formation, and ultimately that of stellar clusters and star formation on larger scales. Here we review the dominant models, and discuss prospects for developing a better understanding of massive star formation in the future.

Recent advances in our understanding of massive star formation have made clear the important role of protostellar disks in mediating accretion. Here we describe a simple, semi-analytic model for young, deeply embedded, massive accretion disks. Our approach enables us to sample a wide parameter space of stellar mass and environmental variables, providing a means to make predictions for a variety of sources that next generation telescopes like ALMA and the EVLA will observe. Moreover we include, at least approximately, multiple mechanisms for angular momentum transport, a comprehensive model for disk heating and cooling, and a realistic estimate for the angular momentum in the gas reservoir. We make predictions for the typical sizes, masses, and temperatures of the disks, and describe the role of gravitational instabilities in determining the binarity fraction and upper mass cut-off.

We present a method for simulating the evolution of HII regions driven by point sources of ionizing radiation in magnetohydrodynamic media, implemented in the three-dimensional Athena MHD code. We compare simulations using our algorithm to analytic solutions and show that the method passes rigorous tests of accuracy and convergence. The tests reveal several conditions that an ionizing radiation-hydrodynamic code must satisfy to reproduce analytic solutions. As a demonstration of our new method, we present the first three-dimensional, global simulation of an HII region expanding into a magnetized gas. The simulation shows that magnetic fields suppress sweeping up of gas perpendicular to magnetic field lines, leading to small density contrasts and extremely weak shocks at the leading edge of the HII region's expanding shell.

We provide a model for how Kennicutt-Schmidt (KS) laws, which describe the correlation between star formation rate and gas surface or volume density, depend on the molecular line chosen to trace the gas. We show that, for lines that can be excited at low temperatures, the KS law depends on how the line critical density compares to the median density in a galaxy's star-forming molecular clouds. High critical density lines trace regions with similar physical properties across galaxy types, and this produces a linear correlation between line luminosity and star formation rate. Low critical density lines probe regions whose properties vary across galaxies, leading to a star formation rate that varies superlinearly with line luminosity. We show that a simple model in which molecular clouds are treated as isothermal and homogenous can quantitatively reproduce the observed correlations between galactic luminosities in far infrared and in the CO(1->0) and HCN(1->0) lines, and naturally explains why these correlations have different slopes. We predict that IR-line luminosity correlations should change slope for galaxies in which the median density is close to the line critical density. This prediction may be tested by observations of lines such as HCO^+(1->0) with intermediate critical densities, or by HCN(1->0) observations of intensely star-forming high redshift galaxies with very high densities. Recent observations by Gao et al. hint at just such a change in slope. We argue that deviations from linearity in the HCN(1->0)-IR correlation at high luminosity are consistent with the assumption of a constant star formation efficiency.

The similarity between the mass and spatial distributions of pre-stellar gas cores in star-forming clouds and young stars in clusters provides strong circumstantial evidence that these gas cores are the direct progenitors of individual stars. Here I describe a physical model for the evolution of massive cores into stars, starting with the initial phases of collapse and fragmentation, through disk formation and fragmentation, the later phases of stellar feedback, and finally interaction of the newly formed stars with their environments. This model shows that a direct mapping from cores to stars is the natural physical outcome of massive core evolution, and thereby allows us to explain many of the properties of young star clusters as direct imprints of their gas-phase progenitors.

We analyze the mixed frame equations of radiation hydrodynamics under the approximations of flux-limited diffusion and a thermal radiation field, and derive the minimal set of evolution equations that includes all terms that are of leading order in any regime of non-relativistic radiation hydrodynamics. Our equations are accurate to first order in v/c in the static diffusion regime. In contrast, we show that previous lower order derivations of these equations omit leading terms in at least some regimes. In comparison to comoving frame formulations of radiation hydrodynamics, our equations have the advantage that they manifestly conserve total energy, making them very well-suited to numerical simulations, particularly with adaptive meshes. For systems in the static diffusion regime, our analysis also suggests an algorithm that is both simpler and faster than earlier comoving frame methods. We implement this algorithm in the Orion adaptive mesh refinement code, and show that it performs well in a range of test problems.

We compute the molecular line emission of massive protostellar disks by solving the equation of radiative transfer through the cores and disks produced by the recent radiation-hydrodynamic simulations of Krumholz, Klein, & McKee. We find that in several representative lines the disks show brightness temperatures of hundreds of Kelvin over velocity channels ~10 km s^-1 wide, extending over regions hundreds of AU in size. We process the computed intensities to model the performance of next-generation radio and submillimeter telescopes. Our calculations show that observations using facilities such as the EVLA and ALMA should be able to detect massive protostellar disks and measure their rotation curves, at least in the nearest massive star-forming regions. They should also detect significant sub-structure and non-axisymmetry in the disks, and in some cases may be able to detect star-disk velocity offsets of a few km s^-1, both of which are the result of strong gravitational instability in massive disks. We use our simulations to explore the strengths and weaknesses of different observational techniques, and we also discuss how observations of massive protostellar disks may be used to distinguish between alternative models of massive star formation.

Rapidly accreting massive protostars undergo a phase of deuterium shell burning during pre-main sequence evolution that causes them to swell to tenths of an AU in radius. During this phase, those with close binary companions will overflow their Roche lobes and begin transferring mass. Since massive stars frequently have companions at distances well under 1 AU, this process may affect the early evolution of a substantial fraction of massive stars. We use a simple protostellar evolution model to determine the range in accretion rates, mass ratios, and orbital separations for which mass transfer will occur, and we compute approximately the stability and final outcome of the transfer process. We discuss how mass transfer affects the demographics of massive binaries, and show that it provides a natural explanation for the heretofore unexplained population of massive "twins", high mass binaries with mass ratios very close to unity.

We simulate the early stages of the evolution of turbulent, virialized, high-mass protostellar cores, with primary attention to how cores fragment, and whether they form a small or large number of protostars. Our simulations use the Orion adaptive mesh refinement code to follow the collapse from ~0.1 pc scales to ~10 AU scales, for durations that cover the main fragmentation phase, using three-dimensional gravito-radiation hydrodynamics. We find that for a wide range of initial conditions radiation feedback from accreting protostars inhibits the formation of fragments, so that the vast majority of the collapsed mass accretes onto one or a few objects. Most of the fragmentation that does occur takes place in massive, self-shielding disks. These are driven to gravitational instability by rapid accretion, producing rapid mass and angular momentum transport that allows most of the gas to accrete onto the central star rather than forming fragments. In contrast, a control run using the same initial conditions but an isothermal equation of state produces much more fragmentation, both in and out of the disk. We conclude that massive cores with observed properties are not likely to fragment into many stars, so that, at least at high masses, the core mass function probably determines the stellar initial mass function. Our results also demonstrate that simulations of massive star forming regions that do not include radiative transfer, and instead rely on a barotropic equation of state or optically thin heating and cooling curves, are likely to produce misleading results.

One of the outstanding puzzles about star formation is why it proceeds so slowly. Giant molecular clouds convert only a few percent of their gas into stars per free-fall time, and recent observations show that this low star formation rate is essentially constant over a range of scales from individual cluster-forming molecular clumps in the Milky Way to entire starburst galaxies. This striking result is perhaps the most basic fact that any theory of star formation must explain. I argue that a model in which star formation occurs in virialized structures at a rate regulated by supersonic turbulence can explain this observation. The turbulence in turn is driven by star formation feedback, which injects energy to offset radiation from isothermal shocks and keeps star-forming structures from wandering too far from virial balance. This model is able to reproduce observational results covering a wide range of scales, from the formation times of young clusters to the extragalactic IR-HCN correlation, and makes additional quantitative predictions that will be testable in the next few years.

It has been known for more than 30 years that star formation in giant molecular clouds (GMCs) is slow, in the sense that only ~1% of the gas forms stars every free-fall time. This result is entirely independent of any particular model of molecular cloud lifetime or evolution. Here we survey observational data on higher density objects in the interstellar medium, including infrared dark clouds and dense molecular clumps, to determine if these objects form stars slowly like GMCs, or rapidly, converting a significant fraction of their mass into stars in one free-fall time. We find no evidence for a transition from slow to rapid star formation in structures covering three orders of magnitude in density. This has important implications for models of star formation, since competing models make differing predictions for the characteristic density at which star formation should transition from slow to rapid. The data are inconsistent with models that predict that star clusters form rapidly and in free-fall collapse. Magnetic- and turbulence-regulated star formation models can reproduce the observations qualitatively, and the turbulence-regulated star formation model of Krumholz & McKee quantitatively reproduces the infrared-HCN luminosity correlation recently reported by Gao & Solomon. Slow star formation also implies that the process of star cluster formation cannot be one of global collapse, but must instead proceed over many free-fall times. This suggests that turbulence in star-forming clumps must be driven, and that the competitive accretion mechanism does not operate in typical cluster-forming molecular clumps.

We present semi-analytic dynamical models for giant molecular clouds evolving under the influence of HII regions launched by newborn star clusters. In contrast to previous work, we neither assume that clouds are in virial or energetic equilibrium, nor do we ignore the effects of star formation feedback. The clouds, which we treat as spherical, can expand and contract homologously. Photoionization drives mass ejection; the recoil of cloud material both stirs turbulent motions and leads to an effective confining pressure. The balance between these effects and the decay of turbulent motions through isothermal shocks determines clouds' dynamical and energetic evolution. We find that for realistic values of the rates of turbulent dissipation, photoevaporation, and energy injection by HII regions, the massive clouds where most molecular gas in the Galaxy resides live for a few crossing times, in good agreement with recent observational estimates that large clouds in local group galaxies survive roughly 20-30 Myr. During this time clouds remain close to equilibrium, with virial parameters of 1-3 and column densities near 10^22 H atoms cm^-2, also in agreement with observed cloud properties. Over their lives they convert 5-10% of their mass into stars, after which point most clouds are destroyed when a large HII region unbinds them. In contrast, small clouds like those found in the solar neighborhood only survive ~1 crossing time before being destroyed.

The current generation of millimeter interferometers have revealed a population of compact (r <~ 0.1 pc), massive (M ~ 100 Msun) gas cores that are the likely progenitors of massive stars. I review models for the evolution of these objects from the observed massive core phase through collapse and into massive star formation, with particular attention to the least well-understood aspects of the problem: fragmentation during collapse, interactions of newborn stars with the gas outside their parent core, and the effects of radiation pressure feedback. Through a combination of observation, analytic argument, and numerical simulation, I develop a model for massive star formation by gravitational collapse in which massive cores collapse to produce single stars or (more commonly) small-multiple systems, and these stars do not gain significant mass from outside their parent core by accretion of either gas or other stars. Collapse is only very slightly inhibited by feedback from the massive star, thanks to beaming of the radiation by a combination of protostellar outflows and radiation-hydrodynamic instabilities. Based on these findings, I argue that many of the observed properties of young star clusters can be understood as direct translations of the properties of their gas phase progenitors. Finally, I discuss unsolved problems in the theory of massive star formation, and directions for future work on them.

We argue that rich star clusters take at least several local dynamical times to form, and so are quasi-equilibrium structures during their assembly. Observations supporting this conclusion include morphologies of star-forming clumps, momentum flux of protostellar outflows from forming clusters, age spreads of stars in the Orion Nebula Cluster (ONC) and other clusters, and the age of a dynamical ejection event from the ONC. We show that these long formation timescales are consistent with the expected star formation rate in turbulent gas, as recently evaluated by Krumholz & McKee. Finally, we discuss the implications of these timescales for star formation efficiencies, the disruption of gas by stellar feedback, mass segregation of stars, and the longevity of turbulence in molecular clumps.

Star formation generally proceeds inside-out, with overdense regions inside protostellar cores collapsing rapidly and progressively less dense regions following later. Consequently, a small protostar will form early in the evolution of a core, and collapsing material will fall to the protostellar surface and radiate away its gravitational potential energy. The resulting accretion luminosity will heat the core and may substantially affect the process of fragmentation. This is of particular interest for massive cores that, at their initial temperatures, have masses much greater than a thermal Jeans mass and thus might be expected to fragment into many stars during collapse. Here I show that accretion luminosity can heat the inner parts of a core to >100 K very early in the star formation process, and that this in turn strongly suppresses fragmentation. This has implications for a number of outstanding problems in star formation, including the mechanism of massive star formation, the origin of the stellar initial mass function and its relationship to the core mass function, the demographics of massive binaries, and the equation of state in star-forming gas.

The physical mechanism that allows massive stars to form is a major unsolved problem in astrophysics. Stars with masses $\gtsim 20$ $\msun$ reach the main sequence while still embedded in their natal clouds, and the immense radiation output they generate once fusion begins can exert a force stronger than gravity on the dust and gas around them. They also produce huge Lyman continuum luminosities, which can ionize and potentially unbind their parent clouds. This makes massive star formation a more daunting problem than the formation of low mass stars. In this review I present the current state of the field, and discuss the two primary approaches to massive star formation. One holds that the most massive stars form by direct collisions between lower mass stars and their disks. The other approach is to see if the radiation barrier can be overcome by improved treatment of the radiation-hydrodynamic accretion process. I discuss the theoretical background to each model, the observational predictions that can be used to test them, and the substantial parts of the problem that neither theory has fully addressed.

Stars with masses of >~ 20 solar masses have short Kelvin times that enable them to reach the main sequence while still accreting from their natal clouds. The resulting nuclear burning produces a huge luminosity and a correspondingly large radiation pressure force on dust grains in the accreting gas. This effect may limit the upper mass of stars that can form by accretion. Indeed, simulations and analytic calculations to date have been unable to resolve the mystery of how stars of 50 solar masses and up form. We present two new ideas to solve the radiation pressure problem. First, we use three-dimensional radiation hydrodynamic adaptive mesh refinement simulations to study the collapse of massive cores. We find that in three dimensions a configuration in which radiation holds up an infalling envelope is Rayleigh-Taylor unstable, leading radiation driven bubbles to collapse and accretion to continue. We also present Monte Carlo radiative transfer calculations showing that the cavities created by protostellar winds provides a valve that allow radiation to escape the accreting envelope, further reducing the ability of radiation pressure to inhibit accretion.

There are now two dominant models of how stars form: gravitational collapse theory holds that star-forming molecular clumps, typically hundreds to thousands of solar masses in mass, fragment into gaseous cores that subsequently collapse to make individual stars or small multiple systems. In contrast, competitive accretion theory suggests that at birth all stars are much smaller than the typical stellar mass (~0.5 solar masses), and that final stellar masses are determined by the subsequent accretion of unbound gas from the clump. Competitive accretion models explain brown dwarfs and free-floating planets as protostars ejected from star-forming clumps before accreting much mass, predicting that they should lack disks, have high velocity dispersions, and form more frequently in denser clumps. They also predict that mean stellar mass should vary within the Galaxy. Here we derive a simple estimate for the rate of competitive accretion as a function of the star-forming environment, based partly on simulations, and determine in what types of environments competitive accretion can occur. We show that no observed star-forming region produces significant competitive accretion, and that simulations that show competitive accretion do so because their properties differ from those determined by observation. Our result shows that stars form by gravitational collapse, and explains why observations have failed to confirm predictions of the competitive accretion scenario.

The Bondi-Hoyle formula gives the approximate accretion rate onto a point particle accreting from a uniform medium. However, in many situations accretion onto point particles occurs from media that are turbulent rather than uniform. In this paper, we give an approximate solution to the problem of a point particle accreting from an ambient medium of supersonically turbulent gas. Accretion in such media is bimodal, at some points resembling classical Bondi-Hoyle flow, and in other cases being closer to the vorticity-dominated accretion flows recently studied by Krumholz, McKee, & Klein. Based on this observation, we develop a theoretical prediction for the accretion rate, and confirm that our predictions are highly consistent with the results of numerical simulations. The distribution of accretion rates is lognormal, and the mean accretion rate in supersonically turbulent gas can be substantially enhanced above the value that would be predicted by a naive application of the Bondi-Hoyle formula. However, it can also be suppressed by the vorticity, just as Krumholz, McKee, & Klein found for non-supersonic vorticity-dominated flows. Magnetic fields, which we have not included in these models, may further inhibit accretion. Our results have significant implications for a number astrophysical problems, ranging from star formation to the black holes in galactic centers. In particular, there are likely to be significant errors in results that assume that accretion from turbulent media occurs at the unmodified Bondi-Hoyle rate, or that are based on simulations that do not resolve the Bondi-Hoyle radius of accreting objects.

We derive an analytic prediction for the star formation rate in environments ranging from normal galactic disks to starbursts and ULIRGs in terms of the observables of those systems. Our calculation is based on three premises: (1) star formation occurs in virialized molecular clouds that are supersonically turbulent; (2) the density distribution within these clouds is lognormal, as expected for supersonic isothermal turbulence; (3) stars form in any sub-region of a cloud that is so overdense that its gravitational potential energy exceeds the energy in turbulent motions. We show that a theory based on this model is consistent with simulations and with the observed star formation rate in the Milky Way. We use our theory to derive the Kennicutt-Schmidt Law from first principles, and make other predictions that can be tested by future observations. We also provide an algorithm for estimating the star formation rate that is suitable for inclusion in numerical simulations.

We consider the effects of an outflow on radiation escaping from the infalling envelope around a massive protostar. Using numerical radiative transfer calculations, we show that outflows with properties comparable to those observed around massive stars lead to significant anisotropy in the stellar radiation field, which greatly reduces the radiation pressure experienced by gas in the infalling envelope. This means that radiation pressure is a much less significant barrier to massive star formation than has previously been thought.

The classical Bondi-Hoyle formula gives the accretion rate onto a point particle of a gas with a uniform density and velocity. However, the Bondi-Hoyle problem considers only gas with no net vorticity, while in a real astrophysical situation accreting gas invariably has at least a small amount of vorticity. We therefore consider the related case of accretion of gas with constant vorticity, for the cases of both small and large vorticity. We confirm the findings of earlier two dimensional simulations that even a small amount of vorticity can substantially change both the accretion rate and the morphology of the gas flow lines. We show that in three dimensions the resulting flow field is non-axisymmetric and time dependent. The reduction in accretion rate is due to an accumulation of circulation near the accreting particle. Using a combination of simulations and analytic treatment, we provide an approximate formula for the accretion rate of gas onto a point particle as a function of the vorticity of the surrounding gas.

We introduce a new computational method for embedding Lagrangian sink particles into an Eulerian calculation. Simulations of gravitational collapse or accretion generally produce regions whose density greatly exceeds the mean density in the simulation. These dense regions require extremely small time steps to maintain numerical stability. Smoothed particle hydrodynamics (SPH) codes approach this problem by introducing non-gaseous, accreting sink particles, and Eulerian codes may introduce fixed sink cells. However, until now there has been no approach that allows Eulerian codes to follow accretion onto multiple, moving objects. We have removed that limitation by extending the sink particle capability to Eulerian hydrodynamics codes. We have tested this new method and found that it produces excellent agreement with analytic solutions. In analyzing our sink particle method, we present a method for evaluating the disk viscosity parameter $\alpha$ due to the numerical viscosity of a hydrodynamics code, and use it to compute $\alpha$ for our Cartesian AMR code. We also present a simple application of this new method: studying the transition from Bondi to Bondi-Hoyle accretion that occurs when a shock hits a particle undergoing Bondi accretion.

The Wavelength-Oriented Microwave Background Analysis Team (WOMBAT) is constructing microwave maps which will be more realistic than previous simulations. Our foreground models represent a considerable improvement: where spatial templates are available for a given foreground, we predict the flux and spectral index of that component at each place on the sky and estimate uncertainties. We will produce maps containing simulated CMB anisotropy combined with expected foregrounds. The simulated maps will be provided to the community as the WOMBAT Challenge, so such maps can be analyzed to extract cosmological parameters by scientists who are unaware of their input values. This will test the efficacy of foreground subtraction, power spectrum analysis, and parameter estimation techniques and help identify the areas most in need of progress. These maps are also part of the FORECAST project, which allows web-based access to the known foreground maps for the planning of CMB missions.

We have tested several models of GRB luminosity and redshift distribution functions for compatibility with the BATSE 4B number versus peak flux relation. Our results disagree with recent claims that current GRB observations can be used to strongly constrain the cosmic star formation history. Instead, we find that relaxing the assumption that GRBs are standard candles renders a very broad range of models consistent with the BATSE number-flux relation. We explicitly construct two sample distributions, one tracing the star formation history and one with a constant comoving density. We show that both distributions are compatible with the observed fluxes and redshifts of the bursts GRB970508, GRB971214, and GRB980703, and we discuss the measurements required to distinguish the two models.