Nathan J. Goldbaum

ENZO: AN ADAPTIVE MESH REFINEMENT CODE FOR ASTROPHYSICS

This paper describes the open-source code Enzo, which uses block-structured adaptive mesh refinement to provide high spatial and temporal resolution for modeling astrophysical fluid flows. The code is Cartesian, can be run in one, two, and three dimensions, and supports a wide variety of physics including hydrodynamics, ideal and non-ideal magnetohydrodynamics, N-body dynamics (and, more broadly, self-gravity of fluids and particles), primordial gas chemistry, optically thin radiative cooling of primordial and metal-enriched plasmas (as well as some optically-thick cooling models), radiation transport, cosmological expansion, and models for star formation and feedback in a cosmological context. In addition to explaining the algorithms implemented, we present solutions for a wide range of test problems, demonstrate the codes parallel performance, and discuss the Enzo collaborations code development methodology.

THE GLOBAL EVOLUTION OF GIANT MOLECULAR CLOUDS. II. THE ROLE OF ACCRETION

We present virial models for the global evolution of giant molecular clouds (GMCs). 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 GMCs. Our model clouds reproduce the scaling relations observed in both galactic and extragalactic clouds. We find that accretion and star formation 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 M ☉ pc–2, in good agreement with observations of GMCs 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 GMCs and associated young star clusters in the Large Magellanic Cloud and find good agreement between our model clouds and the observed relationship between H II regions, young star clusters, and GMCs.

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

Suppression of star formation in dwarf galaxies by photoelectric grain heating feedback

\sim 1

DWARF GALAXIES WITH IONIZING RADIATION FEEDBACK. II: SPATIALLY-RESOLVED STAR FORMATION RELATION

M_\odot

GRACKLE: a chemistry and cooling library for astrophysics

yr

The AGORA High-resolution Galaxy Simulations Comparison Project II: Isolated disk test

^{-1}

GAMER-2: a GPU-accelerated adaptive mesh refinement code — accuracy, performance, and scalability

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.

unyt: Handle, manipulate, and convert data with units in Python

Photoelectric heating—heating of dust grains by far-ultraviolet photons—has long been recognized as the primary source of heating for the neutral interstellar medium. Simulations of spiral galaxies have shown some indication that photoelectric heating could suppress star formation; however, simulations that include photoelectric heating have typically shown that it has little effect on the rate of star formation in either spiral galaxies or dwarf galaxies, which suggests that supernovae are responsible for setting the gas depletion time in galaxies. This result is in contrast with recent work indicating that a star formation law that depends on galaxy metallicity—as is expected with photoelectric heating, but not with supernovae—reproduces the present-day galaxy population better than does a metallicity-independent one. Here we report a series of simulations of dwarf galaxies, the class of galaxy in which the effects of both photoelectric heating and supernovae are expected to be strongest. We simultaneously include space- and time-dependent photoelectric heating in our simulations, and we resolve the energy-conserving phase of every supernova blast wave, which allows us to directly measure the relative importance of feedback by supernovae and photoelectric heating in suppressing star formation. We find that supernovae are unable to account for the observed large gas depletion times in dwarf galaxies. Instead, photoelectric heating is the dominant means by which dwarf galaxies regulate their star formation rate at any given time, suppressing the rate by more than an order of magnitude relative to simulations with only supernovae.