Source author record

Sarah Wellons

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

9works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

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

9 published item(s)

preprint2022arXiv

Running Late: Testing Delayed Supermassive Black Hole Growth Models Against the Quasar Luminosity Function

Observations of massive galaxies at low redshift have revealed approximately linear scaling relations between the mass of a supermassive black hole (SMBH) and properties of its host galaxy. How these scaling relations evolve with redshift and whether they extend to lower-mass galaxies however remain open questions. Recent galaxy formation simulations predict a delayed, or "two-phase", growth of SMBHs: slow, highly intermittent BH growth due to repeated gas ejection by stellar feedback in low-mass galaxies, followed by more sustained gas accretion that eventually brings BHs onto the local scaling relations. The predicted two-phase growth implies a steep increase, or "kink", in BH-galaxy scaling relations at a stellar mass $M_{*}\sim5\times10^{10} M_{\odot}$. We develop a parametric, semi-analytic model to compare different SMBH growth models against observations of the quasar luminosity function (QLF) at $z\sim0.5-4$. We compare models in which the relation between SMBH mass and galaxy mass is purely linear versus two-phase models. The models are anchored to the observed galaxy stellar mass function, and the BH mass functions at different redshifts are consistently connected by the accretion rates contributing to the QLF. The best fits suggest that two-phase evolution is significantly preferred by the QLF data over a purely linear scaling relation. Moreover, when the model parameters are left free, the two-phase model fits imply a transition mass consistent with that predicted by simulations. Our analysis motivates further observational tests, including measurements of BH masses and AGN activity at the low-mass end, which could more directly test two-phase SMBH growth.

preprint2022arXiv

Why do Black Holes Trace Bulges (& Central Surface Densities), Instead of Galaxies as a Whole?

Previous studies of fueling black holes (BHs) in galactic nuclei have argued (on scales ~0.01-1000pc) accretion is dynamical with inflow rates $\dot{M}\simη\,M_{\rm gas}/t_{\rm dyn}$ in terms of gas mass $M_{\rm gas}$, dynamical time $t_{\rm dyn}$, and some $η$. But these models generally neglected expulsion of gas by stellar feedback, or considered extremely high densities where expulsion is inefficient. Studies of star formation, however, have shown on sub-kpc scales the expulsion efficiency $f_{\rm wind}=M_{\rm ejected}/M_{\rm total}$ scales with the gravitational acceleration as $(1-f_{\rm wind})/f_{\rm wind}\sim\bar{a}_{\rm grav}/\langle\dot{p}/m_{\ast}\rangle\sim Σ_{\rm eff}/Σ_{\rm crit}$ where $\bar{a}_{\rm grav}\equiv G\,M_{\rm tot}(<r)/r^{2}$ and $\langle\dot{p}/m_{\ast}\rangle$ is the momentum injection rate from young stars. Adopting this as the simplest correction for stellar feedback, $η\rightarrow η\,(1-f_{\rm wind})$, we show this provides a more accurate description of simulations with stellar feedback at low densities. This has immediate consequences, predicting e.g. the slope and normalization of the $M-σ$ and $M-M_{\rm bulge}$ relation, $L_{\rm AGN}-$SFR relations, and explanations for outliers in compact Es. Most strikingly, because star formation simulations show expulsion is efficient ($f_{\rm wind}\sim1$) below total-mass surface density $M_{\rm tot}/π\,r^{2}<Σ_{\rm crit}\sim3\times10^{9}\,M_{\odot}\,{\rm kpc^{-2}}$ (where $Σ_{\rm crit}=\langle\dot{p}/m_{\ast}\rangle/(π\,G)$), BH mass is predicted to specifically trace host galaxy properties above a critical surface brightness $Σ_{\rm crit}$ (B-band $μ_{\rm B}^{\rm crit}\sim 19\,{\rm mag\,arcsec^{-2}}$). This naturally explains why BH masses preferentially reflect bulge properties or central surface-densities ($Σ_{1\,{\rm kpc}}$), not 'total' galaxy properties.

preprint2020arXiv

Measuring dynamical masses from gas kinematics in simulated high-redshift galaxies

Advances in instrumentation have recently extended detailed measurements of gas kinematics to large samples of high-redshift galaxies. Relative to most nearby, thin disk galaxies, in which gas rotation accurately traces the gravitational potential, the interstellar medium (ISM) of z>1 galaxies is typically more dynamic and exhibits elevated turbulence. If not properly modeled, these effects can strongly bias dynamical mass measurements. We use high-resolution FIRE-2 cosmological zoom-in simulations to analyze the physical effects that must be considered to correctly infer dynamical masses from gas kinematics. Our analysis covers a range of galaxy properties from low-redshift Milky-Way-mass galaxies to massive high-redshift galaxies (M_* > 10^11 M_sun at z=1). Selecting only snapshots where a disk is present, we calculate the rotational profile v_phi(r) of the cool (10^3.5 K < T < 10^4.5 K) gas and compare it to the circular velocity v_c=sqrt(GM/r). In the simulated galaxies, the gas rotation traces the circular velocity at intermediate radii, but the two quantities diverge significantly in the center and in the outer disk. Our simulations appear to over-predict observed rotational velocities in the centers of massive galaxies (likely from a lack of black hole feedback), so we focus on larger radii. Gradients in the turbulent pressure at these radii can provide additional radial support and bias dynamical mass measurements low by up to 40%. In both the interior and exterior, the gas' motion can be significantly non-circular due to e.g. bars, satellites, and inflows/outflows. We discuss the accuracy of commonly-used analytic models for pressure gradients (or "asymmetric drift") in the ISM of high-redshift galaxies.

preprint2020arXiv

Recent Star-formation in a Massive Slowly-Quenched Lensed Quiescent Galaxy at z=1.88

In this letter, we reconstruct the formation pathway of MRG-S0851, a massive, $\log M_*/M_\odot=11.02\pm0.04$, strongly lensed, red, galaxy at $z=1.883\pm0.001$. While the global photometry and spatially-resolved outskirts of MRG-S0851 imply an early-formation scenario with a slowly decreasing or constant star-formation history, a joint fit of 2D grism spectroscopy and photometry reveals a more complex scenario: MRG-S0851 is likely to be experiencing a centrally-concentrated rejuvenation in the inner $\sim$1 kpc in the last $\sim$100 Myr of evolution. We estimate $0.5\pm0.1\%$ of the total stellar mass is formed in this phase. Rejuvenation episodes are suggested to be infrequent for massive galaxies at $z\sim2$, but as our analyses indicate, more examples of complex star-formation histories may yet be hidden within existing data. By adding a FUV color criterion to the standard U-V/V-J diagnostic, thereby heightening our sensitivity to recent star formation, we show that we can select populations of galaxies with similar spectral energy distributions to that of MRG-S0851, but note that deep follow-up spectroscopic observations and/or spatially resolved analyses are necessary to robustly confirm the rejuvenation of these candidates. Using our criteria with MRG-S0851 as a prototype, we estimate that $\sim$1\% of massive quiescent galaxies at $1<z<2$ are potentially rejuvenating.

preprint2015arXiv

An analysis of the evolving comoving number density of galaxies in hydrodynamical simulations

The cumulative comoving number-density of galaxies as a function of stellar mass or central velocity dispersion is commonly used to link galaxy populations across different epochs. By assuming that galaxies preserve their number-density in time, one can infer the evolution of their properties, such as masses, sizes, and morphologies. However, this assumption does not hold in the presence of galaxy mergers or when rank ordering is broken owing to variable stellar growth rates. We present an analysis of the evolving comoving number density of galaxy populations found in the Illustris cosmological hydrodynamical simulation focused on the redshift range $0\leq z \leq 3$. Our primary results are as follows: 1) The inferred average stellar mass evolution obtained via a constant comoving number density assumption is systematically biased compared to the merger tree results at the factor of $\sim$2(4) level when tracking galaxies from redshift $z=0$ out to redshift $z=2(3)$; 2) The median number density evolution for galaxy populations tracked forward in time is shallower than for galaxy populations tracked backward in time; 3) A similar evolution in the median number density of tracked galaxy populations is found regardless of whether number density is assigned via stellar mass, stellar velocity dispersion, or dark matter halo mass; 4) Explicit tracking reveals a large diversity in galaxies' assembly histories that cannot be captured by constant number-density analyses; 5) The significant scatter in galaxy linking methods is only marginally reduced by considering a number of additional physical and observable galaxy properties as realized in our simulation. We provide fits for the forward and backward median evolution in stellar mass and number density and discuss implications of our analysis for interpreting multi-epoch galaxy property observations.

preprint2015arXiv

The diverse evolutionary paths of simulated high-z massive, compact galaxies to z=0

Massive quiescent galaxies have much smaller physical sizes at high redshift than today. The strong evolution of galaxy size may be caused by progenitor bias, major and minor mergers, adiabatic expansion, and/or renewed star formation, but it is difficult to test these theories observationally. Herein, we select a sample of 35 massive, compact galaxies ($M_* = 1-3 \times 10^{11}$ M$_\odot$, $M_*/R^{1.5} > 10^{10.5}$ M$_\odot$/kpc$^{1.5}$) at $z=2$ in the cosmological hydrodynamical simulation Illustris and trace them forward to $z=0$ to uncover their evolution and identify their descendants. By $z=0$, the original factor of 3 difference in stellar mass spreads to a factor of 20. The dark matter halo masses similarly spread from a factor of 5 to 40. The galaxies' evolutionary paths are diverse: about half acquire an ex-situ envelope and are the core of a more massive descendant, a third survive undisturbed and gain very little mass, 15% are consumed in a merger with a more massive galaxy, and a small remainder are thoroughly mixed by major mergers. The galaxies grow in size as well as mass, and only $\sim$10% remain compact by $z=0$. The majority of the size growth is driven by the acquisition of ex-situ mass. The most massive galaxies at $z=0$ are the most likely to have compact progenitors, but this trend possesses significant dispersion which precludes a direct linkage to compact galaxies at $z=2$. The compact galaxies' merger rates are influenced by their $z=2$ environments, so that isolated or satellite compact galaxies (which are protected from mergers) are the most likely to survive to the present day.

preprint2015arXiv

The Formation of Massive, Compact Galaxies at z=2 in the Illustris Simulation

Massive, quiescent galaxies at high redshift have been found to be considerably more compact than galaxies of similar mass in the local universe. How these compact galaxies formed has yet to be determined, though several progenitor populations have been proposed. Here we investigate the formation processes and quantify the assembly histories of such galaxies in Illustris, a suite of hydrodynamical cosmological simulations encompassing a sufficiently large volume to include rare objects, while simultaneously resolving the internal structure of galaxies. We select massive (~10^11 solar masses) and compact (stellar half-mass radius < 2 kpc) galaxies from the simulation at z=2. Within the Illustris suite, we find that these quantities are not perfectly converged, but are reasonably reliable for our purposes. The resulting population is composed primarily of quiescent galaxies, but we also find several star-forming compact galaxies. The simulated compact galaxies are similar to observed galaxies in star formation activity and appearance. We follow their evolution at high redshift in the simulation and find that there are multiple pathways to form these compact galaxies, dominated by two mechanisms: (i) intense, centrally concentrated starbursts generally triggered by gas-rich major mergers between z~2-4, reducing the galaxies' half-mass radii by a factor of a few to below 2 kpc, and (ii) assembly at very early times when the universe was much denser; the galaxies formed compact and remained so until z~2.

preprint2014arXiv

A High-Frequency Doppler Feature in the Power Spectra of Simulated GRMHD Black Hole Accretion Disks

Black hole binaries exhibit a wide range of variability phenomena, from large-scale state changes to broadband noise and quasi-periodic oscillations, but the physical nature of much of this variability is poorly understood. We examine the variability properties of three GRMHD simulations of thin accretion disks around black holes of varying spin, producing light curves and power spectra as would be seen by observers. We find that the simulated power spectra show a broad feature at high frequency, which increases in amplitude with the inclination of the observer. We show that this high-frequency feature is a product of the Doppler effect and that its location is a function of the mass and spin of the black hole. This Doppler feature demonstrates that power spectral properties of the accretion disk can be tied to, and potentially used to determine, physical properties of the black hole.

preprint2012arXiv

Radio Observations Reveal Unusual Circumstellar Environments for Some Type Ibc Supernova Progenitors

We present extensive radio observations of the nearby Type Ibc supernovae 2004cc, 2004dk, and 2004gq spanning 8-1900 days after explosion. Using a dynamical model developed for synchrotron emission from a slightly decelerated shockwave, we estimate the velocity and energy of the fastest ejecta and the density profile of the circumstellar medium. The shockwaves of all three supernovae are characterized by non-relativistic velocities of v ~ (0.1-25)c and associated energies of E ~ (2-10) * 1e47 erg, in line with the expectations for a typical homologous explosion. Smooth circumstellar density profiles are indicated by the early radio data and we estimate the progenitor mass loss rates to be ~ (0.6-13) * 1e-5 M_sun/yr (wind velocity 10^3 km/s). These estimates approach the saturation limit (~1e-4 M_sun/yr) for line-driven winds from Wolf-Rayet stars, the favored progenitors of SNe Ibc including those associated with long-duration GRBs. Intriguingly, at later epochs all three supernovae show evidence for abrupt radio variability that we attribute to large density modulations (factor of ~3-6) at circumstellar radii of r ~ (1-50) * 1e16 cm. If due to variable mass loss, these modulations are associated with progenitor activity on a timescale of ~ 10-100 years before explosion. We consider these results in the context of variable mass loss mechanisms including wind clumping, metallicity-independent continuum-driven ejections, and binary-induced modulations. It may also be possible that the SN shockwaves are dynamically interacting with wind termination shocks, however, this requires the environment to be highly pressurized and/or the progenitor to be rapidly rotating prior to explosion. The proximity of the density modulations to the explosion sites may suggest a synchronization between unusual progenitor mass loss and the SN explosion, reminiscent of Type IIn supernovae. [ABRIDGED]