BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

John C. Forbes

John C. Forbes appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.GAastro-ph.COastro-ph.IMastro-ph.SRastro-ph.EP

Catalog footprint

What is connected

20works
5topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

20 published item(s)

preprint2026arXiv

Pink Dwarfs and the Paths to Stardom: How Brown Dwarfs Pushed Above the Hydrogen Burning Limit Evolve

Brown dwarfs that gain mass through binary interactions may be pushed above the boundary that divides brown dwarfs from low-mass stars: the hydrogen burning limit (HBL). Some of these objects will make their way to the main sequence and may eventually be indistinguishable from ordinary low-mass stars, while others will remain brown dwarf-like, unable to burn hydrogen at a high enough rate to power their surface luminosity. We study the evolution of both types of object to provide a taxonomy and testable observational predictions for these objects depending on their evolutionary path. Using MESA simulations, we find that a subset of the objects that will eventually become stars experience an extended luminosity plateau, where their surface luminosity remains nearly constant on 100 Myr - Gyr timescales. We find that the plateau timescale is set by the amount of energy required to re-heat the cores of these objects to a level sufficient to sustain convection. The timescales required for the cores of these objects to "unfreeze" and arrive at the main sequence is long enough that surveys may be able to find objects in this evolutionary stage. These objects, along with those that never reach the main sequence, occupy a unique space in a mass-luminosity diagram, and would provide a unique constraint on binary mass transfer physics.

0 citations0 reviewsNo code linkedOpen access
preprint2024arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A Bayesian Population Model for the Observed Dust Attenuation in Galaxies

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The physics of gas phase metallicity gradients in galaxies

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Suppression of star formation in dwarf galaxies by grain photoelectric feedback

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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

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

0 citations0 reviewsNo code linkedOpen access