Showing 1–50 of 52 results for author: Forbes, J C

Why do semi-analytic models predict higher scatter in the stellar mass-halo mass relation than cosmological hydrodynamic simulations?

Authors: Antonio J. Porras-Valverde, John C. Forbes, Rachel S. Somerville, Adam R. H. Stevens, Kelly Holley-Bockelmann, Andreas A. Berlind, Shy Genel

Abstract: Semi-analytic models (SAMs) systematically predict higher stellar-mass scatter at a given halo mass than hydrodynamical simulations and most empirical models. Our goal is to investigate the physical origin of this scatter by exploring modifications to the physics in the SAM Dark Sage. We design two black hole formation models that approximate results from the IllustrisTNG 300-1 hydrodynamical simu… ▽ More Semi-analytic models (SAMs) systematically predict higher stellar-mass scatter at a given halo mass than hydrodynamical simulations and most empirical models. Our goal is to investigate the physical origin of this scatter by exploring modifications to the physics in the SAM Dark Sage. We design two black hole formation models that approximate results from the IllustrisTNG 300-1 hydrodynamical simulation. In the first model, we assign a fixed black hole mass of $10^{6}\, \mathrm{M}_{\odot}$ to every halo that reaches $10^{10.5}\, \mathrm{M}_{\odot}$. In the second model, we disregard any black hole growth as implemented in the standard Dark Sage model. Instead, we force all black hole masses to follow the median black hole mass-halo mass relation in IllustrisTNG 300-1 with a fixed scatter. We find that each model on its own does not significantly reduce the scatter in stellar mass. To do this, we replace the native Dark Sage AGN feedback model with a simple model where we turn off cooling for galaxies with black hole masses above $10^{8}\, \mathrm{M}_{\odot}$. With this additional modification, the SMBH seeding and fixed conditional distribution models find a significant reduction in the scatter in stellar mass at halo masses between $10^{11-14}\, \mathrm{M}_{\odot}$. These results suggest that AGN feedback in SAMs acts in a qualitatively different way than feedback implemented in cosmological simulations. Either or both may require substantial modification to match the empirically inferred scatter in the Stellar Mass Halo Mass Relation (SMHMR). △ Less

Submitted 17 October, 2023; originally announced October 2023.

Comments: 21 pages, 16 figures

The Global Structure of Molecular Clouds: I. Trends with Mass and Star Formation Rate

Authors: Nia Imara, John C. Forbes

Abstract: We introduce a model for the large-scale, global 3D structure of molecular clouds. Motivated by the morphological appearance of clouds in surface density maps, we model clouds as cylinders, with the aim of backing out information about the volume density distribution of gas and its relationship to star formation. We test our model by applying it to surface density maps for a sample of nearby cloud… ▽ More We introduce a model for the large-scale, global 3D structure of molecular clouds. Motivated by the morphological appearance of clouds in surface density maps, we model clouds as cylinders, with the aim of backing out information about the volume density distribution of gas and its relationship to star formation. We test our model by applying it to surface density maps for a sample of nearby clouds and find solutions that fit each of the observed radial surface density profiles remarkably well. Our most salient findings are that clouds with higher central volume densities are more compact and also have lower total mass. These same lower-mass clouds tend to have shorter gas depletion times, regardless of whether we consider their total mass or dense mass. Our analyses lead us to conclude that cylindrical clouds can be characterized by a universal structure that sets the timescale on which they form stars. △ Less

Submitted 8 September, 2023; originally announced September 2023.

Comments: Accepted to ApJ. 12 pages, 9 figures

The Galactic Interstellar Object Population: A Framework for Prediction and Inference

Authors: Matthew J. Hopkins, Chris Lintott, Michele T. Bannister, J. Ted Mackereth, John C. Forbes

Abstract: The Milky Way is thought to host a huge population of interstellar objects (ISOs), numbering approximately $10^{15}\mathrm{pc}^{-3}$ around the Sun, which are formed and shaped by a diverse set of processes ranging from planet formation to galactic dynamics. We define a novel framework: firstly to predict the properties of this Galactic ISO population by combining models of processes across planet… ▽ More The Milky Way is thought to host a huge population of interstellar objects (ISOs), numbering approximately $10^{15}\mathrm{pc}^{-3}$ around the Sun, which are formed and shaped by a diverse set of processes ranging from planet formation to galactic dynamics. We define a novel framework: firstly to predict the properties of this Galactic ISO population by combining models of processes across planetary and galactic scales, and secondly to make inferences about the processes modelled, by comparing the predicted population to what is observed. We predict the spatial and compositional distribution of the Galaxy's population of ISOs by modelling the Galactic stellar population with data from the APOGEE survey and combining this with a protoplanetary disk chemistry model. Selecting ISO water mass fraction as an example observable quantity, we evaluate its distribution both at the position of the Sun and averaged over the Galactic disk; our prediction for the Solar neighbourhood is compatible with the inferred water mass fraction of 2I/Borisov. We show that the well-studied Galactic stellar metallicity gradient has a corresponding ISO compositional gradient. We also demonstrate the inference part of the framework by using the current observed ISO composition distribution to constrain the parent star metallicity dependence of the ISO production rate. This constraint, and other inferences made with this framework, will improve dramatically as the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) progresses and more ISOs are observed. Finally, we explore generalisations of this framework to other Galactic populations, such as that of exoplanets. △ Less

Submitted 15 November, 2023; v1 submitted 10 August, 2023; originally announced August 2023.

Comments: Accepted to AJ

Journal ref: AJ 166 241 (2023)

The interplay between feedback, accretion, transport and winds in setting gas-phase metal distribution in galaxies

Authors: Piyush Sharda, Omri Ginzburg, Mark R. Krumholz, John C. Forbes, Emily Wisnioski, Matilde Mingozzi, Henry R. M. Zovaro, Avishai Dekel

Abstract: The recent decade has seen an exponential growth in spatially-resolved metallicity measurements in the interstellar medium (ISM) of galaxies. To first order, these measurements are characterised by the slope of the radial metallicity profile, known as the metallicity gradient. In this work, we model the relative role of star formation feedback, gas transport, cosmic gas accretion, and galactic win… ▽ More The recent decade has seen an exponential growth in spatially-resolved metallicity measurements in the interstellar medium (ISM) of galaxies. To first order, these measurements are characterised by the slope of the radial metallicity profile, known as the metallicity gradient. In this work, we model the relative role of star formation feedback, gas transport, cosmic gas accretion, and galactic winds in driving radial metallicity profiles and setting the mass-metallicity gradient relation (MZGR). We include a comprehensive treatment of these processes by including them as sources that supply mass, metals, and energy to marginally unstable galactic discs in pressure and energy balance. We show that both feedback and accretion that can drive turbulence and enhance metal-mixing via diffusion are crucial to reproduce the observed MZGR in local galaxies. Metal transport also contributes to setting metallicity profiles, but it is sensitive to the strength of radial gas flows in galaxies. While the mass loading of galactic winds is important to reproduce the mass metallicity relation (MZR), we find that metal mass loading is more important to reproducing the MZGR. Specifically, our model predicts preferential metal enrichment of galactic winds in low-mass galaxies. This conclusion is robust against our adopted scaling of the wind mass-loading factor, uncertainties in measured wind metallicities, and systematics due to metallicity calibrations. Overall, we find that at $z \sim 0$, galactic winds and metal transport are more important in setting metallicity gradients in low-mass galaxies whereas star formation feedback and gas accretion dominate setting metallicity gradients in massive galaxies. △ Less

Submitted 11 January, 2024; v1 submitted 28 March, 2023; originally announced March 2023.

Comments: 21 pages (+ appendix), 11 figures, 4 tables. Accepted by MNRAS

Characterizing the Conditional Galaxy Property Distribution using Gaussian Mixture Models

Authors: Yucheng Zhang, Anthony R. Pullen, Rachel S. Somerville, Patrick C. Breysse, John C. Forbes, Shengqi Yang, Yin Li, Abhishek S. Maniyar

Abstract: Line-intensity mapping (LIM) is a promising technique to constrain the global distribution of galaxy properties. To combine LIM experiments probing different tracers with traditional galaxy surveys and fully exploit the scientific potential of these observations, it is necessary to have a physically motivated modeling framework. As part of developing such a framework, in this work we introduce and… ▽ More Line-intensity mapping (LIM) is a promising technique to constrain the global distribution of galaxy properties. To combine LIM experiments probing different tracers with traditional galaxy surveys and fully exploit the scientific potential of these observations, it is necessary to have a physically motivated modeling framework. As part of developing such a framework, in this work we introduce and model the conditional galaxy property distribution (CGPD), i.e. the distribution of galaxy properties conditioned on the host halo mass and redshift. We consider five galaxy properties, including the galaxy stellar mass, molecular gas mass, galaxy radius, gas phase metallicity and star formation rate (SFR), which are important for predicting the emission lines of interest. The CGPD represents the full distribution of galaxies in the five dimensional property space; many important galaxy distribution functions and scaling relations, such as the stellar mass function and SFR main sequence, can be derived from integrating and projecting it. We utilize two different kinds of cosmological galaxy simulations, a semi-analytic model and the IllustrisTNG hydrodynamic simulation, to characterize the CGPD and explore how well it can be represented using a Gaussian mixture model (GMM). We find that with just a few ($\sim 3$) Gaussian components, a GMM can describe the CGPD of the simulated galaxies to high accuracy for both simulations. The CGPD can be mapped to LIM or other observables by constructing the appropriate relationship between galaxy properties and the relevant observable tracers. △ Less

Submitted 22 February, 2023; originally announced February 2023.

Comments: 17 pages, 10 figures

How Low Can Q Go?

Authors: John C. Forbes

Abstract: Gravitational instability plays a substantial role in the evolution of galaxies. Various schemes to include it in galaxy evolution models exist, generally assuming that the Toomre $Q$ parameter is self-regulated to $Q_\mathrm{crit}$, the critical $Q$ dividing stable from unstable conditions in a linear stability analysis. This assumption is in tension with observational estimates of $Q$ that find… ▽ More Gravitational instability plays a substantial role in the evolution of galaxies. Various schemes to include it in galaxy evolution models exist, generally assuming that the Toomre $Q$ parameter is self-regulated to $Q_\mathrm{crit}$, the critical $Q$ dividing stable from unstable conditions in a linear stability analysis. This assumption is in tension with observational estimates of $Q$ that find values far below any plausible value of $Q_\mathrm{crit}$. While the observations are subject to some uncertainty, this tension can more easily be relieved on the theoretical side by relaxing the common assumption that $Q\ge Q_\mathrm{crit}$. Based on observations of both $z\sim 2$ disks and local face-on galaxies, we estimate the effect of gravitational instability necessary to balance out every other physical process that affects $Q$. In particular we find that the disk's response to low $Q$ values can be described by simple functions that depend only on $Q$. These response functions allow galaxies to maintain $Q$ values below $Q_\mathrm{crit}$ in equilibrium over a wide range of parameters. Extremely low values of $Q$ are predicted when the gas surface density is greater than $\sim 10^3$ M$_\odot$ pc$^{-2}$, the rotation curve provides minimal shear, the orbital time becomes long, and/or when the gas is much more unstable than the stellar component. We suggest that these response functions should be used in place of the $Q\ge Q_\mathrm{crit}$ ansatz. △ Less

Submitted 15 February, 2023; originally announced February 2023.

Comments: Submitted to AAS Journals

The Neon Gap: Probing Ionization with Dwarf Galaxies at z~1

Authors: John Pharo, Yicheng Guo, David C. Koo, John C. Forbes, Puragra Guhathakurta

Abstract: We present measurements of [NeIII]λ3869 emission in z~1 low-mass galaxies taken from the Keck/DEIMOS spectroscopic surveys HALO7D and DEEPWinds. We identify 167 individual galaxies with significant [NeIII] emission lines, including 112 "dwarf" galaxies with log(M_{\star}/M_{\odot}) < 9.5, with 0.3 < z < 1.4. We also measure [NeIII] emission from composite spectra derived from all [OII]λλ3727,3729… ▽ More We present measurements of [NeIII]λ3869 emission in z~1 low-mass galaxies taken from the Keck/DEIMOS spectroscopic surveys HALO7D and DEEPWinds. We identify 167 individual galaxies with significant [NeIII] emission lines, including 112 "dwarf" galaxies with log(M_{\star}/M_{\odot}) < 9.5, with 0.3 < z < 1.4. We also measure [NeIII] emission from composite spectra derived from all [OII]λλ3727,3729 line emitters in this range. This provides a unique sample of [NeIII]-emitters in the gap between well-studied emitters at z = 0 and 2 < z < 3. To study evolution in ionization conditions in the ISM over this time, we analyze the log([NeIII]λ3869/[OII]λλ3727,3729) ratio (Ne3O2) as a function of the stellar mass and of the log([OIII]λλ4959,5007/[OII]λλ3727,3729) ratio (O32). We find that the typical star-forming dwarf galaxy at this redshift, as measured from the composite spectra, shares the Ne3O2-M_{\star} relation with local galaxies, but have higher O32 at given Ne3O2. This finding implies that the ionization and metallicity characteristics of the z~1 dwarf population do not evolve substantially from z~1 to z=0, suggesting that the known evolution in those parameter from z~2 has largely taken place by z~1. Individual [NeIII]-detected galaxies have emission characteristics situated between local and z~2 galaxies, with elevated Ne3O2 and O32 emission potentially explained by variations in stellar and nebular metallicity. We also compare our dwarf sample to similarly low-mass z > 7 galaxies identified in JWST Early Release Observations, finding four HALO7D dwarfs with similar size, metallicity, and star formation properties. △ Less

Submitted 18 January, 2023; originally announced January 2023.

Comments: Accepted to ApJL, 16.01.2023

Arkenstone I: A Novel Method for Robustly Capturing High Specific Energy Outflows In Cosmological Simulations

Authors: Matthew C. Smith, Drummond B. Fielding, Greg L. Bryan, Chang-Goo Kim, Eve C. Ostriker, Rachel S. Somerville, Jonathan Stern, Kung-Yi Su, Rainer Weinberger, Chia-Yu Hu, John C. Forbes, Lars Hernquist, Blakesley Burkhart, Yuan Li

Abstract: Arkenstone is a new model for multiphase, stellar feedback driven galactic winds designed for inclusion in coarse resolution cosmological simulations. In this first paper of a series, we describe the features that allow Arkenstone to properly treat high specific energy wind components and demonstrate them using idealised non-cosmological simulations of a galaxy with a realistic CGM, using the Arep… ▽ More Arkenstone is a new model for multiphase, stellar feedback driven galactic winds designed for inclusion in coarse resolution cosmological simulations. In this first paper of a series, we describe the features that allow Arkenstone to properly treat high specific energy wind components and demonstrate them using idealised non-cosmological simulations of a galaxy with a realistic CGM, using the Arepo code. Hot, fast gas phases with low mass loadings are predicted to dominate the energy content of multiphase outflows. In order to treat the huge dynamic range of spatial scales involved in cosmological galaxy formation at feasible computational expense, cosmological volume simulations typically employ a Lagrangian code or else use adaptive mesh refinement with a quasi-Lagrangian refinement strategy. However, it is difficult to inject a high specific energy wind in a Lagrangian scheme without incurring artificial burstiness. Additionally, the low densities inherent to this type of flow result in poor spatial resolution. Arkenstone addresses these issues with a novel scheme for coupling energy into the ISM/CGM transition region which also provides the necessary level of refinement at the base of the wind. In the absence of our improvements, we show that poor spatial resolution near the sonic point of a hot, fast outflow leads to an underestimation of gas acceleration as the wind propagates. We explore the different mechanisms by which low and high specific energy winds can regulate the SFR of galaxies. In future work, we will demonstrate other aspects of the Arkenstone model. △ Less

Submitted 7 March, 2024; v1 submitted 17 January, 2023; originally announced January 2023.

Comments: Published MNRAS, 27 pages, 17 figures; updated for consistency with journal version

Journal ref: MNRAS, 527, 1216 (2024)

A 3D View of Orion: I. Barnard's Loop

Authors: Michael M. Foley, Alyssa Goodman, Catherine Zucker, John C. Forbes, Ralf Konietzka, Cameren Swiggum, João Alves, John Bally, Juan D. Soler, Josefa E. Großschedl, Shmuel Bialy, Michael Y. Grudić, Reimar Leike, Torsten Ensslin

Abstract: Barnard's Loop is a famous arc of H$α$ emission located in the Orion star-forming region. Here, we provide evidence of a possible formation mechanism for Barnard's Loop and compare our results with recent work suggesting a major feedback event occurred in the region around 6 Myr ago. We present a 3D model of the large-scale Orion region, indicating coherent, radial, 3D expansion of the OBP-Near/Br… ▽ More Barnard's Loop is a famous arc of H$α$ emission located in the Orion star-forming region. Here, we provide evidence of a possible formation mechanism for Barnard's Loop and compare our results with recent work suggesting a major feedback event occurred in the region around 6 Myr ago. We present a 3D model of the large-scale Orion region, indicating coherent, radial, 3D expansion of the OBP-Near/Briceño-1 (OBP-B1) cluster in the middle of a large dust cavity. The large-scale gas in the region also appears to be expanding from a central point, originally proposed to be Orion X. OBP-B1 appears to serve as another possible center, and we evaluate whether Orion X or OBP-B1 is more likely to be the cause of the expansion. We find that neither cluster served as the single expansion center, but rather a combination of feedback from both likely propelled the expansion. Recent 3D dust maps are used to characterize the 3D topology of the entire region, which shows Barnard's Loop's correspondence with a large dust cavity around the OPB-B1 cluster. The molecular clouds Orion A, Orion B, and Orion $λ$ reside on the shell of this cavity. Simple estimates of gravitational effects from both stars and gas indicate that the expansion of this asymmetric cavity likely induced anisotropy in the kinematics of OBP-B1. We conclude that feedback from OBP-B1 has affected the structure of the Orion A, Orion B, and Orion $λ$ molecular clouds and may have played a major role in the formation of Barnard's Loop. △ Less

Submitted 2 December, 2022; originally announced December 2022.

Comments: 25 pages, 10 figures. Submitted to ApJ

A unified model for the co-evolution of galaxies and their circumgalactic medium: the relative roles of turbulence and atomic cooling physics

Authors: Viraj Pandya, Drummond B. Fielding, Greg L. Bryan, Christopher Carr, Rachel S. Somerville, Jonathan Stern, Claude-Andre Faucher-Giguere, Zachary Hafen, Daniel Angles-Alcazar, John C. Forbes

Abstract: The circumgalactic medium (CGM) plays a pivotal role in regulating gas flows around galaxies and thus shapes their evolution. However, the details of how galaxies and their CGM co-evolve remain poorly understood. We present a new time-dependent two-zone model that self-consistently tracks not just mass and metal flows between galaxies and their CGM but also the evolution of the global thermal and… ▽ More The circumgalactic medium (CGM) plays a pivotal role in regulating gas flows around galaxies and thus shapes their evolution. However, the details of how galaxies and their CGM co-evolve remain poorly understood. We present a new time-dependent two-zone model that self-consistently tracks not just mass and metal flows between galaxies and their CGM but also the evolution of the global thermal and turbulent kinetic energy of the CGM. Our model accounts for heating and turbulence driven by both supernova winds and cosmic accretion as well as radiative cooling, turbulence dissipation, and halo outflows due to CGM overpressurization. We demonstrate that, depending on parameters, the CGM can undergo a phase transition (``thermalization'') from a cool, turbulence-supported phase to a virial-temperature, thermally-supported phase. This CGM phase transition is largely determined by the ability of radiative cooling to balance heating from supernova winds and turbulence dissipation. We perform an initial calibration of our model to the FIRE-2 cosmological hydrodynamical simulations and show that it can approximately reproduce the baryon cycles of the simulated halos. In particular, we find that, for these parameters, the phase transition occurs at high-redshift in ultrafaint progenitors and at low redshift in classical $M_{\rm vir}\sim10^{11}M_{\odot}$ dwarfs, while Milky Way-mass halos undergo the transition at $z\approx0.5$. We see a similar transition in the simulations though it is more gradual, likely reflecting radial dependence and multi-phase gas not captured by our model. We discuss these and other limitations of the model and possible future extensions. △ Less

Submitted 27 August, 2023; v1 submitted 17 November, 2022; originally announced November 2022.

Comments: Updated to reflect version accepted by ApJ; minor text edits to clarify limitations and caveats of model

Gas Morphology of Milky Way-like Galaxies in the TNG50 Simulation: Signals of Twisting and Stretching

Authors: Thomas K. Waters, Colton Peterson, Razieh Emami, Xuejian Shen, Lars Hernquist, Randall Smith, Mark Vogelsberger, Charles Alcock, Grant Tremblay, Matthew Liska, John C. Forbes, Jorge Moreno

Abstract: We present an in-depth analysis of gas morphologies for a sample of 25 Milky Way-like galaxies from the IllustrisTNG TNG50 simulation. We constrain the morphology of cold, warm, hot gas, and gas particles as a whole using a Local Shell Iterative Method (LSIM) and explore its observational implications by computing the hard-to-soft X-ray ratio, which ranges between $10^{-3}$-$10^{-2}$ in the inner… ▽ More We present an in-depth analysis of gas morphologies for a sample of 25 Milky Way-like galaxies from the IllustrisTNG TNG50 simulation. We constrain the morphology of cold, warm, hot gas, and gas particles as a whole using a Local Shell Iterative Method (LSIM) and explore its observational implications by computing the hard-to-soft X-ray ratio, which ranges between $10^{-3}$-$10^{-2}$ in the inner $\sim 50 \rm kpc$ of the distribution and $10^{-5}$-$10^{-4}$ at the outer portion of the hot gas distribution. We group galaxies into three main categories: simple, stretched, and twisted. These categories are based on the radial reorientation of the principal axes of the reduced inertia tensor. We find that a vast majority ($77\%$) of the galaxies in our sample exhibit twisting patterns in their radial profiles. Additionally, we present detailed comparisons between 1) the gaseous distributions belonging to individual temperature regimes, 2) the cold gas distributions and stellar distributions, and 3) the gaseous distributions and dark matter (DM) halos. We find a strong correlation between the morphological properties of the cold gas and stellar distributions. Furthermore, we find a correlation between gaseous distributions with DM halo that increases with gas temperature, implying that we may use the warm-hot gaseous morphology as a tracer to probe the DM morphology. Finally, we show gaseous distributions exhibit significantly more prolate morphologies than the stellar distributions and DM halos, which we hypothesize is due to stellar and AGN feedback. △ Less

Submitted 5 January, 2024; v1 submitted 3 October, 2022; originally announced October 2022.

Comments: 26 pages, 12 figures

Code Comparison in Galaxy Scale Simulations with Resolved Supernova Feedback: Lagrangian vs. Eulerian Methods

Authors: Chia-Yu Hu, Matthew C. Smith, Romain Teyssier, Greg L. Bryan, Robbert Verbeke, Andrew Emerick, Rachel S. Somerville, Blakesley Burkhart, Yuan Li, John C. Forbes, Tjitske Starkenburg

Abstract: We present a suite of high-resolution simulations of an isolated dwarf galaxy using four different hydrodynamical codes: {\sc Gizmo}, {\sc Arepo}, {\sc Gadget}, and {\sc Ramses}. All codes adopt the same physical model which includes radiative cooling, photoelectric heating, star formation, and supernova (SN) feedback. Individual SN explosions are directly resolved without resorting to sub-grid mo… ▽ More We present a suite of high-resolution simulations of an isolated dwarf galaxy using four different hydrodynamical codes: {\sc Gizmo}, {\sc Arepo}, {\sc Gadget}, and {\sc Ramses}. All codes adopt the same physical model which includes radiative cooling, photoelectric heating, star formation, and supernova (SN) feedback. Individual SN explosions are directly resolved without resorting to sub-grid models, eliminating one of the major uncertainties in cosmological simulations. We find reasonable agreement on the time-averaged star formation rates as well as the joint density-temperature distributions between all codes. However, the Lagrangian codes show significantly burstier star formation, larger supernova-driven bubbles, and stronger galactic outflows compared to the Eulerian code. This is caused by the behavior in the dense, collapsing gas clouds when the Jeans length becomes unresolved: gas in Lagrangian codes collapses to much higher densities than in Eulerian codes, as the latter is stabilized by the minimal cell size. Therefore, more of the gas cloud is converted to stars and SNe are much more clustered in the Lagrangian models, amplifying their dynamical impact. The differences between Lagrangian and Eulerian codes can be reduced by adopting a higher star formation efficiency in Eulerian codes, which significantly enhances SN clustering in the latter. Adopting a zero SN delay time reduces burstiness in all codes, resulting in vanishing outflows as SN clustering is suppressed. △ Less

Submitted 3 May, 2023; v1 submitted 22 August, 2022; originally announced August 2022.

Comments: accepted version by ApJ (including a new simulation in Appendix B suggested by the referee)

Stochastic Modelling of Star Formation Histories III. Constraints from Physically-Motivated Gaussian Processes

Authors: Kartheik G. Iyer, Joshua S. Speagle, Neven Caplar, John C. Forbes, Eric Gawiser, Joel Leja, Sandro Tacchella

Abstract: Galaxy formation and evolution involves a variety of effectively stochastic processes that operate over different timescales. The Extended Regulator model provides an analytic framework for the resulting variability (or `burstiness') in galaxy-wide star formation due to these processes. It does this by relating the variability in Fourier space to the effective timescales of stochastic gas inflow,… ▽ More Galaxy formation and evolution involves a variety of effectively stochastic processes that operate over different timescales. The Extended Regulator model provides an analytic framework for the resulting variability (or `burstiness') in galaxy-wide star formation due to these processes. It does this by relating the variability in Fourier space to the effective timescales of stochastic gas inflow, equilibrium, and dynamical processes influencing GMC creation and destruction using the power spectral density (PSD) formalism. We use the connection between the PSD and auto-covariance function (ACF) for general stochastic processes to reformulate this model as an auto-covariance function, which we use to model variability in galaxy star formation histories (SFHs) using physically-motivated Gaussian Processes in log SFR space. Using stellar population synthesis models, we then explore how changes in model stochasticity can affect spectral signatures across galaxy populations with properties similar to the Milky Way and present-day dwarfs as well as at higher redshifts. We find that, even at fixed scatter, perturbations to the stochasticity model (changing timescales vs overall variability) leave unique spectral signatures across both idealized and more realistic galaxy populations. Distributions of spectral features including H$α$ and UV-based SFR indicators, H$δ$ and Ca-H,K absorption line strengths, D$_n$(4000) and broadband colors provide testable predictions for galaxy populations from present and upcoming surveys with Hubble, Webb \& Roman. The Gaussian process SFH framework provides a fast, flexible implementation of physical covariance models for the next generation of SED modeling tools. Code to reproduce our results can be found at https://github.com/kartheikiyer/GP-SFH △ Less

Submitted 11 August, 2022; originally announced August 2022.

Comments: 33 pages, 18 figures. Submitted to ApJ, comments are welcome! Code available at https://github.com/kartheikiyer/GP-SFH

Empirical Dust Attenuation Model Leads to More Realistic UVJ Diagram for TNG100 Galaxies

Authors: Gautam Nagaraj, John C. Forbes, Joel Leja, Dan Foreman-Mackey, Christopher C. Hayward

Abstract: 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 g… ▽ More 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. △ Less

Submitted 13 April, 2022; originally announced April 2022.

Comments: 12 pages, 5 figures, submitted to ApJ

Gas Accretion Can Drive Turbulence in Galaxies

Authors: John C. Forbes, Razieh Emami, Rachel S. Somerville, Shy Genel, Dylan Nelson, Annalisa Pillepich, Blakesley Burkhart, Greg L. Bryan, Mark R. Krumholz, Lars Hernquist, Stephanie Tonnesen, Paul Torrey, Viraj Pandya, Christopher C. Hayward

Abstract: 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… ▽ More 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$. △ Less

Submitted 11 April, 2022; originally announced April 2022.

Comments: Submitted to AAS Journals, comments welcome

The Dwarf Galaxy Population at $z\sim 0.7$: A Catalog of Emission Lines and Redshifts from Deep Keck Observations

Authors: John Pharo, Yicheng Guo, Guillermo Barro Calvo, Timothy Carleton, S. M. Faber, Puragra Guhathakurta, Susan A. Kassin, David C. Koo, Jack Lonergan, Teja Teppala, Weichen Wang, Hassen M. Yesuf, Fuyan Bian, Romeel Dave, John C. Forbes, Dusan Keres, Pablo Perez-Gonzalez, Alec Martin, A. J. Puleo, Lauryn Williams, Benjamin Winningham

Abstract: 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 spectrosco… ▽ More 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. △ Less

Submitted 25 July, 2022; v1 submitted 17 March, 2022; originally announced March 2022.

Comments: 23 pages, 19 Figures, updated to version accepted by ApJS

On the robustness of the velocity anisotropy parameter in probing the stellar kinematics in Milky Way like galaxies: Take away from TNG50 simulation

Authors: Razieh Emami, Lars Hernquist, Mark Vogelsberger, Xuejian Shen, Joshua S. Speagle, Jorge Moreno, Charles Alcock, Shy Genel, John C. Forbes, Federico Marinacci, Paul Torrey

Abstract: We analyze the velocity anisotropy of stars in real and energy space for a sample of Milky Way-like galaxies in the TNG50 simulation. We employ different selection criteria, including spatial, kinematic and metallicity cuts, and make three halo classes ($\mathcal{A}$-$\mathcal{C}$) which show mild-to-strong sensitivity to different selections. The above classes cover 48%, 16% and 36% of halos, res… ▽ More We analyze the velocity anisotropy of stars in real and energy space for a sample of Milky Way-like galaxies in the TNG50 simulation. We employ different selection criteria, including spatial, kinematic and metallicity cuts, and make three halo classes ($\mathcal{A}$-$\mathcal{C}$) which show mild-to-strong sensitivity to different selections. The above classes cover 48%, 16% and 36% of halos, respectively. We analyze the $β$ radial profiles and divide them into either monotonically increasing radial profiles or ones with peaks and troughs. We demonstrate that halos with monotonically increasing $β$ profiles are mostly from class $\mathcal{A}$, whilst those with peaks/troughs are part of classes $\mathcal{B}$-$\mathcal{C}$. This means that care must be taken as the observationally reported peaks/troughs might be a consequence of different selection criteria. We infer the anisotropy parameter $β$ energy space and compare that against the $β$ radial profile. It is seen that 65% of halos with very mild sensitivity to different selections in real space, are those for which the $β$ radial and energy profiles are closely related. Consequently, we propose that comparing the $β$ radial and energy profiles might be a novel way to examine the sensitivity to different selection criteria and thus examining the robustness of the anisotropy parameter in tracing stellar kinematics. We compare simulated $β$ radial profiles against various observations and demonstrate that, in most cases, the model diversity is comparable with the error bars from different observations, meaning that the TNG50 models are in good overall agreement with observations. △ Less

Submitted 8 August, 2022; v1 submitted 14 February, 2022; originally announced February 2022.

Comments: 21 pages, 16 figures, Accepted for publication in ApJ

Journal ref: ApJ, 2022

A Bayesian Population Model for the Observed Dust Attenuation in Galaxies

Authors: Gautam Nagaraj, John C. Forbes, Joel Leja, Daniel Foreman-Mackey, Christopher C. Hayward

Abstract: 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… ▽ More 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 ($τ$), though less so than in previous studies. $τ$ increases strongly with SFR, though when $\log~{\rm SFR}\lesssim 0$, $τ$ remains roughly constant over a wide range of stellar masses. Edge-on galaxies tend to have larger $τ$ 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. △ Less

Submitted 10 February, 2022; originally announced February 2022.

Comments: 27 pages, 14 figures, 1 table. Submitted to ApJ

Where to find over-massive brown dwarfs: new benchmark systems for binary evolution

Authors: Dorsa Majidi, John C. Forbes, Abraham Loeb

Abstract: 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 w… ▽ More 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. △ Less

Submitted 11 April, 2022; v1 submitted 14 September, 2021; originally announced September 2021.

Comments: Accepted at ApJ, includes a new figure showing predictions for Gaia

A Solar System formation analogue in the Ophiuchus star-forming complex

Authors: John C. Forbes, João Alves, Douglas N. C. Lin

Abstract: Anomalies among the daughter nuclei of the extinct short-lived radionuclides (SLRs) in the calcium-aluminum-rich inclusions (CAIs) indicate that the Solar System must have been born near a source of the SLRs so that they could be incorporated before they decayed away. $γ$-rays from one such living SLR, $^{26}$Al, are detected in only a few nearby star-forming regions. Here we employ multi-waveleng… ▽ More Anomalies among the daughter nuclei of the extinct short-lived radionuclides (SLRs) in the calcium-aluminum-rich inclusions (CAIs) indicate that the Solar System must have been born near a source of the SLRs so that they could be incorporated before they decayed away. $γ$-rays from one such living SLR, $^{26}$Al, are detected in only a few nearby star-forming regions. Here we employ multi-wavelength observations to demonstrate that one such region, Ophiuchus, containing many pre-stellar cores that may serve as analogs for the emerging Solar System, is inundated with $^{26}$Al from the neighboring Upper-Scorpius association, and so may provide concrete guidance for how SLR enrichment proceeded in the Solar System complementary to the meteoritics. We demonstrate via Bayesian forward modeling drawing on a wide range of observational and theoretical results that this $^{26}$Al likely 1) arises from supernova explosions, 2) arises from multiple stars, 3) has enriched the gas prior to the formation of the cores, and 4) gives rise to a broad distribution of core enrichment spanning about two orders of magnitude. This means that if the spread in CAI ages is small, as it is in the Solar System, protoplanetary disks must suffer a global heating event. △ Less

Submitted 20 August, 2021; originally announced August 2021.

Comments: Published in Nature Astronomy, August 16, 2021

Journal ref: Nature Astronomy (2021)

Touching the Stars: Using High-Resolution 3D Printing to Visualize Stellar Nurseries

Authors: Nia Imara, John C. Forbes, James C. Weaver

Abstract: Owing to their intricate variable density architecture, and as a principal site of star formation, molecular clouds represent one of the most functionally significant, yet least understood features of our universe. To unravel the intrinsic structural complexity of molecular clouds, here we leverage the power of high-resolution bitmap-based 3D printing, which provides the opportunity to visualize a… ▽ More Owing to their intricate variable density architecture, and as a principal site of star formation, molecular clouds represent one of the most functionally significant, yet least understood features of our universe. To unravel the intrinsic structural complexity of molecular clouds, here we leverage the power of high-resolution bitmap-based 3D printing, which provides the opportunity to visualize astrophysical structures in a way that uniquely taps into the human brain's ability to recognize patterns suppressed in 2D representations. Using a new suite of nine simulations, each representing different physical extremes in the turbulent interstellar medium, as our source data, our workflow permits the unambiguous visualization of features in the 3D-printed models, such as quasi-planar structures, that are frequently obscured in traditional renderings and animations. Our bitmap-based 3D printing approach thus faithfully reproduces the subtle density gradient distribution within molecular clouds in a tangible, intuitive, and visually stunning manner. While laying the groundwork for the intuitive analysis of other structurally complex astronomical data sets, our 3D-printed models also serve as valuable tools in educational and public outreach endeavors. △ Less

Submitted 30 July, 2021; originally announced August 2021.

Comments: Accepted for publication in ApJ Letters

Characterizing mass, momentum, energy and metal outflow rates of multi-phase galactic winds in the FIRE-2 cosmological simulations

Authors: Viraj Pandya, Drummond B. Fielding, Daniel Anglés-Alcázar, Rachel S. Somerville, Greg L. Bryan, Christopher C. Hayward, Jonathan Stern, Chang-Goo Kim, Eliot Quataert, John C. Forbes, Claude-André Faucher-Giguère, Robert Feldmann, Zachary Hafen, Philip F. Hopkins, Dušan Kereš, Norman Murray, Andrew Wetzel

Abstract: 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… ▽ More 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 ($<10^3$ K) are negligible except in high-redshift dwarfs. Energy, momentum and metal loading factors from the ISM are of order unity in dwarfs and significantly lower in more massive halos. Hot outflows have $2-5\times$ higher specific energy than needed to escape from the gravitational potential of dwarf halos; indeed, in dwarfs, the mass, momentum, and metal outflow rates increase with radius whereas energy is roughly conserved, indicating swept up halo gas. Burst-averaged mass loading factors tend to be larger during more powerful star formation episodes and when the inner halo is not virialized, but we see effectively no trend with the dense ISM gas fraction. We discuss how our results can guide future controlled numerical experiments that aim to elucidate the key parameters governing galactic winds and the resulting associated preventative feedback. △ Less

Submitted 17 September, 2021; v1 submitted 11 March, 2021; originally announced March 2021.

Comments: Accepted to MNRAS with minor revisions, main body is 25 pages with 14 figures

On the origin of the mass-metallicity gradient relation in the local Universe

Authors: Piyush Sharda, Mark R. Krumholz, Emily Wisnioski, Ayan Acharyya, Christoph Federrath, John C. Forbes

Abstract: 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 r… ▽ More 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. △ Less

Submitted 9 September, 2021; v1 submitted 18 February, 2021; originally announced February 2021.

Comments: 3 figures, 1 appendix, 12 pages. Accepted by MNRAS; corrected typo in equation A1

Journal ref: 2021, MNRAS, Vol. 504, Issue 1, pg. 53-64

The physics of gas phase metallicity gradients in galaxies

Authors: Piyush Sharda, Mark R. Krumholz, Emily Wisnioski, John C. Forbes, Christoph Federrath, Ayan Acharyya

Abstract: 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 diffusi… ▽ More 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. △ Less

Submitted 22 February, 2021; v1 submitted 1 February, 2021; originally announced February 2021.

Comments: 16 figures, 3 appendices, 28 pages. MNRAS, in press; v2 - updated references

The $ρ$ Oph region revisited with Gaia EDR3

Authors: Natalie Grasser, Sebastian Ratzenböck, João Alves, Josefa Großschedl, Stefan Meingast, Catherine Zucker, Alvaro Hacar, Charles Lada, Alyssa Goodman, Marco Lombardi, John C. Forbes, Immanuel M. Bomze, Torsten Möller

Abstract: Context. Young and embedded stellar populations are important probes of the star formation process. Paradoxically, we have a better census of nearby embedded young populations than the slightly more evolved optically visible young populations. The high accuracy measurements and all-sky coverage of Gaia data are about to change this situation. Aims. This work aims to construct the most complete sam… ▽ More Context. Young and embedded stellar populations are important probes of the star formation process. Paradoxically, we have a better census of nearby embedded young populations than the slightly more evolved optically visible young populations. The high accuracy measurements and all-sky coverage of Gaia data are about to change this situation. Aims. This work aims to construct the most complete sample to date of YSOs in the $ρ$ Oph region. Methods. We compile a catalog of 1114 Ophiuchus YSOs from the literature and crossmatch it with the Gaia EDR3, Gaia-ESO and APOGEE-2 surveys. We apply a multivariate classification algorithm to this catalog to identify new, co-moving population candidates. Results. We find 191 new high-fidelity YSO candidates in the Gaia EDR3 catalog belonging to the $ρ$ Oph region. The new sources appear to be mainly Class III M-stars and substellar objects and are less extincted than the known members. We find 28 previously unknown sources with disks. The analysis of the proper motion distribution of the entire sample reveals a well-defined bimodality, implying two distinct populations sharing a similar 3D volume. The first population comprises young stars' clusters around the $ρ$ Ophiuchi star and the main Ophiuchus clouds (L1688, L1689, L1709). In contrast, the second population is older ($\sim$ 10 Myr), dispersed, has a distinct proper motion, and is possibly from the Upper Sco group. The two populations are moving away from each other at about 4.1 km/s, and will no longer overlap in about 4 Myr. Finally, we flag 17 sources in the literature as impostors, which are sources that exhibit large deviations from the average distance and proper motion properties of the $ρ$ Oph population. Our results show the importance of accurate 3D space and motion information for improved stellar population analysis. (Abridged) △ Less

Submitted 8 June, 2021; v1 submitted 28 January, 2021; originally announced January 2021.

Comments: Submitted to A&A on January 28th 2021. This is the second (revised) version of the paper. All comments welcome

Journal ref: A&A 652, A2 (2021)

Inferring the Morphology of Stellar Distribution in TNG50: Twisted and Twisted-Stretched shapes

Authors: Razieh Emami, Lars Hernquist, Charles Alcock, Shy Genel, Sownak Bose, Rainer Weinberger, Mark Vogelsberger, Xuejian Shen, Joshua S. Speagle, Federico Marinacci, John C. Forbes, Paul Torrey

Abstract: 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 tenso… ▽ More 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. △ Less

Submitted 5 June, 2021; v1 submitted 22 December, 2020; originally announced December 2020.

Comments: 29 pages, 21 figures, Accepted for publication in ApJ

A Framework for Multiphase Galactic Wind Launching using TIGRESS

Authors: Chang-Goo Kim, Eve C. Ostriker, Drummond B. Fielding, Matthew C. Smith, Greg L. Bryan, Rachel S. Somerville, John C. Forbes, Shy Genel, Lars Hernquist

Abstract: 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 ga… ▽ More 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 (T<10^4 K) and hot (T>10^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 Σ_SFR and the metallicity of the ISM, and fully capture the behavior of the original TIGRESS simulation PDFs over Σ_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. △ Less

Submitted 18 October, 2020; originally announced October 2020.

Comments: ApJL accepted. For associated python package, see https://twind.readthedocs.io/en/latest/

DM halo morphological types of MW-like galaxies in the TNG50 simulation: Simple, Twisted, or Stretched

Authors: Razieh Emami, Shy Genel, Lars Hernquist, Charles Alcock, Sownak Bose, Rainer Weinberger, Mark Vogelsberger, Federico Marinacci, Abraham Loeb, Paul Torrey, John C. Forbes

Abstract: 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 (EVIM), we infer an oblate-to-triaxial shape for the DM halo with the median $T \simeq 0.24 $. We group DM halos in 3 different categories. Simple halos (32% of population) establish principal axes whose ordering in magn… ▽ More 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 (EVIM), we infer an oblate-to-triaxial shape for the DM halo with the median $T \simeq 0.24 $. We group DM halos in 3 different categories. Simple halos (32% of 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% of population), experience levels of gradual rotations throughout their radial profiles. Finally, stretched halos (36% of population) demonstrate a stretching in their principal axes lengths where the ordering of different eigenvalues change with radius. Subsequently, the halo experiences a "rotation" of $\sim$90 deg where the stretching occurs. Visualizing the 3D ellipsoid of each halo, for the first time, we report signs of re-orienting 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 impact of substructure on DM halo shape in both hydro and DMO simulations and confirm that their impacts are subdominant. We study the distribution of satellites in our sample. In simple and twisted halos, the angle of satellites' angular momentum with 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. △ Less

Submitted 23 March, 2021; v1 submitted 19 September, 2020; originally announced September 2020.

Comments: 32 pages, 28 plots, Accepted for publication in ApJ

The Diversity and Variability of Star Formation Histories in Models of Galaxy Evolution

Authors: Kartheik G. Iyer, Sandro Tacchella, Shy Genel, Christopher C. Hayward, Lars Hernquist, Alyson M. Brooks, Neven Caplar, Romeel Davé, Benedikt Diemer, John C. Forbes, Eric Gawiser, Rachel S. Somerville, Tjitske K. Starkenburg

Abstract: 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 (FIR… ▽ More 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. △ Less

Submitted 15 July, 2020; originally announced July 2020.

Comments: 31 pages, 17 figures (+ appendix). Resubmitted to MNRAS after responding to referee's comments. Comments are welcome!

First results from SMAUG: The need for preventative stellar feedback and improved baryon cycling in semi-analytic models of galaxy formation

Authors: Viraj Pandya, Rachel S. Somerville, Daniel Anglés-Alcázar, Christopher C. Hayward, Greg L. Bryan, Drummond B. Fielding, John C. Forbes, Blakesley Burkhart, Shy Genel, Lars Hernquist, Chang-Goo Kim, Stephanie Tonnesen, Tjitske Starkenburg

Abstract: 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… ▽ More 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. △ Less

Submitted 30 October, 2020; v1 submitted 29 June, 2020; originally announced June 2020.

Comments: 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

First results from SMAUG: Uncovering the Origin of the Multiphase Circumgalactic Medium with a Comparative Analysis of Idealized and Cosmological Simulations

Authors: Drummond B. Fielding, Stephanie Tonnesen, Daniel DeFelippis, Miao Li, Kung-Yi Su, Greg L. Bryan, Chang-Goo Kim, John C. Forbes, Rachel S. Somerville, Nicholas Battaglia, Evan E. Schneider, Yuan Li, Ena Choi, Christopher C. Hayward, Lars Hernquist

Abstract: 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, resolu… ▽ More 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. △ Less

Submitted 29 June, 2020; originally announced June 2020.

Comments: 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

First results from SMAUG: Characterization of Multiphase Galactic Outflows from a Suite of Local Star-Forming Galactic Disk Simulations

Authors: Chang-Goo Kim, Eve C. Ostriker, Rachel S. Somerville, Greg L. Bryan, Drummond B. Fielding, John C. Forbes, Christopher C. Hayward, Lars Hernquist, Viraj Pandya

Abstract: 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, temperatu… ▽ More 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 (Σ_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 Σ_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 Σ_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. △ Less

Submitted 27 July, 2020; v1 submitted 29 June, 2020; originally announced June 2020.

Comments: 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

First results from SMAUG: Insights into star formation conditions from spatially-resolved ISM properties in TNG50

Authors: Bhawna Motwani, Shy Genel, Greg L. Bryan, Chang-Goo Kim, Viraj Pandya, Rachel S. Somerville, Matthew C. Smith, Eve C. Ostriker, Dylan Nelson, Annalisa Pillepich, John C. Forbes, Francesco Belfiore, Rüdiger Pakmor, Lars Hernquist

Abstract: 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-… ▽ More 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. △ Less

Submitted 22 October, 2020; v1 submitted 29 June, 2020; originally announced June 2020.

Comments: 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

Stochastic modelling of star-formation histories II: star-formation variability from molecular clouds and gas inflow

Authors: Sandro Tacchella, John C. Forbes, Neven Caplar

Abstract: 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"… ▽ More 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 $β\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. △ Less

Submitted 6 July, 2020; v1 submitted 16 June, 2020; originally announced June 2020.

Comments: Accepted to MNRAS, 26 pages with 12 figures (+ appendix)

A PDF PSA, or Never gonna set_xscale again -- guilty feats with logarithms

Authors: John C. Forbes

Abstract: In the course of doing astronomy, one often encounters plots of densities, for example probability densities, flux densities, and mass functions. Quite frequently the ordinate of these diagrams is plotted logarithmically to accommodate a large dynamic range. In this situation, I argue that it is critical to adjust the density appropriately, rather than simply setting the x-scale to `log' in your f… ▽ More In the course of doing astronomy, one often encounters plots of densities, for example probability densities, flux densities, and mass functions. Quite frequently the ordinate of these diagrams is plotted logarithmically to accommodate a large dynamic range. In this situation, I argue that it is critical to adjust the density appropriately, rather than simply setting the x-scale to `log' in your favorite plotting code. I will demonstrate the basic issue with a pedagogical example, then mention a few common plots where this may arise, and finally some possible exceptions to the rule. △ Less

Submitted 31 March, 2020; originally announced March 2020.

Comments: To be submitted to the Astro-pedantic Journal on April 1, 2020

Turning up the heat on `Oumuamua

Authors: John C. Forbes, Abraham Loeb

Abstract: We explore what may be learned by close encounters between extrasolar minor bodies like `Oumuamua and the Sun. These encounters may yield strong constraints on the bulk composition and possible origin of `Oumuamua-like objects. We find that such objects collide with the Sun once every 30 years, while about 2 pass within the orbit of Mercury each year. We identify preferred orientations for the orb… ▽ More We explore what may be learned by close encounters between extrasolar minor bodies like `Oumuamua and the Sun. These encounters may yield strong constraints on the bulk composition and possible origin of `Oumuamua-like objects. We find that such objects collide with the Sun once every 30 years, while about 2 pass within the orbit of Mercury each year. We identify preferred orientations for the orbits of extrasolar objects and point out known Solar System bodies with these orientations. We conclude using a simple Bayesian analysis that about one of these objects is extrasolar in origin, even if we cannot tell which. △ Less

Submitted 2 January, 2019; originally announced January 2019.

Comments: Submitted to ApJL

Hydrodynamic Shielding and the Survival of Cold Streams

Authors: John C. Forbes, Douglas N. C. Lin

Abstract: 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… ▽ More 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. △ Less

Submitted 30 October, 2018; originally announced October 2018.

Comments: Submitted to AAS Journals, comments welcome

Towards a radially-resolved semi-analytic model for the evolution of disc galaxies tuned with machine learning

Authors: John C. Forbes, Mark R. Krumholz, Joshua S. Speagle

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

Submitted 30 October, 2018; originally announced October 2018.

Comments: Submitted to MNRAS, comments welcome

Modeling the atomic-to-molecular transition in cosmological simulations of galaxy formation

Authors: Benedikt Diemer, Adam R. H. Stevens, John C. Forbes, Federico Marinacci, Lars Hernquist, Claudia del P. Lagos, Amiel Sternberg, Annalisa Pillepich, Dylan Nelson, Gergö Popping, Francisco Villaescusa-Navarro, Paul Torrey, Mark Vogelsberger

Abstract: Large-scale cosmological simulations of galaxy formation currently do not resolve the densities at which molecular hydrogen forms, implying that the atomic-to-molecular transition must be modeled either on the fly or in postprocessing. We present an improved postprocessing framework to estimate the abundance of atomic and molecular hydrogen and apply it to the IllustrisTNG simulations. We compare… ▽ More Large-scale cosmological simulations of galaxy formation currently do not resolve the densities at which molecular hydrogen forms, implying that the atomic-to-molecular transition must be modeled either on the fly or in postprocessing. We present an improved postprocessing framework to estimate the abundance of atomic and molecular hydrogen and apply it to the IllustrisTNG simulations. We compare five different models for the atomic-to-molecular transition, including empirical, simulation-based, and theoretical prescriptions. Most of these models rely on the surface density of neutral hydrogen and the ultraviolet (UV) flux in the Lyman-Werner band as input parameters. Computing these quantities on the kiloparsec scales resolved by the simulations emerges as the main challenge. We show that the commonly used Jeans length approximation to the column density of a system can be biased and exhibits large cell-to-cell scatter. Instead, we propose to compute all surface quantities in face-on projections and perform the modeling in two dimensions. In general, the two methods agree on average, but their predictions diverge for individual galaxies and for models based on the observed midplane pressure of galaxies. We model the UV radiation from young stars by assuming a constant escape fraction and optically thin propagation throughout the galaxy. With these improvements, we find that the five models for the atomic-to-molecular transition roughly agree on average but that the details of the modeling matter for individual galaxies and the spatial distribution of molecular hydrogen. We emphasize that the estimated molecular fractions are approximate due to the significant systematic uncertainties. △ Less

Submitted 23 October, 2018; v1 submitted 6 June, 2018; originally announced June 2018.

Comments: 22 pages, 13 figures

Journal ref: 2018 ApJS 238, 33

On the existence of brown dwarfs more massive than the hydrogen burning limit

Authors: John C. Forbes, Abraham Loeb

Abstract: 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 "ove… ▽ More 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. △ Less

Submitted 12 November, 2018; v1 submitted 30 May, 2018; originally announced May 2018.

Comments: Resubmitted to AAS Journals after revisions

Habitable Evaporated Cores and the Occurrence of Panspermia near the Galactic Center

Authors: Howard Chen, John C. Forbes, Abraham Loeb

Abstract: 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… ▽ More 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. △ Less

Submitted 29 January, 2018; v1 submitted 17 November, 2017; originally announced November 2017.

Comments: 9 pages, 4 figures, accepted for publication in Astrophysical Journal Letters

A Unified Model for Galactic Discs: Star Formation, Turbulence Driving, and Mass Transport

Authors: Mark R. Krumholz, Blakesley Burkhart, John C. Forbes, Roland M. Crocker

Abstract: 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… ▽ More 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. △ Less

Submitted 29 March, 2018; v1 submitted 31 May, 2017; originally announced June 2017.

Comments: 26 pages, 7 figures, MNRAS in press; minor changes from previous version; all plotting scripts and data available from https://bitbucket.org/krumholz/kbfc17

Evaporation of planetary atmospheres due to XUV illumination by quasars

Authors: John C. Forbes, Abraham Loeb

Abstract: 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… ▽ More 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. △ Less

Submitted 18 May, 2017; originally announced May 2017.

Comments: Submitted to MNRAS

Suppression of star formation in dwarf galaxies by grain photoelectric feedback

Authors: John C. Forbes, Mark R. Krumholz, Nathan J. Goldbaum, Avishai Dekel

Abstract: 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 sug… ▽ More 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. △ Less

Submitted 2 May, 2016; originally announced May 2016.

Comments: Nature Accepted

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

Authors: Nathan J. Goldbaum, Mark R. Krumholz, John C. Forbes

Abstract: 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… ▽ More 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. △ Less

Submitted 5 May, 2016; v1 submitted 2 May, 2016; originally announced May 2016.

Comments: 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")

Stellar Mass--Gas-phase Metallicity Relation at $0.5\leq z\leq0.7$: A Power Law with Increasing Scatter toward the Low-mass Regime

Authors: Yicheng Guo, David C. Koo, Yu Lu, John C. Forbes, Marc Rafelski, Jonathan R. Trump, Ricardo Amorín, Guillermo Barro, Romeel Davé, S. M. Faber, Nimish P. Hathi, Hassen Yesuf, Michael C. Cooper, Avishai Dekel, Puragra Guhathakurta, Evan N. Kirby, Anton M. Koekemoer, Pablo G. Pérez-González, Lihwai Lin, Jeffery A. Newman, Joel R. Primack, David J. Rosario, Christopher N. A. Willmer, Renbin Yan

Abstract: 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… ▽ More 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. △ Less

Submitted 14 April, 2016; v1 submitted 15 March, 2016; originally announced March 2016.

Comments: 18 pages, 10 figures. Accepted by ApJ. Typos corrected

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

Authors: Nathan J. Goldbaum, Mark R. Krumholz, John C. Forbes

Abstract: 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… ▽ More 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. △ Less

Submitted 28 October, 2015; originally announced October 2015.

Comments: 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

Curveballs in protoplanetary disks - the effect of the Magnus force on planet formation

Authors: John C. Forbes

Abstract: Spinning planetesimals in a gaseous protoplanetary disk may experience a hydrodynamical force perpendicular to their relative velocities. We examine the effect this force has on the dynamics of these objects using analytical arguments based on a simple laminar disk model and numerical integrations of the equations of motion for individual grains. We focus in particular on meter-sized boulders trad… ▽ More Spinning planetesimals in a gaseous protoplanetary disk may experience a hydrodynamical force perpendicular to their relative velocities. We examine the effect this force has on the dynamics of these objects using analytical arguments based on a simple laminar disk model and numerical integrations of the equations of motion for individual grains. We focus in particular on meter-sized boulders traditionally expected to spiral in to the central star in as little as 100 years from 1 A.U. We find that there are plausible scenarios in which this force extends the lifetime of these solids in the disk by a factor of several. More importantly the velocities induced by the Magnus force can prevent the formation of planetesimals via gravitational instability in the inner disk if the size of the dust particles is larger than of order 10 cm. We find that the fastest growing linear modes of the streaming instability may still grow despite the diffusive effect of the Magnus force, but it remains to be seen how the Magnus force will alter the non-linear evolution of these instabilities. △ Less

Submitted 28 July, 2015; originally announced July 2015.

Comments: MNRAS Accepted

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

Authors: Antoine C. Petit, Mark R. Krumholz, Nathan J. Goldbaum, John C. Forbes

Abstract: 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 gr… ▽ More 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. △ Less

Submitted 4 March, 2015; v1 submitted 27 November, 2014; originally announced November 2014.

Comments: 10 pages, 13 figures

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

Authors: Mark R. Krumholz, John C. Forbes

Abstract: 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… ▽ More 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. △ Less

Submitted 23 February, 2015; v1 submitted 25 June, 2014; originally announced June 2014.

Comments: 58 pages, 13 figures, accepted to Astronomy & Computing; this version includes more discussion, but no other changes; code is available for download from https://bitbucket.org/krumholz/vader/