Search

Dust attenuation varies substantially from galaxy to galaxy and as of yet cannot be reproduced from first principles in theoretical models. In Nagaraj et al. (2022), we developed the first Bayesian population model of dust attenuation as a function of stellar population properties and projected galaxy shape, built on spectral energy distribution (SED) fits of nearly 30,000 galaxies in the 3D-HST grism survey with broadband photometric coverage from the rest-frame UV to IR. In this paper, we apply the model to galaxies from the large-volume cosmological simulation TNG100. We produce a UVJ diagram and compare it with one obtained in previous work by applying approximate radiative transfer to the simulated galaxies. We find that the UVJ diagram based on our empirical model is in better agreement with observations than the previous effort, especially in the number density of dusty star forming galaxies. We also construct the intrinsic dust-free UVJ diagram for TNG and 3D-HST galaxies at z ~ 1, finding qualitative agreement but residual differences at the 10-20% level. These differences can be largely attributed to the finding that TNG galaxies have, on average, 29% younger stellar populations and 0.28 dex higher metallicities than observed galaxies., 12 pages, 5 figures, submitted to ApJ

We present a catalog of spectroscopically measured redshifts over $0 < z < 2$ and emission line fluxes for 1440 galaxies. The majority ($\sim$65\%) of the galaxies come from the HALO7D survey, with the remainder from the DEEPwinds program. This catalog includes redshifts for 646 dwarf galaxies with $\log(M_{\star}/M_{\odot}) < 9.5$. 810 catalog galaxies did not have previously published spectroscopic redshifts, including 454 dwarf galaxies. HALO7D used the DEIMOS spectrograph on the Keck II telescope to take very deep (up to 32 hours exposure, with a median of $\sim$7 hours) optical spectroscopy in the COSMOS, EGS, GOODS-North, and GOODS-South CANDELS fields, and in some areas outside CANDELS. We compare our redshift results to existing spectroscopic and photometric redshifts in these fields, finding only a 1\% rate of discrepancy with other spectroscopic redshifts. We measure a small increase in median photometric redshift error (from 1.0\% to 1.3\%) and catastrophic outlier rate (from 3.5\% to 8\%) with decreasing stellar mass. We obtained successful redshift fits for 75\% of massive galaxies, and demonstrate a similar 70-75\% successful redshift measurement rate in $8.5 < \log(M_{\star}/M_{\odot}) < 9.5$ galaxies, suggesting similar survey sensitivity in this low-mass range. We describe the redshift, mass, and color-magnitude distributions of the catalog galaxies, finding HALO7D galaxies representative of CANDELS galaxies up to \textit{i}-band magnitudes of 25. The catalogs presented will enable studies of star formation (SF), the mass-metallicity relation, SF-morphology relations, and other properties of the $z\sim0.7$ dwarf galaxy population., 23 pages, 19 Figures, updated to version accepted by ApJS

Dust plays a pivotal role in determining the observed spectral energy distribution (SED) of galaxies. Yet our understanding of dust attenuation is limited and our observations suffer from the dust-metallicity-age degeneracy in SED fitting (single galaxies), large individual variances (ensemble measurements), and the difficulty in properly dealing with uncertainties (statistical considerations). In this study, we create a population Bayesian model to rigorously account for correlated variables and non-Gaussian error distributions and demonstrate the improvement over a simple Bayesian model. We employ a flexible 5-D linear interpolation model for the parameters that control dust attenuation curves as a function of stellar mass, star formation rate (SFR), metallicity, redshift, and inclination. Our setup allows us to determine the complex relationships between dust attenuation and these galaxy properties simultaneously. Using Prospector fits of nearly 30,000 3D-HST galaxies, we find that the attenuation slope ($n$) flattens with increasing optical depth ($\tau$), though less so than in previous studies. $\tau$ increases strongly with SFR, though when $\log~{\rm SFR}\lesssim 0$, $\tau$ remains roughly constant over a wide range of stellar masses. Edge-on galaxies tend to have larger $\tau$ than face-on galaxies, but only for $\log~M_*\gtrsim 10$, reflecting the lack of triaxiality for low-mass galaxies. Redshift evolution of dust attenuation is strongest for low-mass, low-SFR galaxies, with higher optical depths but flatter curves at high redshift. Finally, $n$ has a complex relationship with stellar mass, highlighting the intricacies of the star-dust geometry. We have publicly released software (https://github.com/Astropianist/DustE) for users to access our population model., Comment: 27 pages, 14 figures, 1 table. Submitted to ApJ

The driving of turbulence in galaxies is deeply connected with the physics of feedback, star formation, outflows, accretion, and radial transport in disks. The velocity dispersion of gas in galaxies therefore offers a promising observational window into these processes. However, the relative importance of each of these mechanisms remains controversial. In this work we revisit the possibility that turbulence on galactic scales is driven by the direct impact of accreting gaseous material on the disk. We measure this effect in a disk-like star-forming galaxy in IllustrisTNG, using the high-resolution cosmological magnetohydrodynamical simulation TNG50. We employ Lagrangian tracer particles with a high time cadence of only a few Myr to identify accretion and other events, such as star formation, outflows, and movement within the disk. The energies of particles as they arrive in the disk are measured by stacking the events in bins of time before and after the event. The average effect of each event is measured on the galaxy by fitting explicit models for the kinetic and turbulent energies as a function of time in the disk. These measurements are corroborated by measuring the cross-correlation of the turbulent energy in the different annuli of the disk with other time series, and searching for signals of causality, i.e. asymmetries in the cross-correlation across zero time lag. We find that accretion contributes to the large-scale turbulent kinetic energy even if it is not the dominant driver of turbulence in this $\sim 5 \times 10^{9} M_\odot$ stellar mass galaxy. Extrapolating this finding to a range of galaxy masses, we find that there are regimes where energy from direct accretion may dominate the turbulent energy budget, particularly in disk outskirts, galaxies less massive than the Milky Way, and at redshift $\sim 2$., Comment: Submitted to AAS Journals, comments welcome

We characterize mass, momentum, energy and metal outflow rates of multi-phase galactic winds in a suite of FIRE-2 cosmological "zoom-in" simulations from the Feedback in Realistic Environments (FIRE) project. We analyze simulations of low-mass dwarfs, intermediate-mass dwarfs, Milky Way-mass halos, and high-redshift massive halos. Consistent with previous work, we find that dwarfs eject about 100 times more gas from their interstellar medium (ISM) than they form in stars, while this mass "loading factor" drops below one in massive galaxies. Most of the mass is carried by the hot phase ($>10^5$ K) in massive halos and the warm phase ($10^3-10^5$ K) in dwarfs; cold outflows ($, Accepted to MNRAS with minor revisions, main body is 25 pages with 14 figures

Under the right conditions brown dwarfs that gain enough mass late in their lives to cross the hydrogen burning limit will not turn into low-mass stars, but rather remain essentially brown dwarf-like. While these objects, called either beige dwarfs or over-massive brown dwarfs, may exist in principle, it remains unclear exactly how they would form astrophysically. We show that accretion from AGB winds, aided by the wind Roche lobe overflow mechanism, is likely to produce a substantial population of observable overmassive brown dwarfs, though other mechanisms are still plausible. Specifically we predict that sun-like stars born with a massive brown dwarf companion on an orbit with a semi-major axis of order 10 AU will likely produce overmassive brown dwarfs, which may be found today as companions to the donor star's remnant white dwarf. The identification and characterization of such an object would produce unique constraints on binary evolution because there is a solid upper limit on the brown dwarf's initial mass., Accepted at ApJ, includes a new figure showing predictions for Gaia

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., 16 figures, 3 appendices, 28 pages. MNRAS, in press; v2 - updated references

In addition to the well-known gas phase mass-metallicity relation (MZR), recent spatially-resolved observations have shown that local galaxies also obey a mass-metallicity gradient relation (MZGR) whereby metallicity gradients can vary systematically with galaxy mass. In this work, we use our recently-developed analytic model for metallicity distributions in galactic discs, which includes a wide range of physical processes -- radial advection, metal diffusion, cosmological accretion, and metal-enriched outflows -- to simultaneously analyse the MZR and MZGR. We show that the same physical principles govern the shape of both: centrally-peaked metal production favours steeper gradients, and this steepening is diluted by the addition of metal-poor gas, which is supplied by inward advection for low-mass galaxies and by cosmological accretion for massive galaxies. The MZR and the MZGR both bend at galaxy stellar mass $\sim 10^{10} - 10^{10.5}\,\rm{M_{\odot}}$, and we show that this feature corresponds to the transition of galaxies from the advection-dominated to the accretion-dominated regime. We also find that both the MZR and MZGR strongly suggest that low-mass galaxies preferentially lose metals entrained in their galactic winds. While this metal-enrichment of the galactic outflows is crucial for reproducing both the MZR and the MZGR at the low-mass end, we show that the flattening of gradients in massive galaxies is expected regardless of the nature of their winds., Comment: 3 figures, 1 appendix, 12 pages. Accepted by MNRAS; corrected typo in equation A1

We present a comprehensive analysis of the shape of dark matter (DM) halos in a sample of 25 Milky Way-like galaxies in TNG50 simulation. Using an enclosed volume iterative method, we infer an oblate-to-triaxial shape for the DM halo with median T ≃ 0.24. We group DM halos into three different categories. Simple halos (32% of the population) establish principal axes whose ordering in magnitude does not change with radius and whose orientations are almost fixed throughout the halo. Twisted halos (32%) experience levels of gradual rotations throughout their radial profiles. Finally, stretched halos (36%) demonstrate a stretching in the lengths of their principal axes where the ordering of different eigenvalues changes with radius. Subsequently, the halo experiences a “rotation” of ∼90° where the stretching occurs. Visualizing the 3D ellipsoid of each halo, for the first time, we report signs of a reorienting ellipsoid in twisted and stretched halos. We examine the impact of baryonic physics on DM halo shape through a comparison to dark matter only (DMO) simulations. This suggests a triaxial (prolate) halo. We analyze the impacts of substructure on DM halo shape in both hydrodynamical and DMO simulations and confirm that they are subdominant. We study the distribution of satellites in our sample. In simple and twisted halos, the angle between satellites’ angular momentum and the galaxy’s angular momentum grows with radius. However, stretched halos show a flat distribution of angles. Overlaying our theoretical outcome on the observational results presented in the literature establishes a fair agreement.

We investigate the morphology of the stellar distribution in a sample of Milky Way (MW) like galaxies in the TNG50 simulation. Using a local in shell iterative method (LSIM) as the main approach, we explicitly show evidence of twisting (in about 52% of halos) and stretching (in 48% of them) in the real space. This is matched with the re-orientation observed in the eigenvectors of the inertia tensor and gives us a clear picture of having a re-oriented stellar distribution. We make a comparison between the shape profile of dark matter (DM) halo and stellar distribution and quite remarkably see that their radial profiles are fairly close, especially at small galactocentric radii where the stellar disk is located. This implies that the DM halo is somewhat aligned with stars in response to the baryonic potential. The level of alignment mostly decreases away from the center. We study the impact of substructures in the orbital circularity parameter. It is demonstrated that in some cases, far away substructures are counter-rotating compared with the central stars and may flip the sign of total angular momentum and thus the orbital circularity parameter. Truncating them above 150 kpc, however, retains the disky structure of the galaxy as per initial selection. Including the impact of substructures in the shape of stars, we explicitly show that their contribution is subdominant. Overlaying our theoretical results to the observational constraints from previous literature, we establish fair agreement., 29 pages, 21 figures, Accepted for publication in ApJ

Physical and chemical properties of the interstellar medium (ISM) at sub-galactic ($\sim$kpc) scales play an indispensable role in controlling the ability of gas to form stars. As part of the SMAUG (Simulating Multiscale Astrophysics to Understand Galaxies) project, in this paper, we use the TNG50 cosmological simulation to explore the physical parameter space of 8 resolved ISM properties in star-forming regions to constrain the areas of this hyperspace over which most star-forming environments exist. We deconstruct our simulated galaxies spanning a wide range of mass (M$_\star = 10^{7-11}$ M$_\odot$) and redshift ($0 \leq z \leq 3$) into kpc-sized regions, and statistically analyze the gas/stellar surface densities, gas metallicity, vertical stellar velocity dispersion, epicyclic frequency and dark-matter volumetric density representative of each region in the context of their star formation activity and galactic environment (radial galactocentric location). By examining the star formation rate (SFR) weighted distributions of these properties, we show that stars primarily form in two spatially distinct environmental regimes, which are brought about by an underlying bi-component radial SFR surface density profile in galaxies. We examine how the relative prominence of these two regimes depends on host galaxy mass and cosmic time. We also compare our findings with those from integral field spectroscopy observations and achieve a good overall agreement. Further, using dimensionality reduction, we characterise the aforementioned hyperspace to reveal a high-degree of multicollinearity in relationships amongst ISM properties that drive the distribution of star formation at kpc-scales. Based on this, we show that a reduced 3D representation underpinned by a multi-variate radius relationship is sufficient to capture most of the variance in the original 8D space., Updated intro, minor text and abstract modifications. More information on the SMAUG project and the other first result papers can be found here: https://www.simonsfoundation.org/flatiron/center-for-computational-astrophysics/galaxy-formation/smaug/papersplash1

Understanding the variability of galaxy star formation histories (SFHs) across a range of timescales provides insight into the underlying physical processes that regulate star formation within galaxies. We compile the SFHs of galaxies at $z=0$ from an extensive set of models, ranging from cosmological hydrodynamical simulations (Illustris, IllustrisTNG, Mufasa, Simba, EAGLE), zoom simulations (FIRE-2, g14, and Marvel/Justice League), semi-analytic models (Santa Cruz SAM) and empirical models (UniverseMachine), and quantify the variability of these SFHs on different timescales using the power spectral density (PSD) formalism. We find that the PSDs are well described by broken power-laws, and variability on long timescales ($\gtrsim1$ Gyr) accounts for most of the power in galaxy SFHs. Most hydrodynamical models show increased variability on shorter timescales ($\lesssim300$ Myr) with decreasing stellar mass. Quenching can induce $\sim0.4-1$ dex of additional power on timescales $>1$ Gyr. The dark matter accretion histories of galaxies have remarkably self-similar PSDs and are coherent with the in-situ star formation on timescales $>3$ Gyr. There is considerable diversity among the different models in their (i) power due to SFR variability at a given timescale, (ii) amount of correlation with adjacent timescales (PSD slope), (iii) evolution of median PSDs with stellar mass, and (iv) presence and locations of breaks in the PSDs. The PSD framework is a useful space to study the SFHs of galaxies since model predictions vary widely. Observational constraints in this space will help constrain the relative strengths of the physical processes responsible for this variability., 31 pages, 17 figures (+ appendix). Resubmitted to MNRAS after responding to referee's comments. Comments are welcome!

Large scale outflows in star-forming galaxies are observed to be ubiquitous, and are a key aspect of theoretical modeling of galactic evolution in a cosmological context, the focus of the SMAUG (Simulating Multiscale Astrophysics to Understand Galaxies) project. Gas blown out from galactic disks, similar to gas within galaxies, consists of multiple phases with large contrasts of density, temperature, and other properties. To study multiphase outflows as emergent phenomena, we run a suite of ~pc-resolution local galactic disk simulations using the TIGRESS framework. Explicit modeling of the interstellar medium (ISM), including star formation and self-consistent radiative heating plus supernova feedback, regulates ISM properties and drives the outflow. We investigate the scaling of outflow mass, momentum, energy, and metal loading factors with galactic disk properties, including star formation rate (SFR) surface density (\Sigma_SFR~10^{-4}-1 M_sun/kpc^2/yr), gas surface density (~1-100 M_sun/pc^2), and total midplane pressure (or weight) (~10^3-10^6 k_B cm^{-3} K). The main components of outflowing gas are mass-delivering cool gas (T~10^4 K) and energy/metal-delivering hot gas (T~10^6 K). Cool mass outflow rates measured at outflow launch points (one or two scale heights) are 1-100 times the SFR (decreasing with \Sigma_SFR), although in massive galaxies most mass falls back due to insufficient outflow velocity. The hot galactic outflow carries mass comparable to 10% of the SFR, together with 10-20% of the energy and 30-60% of the metal mass injected by SN feedback. The characteristic outflow velocities of both phases scale very weakly with SFR, as v_out \propto \Sigma_SFR^{0.1~0.2}, consistent with observations. Importantly, our analysis demonstrates that in any physically-motivated cosmological wind model, it is crucial to include at least two distinct thermal wind components., Comment: Accepted for publication in ApJ. More information on the SMAUG project and the other first result papers can be found here: https://www.simonsfoundation.org/flatiron/center-for-computational-astrophysics/galaxy-formation/smaug/papersplash1

A key uncertainty in galaxy evolution is the physics regulating star formation, ranging from small-scale processes related to the life-cycle of molecular clouds within galaxies to large-scale processes such as gas accretion onto galaxies. We study the imprint of such processes on the time-variability of star formation with an analytical approach tracking the gas mass of galaxies ("regulator model"). Specifically, we quantify the strength of the fluctuation in the star-formation rate (SFR) on different timescales, i.e. the power spectral density (PSD) of the star-formation history, and connect it to gas inflow and the life-cycle of molecular clouds. We show that in the general case the PSD of the SFR has three breaks, corresponding to the correlation time of the inflow rate, the equilibrium timescale of the gas reservoir of the galaxy, and the average lifetime of individual molecular clouds. On long and intermediate timescales (relative to the dynamical timescale of the galaxy), the PSD is typically set by the variability of the inflow rate and the interplay between outflows and gas depletion. On short timescales, the PSD shows an additional component related to the life-cycle of molecular clouds, which can be described by a damped random walk with a power-law slope of $\beta\approx2$ at high frequencies with a break near the average cloud lifetime. We discuss star-formation "burstiness" in a wide range of galaxy regimes, study the evolution of galaxies about the main sequence ridgeline, and explore the applicability of our method for understanding the star-formation process on cloud-scale from galaxy-integrated measurements., Comment: Accepted to MNRAS, 26 pages with 12 figures (+ appendix)

We examine the properties of the circumgalactic medium (CGM) at low redshift in a range of simulated Milky Way mass halos. The sample is comprised of seven idealized simulations, an adaptive mesh refinement cosmological zoom-in simulation, and two groups of 50 halos with star forming or quiescent galaxies taken from the IllustrisTNG100 simulation. The simulations have very different setups, resolution, and feedback models, but are analyzed in a uniform manner. By comparing median radial profiles and mass distributions of CGM properties, we isolate key similarities and differences. In doing so, we advance the efforts of the SMAUG (Simulating Multiscale Astrophysics to Understand Galaxies) project that aims to understand the inherently multiscale galaxy formation process. In the cosmological simulations, the CGM exhibits nearly flat temperature distributions, and broad pressure and radial velocity distributions. In the idealized simulations, similar distributions are found in the inner CGM ($\lesssim 0.5 \, r_{\rm 200c}$) when strong galactic feedback models are employed, but the outer CGM ($\gtrsim 0.5 \, r_{\rm 200c}$) has a much less prominent cold phase, and narrower pressure and velocity distributions even in models with strong feedback. This comparative analysis demonstrates the dominant role feedback plays in shaping the inner CGM and the increased importance of cosmological effects, such as nonspherical accretion and satellite galaxies, in the outer CGM. Furthermore, our findings highlight that while cosmological simulations are required to capture the multiphase structure of the CGM at large radii, idealized simulations provide a robust framework to study how galactic feedback interacts with the inner CGM and thereby provide a reliable avenue to constrain feedback prescriptions., Comment: More information on the SMAUG project and the other first result papers can be found here: https://www.simonsfoundation.org/flatiron/center-for-computational-astrophysics/galaxy-formation/smaug/papersplash1

Semi-analytic models (SAMs) are a promising means of tracking the physical processes associated with galaxy formation, but many of their approximations have not been rigorously tested. As part of the SMAUG (Simulating Multiscale Astrophysics to Understand Galaxies) project, we compare predictions from the FIRE-2 hydrodynamical "zoom-in" simulations to those from the Santa Cruz SAM run on the same halo merger trees, with an emphasis on the global mass flow cycle. Our study includes 13 halos spanning low-mass dwarfs (M_vir~10^10 M_sun at z=0), intermediate-mass dwarfs (M_vir~10^11 M_sun) and Milky Way-mass galaxies (M_vir~10^12 M_sun). The SAM and FIRE-2 predictions agree relatively well with each other in terms of stellar and interstellar mass, but differ dramatically on circumgalactic mass (the SAM is lower than FIRE-2 by ~3 orders of magnitude for dwarfs). Strikingly, the SAM predicts higher gas accretion rates for dwarfs compared to FIRE-2 by factors of ~10-100, and this is compensated for with higher mass outflow rates in the SAM. We argue that the most severe model discrepancies are caused by the lack of preventative stellar feedback and the assumptions for halo gas cooling and recycling in the SAM. As a first step towards resolving these model tensions, we present a simple yet promising new preventative stellar feedback model in which the energy carried by supernova-driven winds is allowed to heat some fraction of gas outside of halos to at least the virial temperature such that accretion is suppressed., Comment: Accepted to ApJ with minor revisions. More information on the SMAUG project and the other first result papers can be found here: https://www.simonsfoundation.org/flatiron/center-for-computational-astrophysics/galaxy-formation/smaug/papersplash1

Galactic outflows have density, temperature, and velocity variations at least as large as that of the multiphase, turbulent interstellar medium (ISM) from which they originate. We have conducted a suite of parsec-resolution numerical simulations using the TIGRESS framework, in which outflows emerge as a consequence of interaction between supernovae (SNe) and the star-forming ISM. The outflowing gas is characterized by two distinct thermal phases, cool (T10^6 K), with most mass carried by the cool phase and most energy and newly-injected metals carried by the hot phase. Both components have a broad distribution of outflow velocity, and especially for cool gas this implies a varying fraction of escaping material depending on the halo potential. Informed by the TIGRESS results, we develop straightforward analytic formulae for the joint probability density functions (PDFs) of mass, momentum, energy, and metal loading as distributions in outflow velocity and sound speed. The model PDFs have only two parameters, SFR surface density \Sigma_SFR and the metallicity of the ISM, and fully capture the behavior of the original TIGRESS simulation PDFs over \Sigma_SFR~(10^{-4},1)M_sun/kpc^2/yr. Employing PDFs from resolved simulations will enable galaxy formation subgrid model implementations with wind velocity and temperature (as well as total loading factors) that are based on theoretical predictions rather than empirical tuning. This is a critical step to incorporate advances from TIGRESS and other high-resolution simulations in future cosmological hydrodynamics and semi-analytic galaxy formation models. We release a python package to prototype our model and to ease its implementation., Comment: ApJL accepted. For associated python package, see https://twind.readthedocs.io/en/latest/

Cold clouds in hot media are quickly crushed, shredded, and then accelerated as a result of their interaction with the background gas. The persistence of cold clouds moving at substantial velocities in harsh environments is a common yet puzzling feature of many astrophysical systems, from quasar absorption lines probing galactic halos to clouds of dust passing near Sgr $A^*$. Here we run a set of idealized numerical experiments, subjecting a line of cold clouds at a series of mutual separations to a hot background wind. We find that this stream of clouds is able to shield itself from hydrodynamic destruction by accelerating the hot background material, creating a protective layer of co-moving gas. We write down a simple diffusion equation that reproduces the behavior of the simulations, and we discuss the implications for cosmological gas accretion and G2., Submitted to AAS Journals, comments welcome

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

Black holes growing via the accretion of gas emit radiation that can photoevaporate the atmospheres of nearby planets. Here we couple planetary structural evolution models of sub-Neptune mass planets to the growth of the Milky way's central supermassive black-hole, Sgr A$^*$ and investigate how planetary evolution is influenced by quasar activity. We find that, out to ${\sim} 20$ pc from Sgr A$^*$, the XUV flux emitted during its quasar phase can remove several percent of a planet's H/He envelope by mass; in many cases, this removal results in bare rocky cores, many of which situated in the habitable zones (HZs) of G-type stars. The erosion of sub-Neptune sized planets may be one of the most prevalent channels by which terrestrial super-Earths are created near the Galactic Center. As such, the planet population demographics may be quite different close to Sgr A$^*$ than in the Galaxy's outskirts. The high stellar densities in this region (about seven orders of magnitude greater than the solar neighborhood) imply that the distance between neighboring rocky worlds is short ($500-5000$~AU). The proximity between potentially habitable terrestrial planets may enable the onset of widespread interstellar panspermia near the nuclei of galaxies. More generally, we predict these phenomena to be ubiquitous for planets in nuclear star clusters and ultra-compact dwarfs. Globular clusters, on the other hand, are less affected by the black holes., Comment: 9 pages, 4 figures, accepted for publication in Astrophysical Journal Letters

Planetary atmospheres are subject to mass loss through a variety of mechanisms including irradiation by XUV photons from their host star. Here we explore the consequences of XUV irradiation by supermassive black holes as they grow by the accretion of gas in galactic nuclei. Based on the mass distribution of stars in galactic bulges and disks and the luminosity history of individual black holes, we estimate the probability distribution function of XUV fluences as a function of galaxy halo mass, redshift, and stellar component. We find that about 50% of all planets in the universe may lose the equivalent of a Martian atmosphere, 10% may lose an Earth's atmosphere, and 0.2% may lose the mass of Earth's oceans. The fractions are appreciably higher in the spheroidal components of galaxies, and depend strongly on galaxy mass, but only weakly on redshift., Comment: Submitted to MNRAS

Almost by definition brown dwarfs are objects with masses below the hydrogen burning limit, around $0.07\ M_\odot$. Below this mass, objects never reach a steady state where they can fuse hydrogen. Here we demonstrate, in contrast to this traditional view, that brown dwarfs with masses greater than the hydrogen burning limit may in principle exist in the universe. These objects, which we term "overmassive brown dwarfs" form a continuous sequence with traditional brown dwarfs in any property (mass, effective temperature, radius, luminosity). To form an overmassive brown dwarf, mass must be added sufficiently slowly to a sufficiently old traditional brown dwarf below the hydrogen burning limit. We identify mass transfer in binary brown dwarf systems via Roche lobe overflow driven by gravitational waves to be the most plausible mechanism to produce the bulk of the putative overmassive brown dwarf population., Comment: Resubmitted to AAS Journals after revisions

We present the stellar mass ($M_{*}$)--gas-phase metallicity relation (MZR) and its scatter at intermediate redshifts ($0.5\leq z\leq0.7$) for 1381 field galaxies collected from deep spectroscopic surveys. The star formation rate (SFR) and color at a given $M_{*}$ of this magnitude-limited ($R\lesssim24$ AB) sample are representative of normal star-forming galaxies. For masses below $10^9 M_\odot$, our sample of 237 galaxies is $\sim$10 times larger than those in previous studies beyond the local universe. This huge gain in sample size enables superior constraints on the MZR and its scatter in the low-mass regime. We find a power-law MZR at $10^{8} M_\odot < M_{*} < 10^{11} M_\odot$: ${12+log(O/H) = (5.83\pm0.19) + (0.30\pm0.02)log(M_{*}/M_\odot)}$. Our MZR shows good agreement with others measured at similar redshifts in the literature in the intermediate and massive regimes, but is shallower than the extrapolation of the MZRs of others to masses below $10^{9} M_\odot$. The SFR dependence of the MZR in our sample is weaker than that found for local galaxies (known as the Fundamental Metallicity Relation). Compared to a variety of theoretical models, the slope of our MZR for low-mass galaxies agrees well with predictions incorporating supernova energy-driven winds. Being robust against currently uncertain metallicity calibrations, the scatter of the MZR serves as a powerful diagnostic of the stochastic history of gas accretion, gas recycling, and star formation of low-mass galaxies. Our major result is that the scatter of our MZR increases as $M_{*}$ decreases. Our result implies that either the scatter of the baryonic accretion rate or the scatter of the $M_{*}$--$M_{halo}$ relation increases as $M_{*}$ decreases. Moreover, our measures of scatter at $z=0.7$ appears consistent with that found for local galaxies., 18 pages, 10 figures. Accepted by ApJ. Typos corrected

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.

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., Comment: 18 pages, 9 figures. Accepted for publication in ApJ. 6.5 TB of simulation data and processed derived data are available at http://dx.doi.org/10.13012/J8F769GV

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., Comment: 28 pages, 17 figures, submitted to MNRAS (comments welcome) - conspicuous typos corrected

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., Comment: 16 pages, 13 figures. Accepted for publication in the Astrophysical Journal. 6 TB of simulation data and processed derived data are available at http://dx.doi.org/10.13012/J85Q4T1T and an interactive data exploration widget is available at https://demo.use.yt (click "galaxy_visualization.ipynb"; once the Jupyter notebook launches, click "Cell -> Run All")

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., Comment: 18 pages, 6 figures, accepted for publication in ApJ. Code available at http://www.ucolick.org/~jforbes/gidget.html