Source author record

Juna A. Kollmeier

Juna A. Kollmeier 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

44works
8topics
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

44 published item(s)

preprint2024arXiv

Spectacular nucleosynthesis from early massive stars

Stars formed with initial mass over 50 Msun are very rare today, but they are thought to be more common in the early universe. The fates of those early, metal-poor, massive stars are highly uncertain. Most are expected to directly collapse to black holes, while some may explode as a result of rotationally powered engines or the pair-creation instability. We present the chemical abundances of J0931+0038, a nearby low-mass star identified in early followup of SDSS-V Milky Way Mapper, which preserves the signature of unusual nucleosynthesis from a massive star in the early universe. J0931+0038 has relatively high metallicity ([Fe/H] = -1.76 +/- 0.13) but an extreme odd-even abundance pattern, with some of the lowest known abundance ratios of [N/Fe], [Na/Fe], [K/Fe], [Sc/Fe], and [Ba/Fe]. The implication is that a majority of its metals originated in a single extremely metal-poor nucleosynthetic source. An extensive search through nucleosynthesis predictions finds a clear preference for progenitors with initial mass > 50 Msun, making J0931+0038 one of the first observational constraints on nucleosynthesis in this mass range. However the full abundance pattern is not matched by any models in the literature. J0931+0038 thus presents a challenge for the next generation of nucleosynthesis models and motivates study of high-mass progenitor stars impacted by convection, rotation, jets, and/or binary companions. Though rare, more examples of unusual early nucleosynthesis in metal-poor stars should be found in upcoming large spectroscopic surveys.

preprint2022arXiv

A Spectroscopic Road Map for Cosmic Frontier: DESI, DESI-II, Stage-5

In this white paper, we present an experimental road map for spectroscopic experiments beyond DESI. DESI will be a transformative cosmological survey in the 2020s, mapping 40 million galaxies and quasars and capturing a significant fraction of the available linear modes up to z=1.2. DESI-II will pilot observations of galaxies both at much higher densities and extending to higher redshifts. A Stage-5 experiment would build out those high-density and high-redshift observations, mapping hundreds of millions of stars and galaxies in three dimensions, to address the problems of inflation, dark energy, light relativistic species, and dark matter. These spectroscopic data will also complement the next generation of weak lensing, line intensity mapping and CMB experiments and allow them to reach their full potential.

preprint2022arXiv

The MegaMapper: A Stage-5 Spectroscopic Instrument Concept for the Study of Inflation and Dark Energy

In this white paper, we present the MegaMapper concept. The MegaMapper is a proposed ground-based experiment to measure Inflation parameters and Dark Energy from galaxy redshifts at $2<z<5$. In order to achieve path-breaking results with a mid-scale investment, the MegaMapper combines existing technologies for critical path elements and pushes innovative development in other design areas. To this aim, we envision a 6.5-m Magellan-like telescope, with a newly designed wide field, coupled with DESI spectrographs, and small-pitch robots to achieve multiplexing of at least 26,000. This will match the expected achievable target density in the redshift range of interest and provide a 10x capability over the existing state-of the art, without a 10x increase in project budget.

preprint2022arXiv

Zeta-Payne: a fully automated spectrum analysis algorithm for the Milky Way Mapper program of the SDSS-V survey

The Sloan Digital Sky Survey has recently initiated its 5th survey generation (SDSS-V), with a central focus on stellar spectroscopy. In particular, SDSS-V Milky Way Mapper program will deliver multi-epoch optical and near-infrared spectra for more than 5 million stars across the entire sky, covering a large range in stellar mass, surface temperature, evolutionary stage, and age. About 10% of those spectra will be of hot stars of OBAF spectral types, for whose analysis no established survey pipelines exist. Here we present the spectral analysis algorithm, Zeta-Payne, developed specifically to obtain stellar labels from SDSS-V spectra of stars with these spectral types and drawing on machine learning tools. We provide details of the algorithm training, its test on artificial spectra, and its validation on two control samples of real stars. Analysis with Zeta-Payne leads to only modest internal uncertainties in the near-IR with APOGEE (optical with BOSS): 3-10% (1-2%) for Teff, 5-30% (5-25%) for v*sin(i), 1.7-6.3 km/s(0.7-2.2 km/s) for RV, $<0.1$ dex ($<0.05$ dex) for log(g), and 0.4-0.5 dex (0.1 dex) for [M/H] of the star, respectively. We find a good agreement between atmospheric parameters of OBAF-type stars when inferred from their high- and low-resolution optical spectra. For most stellar labels the APOGEE spectra are (far) less informative than the BOSS spectra of these stars, while log(g), v*sin(i), and [M/H] are in most cases too uncertain for meaningful astrophysical interpretation. This makes BOSS low-resolution optical spectra better for stellar labels of OBAF-type stars, unless the latter are subject to high levels of extinction.

preprint2021arXiv

Final Targeting Strategy for the SDSS-IV APOGEE-2N Survey

APOGEE-2 is a dual-hemisphere, near-infrared (NIR), spectroscopic survey with the goal of producing a chemo-dynamical mapping of the Milky Way Galaxy. The targeting for APOGEE-2 is complex and has evolved with time. In this paper, we present the updates and additions to the initial targeting strategy for APOGEE-2N presented in Zasowski et al. (2017). These modifications come in two implementation modes: (i) "Ancillary Science Programs" competitively awarded to SDSS-IV PIs through proposal calls in 2015 and 2017 for the pursuit of new scientific avenues outside the main survey, and (ii) an effective 1.5-year expansion of the survey, known as the Bright Time Extension, made possible through accrued efficiency gains over the first years of the APOGEE-2N project. For the 23 distinct ancillary programs, we provide descriptions of the scientific aims, target selection, and how to identify these targets within the APOGEE-2 sample. The Bright Time Extension permitted changes to the main survey strategy, the inclusion of new programs in response to scientific discoveries or to exploit major new datasets not available at the outset of the survey design, and expansions of existing programs to enhance their scientific success and reach. After describing the motivations, implementation, and assessment of these programs, we also leave a summary of lessons learned from nearly a decade of APOGEE-1 and APOGEE-2 survey operations. A companion paper, Santana et al. (submitted), provides a complementary presentation of targeting modifications relevant to APOGEE-2 operations in the Southern Hemisphere.

preprint2020arXiv

LATIS: The Ly$α$ Tomography IMACS Survey

We introduce LATIS, the Ly$α$ Tomography IMACS Survey, a spectroscopic survey at Magellan designed to map the z=2.2-2.8 intergalactic medium (IGM) in three dimensions by observing the Ly$α$ forest in the spectra of galaxies and QSOs. Within an area of 1.7 deg${}^2$, we will observe approximately half of $\gtrsim L^*$ galaxies at z=2.2-3.2 for typically 12 hours, providing a dense network of sightlines piercing the IGM with an average transverse separation of 2.5 $h^{-1}$ comoving Mpc (1 physical Mpc). At these scales, the opacity of the IGM is expected to be closely related to the dark matter density, and LATIS will therefore map the density field in the $z \sim 2.5$ universe at $\sim$Mpc resolution over the largest volume to date. Ultimately LATIS will produce approximately 3800 spectra of z=2.2-3.2 galaxies that probe the IGM within a volume of $4 \times 10^6 h^{-3}$ Mpc${}^3$, large enough to contain a representative sample of structures from protoclusters to large voids. Observations are already complete over one-third of the survey area. In this paper, we describe the survey design and execution. We present the largest IGM tomographic maps at comparable resolution yet made. We show that the recovered matter overdensities are broadly consistent with cosmological expectations based on realistic mock surveys, that they correspond to galaxy overdensities, and that we can recover structures identified using other tracers. LATIS is conducted in Canada-France-Hawaii Telescope Legacy Survey fields, including COSMOS. Coupling the LATIS tomographic maps with the rich data sets collected in these fields will enable novel studies of environment-dependent galaxy evolution and the galaxy-IGM connection at cosmic noon.

preprint2020arXiv

SDSS-V Algorithms: Fast, Collision-Free Trajectory Planning for Heavily Overlapping Robotic Fiber Positioners

Robotic fiber positioner (RFP) arrays are becoming heavily adopted in wide field massively multiplexed spectroscopic survey instruments. RFP arrays decrease nightly operational overheads through rapid reconfiguration between fields and exposures. In comparison to similar instruments, SDSS-V has selected a very dense RFP packing scheme where any point in a field is typically accessible to three or more robots. This design provides flexibility in target assignment. However, the task of collision-less trajectory planning is especially challenging. We present two multi-agent distributed control strategies that are highly efficient and computationally inexpensive for determining collision-free paths for RFPs in heavily overlapping workspaces. We demonstrate that a reconfiguration path between two arbitrary robot configurations can be efficiently found if "folded" state, in which all robot arms are retracted and aligned in a lattice-like orientation, is inserted between the initial and final states. Although developed for SDSS-V, the approach we describe is generic and so applicable to a wide range of RFP designs and layouts. Robotic fiber positioner technology continues to advance rapidly, and in the near future ultra-densely packed RFP designs may be feasible. Our algorithms are especially capable in routing paths in very crowded environments, where we see efficient results even in regimes significantly more crowded than the SDSS-V RFP design.

preprint2020arXiv

The Impact of Wind Scalings on Stellar Growth and the Baryon Cycle in Cosmological Simulations

Many phenomenologically successful cosmological galaxy formation simulations employ kinetic winds to model galactic outflows, a crucial ingredient in obtaining predictions that agree with various observations. Yet systematic studies of how variations in kinetic wind scalings might alter observable galaxy properties are rare. Here we employ GADGET-3 simulations to study how the baryon cycle, stellar mass function, and other galaxy and CGM predictions vary as a function of the assumed outflow speed $v_w$ and the scaling of the mass loading factor $η$ with velocity dispersion $σ$. We design our fiducial model to reproduce the measured wind properties at 25% of the virial radius from the Feedback In Realistic Environments (FIRE) simulations. We find that a strong dependence of $η\sim σ^5$ in low mass haloes with $σ< 106\ \mathrm{km\ s^{-1}}$ is required to match the faint end of the stellar mass functions at $z > 1$. The wind speed also has a major impact, with faster winds significantly reducing wind recycling and heating more halo gas. Both effects result in less stellar mass growth in massive haloes and impact high ionization absorption in halo gas. We cannot simultaneously match the stellar content at $z=2$ and $z=0$ within a single model, suggesting that an additional feedback source such as AGN might be required in massive galaxies at lower redshifts, but the amount needed depends strongly on assumptions regarding the outflow properties. We run a 50 $\mathrm{Mpc/h}$, $2\times576^3$ simulation with our fiducial parameters and show that it matches a range of star-forming galaxy properties at $z\sim0-2$. In closing, the results from simulations of galaxy formation are much more sensitive to small changes in the feedback implementation than to the hydrodynamic technique.

preprint2020arXiv

The Sixteenth Data Release of the Sloan Digital Sky Surveys: First Release from the APOGEE-2 Southern Survey and Full Release of eBOSS Spectra

This paper documents the sixteenth data release (DR16) from the Sloan Digital Sky Surveys; the fourth and penultimate from the fourth phase (SDSS-IV). This is the first release of data from the southern hemisphere survey of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2); new data from APOGEE-2 North are also included. DR16 is also notable as the final data release for the main cosmological program of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), and all raw and reduced spectra from that project are released here. DR16 also includes all the data from the Time Domain Spectroscopic Survey (TDSS) and new data from the SPectroscopic IDentification of ERosita Survey (SPIDERS) programs, both of which were co-observed on eBOSS plates. DR16 has no new data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey (or the MaNGA Stellar Library "MaStar"). We also preview future SDSS-V operations (due to start in 2020), and summarize plans for the final SDSS-IV data release (DR17).

preprint2016arXiv

GW150914: Spin based constraints on the merger time of the progenitor system

We explore the implications of the observed low spin of GW150914 within the context of stellar astrophysics and progenitor models. We conclude that many of the recently proposed scenarios are in marked tension with this observation. We derive a simple model for the observed spin in the case that the progenitor system was a field binary composed of a black hole (BH) and a Wolf-Rayet star and explore the implications of the observed spin for this model. The spin observation allows us to place a lower limit for the delay time between the formation of the BH+BH binary and the actual merger, $t_{\textrm{merge}}$. We use typical values for these systems to derive $t_{\textrm{merge}}\gtrsim10^{8}\,\textrm{yr}$, which proves to be an important diagnostic for different progenitor models. We anticipate the next series of events, and the associated spin parameters, will ultimately yield critical constraints on formation scenarios and on stellar parameters describing the late-stage evolution of massive stars.

preprint2016arXiv

Predicted Information Content of Gaia-Hipparcos

The Gaia-Tycho release, scheduled for 14 September, is forecast to yield parallax errors of $σ(π)\sim 300\,μ$as for about 2 million Tycho stars. We show analytically that the actual performance should be $$ σ(π) = {\rm max}(σ_{1991}/96,20\,μ{\rm as}) $$ where $σ_{1991}$ is the positional error from the Hipparcos mission. For typical Tycho stars, $σ_{1991}\sim 30\,$mas, so this reproduces the usual claims. However, for the 100,000 star Hipparcos subset of this sample, $σ_{1991}$ is a factor 15 or more smaller. These much lower Hipparcos positional errors apply even to stars at the Hipparcos-Tycho limit, $V\sim 12$. This is especially important for RR Lyrae stars, as well as other special classes, that were systematically included in the Hipparcos catalog down to this limit because of their exceptional scientific importance. This predicted performance will provide an early test of the Gaia algorithms.

preprint2016arXiv

TGAS Error Renormalization from the RR Lyrae Period-Luminosity Relation

The Gaia team has applied a renormalization to their internally-derived parallax errors $σ_{\rm int}(π)$ $$ σ_{tgas}(π) = \sqrt{[Aσ_{int}(π)]^2 + σ_0^2}; \ \ \ \ (A,σ_0) = (1.4,0.20\ \rm mas) $$ based on comparison to Hipparcos astrometry. We use a completely independent method based on the RR Lyrae $K$-band period-luminosity relation to derive a substantially different result, with smaller ultimate errors $$ (A,σ_0) = (1.1,0.12\ \rm mas) \ \ \ \ (this\ paper). $$ We argue that our estimate is likely to be more accurate and therefore that the reported TGAS parallax errors should be reduced according to the prescription: $$ σ_{true}(π) = \sqrt{(0.79σ_{tgas}(π))^2 - (0.10\ \rm mas)^2}. $$

preprint2016arXiv

The Carnegie-Chicago Hubble Program. I. An Independent Approach to the Extragalactic Distance Scale Using only Population II Distance Indicators

We present an overview of the Carnegie-Chicago Hubble Program, an ongoing program to obtain a 3 per cent measurement of the Hubble constant using alternative methods to the traditional Cepheid distance scale. We aim to establish a completely independent route to the Hubble constant using RR Lyrae variables, the tip of the red giant branch (TRGB), and Type Ia supernovae (SNe Ia). This alternative distance ladder can be applied to galaxies of any Hubble Type, of any inclination, and, utilizing old stars in low density environments, is robust to the degenerate effects of metallicity and interstellar extinction. Given the relatively small number of SNe Ia host galaxies with independently measured distances, these properties provide a great systematic advantage in the measurement of the Hubble constant via the distance ladder. Initially, the accuracy of our value of the Hubble constant will be set by the five Galactic RR Lyrae calibrators with Hubble Space Telescope Fine-Guidance Sensor parallaxes. With Gaia, both the RR Lyrae zero point and TRGB method will be independently calibrated, the former with at least an order of magnitude more calibrators and the latter directly through parallax measurement of tip red giants. As the first end-to-end "distance ladder" completely independent of both Cepheid variables and the Large Magellanic Cloud, this path to the Hubble constant will allow for the high precision comparison at each rung of the traditional distance ladder that is necessary to understand tensions between this and other routes to the Hubble constant.

preprint2016arXiv

Ultrahigh-Energy Cosmic Rays from the "En Caul" Birth of Magnetars

Rapidly-spinning magnetars can potentially form by the accretion induced collapse of a white dwarf or by neutron star mergers if the equation of state of nuclear density matter is such that two low mass neutron stars can sometimes form a massive neutron star rather than a black hole. In either case, the newly born magnetar is an attractive site for producing ultrahigh-energy cosmic rays (particles with individual energies exceeding $10^{18}\,{\rm eV}$; UHECRs). The short-period spin and strong magnetic field are able to accelerate particles up to the appropriate energies, and the composition of material on and around the magnetar may naturally explain recent inferences of heavy elements in UHECRs. We explore whether the small amount of natal debris surrounding these magnetars allows the UHECRs to easily escape. We also investigate the impact on the UHECRs of the unique environment around the magnetar, which consists of a bubble of relativistic particles and magnetic field within the debris. Rates and energetics of UHECRs are consistent with such an origin even though the rates of events that produce rapidly-spinning magnetars remain very uncertain. The low ejecta mass also helps limit the high-energy neutrino background associated with this scenario to be below current IceCube constraints over most of the magnetar parameter space. A unique prediction is that UHECRs may be generated in old stellar environments without strong star formation in contrast to what would be expected for other UHECR scenarios, such as active galactic nuclei or long gamma-ray bursts.

preprint2015arXiv

Baryon Cycling in the Low-Redshift Circumgalactic Medium: A Comparison of Simulations to the COS-Halos Survey

We analyze the low-redshift (z~0.2) circumgalactic medium by comparing absorption-line data from the COS-Halos Survey to absorption around a matched galaxy sample from two cosmological hydrodynamic simulations. The models include different prescriptions for galactic outflows, namely hybrid energy/momentum driven wind (ezw), and constant winds (cw). We extract for comparison direct observables including equivalent widths, covering factors, ion ratios, and kinematics. Both wind models are generally in good agreement with these observations for HI and certain low ionization metal lines, but show poorer agreement with higher ionization metal lines including SiIII and OVI that are well-observed by COS-Halos. These discrepancies suggest that both wind models predict too much cool, metal-enriched gas and not enough hot gas, and/or that the metals are not sufficiently well-mixed. This may reflect our model assumption of ejecting outflows as cool and unmixing gas. Our ezw simulation includes a heuristic prescription to quench massive galaxies by super-heating its ISM gas, which we show yields sufficient low ionisation absorption to be broadly consistent with observations, but also substantial OVI absorption that is inconsistent with data, suggesting that gas around quenched galaxies in the real Universe does not cool. At impact parameters of <50 kpc, recycling winds dominate the absorption of low ions and even HI, while OVI almost always arises from metals ejected longer than 1 Gyr ago. The similarity between the wind models is surprising, since we show that they differ substantially in their predicted amount and phase distribution of halo gas. We show that this similarity owes mainly to our comparison here at fixed stellar mass rather than at fixed halo mass in our previous works, which suggests that CGM properties are more closely tied to the stellar mass of galaxies rather than halo mass.

preprint2015arXiv

Rest-frame UV single-epoch black hole mass estimates of low-luminosity AGN at intermediate redshifts

The ability to accurately derive black hole (BH) masses at progressively higher redshifts and over a wide range of continuum luminosities has become indispensable in the era of large-area extragalactic spectroscopic surveys. In this paper we present an extension of existing comparisons between rest-frame UV and optical virial BH mass estimators to intermediate redshifts and luminosities comparable to the local H$β$ reverberation mapped active galactic nuclei (AGN). We focus on the MgII, CIV, and CIII] broad emission lines and compare them to both H$α$ and H$β$. We use newly acquired near-infrared spectra from the FMOS instrument on the Subaru telescope for 89 broad-lined AGN at redshifts between 0.3 and 3.5, complemented by data from the AGES survey. We employ two different prescriptions for measuring the emission line widths and compare the results. We confirm that MgII shows a tight correlation with H$α$ and H$β$, with a scatter of ~0.25 dex. The CIV and CIII] estimators, while showing larger scatter, are viable virial mass estimators after accounting for a trend with the UV-to-optical luminosity ratio. We find an intrinsic scatter of ~0.37 dex between Balmer and carbon virial estimators by combining our dataset with previous high redshift measurements. This updated comparison spans a total of 3 decades in BH mass. We calculate a virial factor for CIV/CIII] logf(CIV/CIII])=0.87 with an estimated systematic uncertainty of ~0.4 dex and find excellent agreement between the local reverberation mapped AGN sample and our high-z sample.

preprint2015arXiv

The impact of environment and mergers on the HI content of galaxies in hydrodynamic simulations

We quantitatively examine the effects of accretion and environment on the HI content of galaxies within a cosmological hydrodynamic simulation that reproduces basic observed trends of HI in galaxies. We show that our model broadly reproduces the observed scatter in HI at a given stellar mass as quantified by the HI mass function in bins of stellar mass, as well as the HI richness versus local galaxy density. This shows that the predicted HI fluctuations and environmental effects are roughly consistent with data with few minor discrepancies. For satellite galaxies in >= 10^12M_* halos, the HI richness distribution is bimodal and drops towards the largest halo masses. The depletion rate of HI once a galaxy enters a more massive halo is more rapid at higher halo mass, in contrast to the specific star formation rate which shows much less variation in the attenuation rate versus halo mass. This suggests that, up to halo mass scales probed here (<= 10^14M_*), star formation is mainly attenuated by starvation, but HI is additionally removed by stripping once a hot gaseous halo is present. In low mass halos, the HI richness of satellites is independent of radius, while in high mass halos they become gas-poor towards the center. We show that the gas fraction of satellite and central galaxies decreases from z=5 to z=0,tracking each other until z~1 after which the satellites' HI content drops much more quickly, particularly for the highest halo masses. Mergers somewhat increase the HI richness and its scatter about the mean relation, tracking the metallicity in a way consistent with arising from inflow fluctuations, while star formation is significantly boosted relative to HI.

preprint2015arXiv

Torque-Limited Growth of Massive Black Holes in Galaxies Across Cosmic Time

We combine cosmological hydrodynamic simulations with analytic models to evaluate the role of galaxy-scale gravitational torques on the evolution of massive black holes at the centers of star-forming galaxies. We confirm and extend our earlier results to show that torque-limited growth yields black holes and host galaxies evolving on average along the Mbh-Mbulge relation from early times down to z = 0 and that convergence onto the scaling relation occurs independent of the initial conditions and with no need for mass averaging through mergers or additional self-regulation processes. Smooth accretion dominates the long-term evolution, with black hole mergers with mass ratios >1:5 representing typically a small fraction of the total growth. Winds from the accretion disk are required to eject significant mass to suppress black hole growth, but there is no need for coupling this wind to galactic-scale gas to regulate black holes in a non-linear feedback loop. Torque-limited growth yields a close-to-linear relation for the star formation rate and the black hole accretion rate averaged over galaxy evolution time scales. However, the SFR-AGN connection has significant scatter owing to strong variability of black hole accretion at all resolved time scales. Eddington ratios can be described by a broad lognormal distribution with median value evolving roughly as (1 + z)^1.9, suggesting a main sequence for black hole growth similar to the cosmic evolution of specific SFRs. Our results offer an attractive scenario consistent with available observations in which cosmological gas infall and transport of angular momentum in the galaxy by gravitational instabilities regulate the long-term co-evolution of black holes and star-forming galaxies.

preprint2014arXiv

Calcium-rich Gap Transients: Solving the Calcium Conundrum in the Intracluster Medium

X-ray measurements suggest the abundance of Calcium in the intracluster medium is higher than can be explained using favored models for core-collapse and Type Ia supernovae alone. We investigate whether the Calcium conundrum in the intracluster medium can be alleviated by including a contribution from the recently discovered subclass of supernovae known as Calcium-rich gap transients. Although the Calcium-rich gap transients make up only a small fraction of all supernovae events, we find that their high Calcium yields are sufficient to reproduce the X-ray measurements found for nearby rich clusters. We find the $χ^{2}$ goodness-of-fit metric improves from 84 to 2 by including this new class. Moreover, Calcium-rich supernovae preferentially occur in the outskirts of galaxies making it easier for the nucleosynthesis products of these events to be incorporated in the intracluster medium via ram-pressure stripping. The discovery of a Calcium-rich gap transients in clusters and groups far from any individual galaxy suggests supernovae associated with intracluster stars may play an important role in enriching the intracluster medium. Calcium-rich gap transients may also help explain anomalous Calcium abundances in many other astrophysical systems including individual stars in the Milky Way, the halos of nearby galaxies and the circumgalactic medium. Our work highlights the importance of considering the diversity of supernovae types and corresponding yields when modeling the abundance of the intracluster medium and other gas reservoirs.

preprint2014arXiv

The Circular Velocity Function of Group Galaxies

A robust prediction of $Λ{\rm CDM}$ cosmology is the halo circular velocity function (CVF), a dynamical cousin of the halo mass function. The correspondence between theoretical and observed CVFs is uncertain, however: cluster galaxies are reported to exhibit a power-law CVF consistent with $N$-body simulations, but that of the field is distinctly Schechter-like, flattened relative to $Λ{\rm CDM}$ expectations at circular velocities $v_c \lesssim 200\, {\rm km\, s^{-1}}$. Groups offer a powerful probe of the role environment plays in this discrepancy as they bridge the field and clusters. Here, we construct the CVF for a large, mass- and multiplicity-complete sample of group galaxies from the Sloan Digital Sky Survey. Using independent photometric $v_c$ estimators, we find no transition from a field- to $Λ{\rm CDM}$-shaped CVF above $v_c = 50\, {\rm km\, s^{-1}}$ as a function of group halo mass. All groups with $12.4 \lesssim \log M_{\rm halo} / M_{\odot} \lesssim 15.1$ (Local Group analogs to rich clusters) display similar Schechter-like CVFs marginally suppressed at low-$v_c$ compared to that of the field. Conversely, some agreement with $N$-body results emerges for samples saturated with late-type galaxies, with isolated late-types displaying a CVF similar in shape to $Λ{\rm CDM}$ predictions. We conclude that the flattening of the low-$v_c$ slope in groups is due to their depressed late-type fractions -- environment affecting the CVF only to the extent that it correlates with this quantity -- and that previous cluster analyses may suffer from interloper contamination. These results serve as useful benchmarks for cosmological simulations of galaxy formation.

preprint2014arXiv

The COS-Dwarfs Survey: The Carbon Reservoir Around sub-L* Galaxies

We report new observations of circumgalactic gas from the COS-Dwarfs survey, a systematic investigation of the gaseous halos around 43 low-mass z $\leq$ 0.1 galaxies using background QSOs observed with the Cosmic Origins Spectrograph. From the projected 1D and 2D distribution of C IV absorption, we find that C IV absorption is detected out to ~ 0.5 R$_{vir}$ of the host galaxies. The C IV absorption strength falls off radially as a power law and beyond 0.5 R$_{vir}$, no C IV absorption is detected above our sensitivity limit of ~ 50-100 m$Å$. We find a tentative correlation between detected C IV absorption strength and star formation, paralleling the strong correlation seen in highly ionized oxygen for L~L* galaxies by the COS-Halos survey. The data imply a large carbon reservoir in the CGM of these galaxies, corresponding to a minimum carbon mass of $\gtrsim$ 1.2$\times 10^6$ $M_\odot$ out to ~ 110 kpc. This mass is comparable to the carbon mass in the ISM and more than the carbon mass currently in stars of these galaxies. The C IV absorption seen around these sub-L* galaxies can account for almost two-thirds of all $W_r$> 100 m$Å$ C IV absorption detected at low z. Comparing the C IV covering fraction with hydrodynamical simulations, we find that an energy-driven wind model is consistent with the observations whereas a wind model of constant velocity fails to reproduce the CGM or the galaxy properties.

preprint2014arXiv

The Photon Underproduction Crisis

We examine the statistics of the low-redshift Lyman-alpha forest from smoothed particle hydrodynamic simulations in light of recent improvements in the estimated evolution of the cosmic ultraviolet background (UVB) and recent observations from the Cosmic Origins Spectrograph (COS). We find that the value of the metagalactic photoionization rate required by our simulations to match the observed properties of the low-redshift Lyman-alpha forest is a factor of 5 larger than the value predicted by state-of-the art models for the evolution of this quantity. This mismatch results in the mean flux decrement of the Lyman-alpha forest being underpredicted by at least a factor of 2 (a 10-sigma discrepancy with observations) and a column density distribution of Lyman-alpha forest absorbers systematically and significantly elevated compared to observations over nearly two decades in column density. We examine potential resolutions to this mismatch and find that either conventional sources of ionizing photons (galaxies and quasars) must be significantly elevated relative to current observational estimates or our theoretical understanding of the low-redshift universe is in need of substantial revision.

preprint2013arXiv

A Preliminary Calibration of the RR Lyrae Period-Luminosity Relation at Mid-Infrared Wavelengths: WISE Data

Using time-resolved, mid-infrared data from WISE and geometric parallaxes from HST for four Galactic RR Lyrae variables, we derive the following Population II Period-Luminosity (PL) relations for the WISE [W1], [W2] and [W3] bands at 3.4, 4.6 & 12 um, respectively: M[W1] = -2.44 (+/- 0.95) x logP - 1.26 (+/- 0.25) sigma = 0.10 M[W2] = -2.55 (+/- 0.89) x logP - 1.29 (+/- 0.23) sigma = 0.10 M[W3] = -2.58 (+/- 0.97) x logP - 1.32 (+/- 0.25) sigma = 0.10 The slopes and the scatter around the fits are consistent with a smooth extrapolation of those same quantities from previously-published K-band observations at 2.2 um, where the asymptotic (long-wavelength) behavior is consistent with a Period-Radius relation having a slope of 0.5. No obvious correlation with metallicity (spanning 0.4 dex in [Fe/H]) is found in the residuals of the four calibrating RR Lyrae stars about the mean PL regression line.

preprint2013arXiv

Hydrogen and Metal Line Absorption Around Low-Redshift Galaxies in Cosmological Hydrodynamic Simulations

We study the physical conditions of the circum-galactic medium (CGM) around z=0.25 galaxies as traced by HI and metal line absorption, using cosmological hydrodynamic simulations that include galactic outflows. Using lines of sight targeted at impact parameters from 10 kpc to 1 Mpc around galaxies with halo masses from 10^11-10^13 M_solar, we study the physical conditions and their variation with impact parameter b and line-of-sight velocity delta v in the CGM as traced by HI, MgII, SiIV, CIV, OVI, and NeVIII absorbers. All ions show a strong excess of absorption near galaxies compared to random lines of sight. The excess continues beyond 1 Mpc, reflecting the correlation of metal absorption with large-scale structure. Absorption is particularly enhanced within about v<300 km/sec and roughly 300 kpc of galaxies (with distances somewhat larger for the highest ion), approximately delineating the CGM; this range contains the majority of global metal absorption. Low ions like MgII and SiIV predominantly arise in denser gas closer to galaxies and drop more rapidly with b, while high ions OVI and NeVIII trace more diffusely distributed gas with a comparatively flat radial profile; CIV is intermediate. All ions predominantly trace T~10^4-4.5 K photo-ionised gas at all b, but when hot CGM gas is present (mostly in larger halos), we see strong collisionally-ionised OVI and NeVIII at b <= 100 kpc. Larger halo masses generally produce more absorption, though overall the trends are not as strong as that with impact parameter. These findings arise using our favoured outflow scalings as expected for momentum-driven winds; with no winds, the CGM gas remains mostly unenriched, while our outflow model with a constant velocity and mass loading factor produce hotter, more widely dispersed metals.

preprint2013arXiv

The Neutral Hydrogen Content of Galaxies in Cosmological Hydrodynamic Simulations

We examine the global HI properties of galaxies in quarter-billion particle cosmological simulations using Gadget-2, focusing on how galactic outflows impact HI content. We consider four outflow models, including a new one (ezw) motivated by recent interstellar medium simulations in which the wind speed and mass loading factor scale as expected for momentum-driven outflows for larger galaxies and energy-driven outflows for dwarfs (sigma<75 km/s). To obtain predicted HI masses, we employ a simple but effective local correction for particle self-shielding, and an observationally-constrained transition from neutral to molecular hydrogen. Our ezw simulation produces an HI mass function whose faint-end slope of -1.3 agrees well with observations from the ALFALFA survey; other models agree less well. Satellite galaxies have a bimodal distribution in HI fraction versus halo mass, with smaller satellites and/or those in larger halos more often being HI-deficient. At a given stellar mass, HI content correlates with star formation rate and inversely correlates with metallicity, as expected if driven by stochasticity in the accretion rate. To higher redshifts, massive HI galaxies disappear and the mass function steepens. The global cosmic HI density conspires to remain fairly constant from z~5-0, but the relative contribution from smaller galaxies increases with redshift.

preprint2013arXiv

The Spectrally Resolved Lyman-alpha Emission of Three Lyman-alpha Selected Field Galaxies at z~2.4 from the HETDEX Pilot Survey

We present new results on the spectrally resolved Lyman-alpha (LyA) emission of three LyA emitting field galaxies at z~2.4 with high LyA equivalent width (>100 Angstroms) and LyA luminosity (~10^43 erg/s). At 120 km/s (FWHM) spectral resolution, the prominent double-peaked LyA profile straddles the systemic velocity, where the velocity zero-point is determined from spectroscopy of the galaxies' rest-frame optical nebular emission lines. The average velocity offset from systemic of the stronger redshifted emission component for our sample is 176 km/s while the average total separation between the redshifted and main blueshifted emission components is 380 km/s. These measurements are a factor of ~2 smaller than for UV continuum-selected galaxies that show LyA in emission with lower LyA equivalent width. We compare our LyA spectra to the predicted line profiles of a spherical "expanding shell" LyA radiative transfer grid that models large-scale galaxy outflows. Specifically blueward of the systemic velocity where two galaxies show a weak, highly blueshifted (by ~1000 km/s) tertiary emission peak, the model line profiles are a relatively poor representation of the observed spectra. Since the neutral gas column density has a dominant influence over the shape of the LyA line profile, we caution against equating the observed LyA velocity offset with a physical outflow velocity, especially at lower spectral resolution where the unresolved LyA velocity offset is a convoluted function of several degenerate parameters. Referring to rest-frame ultraviolet and optical Hubble Space Telescope imaging, we find that galaxy-galaxy interactions may play an important role in inducing a starburst that results in copious LyA emission, as well as perturbing the gas distribution and velocity field which have strong influence over the LyA emission line profile.

preprint2013arXiv

Tracing Inflows and Outflows with Absorption Lines in Circumgalactic Gas

We examine how HI and metal absorption lines within low-redshift galaxy halos trace the dynamical state of circumgalactic gas, using cosmological hydrodynamic simulations that include a well-vetted heuristic model for galactic outflows. We categorize inflowing, outflowing, and ambient gas based on its history and fate as tracked in our simulation. Following our earlier work showing that the ionisation level of absorbers was a primary factor in determining the physical conditions of absorbing gas, we show here that it is also a governing factor for its dynamical state. Low-ionisation metal absorbers (e.g. MgII) tend to arise in gas that will fall onto galaxies within several Gyr, while high-ionisation metal absorbers (e.g. OVI) generally trace material that was deposited by outflows many Gyr ago. Inflowing gas is dominated by enriched material that was previously ejected in an outflow, hence accretion at low redshifts is typically substantially enriched. Recycling wind material is preferentially found closer to galaxies, and is more dominant in lower-mass halos since high-mass halos have more hot gas that is able to support itself against infall. Low-mass halos also tend to re-eject more of their accreted material, owing to our outflow prescription that employs higher mass loading factors for lower-mass galaxies. Typical HI absorbers trace unenriched ambient material that is not participating in the baryon cycle, but stronger HI absorbers arise in cool, enriched inflowing gas. Instantaneous radial velocity measures of absorbers are generally poor at distinguishing between inflowing and outflowing gas, except in the case of very recent outflows. These results suggest that probing halo gas using a range of absorbers can provide detailed information about the amount and physical conditions of material that is participating in the baryon cycle.

preprint2012arXiv

The Absolute Magnitude of RRc Variables From Statistical Parallax

We present the first definitive measurement of the absolute magnitude of RR Lyrae c-type variable stars (RRc) determined purely from statistical parallax. We use a sample of 247 RRc selected from the All Sky Automated Survey (ASAS) for which high-quality light curves, photometry and proper motions are available. We obtain high-resolution echelle spectra for these objects to determine radial velocities and abundances as part of the Carnegie RR Lyrae Survey (CARRS). We find that M_(V,RRc) = 0.52 +/- 0.11 at a mean metallicity of [Fe/H] = -1.59. This is to be compared with previous estimates for RRab stars (M_(V,RRab) = 0.75 +/- 0.13 and the only direct measurement of an RRc absolute magnitude (RZ Cephei, M_(V, RRc) = 0.27 +/- 0.17). We find the bulk velocity of the halo to be (W_pi, W_theta, W_z) = (10.9,34.9,7.2) km/s in the radial, rotational and vertical directions with dispersions (sigma_(W_pi), sigma_(W_theta), sigma_(W_z)) = (154.7, 103.6, 93.8) km/s. For the disk, we find (W_pi, W_theta, W_z) = (8.5, 213.2, -22.1) km/s with dispersions (sigma_(W_pi), sigma_(W_theta), sigma_(W_z)) = (63.5, 49.6, 51.3) km/s. Finally, we suggest that UCAC2 proper motion errors may be overestimated by about 25%

preprint2012arXiv

Warm-hot gas in groups and galaxies toward H2356-309

We present a detailed analysis of the galaxy and group distributions around three reported X-ray absorption line systems in the spectrum of the quasar H2356-309. Previous studies associated these absorbers with known large-scale galaxy structures (i.e., walls and filaments) along the line of sight. Such absorption lines typically trace 10^{5-7} K gas, and may be evidence of the elusive warm-hot intergalactic medium (WHIM) thought to harbor the bulk of the low-redshift "missing baryons;" alternatively, they may be linked to individual galaxies or groups in the filaments. Here we combine existing galaxy survey data with new, multi-object Magellan spectroscopy to investigate the detailed galaxy distribution near each absorber. All of these three absorption systems are within the projected virial radii of nearby galaxies and/or groups, and could therefore arise in these virialized structures rather than (or in addition to) the WHIM. However, we find no additional galaxies near a fourth "void" absorber recently found in the spectrum, suggesting that this system may indeed trace gas unassociated with any individual halo. Though the number of known systems is still small, spatial coincidences suggest that some X-ray absorbers lie in galaxy and/or group environments, though others could still trace the large-scale filamentary WHIM gas predicted by simulations.

preprint2011arXiv

SDSS-III: Massive Spectroscopic Surveys of the Distant Universe, the Milky Way Galaxy, and Extra-Solar Planetary Systems

Building on the legacy of the Sloan Digital Sky Survey (SDSS-I and II), SDSS-III is a program of four spectroscopic surveys on three scientific themes: dark energy and cosmological parameters, the history and structure of the Milky Way, and the population of giant planets around other stars. In keeping with SDSS tradition, SDSS-III will provide regular public releases of all its data, beginning with SDSS DR8 (which occurred in Jan 2011). This paper presents an overview of the four SDSS-III surveys. BOSS will measure redshifts of 1.5 million massive galaxies and Lya forest spectra of 150,000 quasars, using the BAO feature of large scale structure to obtain percent-level determinations of the distance scale and Hubble expansion rate at z<0.7 and at z~2.5. SEGUE-2, which is now completed, measured medium-resolution (R=1800) optical spectra of 118,000 stars in a variety of target categories, probing chemical evolution, stellar kinematics and substructure, and the mass profile of the dark matter halo from the solar neighborhood to distances of 100 kpc. APOGEE will obtain high-resolution (R~30,000), high signal-to-noise (S/N>100 per resolution element), H-band (1.51-1.70 micron) spectra of 10^5 evolved, late-type stars, measuring separate abundances for ~15 elements per star and creating the first high-precision spectroscopic survey of all Galactic stellar populations (bulge, bar, disks, halo) with a uniform set of stellar tracers and spectral diagnostics. MARVELS will monitor radial velocities of more than 8000 FGK stars with the sensitivity and cadence (10-40 m/s, ~24 visits per star) needed to detect giant planets with periods up to two years, providing an unprecedented data set for understanding the formation and dynamical evolution of giant planet systems. (Abridged)

preprint2011arXiv

The Eighth Data Release of the Sloan Digital Sky Survey: First Data from SDSS-III

The Sloan Digital Sky Survey (SDSS) started a new phase in August 2008, with new instrumentation and new surveys focused on Galactic structure and chemical evolution, measurements of the baryon oscillation feature in the clustering of galaxies and the quasar Ly alpha forest, and a radial velocity search for planets around ~8000 stars. This paper describes the first data release of SDSS-III (and the eighth counting from the beginning of the SDSS). The release includes five-band imaging of roughly 5200 deg^2 in the Southern Galactic Cap, bringing the total footprint of the SDSS imaging to 14,555 deg^2, or over a third of the Celestial Sphere. All the imaging data have been reprocessed with an improved sky-subtraction algorithm and a final, self-consistent photometric recalibration and flat-field determination. This release also includes all data from the second phase of the Sloan Extension for Galactic Understanding and Evolution (SEGUE-2), consisting of spectroscopy of approximately 118,000 stars at both high and low Galactic latitudes. All the more than half a million stellar spectra obtained with the SDSS spectrograph have been reprocessed through an improved stellar parameters pipeline, which has better determination of metallicity for high metallicity stars.

preprint2011arXiv

The intergalactic medium over the last 10 billion years II: Metal-line absorption and physical conditions

We investigate the metallicity evolution and content of the intergalactic medium (IGM) and galactic halo gas from z=2->0 using 110-million particle cosmological hydrodynamic simulations. We focus on the detectability and physical properties of UV resonance metal-line absorbers observable with Hubble's Cosmic Origins Spectrograph (COS). We confirm that galactic superwind outflows are required to enrich the IGM to observed levels down to z=0 using three wind prescriptions contrasted to a no-wind simulation. Our favoured momentum-conserved wind prescription deposits metals closer to galaxies owing to its moderate energy input, while the more energetic constant wind model enriches the warm-hot IGM 6.4x more. Despite these significant differences, all wind models produce metal-line statistics within a factor of two of existing observations. This is because OVI, CIV, SiIV, and NeVIII absorbers primarily arise from T<10^5 K, photo-ionised gas that is enriched to similar levels in the three feedback schemes. OVI absorbers trace the diffuse phase with delta<100, which is enriched to ~1/50 Zsol at z=0, although the absorbers themselves usually exceed 0.3 Zsol and arise from inhomogeneously distributed, un-mixed winds. CIV and SiIV absorbers trace primarily T~10^4 K gas inside haloes. We predict COS will observe a population of NeVIII photo-ionised absorbers tracing T<10^5 K, delta~10 gas. MgX and SiXII are rarely detected in COS S/N=30 simulated sight lines although detections trace T=10^(6-7) K halo gas. In general, the IGM is enriched in an outside-in manner, where wind-blown metals released at higher redshift reach lower overdensities, resulting in higher ionisation species tracing lower-density, older metals. At z=0, the 90% of baryons outside of galaxies are enriched to 0.096 Zsol, but the 65% of unbound baryons in the IGM have 0.018 Zsol and contain only 4% of all metals.

preprint2011arXiv

The Kinematics of Multiple-Peaked Ly-alpha Emission in Star-Forming Galaxies at z~2-3

We present new results on the Ly-alpha emission-line kinematics of 18 z~2-3 star-forming galaxies with multiple-peaked Ly-alpha profiles. With our large spectroscopic database of UV-selected star-forming galaxies at these redshifts, we have determined that ~30% of such objects with detectable Ly-alpha emission display multiple-peaked emission profiles. These profiles provide additional constraints on the escape of Ly-alpha photons due to the rich velocity structure in the emergent line. Despite recent advances in modeling the escape of Ly-alpha from star-forming galaxies at high redshifts, comparisons between models and data are often missing crucial observational information. Using Keck II NIRSPEC spectra of H-alpha (z~2) and [OIII] 5007 (z~3), we have measured accurate systemic redshifts, rest-frame optical nebular velocity dispersions and emission-line fluxes for the objects in the sample. Accurate systemic redshifts allow us to translate the multiple-peaked Ly-alpha profiles into velocity space, revealing that the majority (11/18) display double-peaked emission straddling the velocity-field zeropoint with stronger red-side emission. We compare our data with the predictions of simple models for outflowing and infalling gas distributions around high-redshift galaxies. While popular "shell" models provide a qualitative match with many of the observations of Ly-alpha emission, we find that in detail there are important discrepancies between the models and data, as well as problems with applying the framework of an expanding thin shell of gas to explain high-redshift galaxy spectra. Our data highlight these inconsistencies, as well as illuminating critical elements for success in future models of outflow and infall in high-redshift galaxies. [Abridged]

preprint2010arXiv

Discovery of a large-scale galaxy filament near a candidate intergalactic X-ray absorption system

We present an analysis of the large-scale galaxy distribution around two possible warm-hot intergalactic medium (WHIM) absorption systems reported along the Markarian 421 sightline. Using the Sloan Digital Sky Survey, we find a prominent galaxy filament at the redshift of the z=0.027 X-ray absorption line system. The filament exhibits a width of 3.2 Mpc and length of at least 20 Mpc, comparable to the size of WHIM filaments seen in cosmological simulations. No individual galaxies fall within 350 projected kpc so it is unlikely that the absorption is associated with gas in a galaxy halo or outflow. Another, lower-significance X-ray absorption system was reported in the same Chandra spectrum at z=0.011, but the large-scale structure in its vicinity is far weaker and may be a spurious alignment. By searching for similar galaxy structures in 140 random smoothed SDSS fields, we estimate a ~5-10% probability of the z=0.027 absorber-filament alignment occurring by chance. If these two systems are indeed physically associated, this would represent the first known coincidence between large-scale galaxy structure and a blind X-ray WHIM detection.

preprint2010arXiv

Feedback and Recycled Wind Accretion: Assembling the z=0 Galaxy Mass Function

We analyse cosmological hydrodynamic simulations that include observationally-constrained prescriptions for galactic outflows. If these simulated winds accurately represent winds in the real Universe, then material previously ejected in winds provides the dominant source of gas infall for new star formation at redshifts z<1. This recycled wind accretion, or wind mode, provides a third physically distinct accretion channel in addition to the "hot" and "cold" modes emphasised in recent theoretical studies. Because of the interaction between outflows and gas in and around halos, the recycling timescale of wind material (t_rec) is shorter in higher-mass systems, which reside in denser gaseous environments. In these simulations, this differential recycling plays a central role in shaping the present-day galaxy stellar mass function (GSMF). If we remove all particles that were ever ejected in a wind, then the predicted GSMFs are much steeper than observed; galaxy masses are suppressed both by the direct removal of gas and by the hydrodynamic heating of their surroundings, which reduces subsequent infall. With wind recycling included, the simulation that incorporates our favoured momentum-driven wind scalings reproduces the observed GSMF for stellar masses 10^9 < M < 5x10^10 Msolar. At higher masses, wind recycling leads to excessive galaxy masses and excessive star formation rates relative to observations. In these massive systems, some quenching mechanism must suppress the re-accretion of gas ejected from star-forming galaxies. In short, as has long been anticipated, the form of the GSMF is governed by outflows; the unexpected twist here for our simulated winds is that it is not primarily the ejection of material but how the ejected material is re-accreted that governs the GSMF.

preprint2010arXiv

The intergalactic medium over the last 10 billion years I: Lyman alpha absorption and physical conditions

The intergalactic medium (IGM) is the dominant reservoir of baryons at all cosmic epochs. We investigate the evolution of the IGM from z=2-0 in 48 Mpc/h, 110-million particle cosmological hydrodynamic simulations using three prescriptions for galactic outflows. We focus on the evolution of IGM physical properties, and how such properties are traced by Ly-alpha absorption as detectable using HST/COS. Our results broadly confirm the canonical picture that most Ly-alpha absorbers arise from highly ionized gas tracing filamentary large-scale structure. Growth of structure causes gas to move from the diffuse photoionized IGM into other cosmic phases, namely stars, cold and hot gas within galaxy halos, and the unbound and shock-heated warm-hot intergalactic medium (WHIM). By today, baryons are roughly equally divided between bound phases (35%), the diffuse IGM (41%), and the WHIM (24%). Here we (re)define the WHIM as gas with overdensities lower than that in halos and temperatures >10^5 K, in order to more closely align it with "missing baryons". When we tune our photoionizing background to match the observed evolution of the Ly-alpha mean flux decrement, we obtain a line count evolution that broadly agrees with available data. We predict a column density distribution slope of -1.70 for our favored momentum-driven wind model, in agreement with recent observations, and it becomes shallower with redshift. With improved statistics, the frequency of strong lines can be a valuable diagnostic of outflows, and our favored wind model matches existing data best among our models. The relationship between column density and physical density is fairly tight from z=2-0, and evolves as rho N_HI^0.74 10^(-0.37z) for diffuse absorbers. Linewidths only loosely reflect the temperature of the absorbing gas, which will hamper attempts to quantify the WHIM using broad Ly-alpha absorbers. [Abridged]

preprint2010arXiv

The star capture model for fueling quasar accretion disks

Although the powering mechanism for quasars is now widely recognized to be the accretion of matter in a geometrically thin disk, the transport of matter to the inner region of the disk where luminosity is emitted remains an unsolved question. Miralda-Escudé & Kollmeier (2005) proposed a model whereby quasars are fuelled when stars are captured by the accretion disk as they plunge through the gas. Such plunging stars can then be destroyed and deliver their mass to the accretion disk. Here we present the first detailed calculations for the capture of stars originating far from the accretion disk near the zone of influence of the central black hole. In particular we examine the effect of adding a perturbing mass to a fixed stellar cusp potential on bringing stars into the accretion disk where they can be captured. The work presented here will be discussed in detail in an upcoming publication Kennedy et al. (2010).

preprint2009arXiv

Lyman-alpha Emission From Cosmic Structure I: Fluorescence

We present predictions for the fluorescent Lyman-alpha emission signature arising from photoionized, optically thick structures in Smoothed Particle Hydrodynamic (SPH) cosmological simulations of a Lambda-CDM universe using a Monte Carlo Lyman-alpha radiative transfer code. We calculate the expected Lyman-alpha image and 2-dimensional spectra for gas exposed to a uniform ultraviolet ionizing background as well as gas exposed additionally to the photoionizing radiation from a local quasar, after correcting for the self-shielding of hydrogen. As a test of our numerical methods and for application to current observations, we examine simplified analytic structures that are uniformly or anisotropically illuminated. We compare these results with recent observations. We discuss future observing campaigns on large telescopes and realistic strategies for detecting fluorescence owing to the ambient metagalactic ionization and in regions close to bright quasars. While it will take hundreds of hours on the current generation of telescopes to detect fluorescence caused by the ultraviolet background (UVB) alone, our calculations suggest that of order ten sources of quasar-induced fluorescent Lyman-alpha emission should be detectable after a 10 hour exposure in a 10 arcmin^2 field around a bright quasar. These observations will help probe the physical conditions in the densest regions of the intergalactic medium as well as the temporal light curves and isotropy of quasar radiation.

preprint2009arXiv

Spectroscopic Confirmation of the Pisces Overdensity

We present spectroscopic confirmation of the "Pisces Overdensity", also known as "Structure J", a photometric overdensity of RR Lyrae stars discovered by the Sloan Digital Sky Survey (SDSS) at an estimated photometric distance of ~85kpc. We measure radial velocities for 8 RR Lyrae stars within Pisces. We find that 5 of the 8 stars have heliocentric radial velocities within a narrow range of -87 km/s < v < -67 km/s, suggesting that the photometric overdensity is mainly due to a physically associated system, probably a dwarf galaxy or a disrupted galaxy. Two of the remaining 3 stars differ from one another by only 9 km/s, but it would be premature to identify them as a second system.

preprint2008arXiv

Old-Population Hypervelocity Stars from the Galactic Center: Limits from the SDSS

We present limits on the ejection of old-population HVS from a sample of over 290,000 stars selected from the SDSS. We derive the speed at the solar circle from the measured positions and radial velocities by assuming a radial orbit and adopting a simple isothermal model of the Galactic halo, which enables us to identify candidate bound and unbound ejectees. We find 4 candidate bound F-stars from this sample, all with negative Galactocentric radial velocity (i.e., returning toward the GC). We additionally find 2 candidate unbound stars (one F and one G), however, existing proper motion measurements make these unlikely to be emerging from the GC. These data place an upper limit on the rate of ejection of old-population stars from the GC of ~45/Myr. Comparing to the rate for more massive B-star ejectees of ~0.5/Myr, our limit on the rate of ejection of old-population HVS shows that the mass function at the GC is not bottom-heavy and is consistent with being normal. Future targeted surveys of old-population HVS could determine if it is indeed top-heavy.

preprint2007arXiv

Where are the Old-Population High Velocity Stars?

To date, all of the reported high velocity stars (HVSs), which are believed to be ejected from the Galactic center, are blue and therefore almost certainly young. Old-population HVSs could be much more numerous than the young ones that have been discovered, but still have escaped detection because they are hidden in a much denser background of Galactic halo stars. Discovery of these stars would shed light on star formation at the Galactic center, the mechanism by which they are ejected from it, and, if they prove numerous, enable detailed studies of the structure of the dark halo. We analyze the problem of finding these stars and show that the search should be concentrated around the main-sequence turnoff (0.3<g-i<1.1) at relatively faint magnitudes (19.5<g<21.5). If the ratio of turnoff stars to B stars is the same for HVSs as it is in the local disk, such a search would yield about 1 old-population HVS per 10 deg^2. A telescope similar to the Sloan 2.5m could search about 20 deg^2 per night, implying that in short order such a population, should it exist, would show up in interesting numbers.

preprint2006arXiv

Black Hole Masses and Eddington Ratios at 0.3<z<4

We study the distribution of Eddington luminosity ratios, L_bol/L_edd, of active galactic nuclei (AGNs) discovered in the AGN and Galaxy Evolution Survey (AGES). We combine H-beta, MgII, and CIV line widths with continuum luminosities to estimate black hole (BH) masses in 407 AGNs, covering the redshift range z~0.3-4 and the bolometric luminosity range L_bol~10^45-10^47 erg/s. The sample consists of X-ray or mid-infrared (24 micron) point sources with optical magnitude R<=21.5 mag and optical emission line spectra characteristic of AGNs. For the range of luminosity and redshift probed by AGES, the distribution of estimated Eddington ratios is well described as log-normal with a peak at L_bol/L_edd ~ 1/4 and a dispersion of 0.3 dex. Since additional sources of scatter are minimal, this dispersion must account for contributions from the scatter between estimated and true BH mass and the scatter between estimated and true bolometric luminosity. Therefore, we conclude that: (1) neither of these sources of error can contribute more than \~0.3 dex rms; and (2) the true Eddington ratios of optically luminous AGNs are even more sharply peaked. Because the mass estimation errors must be smaller than ~0.3 dex, we can also investigate the distribution of Eddington ratios at fixed BH mass. We show for the first time that the distribution of Eddington ratios at fixed BH mass is peaked, and that the dearth of AGNs at a factor ~10 below Eddington is real and not an artifact of sample selection. These results provide strong evidence that supermassive BHs gain most of their mass while radiating close to the Eddington limit, and they suggest that the fueling rates in luminous AGNs are ultimately determined by BH self-regulation of the accretion flow rather than galactic scale dynamical disturbances.

preprint2004arXiv

Star Captures by Quasar Accretion Disks: A Possible Explanation of the M-sigma Relation

A new theory of quasars is presented in which the matter of thin accretion disks around black holes is supplied by stars that plunge through the disk. Stars in the central part of the host galaxy are randomly perturbed to highly radial orbits, and as they repeatedly cross the disk they lose orbital energy by drag, eventually merging into the disk. Requiring the rate of stellar mass capture to equal the mass accretion rate into the black hole, a relation between the black hole mass and the stellar velocity dispersion is predicted of the form M_{BH} \propto sigma_*^{30/7}. The normalization depends on various uncertain parameters such as the disk viscosity, but is consistent with observation for reasonable assumptions. We show that a seed central black hole in a newly formed stellar system can grow at the Eddington rate up to this predicted mass via stellar captures by the accretion disk. Once this mass is reached, star captures are insufficient to maintain an Eddington accretion rate, and the quasar may naturally turn off as the accretion switches to a low-efficiency advection mode. The model provides a mechanism to deliver mass to the accretion disk at small radius, probably solving the problem of gravitational instability to star formation in the disk at large radius. We note that the matter from stars that is incorporated to the disk has an average specific angular momentum that is very small or opposite to that of the disk, and discuss how a rotating disk may be maintained as it captures this matter if a small fraction of the accreted mass comes from stellar winds that form a disk extending to larger radius. We propose several observational tests and consequences of this theory.

preprint2003arXiv

Proper Motion Catalog From SDSS \intersect USNO-B

We present a catalog of 345,000 stars with proper motions μ> 20 mas/yr and magnitudes r' <20, drawn from the intersection of the SDSS Data Release One and USNO-B. We first use SDSS quasars to correct the proper motions in each of the two source catalogs, finding that these corrections for systematics are comparable to the statistical errors. The combined proper motions have errors of σ_μ=4.0 mas/yr, significantly better than any previous catalog at these magnitudes. While cross-correlating SDSS and USNO-B removes the vast majority of the very large number of spurious proper-motion stars in each, some contamination remains, especially at faint magnitudes, r' > 19.5, and high proper motions. We present a diagrammatic method to estimate the level of contamination in selected subsets of the catalog, which should be useful for choosing appropriate selection criteria for various applications.