Source author record

Mark R. Krumholz

Mark R. Krumholz appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.GA astro-ph.SR astro-ph.CO astro-ph.IM astro-ph.HE astro-ph.EP astro-ph physics.flu-dyn physics.comp-ph physics.plasm-ph

Catalog footprint

What is connected

112works

10topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2026arXiv

A survey of molecular clouds in the Galactic center's outflow

The nucleus of the Milky Way is known to drive a large-scale, multiphase galactic outflow, with gas phases ranging from the hot highly-ionized to the cold molecular component. In this work, we present the first systematic search for molecules in the Milky Way wind. We use the Atacama Pathfinder EXperiment (APEX) to observe the 12CO(2-1) emission line in 19 fields centered on previously known high-velocity atomic hydrogen (HI) clouds associated with the outflow. Over 200 CO clumps are detected within 16 different HI clouds. These clumps have typical radii of 1 - 3 parsec, high velocity dispersions of 1 - 6 km/s and molecular gas masses ranging from a few to several hundred solar masses. Molecular clumps in the wind sit on the low-mass end of the mass - size relation of regular molecular clouds, but are far displaced from the mass (or size) - linewidth relation, being generally more turbulent and showing high internal pressures. Nearly 90% of the clumps are gravitationally unbound with virial parameters >> 10 - 100, indicating that these structures are either being disrupted or they must be confined by external pressure from the surrounding hot medium. While the observed properties of CO clumps do not seem to evolve clearly with latitude, we find that molecular gas is not detected in any of the 6 HI clouds with projected distances over 1 kpc from the Galactic Center, suggesting the existence of a maximum timescale of ~ 3 Myr for the dissociation of molecular gas within the wind. Overall, current observations in the Galactic center support a scenario in which a hot wind entrains cold gas clouds from the disk, driving their progressive transformation from molecular to atomic and ultimately ionized gas through stripping, turbulence, and dissociation.

preprint2026arXiv

QED V: Variations in metal loading of galactic winds with element nucleosynthetic origin

Type Ia supernovae, type II supernovae, and asymptotic giant branch (AGB) stars are important sites of stellar nucleosynthesis, but they differ greatly in their rates, their location within a galaxy, and the mean thermal energy and abundance distribution of their ejecta. In earlier papers in this series we have shown that a significant fraction of metals newly synthesized by type II supernovae are promptly lost to galactic winds -- i.e., galactic winds are metal loaded. Here we investigate whether the elements returned by type Ia supernovae and AGB stars are similarly metal loaded, or whether metal loading varies significantly with nucleosynthetic site. We use a series of high-resolution ``tall box'' simulations of the interstellar medium with the \quokka~GPU-accelerated code, within which we systematically vary the galaxy gas surface density, metallicity, and the scale heights and relative rates of the different nucleosynthetic sources. We show that the metal loadings of galactic winds differ substantially between metals produced by different sources, with typical variations at the level of $\approx 0.3$ dex, a phenomenon we term differential metal loading. Which set of metals suffers preferential loss from this phenomenon varies depending on the galactic environment, and is not easily predictable \textit{a priori}. Our findings call into question the the interpretation of diagnostics of galaxy formation, for example star formation timescales and initial mass functions, based on abundance diagnostics, since the abundance variations upon which these techniques rely are often at levels comparable to those we show can be induced by differential metal loading.

preprint2024arXiv

ALMA-LEGUS I: The Influence of Galaxy Morphology on Molecular Cloud Properties

We present a comparative study of the molecular gas in two galaxies from the LEGUS sample: barred spiral NGC 1313 and flocculent spiral NGC 7793. These two galaxies have similar masses, metallicities, and star formation rates, but NGC 1313 is forming significantly more massive star clusters than NGC 7793, especially young massive clusters (<10 Myr, >10^4 Msol). Using ALMA CO(2-1) observations of the two galaxies with the same sensitivities and resolutions of 13 pc, we directly compare the molecular gas in these two similar galaxies to determine the physical conditions responsible for their large disparity in cluster formation. By fitting size-linewidth relations for the clouds in each galaxy, we find that NGC 1313 has a higher intercept than NGC 7793, implying that its clouds have higher kinetic energies at a given size scale. NGC 1313 also has more clouds near virial equilibrium than NGC 7793, which may be connected to its higher rate of massive cluster formation. However, these virially bound clouds do not show a stronger correlation with young clusters than that of the general cloud population. We find surprisingly small differences between the distributions of molecular cloud populations in the two galaxies, though the largest of those differences are that NGC 1313 has higher surface densities and lower free-fall times.

preprint2024arXiv

ALMA-LEGUS II: The Influence of Sub-Galactic Environment on Molecular Cloud Properties

We compare the molecular cloud properties in sub-galactic regions of two galaxies, barred spiral NGC 1313, which is forming many massive clusters, and flocculent spiral NGC 7793, which is forming significantly fewer massive clusters despite having a similar star formation rate to NGC 1313. We find that there are larger variations in cloud properties between different regions within each galaxy than there are between the galaxies on a global scale, especially for NGC 1313. There are higher masses, linewidths, pressures, and virial parameters in the arms of NGC 1313 and center of NGC 7793 than in the interarm and outer regions of the galaxies. The massive cluster formation of NGC 1313 may be driven by its greater variation in environments, allowing more clouds with the necessary conditions to arise, although no one parameter seems primarily responsible for the difference in star formation. Meanwhile NGC 7793 has clouds that are as massive and have as much kinetic energy as clouds in the arms of NGC 1313, but have densities and pressures more similar to the interarm regions and so are less inclined to collapse and form stars. The cloud properties in NGC 1313 and NGC 7793 suggest that spiral arms, bars, interarm regions, and flocculent spirals each represent distinct environments with regard to molecular cloud populations. We see surprisingly little difference in surface densities between the regions, suggesting that the differences in surface densities frequently seen between arm and interarm regions of lower-resolution studies are indicative of the sparsity of molecular clouds, rather than differences in their true surface density.

preprint2022arXiv

Energy balance and Alfvén Mach numbers in compressible magnetohydrodynamic turbulence with a large-scale magnetic field

Energy equipartition is a powerful theoretical tool for understanding astrophysical plasmas. It is invoked, for example, to measure magnetic fields in the interstellar medium (ISM), as evidence for small-scale turbulent dynamo action, and, in general, to estimate the energy budget of star-forming molecular clouds. In this study we motivate and explore the role of the volume-averaged root-mean-squared (rms) magnetic coupling term between the turbulent, $δ\mathbf{B}$ and large-scale, $\mathbf{B}_0$ fields, $\left< (δ\mathbf{B}\cdot\mathbf{B}_0)^{2} \right>^{1/2}_{\mathcal{V}}$. By considering the second moments of the energy balance equations we show that the rms coupling term is in energy equipartition with the volume-averaged turbulent kinetic energy for turbulence with a sub-Alfvénic large-scale field. Under the assumption of exact energy equipartition between these terms, we derive relations for the magnetic and coupling term fluctuations, which provide excellent, parameter-free agreement with time-averaged data from 280 numerical simulations of compressible MHD turbulence. Furthermore, we explore the relation between the turbulent, mean-field and total Alfvén Mach numbers, and demonstrate that sub-Alfvénic turbulence can only be developed through a strong, large-scale magnetic field, which supports an extremely super-Alfvénic turbulent magnetic field. This means that the magnetic field fluctuations are significantly subdominant to the velocity fluctuations in the sub-Alfvénic large-scale field regime. Throughout our study, we broadly discuss the implications for observations of magnetic fields and understanding the dynamics in the magnetised ISM.

preprint2022arXiv

Gamma-ray emission from the Sagittarius Dwarf Spheroidal galaxy due to millisecond pulsars

The Fermi Bubbles are giant, gamma-ray emitting lobes emanating from the nucleus of the Milky Way discovered in ~1-100 GeV data collected by the Large Area Telescope on board the Fermi Gamma-Ray Space Telescope. Previous work has revealed substructure within the Fermi Bubbles that has been interpreted as a signature of collimated outflows from the Galaxy's super-massive black hole. Here we show via a spatial template analysis that much of the gamma-ray emission associated to the brightest region of substructure -- the so-called cocoon -- is likely due to the Sagittarius dwarf spheroidal (Sgr dSph) galaxy. This large Milky Way satellite is viewed through the Fermi Bubbles from the position of the Solar System. As a tidally and ram-pressure stripped remnant, the Sgr dSph has no on-going star formation, but we nevertheless demonstrate that the dwarf's millisecond pulsar (MSP) population can plausibly supply the gamma-ray signal that our analysis associates to its stellar template. The measured spectrum is naturally explained by inverse Compton scattering of cosmic microwave background photons by high-energy electron-positron pairs injected by MSPs belonging to the Sgr dSph, combined with these objects' magnetospheric emission. This finding plausibly suggests that MSPs produce significant gamma-ray emission amongst old stellar populations, potentially confounding indirect dark matter searches in regions such as the Galactic Centre, the Andromeda galaxy, and other massive Milky Way dwarf spheroidals.

preprint2022arXiv

Ion Alfvén velocity fluctuations and implications for the diffusion of streaming cosmic rays

The interstellar medium (ISM) of star-forming galaxies is magnetized and turbulent. Cosmic rays (CRs) propagate through it, and those with energies from $\sim\,\rm{GeV} - \rm{TeV}$ are likely subject to the streaming instability, whereby the wave damping processes balances excitation of resonant ionic Alfvén waves by the CRs, reaching an equilibrium in which the propagation speed of the CRs is very close to the local ion Alfvén velocity. The transport of streaming CRs is therefore sensitive to ionic Alfvén velocity fluctuations. In this paper we systematically study these fluctuations using a large ensemble of compressible MHD turbulence simulations. We show that for sub-Alfvénic turbulence, as applies for a strongly magnetized ISM, the ionic Alfvén velocity probability density function (PDF) is determined solely by the density fluctuations from shocked gas forming parallel to the magnetic field, and we develop analytical models for the ionic Alfvén velocity PDF up to second moments. For super-Alfvénic turbulence, magnetic and density fluctuations are correlated in complex ways, and these correlations as well as contributions from the magnetic fluctuations sets the ionic Alfvén velocity PDF. We discuss the implications of these findings for underlying "macroscopic" diffusion mechanisms in CRs undergoing the streaming instability, including modeling the macroscopic diffusion coefficient for the parallel transport in sub-Alfvénic plasmas. We also describe how, for highly-magnetized turbulent gas, the gas density PDF, and hence column density PDF, can be used to access information about ionic Alfvén velocity structure from observations of the magnetized ISM.

preprint2022arXiv

Quokka: A code for two-moment AMR radiation hydrodynamics on GPUs

We present Quokka, a new subcycling-in-time, block-structured adaptive mesh refinement (AMR) radiation hydrodynamics code optimised for graphics processing units (GPUs). Quokka solves the equations of hydrodynamics with the piecewise parabolic method (PPM) in a method-of-lines formulation, and handles radiative transfer via the variable Eddington tensor (VET) radiation moment equations with a local closure. We use the AMReX library to handle the adaptive mesh management. In order to maximise GPU performance, we combine explicit-in-time evolution of the radiation moment equations with the reduced speed-of-light approximation. We show results for a wide range of test problems for hydrodynamics, radiation, and coupled radiation hydrodynamics. On uniform grids in 3D on a single GPU, our code achieves > 250 million hydrodynamic updates per second and almost 40 million radiation hydrodynamic updates per second. For radiation hydrodynamics problems on uniform grids in 3D, our code scales from 4 GPUs to 256 GPUs with an efficiency of 76 per cent. The code is publicly released under an open-source license on GitHub.

preprint2022arXiv

The evolution of turbulent galactic discs: gravitational instability, feedback and accretion

We study the driving of turbulence in star-froming disc galaxies of different masses at different epochs, using an analytic "bathtub" model. The disc of gas and stars is assumed to be in marginal Toomre instability. Turbulence is assumed to be sustained via an energy balance between its dissipation and three simultaneous energy sources. These are stellar feedback, inward transport due to disc instability and clumpy accretion via streams. The transport rate is computed with two different formalisms, with similar results. To achieve the energy balance, the disc self-regulates either the mass fraction in clumps or the turbulent viscous torque parameter. In this version of the model, the efficiency by which the stream kinetic energy is converted into turbulence is a free parameter, $ξ_a$. We find that the contributions of the three energy sources are in the same ball park, within a factor of $\sim\!2$ in all discs at all times. In haloes that evolve to a mass $\leq 10^{12}\,\Msun$ by $z=0$ ($\leq 10^{11.5}\,\Msun$ at $z\!\sim\!2$), feedback is the main driver throughout their lifetimes. Above this mass, the main driver is either transport or accretion for very low or very high values of $ξ_a$, respectively. For an assumed $ξ_a(t)$ that declines in time, galaxies in halos with present-day mass $>\!10^{12}$ M$_\odot$ make a transition from accretion to transport dominance at intermediate redshifts, $z\! \sim\!3$, when their mass was $\geq\!10^{11.5}\,\Msun$. The predicted relation between star-formation rate and gas velocity dispersion is consistent with observations.

preprint2022arXiv

The Life and Times of Giant Molecular Clouds

Giant molecular clouds (GMCs) are the sites of star formation and stellar feedback in galaxies. Their properties set the initial conditions for star formation and their lifecycles determine how feedback regulates galaxy evolution. In recent years, the advent of high-resolution telescopes has enabled systematic GMC-scale studies of the molecular interstellar medium in nearby galaxies, now covering a wide range of physical conditions and allowing for the first studies of how GMC properties depend on galactic environment. These observational developments have been accompanied by numerical simulations of improving resolution that are increasingly accurately accounting for the effects of the galactic-scale environment on GMCs, while simultaneously improving the treatment of the small-scale processes of star-formation and stellar feedback within them. The combination of these recent developments has greatly improved our understanding of the formation, evolution, and destruction of GMCs. We review the current state of the field, highlight current open questions, and discuss promising avenues for future studies.

preprint2022arXiv

The Volatile Carbon to Oxygen Ratio as a Tracer for the Formation Locations of Interstellar Comets

Based on the occurrence rates implied by the discoveries of 1I/`Oumuamua and 2I/Borisov, the forthcoming Rubin Observatory Legacy Survey of Space and Time (LSST) should detect $\ge1$ interstellar objects every year (Hoover et al. 2021). We advocate for future measurements of the production rates of H$_2$O, CO$_2$ and CO in these objects to estimate their carbon to oxygen ratios, which traces formation locations within their original protoplanetary disks. We review similar measurements for Solar System comets, which indicate formation interior to the CO snowline. By quantifying the relative processing in the interstellar medium and Solar System, we estimate that production rates will not be representative of primordial compositions for the majority of interstellar comets. Preferential desorption of CO and CO$_2$ relative to H$_2$O in the interstellar medium implies that measured C/O ratios represent lower limits on the primordial ratios. Specifically, production rate ratios of ${\rm Q}({\rm CO})/{\rm Q}({\rm H_2O})<.2$ and ${\rm Q}({\rm CO})/{\rm Q}({\rm H_2O})>1$ likely indicate formation interior and exterior to the CO snowline, respectively. The high C/O ratio of 2I/Borisov implies that it formed exterior to the CO snowline. We provide an overview of the currently operational facilities capable of obtaining these measurements that will constrain the fraction of ejected comets that formed exterior to the CO snowline. This fraction will provide key insights into the efficiency of and mechanisms for cometary ejection in exoplanetary systems.

preprint2022arXiv

Theoretical and Observational Evidence for Coriolis Effects in Coronal Magnetic Fields Via Direct Current Driven Flaring Events

All stars produce explosive surface events such as flares and coronal mass ejections. These events are driven by the release of energy stored in coronal magnetic fields, generated by the stellar dynamo. However, it remains unclear if the energy deposition in the magnetic fields is driven by direct or alternating currents. Recently, we presented observational measurements of the flare intensity distributions for a sample of $\sim10^5$ stars across the main sequence observed by $\textit{TESS}$, all of which exhibited power-law distributions similar to those observed in the Sun, albeit with varying slopes. Here we investigate the mechanisms required to produce such a distribution of flaring events via direct current energy deposition, in which coronal magnetic fields braid, reconnect, and produce flares. We adopt a topological model for this process which produces a power-law distribution of energetic flaring events. We expand this model to include the Coriolis effect, which we demonstrate produces a shallower distribution of flare energies in stars that rotate more rapidly (corresponding to a weaker decline in occurrence rates toward increasing flare energies). We present tentative evidence for the predicted rotation-power-law index correlation in the observations. We advocate for future observations of stellar flares that would improve our measurements of the power-law exponents, and yield key insights into the underlying dynamo mechanisms that underpin the self-similar flare intensity distributions.

preprint2022arXiv

VETTAM: A scheme for radiation hydrodynamics with adaptive mesh refinement using the variable Eddington tensor method

We present Variable Eddington Tensor-closed Transport on Adaptive Meshes (\texttt{VETTAM}), a new algorithm to solve the equations of radiation hydrodynamics (RHD) with support for adaptive mesh refinement (AMR) in a frequency-integrated, two-moment formulation. The method is based on a non-local Variable Eddington Tensor (VET) closure computed with a hybrid characteristics scheme for ray tracing. We use a Godunov method for the hyperbolic transport of radiation with an implicit backwards-Euler temporal update to avoid the explicit timestep constraint imposed by the light-crossing time, and a fixed-point Picard iteration scheme to handle the nonlinear gas-radiation exchange term, with the two implicit update stages jointly iterated to convergence. We also develop a modified wave-speed correction method for AMR, which we find to be crucial for obtaining accurate results in the diffusion regime. We demonstrate the robustness of our scheme with a suite of pure radiation and RHD tests, and show that it successfully captures the streaming, static diffusion, and dynamic diffusion regimes and the spatial transitions between them, casts sharp shadows, and yields accurate results for rates of momentum and energy exchange between radiation and gas. A comparison between different closures for the radiation moment equations, with the Eddington approximation (0th-moment closure) and the $M_1$ approximation (1st-moment closure), demonstrates the advantages of the VET method (2nd-moment closure) over the simpler closure schemes. \texttt{VETTAM} has been coupled to the AMR \texttt{FLASH} (magneto-)hydrodynamics code and we summarize by reporting performance features and bottlenecks of our implementation.

preprint2021arXiv

Cosmic rays across the star-forming galaxy sequence. I: Cosmic ray pressures and calorimetry

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.

preprint2021arXiv

High-precision star formation efficiency measurements in nearby clouds

On average molecular clouds convert only a small fraction epsilon_ff of their mass into stars per free-fall time, but differing star formation theories make contrasting claims for how this low mean efficiency is achieved. To test these theories, we need precise measurements of both the mean value and the scatter of epsilon_ff, but high-precision measurements have been difficult because they require determining cloud volume densities, from which we can calculate free-fall times. Until recently, most density estimates assume clouds as uniform spheres, while their real structures are often filamentary and highly non-uniform, yielding systematic errors in epsilon_ff estimates and smearing real cloud-to-cloud variations. We recently developed a theoretical model to reduce this error by using column density distributions in clouds to produce more accurate volume density estimates. In this letter, we apply this model to recent observations of 12 nearby molecular clouds. Compared to earlier analyses, our method reduces the typical dispersion of epsilon_ff within individual clouds from 0.35 dex to 0.31 dex, and decreases the median value of epsilon_ff over all clouds from ~ 0.02 to ~ 0.01. However, we find no significant change in the ~ 0.2 dex cloud-to-cloud dispersion of epsilon_ff, suggesting the measured dispersions reflect real structural differences between clouds.

preprint2021arXiv

Magnetic field amplification in accretion discs around the first stars: implications for the primordial IMF

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 first stars, and thus affect the primordial initial mass function (IMF). We perform a suite of 3D ideal 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 'small-scale turbulent dynamo' and the '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.

preprint2021arXiv

On the origin of nitrogen at low metallicity

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.

preprint2021arXiv

Reconstructing three-dimensional densities from two-dimensional observations of molecular gas

Star formation has long been known to be an inefficient process, in the sense that only a small fraction $ε_{\rm ff}$ of the mass of any given gas cloud is converted to stars per cloud free-fall time. However, developing a successful theory of star formation will require measurements of both the mean value of $ε_{\rm ff}$ and its scatter from one molecular cloud to another. Because $ε_{\rm ff}$ is measured relative to the free-fall time, such measurements require accurate determinations of cloud volume densities. Efforts to measure the volume density from two-dimensional projected data, however, have thus far relied on treating molecular clouds as simple uniform spheres, while their real shapes are likely filamentary and their density distributions far from uniform. The resulting uncertainty in the true volume density is likely one of the major sources of error in observational estimates of $ε_{\rm ff}$. In this paper, we use a suite of simulations of turbulent, magnetized, radiative, self-gravitating star-forming clouds to examine whether it is possible to obtain more accurate volume density estimates and thereby reduce this error. We create mock observations from simulations, and show that current analysis methods relying on the spherical assumption likely yield ~ 0.26 dex underestimations and ~ 0.51 dex errors in volume density estimates, corresponding to a ~ 0.13 dex overestimation and a ~ 0.25 dex scatter in $ε_{\rm ff}$, comparable to the scatter in observed cloud samples. We build a predictive model that uses information accessible in two-dimensional measurements -- most significantly the Gini coefficient of the surface density distribution -- to estimate volume density with ~ 0.3 dex less scatter. We test our method on a recent observation of the Ophiuchus cloud, and show that it successfully reduces the $ε_{\rm ff}$ scatter.

preprint2021arXiv

The physics of gas phase metallicity gradients in galaxies

We present a new model for the evolution of gas phase metallicity gradients in galaxies from first principles. We show that metallicity gradients depend on four ratios that collectively describe the metal equilibration timescale, production, transport, consumption, and loss. Our model finds that most galaxy metallicity gradients are in equilibrium at all redshifts. When normalized by metal diffusion, metallicity gradients are governed by the competition between radial advection, metal production, and accretion of metal-poor gas from the cosmic web. The model naturally explains the varying gradients measured in local spirals, local dwarfs, and high-redshift star-forming galaxies. We use the model to study the cosmic evolution of gradients across redshift, showing that the gradient in Milky Way-like galaxies has steepened over time, in good agreement with both observations and simulations. We also predict the evolution of metallicity gradients with redshift in galaxy samples constructed using both matched stellar masses and matched abundances. Our model shows that massive galaxies transition from the advection-dominated to the accretion-dominated regime from high to low redshifts, which mirrors the transition from gravity-driven to star formation feedback-driven turbulence. Lastly, we show that gradients in local ultraluminous infrared galaxies (major mergers) and inverted gradients seen both in the local and high-redshift galaxies may not be in equilibrium. In subsequent papers in this series, we show that the model also explains the observed relationship between galaxy mass and metallicity gradients, and between metallicity gradients and galaxy kinematics.

preprint2020arXiv

An uncertainty principle for star formation -- V. The influence of dust extinction on star formation rate tracer lifetimes and the inferred molecular cloud lifecycle

Recent observational studies aiming to quantify the molecular cloud lifecycle require the use of known 'reference time-scales' to turn the relative durations of different phases of the star formation process into absolute time-scales. We previously constrained the characteristic emission time-scales of different star formation rate (SFR) tracers, as a function of the SFR surface density and metallicity. However, we omitted the effects of dust extinction. Here, we extend our suite of SFR tracer emission time-scales by accounting for extinction, using synthetic emission maps of a high-resolution hydrodynamical simulation of an isolated, Milky-Way-like disc galaxy. The stellar feedback included in the simulation is inefficient compared to observations, implying that it represents a limiting case in which the duration of embedded star formation (and the corresponding effect of extinction) is overestimated. Across our experiments, we find that extinction mostly decreases the SFR tracer emission time-scale, changing the time-scales by factors of 0.04-1.74, depending on the gas column density. UV filters are more strongly affected than H$α$ filters. We provide the limiting correction factors as a function of the gas column density and flux sensitivity limit for a wide variety of SFR tracers. Applying these factors to observational characterisations of the molecular cloud lifecycle produces changes that broadly fall within the quoted uncertainties, except at high kpc-scale gas surface densities ($Σ_{\rm g}\gtrsim20~{\mathrm{M_{\odot}\,pc^{-2}}}$). Under those conditions, correcting for extinction may decrease the measured molecular cloud lifetimes and feedback time-scales, which further strengthens previous conclusions that molecular clouds live for a dynamical time and are dispersed by early, pre-supernova feedback.

preprint2020arXiv

Cosmic ray transport in starburst galaxies

Starburst galaxies are efficient $γ$-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 $γ$-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 $γ$-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 $γ$-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.

preprint2020arXiv

Cosmic rays across the star-forming galaxy sequence. II: Stability limits and the onset of cosmic ray-driven outflows

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.

preprint2020arXiv

Distribution and kinematics of 26Al in the Galactic disc

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.

preprint2020arXiv

Dynamics of small grains in transitional discs

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, and where the ratio of accretion rate to stellar luminosity is low ($\dot{M}/L \lesssim 10^{-10}$ $M_\odot$ yr$^{-1}$ $L_\odot^{-1}$), radiative clearing of any remaining $\sim 0.5$ $μ$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.

preprint2020arXiv

Evolution of Stellar Feedback in HII Regions

Stellar feedback is needed to produce realistic giant molecular clouds (GMCs) and galaxies in simulations, but due to limited numerical resolution, feedback must be implemented using subgrid models. Observational work is an important means to test and anchor these models, but limited studies have assessed the relative dynamical role of multiple feedback modes, particularly at the earliest stages of expansion when HII regions are still deeply embedded. In this paper, we use multiwavelength (radio, infrared, and X-ray) data to measure the pressures associated with direct radiation ($P_{\rm dir}$), dust-processed radiation ($P_{\rm IR}$), photoionization heating ($P_{\rm HII}$), and shock-heating from stellar winds ($P_{\rm X}$) in a sample of 106 young, resolved HII regions with radii $\lesssim$0.5 pc to determine how stellar feedback drives their expansion. We find that the $P_{\rm IR}$ dominates in 84% of the regions and that the median $P_{\rm dir}$ and $P_{\rm HII}$ are smaller than the median $P_{\rm IR}$ by factors of $\approx 6$ and $\approx 9$, respectively. Based on the radial dependences of the pressure terms, we show that HII regions transition from $P_{\rm IR}$-dominated to $P_{\rm HII}$-dominated at radii of $\sim$3 pc. We find a median trapping factor of $f_{\rm trap} \sim$ 8 without any radial dependence for the sample, suggesting this value can be adopted in sub-grid feedback models. Moreover, we show that the total pressure is greater than the gravitational pressure in the majority of our sample, indicating that the feedback is sufficient to expel gas from the regions.

preprint2020arXiv

How do bound star clusters form?

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.

preprint2020arXiv

Quantifying the effects of spatial resolution and noise on galaxy metallicity gradients

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.

preprint2020arXiv

The importance of magnetic fields for the initial mass function of the first stars

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.

preprint2020arXiv

The Life Cycle of the Central Molecular Zone. II: Distribution of atomic and molecular gas tracers

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.

preprint2020arXiv

The Role of Outflows, Radiation Pressure, and Magnetic Fields in Massive Star Formation

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.

preprint2020arXiv

The SAMI Galaxy Survey: Gas velocity dispersions in low-$z$ star-forming galaxies and the drivers of turbulence

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 ($σ_{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 ($Σ_\text{SFR}$). Using the total sample, we find $σ_{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 $σ_{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 $σ_{v, z}$ floor driven by star-formation feedback processes with an upturn in $σ_{v, z}$ at higher SFR driven by gravitational transport of gas through the disc.

preprint2019arXiv

Momentum Injection by Clustered Supernovae: Testing Subgrid Feedback Prescriptions

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.

preprint2016arXiv

An Unstable Truth: How Massive Stars get their Mass

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.

preprint2016arXiv

Enhanced Momentum Feedback from Clustered Supernovae

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.

preprint2016arXiv

Hybrid Adaptive Ray-Moment Method (HARM$^2$): A Highly Parallel Method for Radiation Hydrodynamics on Adaptive Grids

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.

preprint2016arXiv

Instability of Supersonic Cold Streams Feeding Galaxies I: Linear Kelvin-Helmholtz Instability with Body Modes

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.

preprint2016arXiv

Is Turbulence in the Interstellar Medium Driven by Feedback or Gravity? An Observational Test

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 σ^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 σ$). 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.

preprint2016arXiv

Mass Transport and Turbulence in Gravitationally Unstable Disk Galaxies II: The Effects of Star Formation Feedback

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.

preprint2016arXiv

Suppression of star formation in dwarf galaxies by grain photoelectric feedback

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.

preprint2016arXiv

The Bursty Star Formation Histories of Low-mass Galaxies at $0.4<z<1$ Revealed by Star Formation Rates Measured from H$β$ and FUV

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$β$ and FUV (1500 Å) (H$β$--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$β$- and FUV-derived SFRs is measured, dust extinction correction is negligible (except for very dusty galaxies) with the Calzetti attenuation curve. The H$β$--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$β$--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$β$--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.

preprint2016arXiv

What Physics Determines the Peak of the IMF? Insights from the Structure of Cores in Radiation-Magnetohydrodynamic Simulations

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$.

preprint2015arXiv

A Dynamical Model for the Formation of Gas Rings and Episodic Starbursts Near Galactic Centres

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.

preprint2015arXiv

Mass Transport and Turbulence in Gravitationally Unstable Disk Galaxies. I: The Case of Pure Self-Gravity

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.

preprint2015arXiv

Mixing and transport of metals by gravitational instability-driven turbulence in galactic discs

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.

preprint2015arXiv

Simulated Photoevaporative Mass Loss from Hot Jupiters in 3D

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.

preprint2015arXiv

SLUG -- Stochastically Lighting Up Galaxies. III: A Suite of Tools for Simulated Photometry, Spectroscopy, and Bayesian Inference with Stochastic Stellar Populations

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.

preprint2015arXiv

Star Cluster Properties in Two LEGUS Galaxies Computed with Stochastic Stellar Population Synthesis Models

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.

preprint2015arXiv

Sub-Eddington Star-Forming Regions are Super-Eddington: Momentum Driven Outflows from Supersonic Turbulence

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.

preprint2015arXiv

VADER: A Flexible, Robust, Open-Source Code for Simulating Viscous Thin Accretion Disks

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.

preprint2014arXiv

Early turbulent mixing as the origin of chemical homogeneity in open star clusters

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.

preprint2014arXiv

Gone With the Wind: Where is the Missing Stellar Wind Energy from Massive Star Clusters?

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.

preprint2014arXiv

Massive Star Formation

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.

preprint2014arXiv

SLUG - Stochastically Lighting Up Galaxies. II: Quantifying the Effects of Stochasticity on Star Formation Rate Indicators

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.

preprint2014arXiv

Star Cluster Formation and Feedback

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.

preprint2014arXiv

The Big Problems in Star Formation: the Star Formation Rate, Stellar Clustering, and the Initial Mass Function

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.

preprint2014arXiv

The Formation of Very Massive Stars

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.

preprint2014arXiv

The Role of Stellar Feedback in the Dynamics of HII Regions

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.

preprint2013arXiv

An IMF Study of the Dwarf Starburst Galaxy NGC 4214

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.

preprint2013arXiv

Balance among gravitational instability, star formation, and accretion determines the structure and evolution of disk galaxies

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.

preprint2013arXiv

Close Stellar Encounters in Young, Substructured, Dissolving Star Clusters: Statistics and Effects on Planetary Systems

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.

preprint2013arXiv

DESPOTIC -- A New Software Library to Derive the Energetics and SPectra of Optically Thick Interstellar Clouds

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.

preprint2013arXiv

Formation of Molecular Clouds and Global Conditions for Star Formation

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.

preprint2013arXiv

Line Emission from Radiation-Pressurized HII Region II: Dynamics and Population Synthesis

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.

preprint2013arXiv

Line Emission from Radiation-Pressurized HII Regions I: Internal Structure and Line Ratios

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.

preprint2013arXiv

Numerical Simulations of Radiatively-Driven Dusty Winds

[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.

preprint2013arXiv

On the origin of the fundamental metallicity relation and the scatter in galaxy scaling relations

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.

preprint2013arXiv

Protostellar Disk Formation Enabled by Weak, Misaligned Magnetic Fields

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.

preprint2013arXiv

Steady Outflows in Giant Clumps of High-z Disk Galaxies During Migration and Growth by Accretion

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.

preprint2013arXiv

The Fragmentation of Magnetized, Massive Star-Forming Cores with Radiative Feedback

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.

preprint2013arXiv

The Star Formation Law in Molecule-Poor Galaxies

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.

preprint2013arXiv

Very Massive Stars (VMS) in the Local Universe

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.

preprint2012arXiv

A General Model for the CO-H2 Conversion Factor in Galaxies with Applications to the Star Formation Law

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.

preprint2012arXiv

Direct Numerical Simulation of Radiation Pressure-Driven Turbulence and Winds in Star Clusters and Galactic Disks

[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.

preprint2012arXiv

Evolving Gravitationally Unstable Disks Over Cosmic Time: Implications For Thick Disk Formation

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.

preprint2012arXiv

Metallicity-Dependent quenching of Star Formation at High Redshift in Small Galaxies

[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.

preprint2012arXiv

Radiation-Hydrodynamic Simulations of the Formation of Orion-Like Star Clusters II. The Initial Mass Function from Winds, Turbulence, and Radiation

[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.

preprint2012arXiv

Radiatively Efficient Magnetized Bondi Accretion

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 β^{1/2} compared to the Bondi value, where β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.

preprint2012arXiv

Star Formation in Atomic Gas

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.

preprint2012arXiv

Was The Sun Born In A Massive Cluster?

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.

preprint2012arXiv

What Sets the Initial Rotation Rates of Massive Stars?

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.

preprint2011arXiv

A Comparison of Methods for Determining the Molecular Content of Model Galaxies

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.

preprint2011arXiv

A Universal, Local Star Formation Law in Galactic Clouds, Nearby Galaxies, High-Redshift Disks, and Starbursts

[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.

preprint2011arXiv

COLD GASS, an IRAM legacy survey of molecular gas in massive galaxies: I. Relations between H2, HI, stellar content and structural properties

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.

preprint2011arXiv

COLD GASS, an IRAM Legacy Survey of Molecular Gas in Massive Galaxies: II. The non-universality of the Molecular Gas Depletion Timescale

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.

preprint2011arXiv

Evolution of Blister-Type HII Regions in a Magnetized Medium

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.

preprint2011arXiv

Metallicity and the Universality of the IMF

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.

preprint2011arXiv

On the Origin of Stellar Masses

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.

preprint2011arXiv

On the Reliability of Stellar Ages and Age Spreads Inferred from Pre-Main Sequence Evolutionary Models

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.

preprint2011arXiv

Radiation-Hydrodynamic Simulations of Massive Star Formation with Protostellar Outflows

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.

preprint2011arXiv

Radiation-Hydrodynamic Simulations of the Formation of Orion-Like Star Clusters I. Implications for the Origin of the Initial Mass Function

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.

preprint2011arXiv

Spin down of protostars through gravitational torques

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.

preprint2011arXiv

Star Formation in Molecular Clouds

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.

preprint2011arXiv

Stochastic star formation and a (nearly) uniform stellar initial mass function

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.

preprint2011arXiv

The Global Evolution of Giant Molecular Clouds II: The Role of Accretion

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.

preprint2011arXiv

What Drives the Expansion of Giant HII Regions?: A Study of Stellar Feedback in 30 Doradus

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.

preprint2011arXiv

What Phase of the Interstellar Medium Correlates with the Star Formation Rate?

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.

preprint2010arXiv

On the Dynamics and Evolution of Gravitational Instability-Dominated Disks

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.

preprint2010arXiv

Radiation Feedback, Fragmentation, and the Environmental Dependence of the Initial Mass Function

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.

preprint2010arXiv

Star Formation with Adaptive Mesh Refinement Radiation Hydrodynamics

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.

preprint2010arXiv

Stellar Feedback in Molecular Clouds and its Influence on the Mass Function of Young Star Clusters

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.

preprint2010arXiv

Survival of Star-Forming Giant Clumps in High-Redshift Galaxies

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.

preprint2010arXiv

Testing models for molecular gas formation in galaxies: hydrostatic pressure or gas and dust shielding?

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.

preprint2010arXiv

The Atomic to Molecular Transition and its Relation to the Scaling Properties of Galaxy Disks in the Local Universe

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 $λ$. 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 $Σ_{SFR}\proptoΣ_{H_2}$ in the regime where the molecular gas dominates the total gas surface density, and $Σ_{SFR}\propto Σ_{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.

preprint2009arXiv

Molecular hydrogen deficiency in HI-poor galaxies and its implications for star formation

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.

preprint2009arXiv

Numerical Star-Formation Studies -- A Status Report

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.

preprint2009arXiv

On the Role of Disks in the Formation of Stellar Systems: A Numerical Parameter Study of Rapid Accretion

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.

preprint2009arXiv

The Atomic-to-Molecular Transition in Galaxies. III. A New Method for Determining the Molecular Content of Primordial and Dusty Clouds

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.

preprint2009arXiv

The Formation of Massive Star Systems by Accretion

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.

preprint2009arXiv

The Star Formation Law in Atomic and Molecular Gas

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.

preprint2007arXiv

Mass Transfer in Close, Rapidly Accreting Protobinaries: An Origin for Massive Twins?

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.

preprint2006arXiv

Turbulence, Feedback, and Slow Star Formation

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.

preprint2005arXiv

A General Theory of Turbulence-Regulated Star Formation, From Spirals to ULIRGs

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.

Mark R. Krumholz

What is connected

Connect this record

See the researcher in context

Building this map preview

112 published item(s)

A survey of molecular clouds in the Galactic center's outflow

QED V: Variations in metal loading of galactic winds with element nucleosynthetic origin

ALMA-LEGUS I: The Influence of Galaxy Morphology on Molecular Cloud Properties

ALMA-LEGUS II: The Influence of Sub-Galactic Environment on Molecular Cloud Properties

Energy balance and Alfvén Mach numbers in compressible magnetohydrodynamic turbulence with a large-scale magnetic field

Gamma-ray emission from the Sagittarius Dwarf Spheroidal galaxy due to millisecond pulsars

Ion Alfvén velocity fluctuations and implications for the diffusion of streaming cosmic rays

Quokka: A code for two-moment AMR radiation hydrodynamics on GPUs

The evolution of turbulent galactic discs: gravitational instability, feedback and accretion

The Life and Times of Giant Molecular Clouds

The Volatile Carbon to Oxygen Ratio as a Tracer for the Formation Locations of Interstellar Comets

Theoretical and Observational Evidence for Coriolis Effects in Coronal Magnetic Fields Via Direct Current Driven Flaring Events

VETTAM: A scheme for radiation hydrodynamics with adaptive mesh refinement using the variable Eddington tensor method

Cosmic rays across the star-forming galaxy sequence. I: Cosmic ray pressures and calorimetry

High-precision star formation efficiency measurements in nearby clouds

Magnetic field amplification in accretion discs around the first stars: implications for the primordial IMF

On the origin of nitrogen at low metallicity

Reconstructing three-dimensional densities from two-dimensional observations of molecular gas

The physics of gas phase metallicity gradients in galaxies

An uncertainty principle for star formation -- V. The influence of dust extinction on star formation rate tracer lifetimes and the inferred molecular cloud lifecycle

Cosmic ray transport in starburst galaxies

Cosmic rays across the star-forming galaxy sequence. II: Stability limits and the onset of cosmic ray-driven outflows

Distribution and kinematics of 26Al in the Galactic disc

Dynamics of small grains in transitional discs

Evolution of Stellar Feedback in HII Regions

How do bound star clusters form?

Quantifying the effects of spatial resolution and noise on galaxy metallicity gradients

The importance of magnetic fields for the initial mass function of the first stars

The Life Cycle of the Central Molecular Zone. II: Distribution of atomic and molecular gas tracers

The Role of Outflows, Radiation Pressure, and Magnetic Fields in Massive Star Formation

The SAMI Galaxy Survey: Gas velocity dispersions in low-$z$ star-forming galaxies and the drivers of turbulence

Momentum Injection by Clustered Supernovae: Testing Subgrid Feedback Prescriptions

An Unstable Truth: How Massive Stars get their Mass

Enhanced Momentum Feedback from Clustered Supernovae

Hybrid Adaptive Ray-Moment Method (HARM$^2$): A Highly Parallel Method for Radiation Hydrodynamics on Adaptive Grids

Instability of Supersonic Cold Streams Feeding Galaxies I: Linear Kelvin-Helmholtz Instability with Body Modes

Is Turbulence in the Interstellar Medium Driven by Feedback or Gravity? An Observational Test

Mass Transport and Turbulence in Gravitationally Unstable Disk Galaxies II: The Effects of Star Formation Feedback

Suppression of star formation in dwarf galaxies by grain photoelectric feedback

The Bursty Star Formation Histories of Low-mass Galaxies at $0.4<z<1$ Revealed by Star Formation Rates Measured from H$β$ and FUV

What Physics Determines the Peak of the IMF? Insights from the Structure of Cores in Radiation-Magnetohydrodynamic Simulations

A Dynamical Model for the Formation of Gas Rings and Episodic Starbursts Near Galactic Centres

Mass Transport and Turbulence in Gravitationally Unstable Disk Galaxies. I: The Case of Pure Self-Gravity

Mixing and transport of metals by gravitational instability-driven turbulence in galactic discs

Simulated Photoevaporative Mass Loss from Hot Jupiters in 3D

SLUG -- Stochastically Lighting Up Galaxies. III: A Suite of Tools for Simulated Photometry, Spectroscopy, and Bayesian Inference with Stochastic Stellar Populations

Star Cluster Properties in Two LEGUS Galaxies Computed with Stochastic Stellar Population Synthesis Models

Sub-Eddington Star-Forming Regions are Super-Eddington: Momentum Driven Outflows from Supersonic Turbulence

VADER: A Flexible, Robust, Open-Source Code for Simulating Viscous Thin Accretion Disks

Early turbulent mixing as the origin of chemical homogeneity in open star clusters

Gone With the Wind: Where is the Missing Stellar Wind Energy from Massive Star Clusters?

Massive Star Formation

SLUG - Stochastically Lighting Up Galaxies. II: Quantifying the Effects of Stochasticity on Star Formation Rate Indicators

Star Cluster Formation and Feedback

The Big Problems in Star Formation: the Star Formation Rate, Stellar Clustering, and the Initial Mass Function

The Formation of Very Massive Stars

The Role of Stellar Feedback in the Dynamics of HII Regions

An IMF Study of the Dwarf Starburst Galaxy NGC 4214

Balance among gravitational instability, star formation, and accretion determines the structure and evolution of disk galaxies

Close Stellar Encounters in Young, Substructured, Dissolving Star Clusters: Statistics and Effects on Planetary Systems

DESPOTIC -- A New Software Library to Derive the Energetics and SPectra of Optically Thick Interstellar Clouds

Formation of Molecular Clouds and Global Conditions for Star Formation

Line Emission from Radiation-Pressurized HII Region II: Dynamics and Population Synthesis

Line Emission from Radiation-Pressurized HII Regions I: Internal Structure and Line Ratios

Numerical Simulations of Radiatively-Driven Dusty Winds

On the origin of the fundamental metallicity relation and the scatter in galaxy scaling relations

Protostellar Disk Formation Enabled by Weak, Misaligned Magnetic Fields

Steady Outflows in Giant Clumps of High-z Disk Galaxies During Migration and Growth by Accretion

The Fragmentation of Magnetized, Massive Star-Forming Cores with Radiative Feedback

The Star Formation Law in Molecule-Poor Galaxies

Very Massive Stars (VMS) in the Local Universe

A General Model for the CO-H2 Conversion Factor in Galaxies with Applications to the Star Formation Law

Direct Numerical Simulation of Radiation Pressure-Driven Turbulence and Winds in Star Clusters and Galactic Disks

Evolving Gravitationally Unstable Disks Over Cosmic Time: Implications For Thick Disk Formation