We present DES16C3cje, a low-luminosity, long-lived type II supernova (SN II) at redshift 0.0618, detected by the Dark Energy Survey (DES). DES16C3cje is a unique SN. The spectra are characterized by extremely narrow photospheric lines corresponding to very low expansion velocities of $\lesssim1500$ km s$^{-1}$, and the light curve shows an initial peak that fades after 50 days before slowly rebrightening over a further 100 days to reach an absolute brightness of M$_r\sim -15.5$ mag. The decline rate of the late-time light curve is then slower than that expected from the powering by radioactive decay of $^{56}$Co but is comparable to that expected from accretion power. Comparing the bolometric light curve with hydrodynamical models, we find that DES16C3cje can be explained by either i) a low explosion energy (0.11 foe) and relatively large $^{56}$Ni production of 0.075 M$_{\odot}$ from a $\sim15$ M$_{\odot}$ red supergiant progenitor typical of other SNe II, or ii) a relatively compact $\sim40$ M$_{\odot}$ star, explosion energy of 1 foe, and 0.08 M$_{\odot}$ of $^{56}$Ni. Both scenarios require additional energy input to explain the late-time light curve, which is consistent with fallback accretion at a rate of $\sim0.5\times{10^{-8}}$ M$_{\odot}$ s$^{-1}$.

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).

We introduce LATIS, the Ly$\alpha$ 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$\alpha$ 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.

We test and improve the numerical schemes in our smoothed particle hydrodynamics (SPH) code for cosmological simulations, including the pressure-entropy formulation (PESPH), a time-dependent artificial viscosity, a refined timestep criterion, and metal-line cooling that accounts for photoionisation in the presence of a recently refined Haardt \& Madau (2012) model of the ionising background. The PESPH algorithm effectively removes the artificial surface tension present in the traditional SPH formulation, and in our test simulations it produces better qualitative agreement with mesh-code results for Kelvin-Helmholtz instability and cold cloud disruption. Using a set of cosmological simulations, we examine many of the quantities we have studied in previous work. Results for galaxy stellar and HI mass functions, star formation histories, galaxy scaling relations, and statistics of the Ly$\alpha$ forest are robust to the changes in numerics and microphysics. As in our previous simulations, cold gas accretion dominates the growth of high-redshift galaxies and of low mass galaxies at low redshift, and recycling of winds dominates the growth of massive galaxies at low redshift. However, the PESPH simulation removes spurious cold clumps seen in our earlier simulations, and the accretion rate of hot gas increases by up to an order of magnitude at some redshifts. The new numerical model also influences the distribution of metals among gas phases, leading to considerable differences in the statistics of some metal absorption lines, most notably NeVIII.

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 $\eta$ with velocity dispersion $\sigma$. 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 $\eta \sim \sigma^5$ in low mass haloes with $\sigma < 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.

The past two decades have seen a tremendous investment in observational facilities that promise to reveal new and unprecedented discoveries about the universe. In comparison, the investment in theoretical work is completely dwarfed, even though theory plays a crucial role in the interpretation of these observations, predicting new types of phenomena, and informing observing strategies. In this white paper, we argue that in order to reach the promised critical breakthroughs in astrophysics over the next decade and well beyond, the national agencies must take a serious approach to investment in theoretical astrophysics research. We discuss the role of theory in shaping our understanding of the universe, and then we provide a multi-level strategy, from the grassroots to the national, to address the current underinvestment in theory relative to observational work.

We present optical photometry and spectroscopy of the Type Ia supernova SN2018cqj/ATLAS18qtd. The supernova exploded in an isolated region at $\sim 65$~kpc from the S0 galaxy IC~550 at $z=0.0165$ ($D\approx 74$~Mpc) and has a redshift consistent with a physical association to this galaxy. Multicolor photometry show that SN2018cqj/ATLAS18qtd is a low-luminosity ($M_{B_{max}}\approx -17.9$ mag), fast-declining Type Ia with color stretch $s_{BV} \approx 0.6$ and $B$-band decline rate $\Delta m_{15}(B) \approx 1.77$ mag. Two nebular-phase spectra obtained as part of the 100IAS survey at +193 and +307 days after peak show the clear detection of a narrow H$\alpha$ line in emission that is resolved in the first spectrum with $\rm FWHM \approx 1200$ km s$^{-1}$ and $L_{H\alpha} \approx 3.8\times 10^{37}$ erg s$^{-1}$. The detection of a resolved H$\alpha$ line with a declining luminosity is broadly consistent with recent models where hydrogen is stripped from the non-degenerate companion in a single-degenerate progenitor system. However, the amount of hydrogen consistent with the luminosities of the H$\alpha$ line would be $\sim 10^{-3}$ M$_{\odot}$, significantly less than theoretical model predictions in the classical single-degenerate progenitor systems. SN2018cqj/ATLAS18qtd is the second low-luminosity, fast-declining Type Ia SN after SN2018fhw/ASASSN-18tb that shows narrow H$\alpha$ in emission in its nebular-phase spectra.

MegaMapper is a proposed ground-based experiment to measure Inflation parameters and Dark Energy from galaxy redshifts at 2<z<5. A 6.5-m Magellan telescope will be coupled with DESI spectrographs to achieve multiplexing of 20,000. MegaMapper would be located at Las Campanas Observatory to fully access LSST imaging for target selection.

We report comprehensive multi-wavelength observations of a peculiar Type Ia-like supernova ("SN Ia-pec") ASASSN-15pz. ASASSN-15pz is a spectroscopic "twin" of SN 2009dc, a so-called "Super-Chandrasekhar-mass" SN, throughout its evolution, but it has a peak luminosity M_B,peak = -19.69 +/- 0.12 mag that is \approx 0.6 mag dimmer and comparable to the SN 1991T sub-class of SNe Ia at the luminous end of the normal width-luminosity relation. The synthesized Ni56 mass of M_Ni56 = 1.13 +/- 0.14 M_sun is also substantially less than that found for several 2009dc-like SNe. Previous well-studied 2009dc-like SNe have generally suffered from large and uncertain amounts of host-galaxy extinction, which is negligible for ASASSN-15pz. Based on the color of ASASSN-15pz, we estimate a host extinction for SN 2009dc of E(B-V)_host=0.12 mag and confirm its high luminosity (M_B, peak[2009dc] \approx -20.3 mag). The 2009dc-like SN population, which represents ~1% of SNe Ia, exhibits a range of peak luminosities, and do not fit onto the tight width-luminosity relation. Their optical light curves also show significant diversity of late-time (>~ 50 days) decline rates. The nebular-phase spectra provide powerful diagnostics to identify the 2009dc-like events as a distinct class of SNe Ia. We suggest referring to these sources using the phenomenology-based "2009dc-like SN Ia-pec" instead of "Super-Chandrasekhar SN Ia," which is based on an uncertain theoretical interpretation.

The next decade affords tremendous opportunity to achieve the goals of Galactic archaeology. That is, to reconstruct the evolutionary narrative of the Milky Way, based on the empirical data that describes its current morphological, dynamical, temporal and chemical structures. Here, we describe a path to achieving this goal. The critical observational objective is a Galaxy-scale, contiguous, comprehensive mapping of the disk's phase space, tracing where the majority of the stellar mass resides. An ensemble of recent, ongoing, and imminent surveys are working to deliver such a transformative stellar map. Once this empirical description of the dust-obscured disk is assembled, we will no longer be operationally limited by the observational data. The primary and significant challenge within stellar astronomy and Galactic archaeology will then be in fully utilizing these data. We outline the next-decade framework for obtaining and then realizing the potential of the data to chart the Galactic disk via its stars. One way to support the investment in the massive data assemblage will be to establish a Galactic Archaeology Consortium across the ensemble of stellar missions. This would reflect a long-term commitment to build and support a network of personnel in a dedicated effort to aggregate, engineer, and transform stellar measurements into a comprehensive perspective of our Galaxy.

The problem of the origin of the elements is a fundamental one in astronomy and one that has many open questions. Prominent examples include (1) the nature of Type Ia supernovae and the timescale of their contributions; (2) the observational identification of elements such as titanium and potassium with the $\alpha$-elements in conflict with core-collapse supernova predictions; (3) the number and relative importance of r-process sites; (4) the origin of carbon and nitrogen and the influence of mixing and mass loss in winds; and (5) the origin of the intermediate elements, such as Cu, Ge, As, and Se, that bridge the region between charged-particle and neutron-capture reactions. The next decade will bring to maturity many of the new tools that have recently made their mark, such as large-scale chemical cartography of the Milky Way and its satellites, the addition of astrometric and asteroseismic information, the detection and characterization of gravitational wave events, 3-D simulations of convection and model atmospheres, and improved laboratory measurements for transition probabilities and nuclear masses. All of these areas are key for continued improvement, and such improvement will benefit many areas of astrophysics.

We report the discovery of the microlensing planet OGLE-2018-BLG-0740Lb. The planet is detected with a very strong signal of $\Delta\chi^2\sim 4630$, but the interpretation of the signal suffers from two types of degeneracies. One type is caused by the previously known close/wide degeneracy, and the other is caused by an ambiguity between two solutions, in which one solution requires to incorporate finite-source effects, while the other solution is consistent with a point-source interpretation. Although difficult to be firmly resolved based on only the photometric data, the degeneracy is resolved in strong favor of the point-source solution with the additional external information obtained from astrometric and spectroscopic observations. The small astrometric offset between the source and baseline object supports that the blend is the lens and this interpretation is further secured by the consistency of the spectroscopic distance estimate of the blend with the lensing parameters of the point-source solution. The estimated mass of the host is $1.0\pm 0.1~M_\odot$ and the mass of the planet is $4.5\pm 0.6~M_{\rm J}$ (close solution) or $4.8\pm 0.6~M_{\rm J}$ (wide solution) and the lens is located at a distance of $3.2\pm 0.5$~kpc. The bright nature of the lens, with $I\sim 17.1$ ($V\sim 18.2$), combined with its dominance of the observed flux suggest that radial-velocity (RV) follow-up observations of the lens can be done using high-resolution spectrometers mounted on large telescopes, e.g., VLT/ESPRESSO, and this can potentially not only measure the period and eccentricity of the planet but also probe for close-in planets. We estimate that the expected RV amplitude would be $\sim 60\sin i ~{\rm m~s}^{-1}$.

Theoretical models of protoplanetary disc dispersal predict a phase where photoevaporation has truncated the disc at several AU, creating a pressure trap which is dust-rich. Previous models predicted this phase could be long-lived (~Myr), contrary to the observational constraints. We show that dust in the pressure trap can be removed from the disc by radiation pressure exerting a significant acceleration, and hence radial velocity, on small dust particles that reside in the surface layers of the disc. The dust in the pressure trap is not subject to radial drift so it can grow to reach sizes large enough to fragment. Hence small particles removed from the surface layers are replaced by the fragments of larger particles. This link means radiation pressure can deplete the dust at all particle sizes. Through a combination of 1D and 2D models, along with secular models that follow the disc's long-term evolution, we show that radiation pressure can deplete dust from pressure traps created by photoevaporation in ~1e5 years, while the photoevaporation created cavity still resides at 10s of AU. After this phase of radiation pressure removal of dust, the disc is gas-rich and dust depleted and radially optically thin to stellar light, having observational signatures similar to a gas-rich, young debris disc. Indeed many of the young stars (~<10 Myr old) classified as hosting a debris disc may rather be discs that have undergone this process.

We report the detection of a hot Jupiter ($M_{p}=1.75_{-0.17}^{+0.14}\ M_{J}$, $R_{p}=1.38\pm0.04\ R_{J}$) orbiting a middle-aged star ($\log g=4.152^{+0.030}_{-0.043}$) in the Transiting Exoplanet Survey Satellite (TESS) southern continuous viewing zone ($\beta=-79.59^{\circ}$). We confirm the planetary nature of the candidate TOI-150.01 using radial velocity observations from the APOGEE-2 South spectrograph and the Carnegie Planet Finder Spectrograph, ground-based photometric observations from the robotic Three-hundred MilliMeter Telescope at Las Campanas Observatory, and Gaia distance estimates. Large-scale spectroscopic surveys, such as APOGEE/APOGEE-2, now have sufficient radial velocity precision to directly confirm the signature of giant exoplanets, making such data sets valuable tools in the TESS era. Continual monitoring of TOI-150 by TESS can reveal additional planets and subsequent observations can provide insights into planetary system architectures involving a hot Jupiter around a star about halfway through its main-sequence life.

Crucial progress in our understanding of star formation and feedback will depend on the ability to obtain spatially resolved spectroscopic observations of \ion{H}{ii} regions, from which reliable instantaneous measurements of their physical conditions can be obtained. Acquiring these datasets across full galactic systems will prove crucial for obtaining population samples that enable us to understand the time evolution of similar regions, and the variability of conditions among coeval regions. Separating the spatial and temporal dependencies in such way for different physical processes involved in star formation and the injection of feedback is crucial to overcome the inherit degeneracies associated with observing instantaneous snapshots of a dynamic ISM at any given time. Emission line diagnostics are at the core of measuring the physical condition in \ion{H}{ii} regions (e.g. dynamics, SFR, chemical abundances, dust extinction, ionization and excitation, etc.). These measurements require high spatial resolution, contiguous coverage across full galactic systems, and sensitivities significantly deeper than past efforts. The spatial scale required to resolve the \ion{H}{ii} regions of a few pc is only attainable in the Local Group where very large sky coverage is necessary.

Practically all known planet hosts will evolve into white dwarfs, and large parts of their planetary systems will survive this transition - the same is true for the solar system beyond the orbit of Mars. Spectroscopy of white dwarfs accreting planetary debris provides the most accurate insight into the bulk composition of exo-planets. Ground-based spectroscopic surveys of ~260, 000 white dwarfs detected with Gaia will identify >1000 evolved planetary systems, and high-throughput high-resolution space-based ultraviolet spectroscopy is essential to measure in detail their abundances. So far, evidence for two planetesimals orbiting closely around white dwarfs has been obtained, and their study provides important constraints on the composition and internal structure of these bodies. Major photometric and spectroscopic efforts will be necessary to assemble a sample of such close-in planetesimals that is sufficiently large to establish their properties as a population, and to deduce the architectures of the outer planetary systems from where they originated. Mid-infrared spectroscopy of the dusty disks will provide detailed mineralogical information of the debris, which, in combination with the elemental abundances measured from the white dwarf spectroscopy, will enable detailed physical modelling of the chemical, thermodynamic, and physical history of the accreted material. Flexible multi-epoch infrared observations are essential to determine the physical nature, and origin of the variability observed in many of the dusty disks. Finally, the direct detection of the outer reservoirs feeding material to the white dwarfs will require sensitive mid- and far-infrared capabilities.

A large fraction of the baryons and most of the metals in the Universe are unaccounted for. They likely lie in extended galaxy halos, galaxy groups, and the cosmic web, and measuring their nature is essential to understanding galaxy formation. These environments have virial temperatures >10^5.5 K, so the gas should be visible in X-rays. Here we show the breakthrough capabilities of grating spectrometers to 1) detect these reservoirs of hidden metals and mass, and 2) quantify hot gas flows, turbulence, and rotation around the Milky Way and external galaxies. Grating spectrometers are essential instruments for future X-ray missions, and existing technologies provide 50-1500-fold higher throughput compared to current orbiting instruments.

The majority of baryons reside beyond the optical extent of a galaxy in the circumgalactic and intergalactic media (CGM/IGM). Gaseous halos are inextricably linked to the appearance of their host galaxies through a complex story of accretion, feedback, and continual recycling. The energetic processes, which define the state of gas in the CGM, are the same ones that 1) regulate stellar growth so that it is not over-efficient, and 2) create the diversity of today's galaxy colors, SFRs, and morphologies spanning Hubble's Tuning Fork Diagram. They work in concert to set the speed of growth on the star-forming Main Sequence, transform a galaxy across the Green Valley, and maintain a galaxy's quenched appearance on the Red Sequence. Most baryons in halos more massive than 10^12 Msolar along with their high-energy physics and dynamics remain invisible because that gas is heated above the UV ionization states. We argue that information on many of the essential drivers of galaxy evolution is primarily contained in this "missing" hot gas phase. Completing the picture of galaxy formation requires uncovering the physical mechanisms behind stellar and SMBH feedback driving mass, metals, and energy into the CGM. By opening galactic hot halos to new wavebands, we not only obtain fossil imprints of >13 Gyrs of evolution, but observe on-going hot-mode accretion, the deposition of superwind outflows into the CGM, and the re-arrangement of baryons by SMBH feedback. A description of the flows of mass, metals, and energy will only be complete by observing the thermodynamic states, chemical compositions, structure, and dynamics of T>=10^6 K halos. These measurements are uniquely possible with a next-generation X-ray observatory if it provides the sensitivity to detect faint CGM emission, spectroscopic power to measure absorption lines and gas motions, and high spatial resolution to resolve structures.

The distribution of gas-phase abundances in galaxy disks encodes the history of nucleosynthesis and transport through the interstellar medium (ISM) over cosmic time. Multi-object and high resolution integral-field spectroscopy have started to measure these distributions across hundreds of HII regions individually resolved at $\lesssim 100$ pc scales in a handful of objects, but in the coming decade these studies will expand to larger samples of galaxies. This will allow us to understand the role of feedback and turbulence in driving the mixing and diffusion of metals in the ISM, and statistically assess the role of galaxy environment and disk dynamics in modifying how mixing proceeds. Detailed searches for over- and under-enriched regions can address to what extent star formation is triggered by previous generations of star formation and by pristine and recycled gas flows. Local galaxies, for which these detailed measurements will be possible within the next decade, will inform the interpretation of integrated measurements at high-z, where very different dynamical gas-rich environments are found in early disk galaxies. Currently, progress in the field is severely hampered by the 0.2-0.3 dex level systematic uncertainties plaguing nebular abundance diagnostics. Improving our detailed understanding of ionized nebulae at $<$20 pc scales will help us find a solution to this problem, which will prove key to the study of metal enrichment and mixing across the galaxy population in the next decade.

The expansion of the Universe is understood to have accelerated during two epochs: in its very first moments during a period of Inflation and much more recently, at $z < 1$, when Dark Energy is hypothesized to drive cosmic acceleration. The undiscovered mechanisms behind these two epochs represent some of the most important open problems in fundamental physics. The large cosmological volume at $2 < z < 5$, together with the ability to efficiently target high-$z$ galaxies with known techniques, enables large gains in the study of Inflation and Dark Energy. A future spectroscopic survey can test the Gaussianity of the initial conditions up to a factor of ~50 better than our current bounds, crossing the crucial theoretical threshold of $\sigma(f_{NL}^{\rm local})$ of order unity that separates single field and multi-field models. Simultaneously, it can measure the fraction of Dark Energy at the percent level up to $z = 5$, thus serving as an unprecedented test of the standard model and opening up a tremendous discovery space.

The evolution of a star is driven by the physical processes in its interior making the theory of stellar structure and evolution the most crucial ingredient for not only stellar evolution studies, but any field of astronomy which relies on the yields along stellar evolution. High-precision time-series photometric data assembled by recent space missions revealed that current models of stellar structure and evolution show major shortcomings already in the two earliest nuclear burning phases, impacting all subsequent phases prior to the formation of the end-of-life remnant. This white paper focuses specifically on the transport of chemical elements and of angular momentum in the stellar structure and evolution models of stars born with convective core, as revealed by their gravity-mode oscillations.

Interacting binaries containing white dwarfs can lead to a variety of outcomes that range from powerful thermonuclear explosions, which are important in the chemical evolution of galaxies and as cosmological distance estimators, to strong sources of low frequency gravitational wave radiation, which makes them ideal calibrators for the gravitational low-frequency wave detector LISA mission. However, current theoretical evolution models still fail to explain the observed properties of the known populations of white dwarfs in both interacting and detached binaries. Major limitations are that the existing population models have generally been developed to explain the properties of sub-samples of these systems, occupying small volumes of the vast parameter space, and that the observed samples are severely biased. The overarching goal for the next decade is to assemble a large and homogeneous sample of white dwarf binaries that spans the entire range of evolutionary states, to obtain precise measurements of their physical properties, and to further develop the theory to satisfactorily reproduce the properties of the entire population. While ongoing and future all-sky high- and low-resolution optical spectroscopic surveys allow us to enlarge the sample of these systems, high-resolution ultraviolet spectroscopy is absolutely essential for the characterization of the white dwarfs in these binaries. The Hubble Space Telescope is currently the only facility that provides ultraviolet spectroscopy, and with its foreseeable demise, planning the next ultraviolet mission is of utmost urgency.

Galactic interstellar dust has a profound impact not only on our observations of objects throughout the Universe, but also on the morphology, star formation, and chemical evolution of the Galaxy. The advent of massive imaging and spectroscopic surveys (particularly in the infrared) places us on the threshold of being able to map the properties and dynamics of dust and the interstellar medium (ISM) in three dimensions throughout the Milky Way disk and bulge. These developments will enable a fundamentally new understanding of dust properties, including how grains respond to their local environment and how those environments affect dust attenuation of background objects of interest. Distance-resolved maps of dust motion also hold great promise for tracing the flow of interstellar material throughout the Galaxy on a variety of scales, from bar-streaming motions to the collapse and dissolution of individual molecular clouds. These advances require optical and infrared imaging of stars throughout the Galactic midplane, stretching many kiloparsecs from the Sun, matched with very dense spectroscopic coverage to probe the ISM's fine-grained structure.

The ubiquitous variability of quasars across a wide range of wavelengths and timescales encodes critical information about the structure and dynamics of the circumnuclear emitting regions that are too small to be directly resolved, as well as detailed underlying physics of accretion and feedback processes in these active supermassive black holes. We emphasize the importance of studying quasar variability with time-domain spectroscopy, focusing on two science cases: (1) reverberation mapping (RM) to measure the broad-line region sizes and black hole masses in distant quasars; (2) spectroscopic follow-up of extreme variability quasars that dramatically change their continuum and broad-line flux within several years. We highlight the need for dedicated optical-infrared spectroscopic survey facilities in the coming decades to accompany wide-area time-domain imaging surveys, including: (1) the next phase of the Sloan Digital Sky Survey (SDSS-V; ~2020-2025), an all-sky, time-domain multi-object spectroscopic survey with 2.5m-class telescopes; (2) the planned Maunakea Spectroscopic Explorer, a dedicated 10m-class spectroscopic survey telescope with a 1.5 sq. deg field-of-view and multiplex of thousands of fibers in both optical and near-IR (J+H) to begin operations in 2029; (3) the Time-domain Spectroscopic Observatory (TSO), a proposed Probe-class ~1.3m telescope at L2, with imaging and spectroscopy (R=200, 1800) in 4 bands (0.3 - 5 micron) and rapid slew capability to 90% of sky, which will extend the coverage of Hbeta to z=8.

As part of the 100IAS survey, a program aimed to obtain nebular-phase spectra for a volume-limited and homogeneous sample of Type Ia supernovae (SNe Ia), we observed ASASSN-18tb (SN 2018fhw) at 139 days past maximum light. ASASSN-18tb was a fast-declining, sub-luminous event that fits well within the observed photometric and spectroscopic distributions of the SN Ia population. We detect a prominent H$\alpha$ emission line of $L_{{\rm H}\alpha}=2.2\pm0.2\times10^{38}$ ergs s$^{-1}$ with FWHM $\approx1100$ km s$^{-1}$ in the nebular-phase spectrum of this SN Ia. High luminosity H$\alpha$ emission ($L_{{\rm H}\alpha}\gtrsim 10^{40}$ ergs~s$^{-1}$) has previously been discovered in a rare class of SNe Ia-like objects showing CSM interactions (SNe Ia-CSM). They predominantly belong to over-luminous ($M_{\rm max}<-19$ mag in optical) 1991T-like SNe Ia and are exclusively found in star-forming galaxies. By contrast, ASASSN-18tb is a sub-luminous SN Ia ($M_{B, {\rm max}}\sim -17.7$ mag) found in an early-type galaxy dominated by old stellar populations. We discuss possible origins for the observed hydrogen. Out of 75 SNe Ia for which we have so far obtained nebular spectra in 100IAS, no other SN shows a $\sim 1000 \,{\rm km s^{-1}}$ H$\alpha$ emission line with comparable line luminosity as ASASSN-18tb, emphasizing the rarity of such emission in the nebular phase. Based on preliminary results from our survey, the rate for ASASSN-18tb-like nebular H$\alpha$ emission could be as high as $\sim 10\%$ level among sub-luminous SNe Ia.

Each of the giant planets within the Solar System has large moons but none of these moons have their own moons (which we call ${\it submoons}$). By analogy with studies of moons around short-period exoplanets, we investigate the tidal-dynamical stability of submoons. We find that 10 km-scale submoons can only survive around large (1000 km-scale) moons on wide-separation orbits. Tidal dissipation destabilizes the orbits of submoons around moons that are small or too close to their host planet; this is the case for most of the Solar System's moons. A handful of known moons are, however, capable of hosting long-lived submoons: Saturn's moons Titan and Iapetus, Jupiter's moon Callisto, and Earth's Moon. Based on its inferred mass and orbital separation, the newly-discovered exomoon candidate Kepler-1625b-I can in principle host a large submoon, although its stability depends on a number of unknown parameters. We discuss the possible habitability of submoons and the potential for subsubmoons. The existence, or lack thereof, of submoons, may yield important constraints on satellite formation and evolution in planetary systems.

We present photometry of 30 Galactic RR Lyrae variables taken with HST WFC3/IR for the Carnegie-Chicago Hubble Program. These measurements form the base of the distance ladder measurements that comprise a pure Population II base to a measurement of Ho at an accuracy of 3%. These data are taken with the same instrument and filter (F160W) as our observations of RR Lyrae stars in external galaxies so as to to minimize sources of systematic error in our calibration of the extragalactic distance scale. We calculate mean magnitudes based on one to three measurements for each RR Lyrae star using star-by-star templates generated from densely time-sampled data at optical and mid-infrared wavelengths. We use four RR Lyrae stars from our sample with well-measured HST parallaxes to determine a zero point. This zero point will soon be improved with the large number of precise parallaxes to be provided by Gaia. We also provide preliminary calibration with the TGAS & Gaia DR2 data, and all three zero points are in agreement, to within their uncertainties.

Twenty years have passed since first light for the Sloan Digital Sky Survey (SDSS). Here, we release data taken by the fourth phase of SDSS (SDSS-IV) across its first three years of operation (July 2014-July 2017). This is the third data release for SDSS-IV, and the fifteenth from SDSS (Data Release Fifteen; DR15). New data come from MaNGA - we release 4824 datacubes, as well as the first stellar spectra in the MaNGA Stellar Library (MaStar), the first set of survey-supported analysis products (e.g. stellar and gas kinematics, emission line, and other maps) from the MaNGA Data Analysis Pipeline (DAP), and a new data visualisation and access tool we call "Marvin". The next data release, DR16, will include new data from both APOGEE-2 and eBOSS; those surveys release no new data here, but we document updates and corrections to their data processing pipelines. The release is cumulative; it also includes the most recent reductions and calibrations of all data taken by SDSS since first light. In this paper we describe the location and format of the data and tools and cite technical references describing how it was obtained and processed. The SDSS website (www.sdss.org) has also been updated, providing links to data downloads, tutorials and examples of data use. While SDSS-IV will continue to collect astronomical data until 2020, and will be followed by SDSS-V (2020-2025), we end this paper by describing plans to ensure the sustainability of the SDSS data archive for many years beyond the collection of data.

We report a luminous Type II supernova, ASASSN-15nx, with a peak luminosity of M_V=-20 mag, that is between typical core-collapse supernovae and super-luminous supernovae. The post-peak optical light curves show a long, linear decline with a steep slope of 2.5 mag/100 days (i.e., an exponential decline in flux), through the end of observations at phase ~260 days. In contrast, the light curves of hydrogen rich supernovae (SNe II-P/L) always show breaks in their light curves at phase ~100 days, before settling onto Co56 radioactive decay tails with a decline rate of about 1 mag/100 days. The spectra of ASASSN-15nx do not exhibit the narrow emission-line features characteristic of Type IIn SNe, which can have a wide variety of light-curve shapes usually attributed to strong interactions with a dense circumstellar medium (CSM). ASASSN-15nx has a number of spectroscopic peculiarities, including a relatively weak and triangularly-shaped H-alpha emission profile with no absorption component. The physical origin of these peculiarities is unclear, but the long and linear post-peak light curve without a break suggests a single dominant powering mechanism. Decay of a large amount of Ni56 (M_Ni56 = 1.6 +/- 0.2 M_sun) can power the light curve of ASASSN-15nx, and the steep light-curve slope requires substantial gamma-ray escape from the ejecta, which is possible given a low-mass hydrogen envelope for the progenitor. Another possibility is strong CSM interactions powering the light curve, but the CSM needs to be sculpted to produce the unique light-curve shape and to avoid producing SN IIn-like narrow emission lines.

We present nebular-phase spectra of the Type Ia supernova (SN Ia) 2016brx, a member of the 1991bg-like subclass that lies at the faint end of the SN Ia luminosity function. Nebular spectra are available for only three other 1991bg-like SNe, and their Co line centers are all within <~ 500 km/s of each other. In contrast, the nebular Co line center of SN 2016brx is blue-shifted by >1500 km/s compared to them and by ~1200 km/s compared to the rest frame. This is a significant shift relative to the narrow nebular line velocity dispersion of <~ 2000 km/s of these SNe. The large range of nebular line shifts implies that the Ni56 in the ejecta of SN 1991bg-like events is off-center by ~1000 km/s rather than universally centrally confined as previously suggested. With the addition of SN 2016brx, the Co nebular line shapes of 1991bg-like objects appear to connect with the brighter SNe Ia that show double-peak profiles, hinting at a continuous distribution of line profiles among SNe Ia. One class of models to produce both off-center and bi-modal Ni56 distributions is collisions of white dwarfs with unequal and equal masses.

We describe test observations made with a customized 640 x 512 pixel Indium Gallium Arsenide (InGaAs) prototype astronomical camera on the 100" DuPont telescope. This is the first test of InGaAs as a cost-effective alternative to HgCdTe for research-grade astronomical observations. The camera exhibits an instrument background of 113 e-/sec/pixel (dark + thermal) at an operating temperature of -40C for the sensor, maintained by a simple thermo-electric cooler. The optical train and mechanical structure float at ambient temperature with no cold stop, in contrast to most IR instruments which must be cooled to mitigate thermal backgrounds. Measurements of the night sky using a reimager with plate scale of 0.4 arc seconds / pixel show that the sky flux in Y is comparable to the dark current. At J the sky brightness exceeds dark current by a factor of four, and hence dominates the noise budget. The sensor read noise of ~43e- falls below sky+dark noise for exposures of t>7 seconds in Y and 3.5 seconds in J. We present test observations of several selected science targets, including high-significance detections of a lensed Type Ia supernova, a type IIb supernova, and a z=6.3 quasar. Deeper images are obtained for two local galaxies monitored for IR transients, and a galaxy cluster at z=0.87. Finally, we observe a partial transit of the hot JupiterHATS34b, demonstrating the photometric stability required over several hours to detect a 1.2% transit depth at high significance. A tiling of available larger-format sensors would produce an IR survey instrument with significant cost savings relative to HgCdTe-based cameras, if one is willing to forego the K band. Such a camera would be sensitive for a week or more to isotropic emission from r-process kilonova ejecta similar to that observed in GW170817, over the full 190 Mpc horizon of Advanced LIGO's design sensitivity for neutron star mergers.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) has been in operation since July 2014. This paper describes the second data release from this phase, and the fourteenth from SDSS overall (making this, Data Release Fourteen or DR14). This release makes public data taken by SDSS-IV in its first two years of operation (July 2014-2016). Like all previous SDSS releases, DR14 is cumulative, including the most recent reductions and calibrations of all data taken by SDSS since the first phase began operations in 2000. New in DR14 is the first public release of data from the extended Baryon Oscillation Spectroscopic Survey (eBOSS); the first data from the second phase of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE-2), including stellar parameter estimates from an innovative data driven machine learning algorithm known as "The Cannon"; and almost twice as many data cubes from the Mapping Nearby Galaxies at APO (MaNGA) survey as were in the previous release (N = 2812 in total). This paper describes the location and format of the publicly available data from SDSS-IV surveys. We provide references to the important technical papers describing how these data have been taken (both targeting and observation details) and processed for scientific use. The SDSS website (www.sdss.org) has been updated for this release, and provides links to data downloads, as well as tutorials and examples of data use. SDSS-IV is planning to continue to collect astronomical data until 2020, and will be followed by SDSS-V.

Swope Supernova Survey 2017a (SSS17a) was discovered as the first optical counterpart to the gravitational wave event GW170817. Although its light curve on the timescale of weeks roughly matches the expected luminosity and red color of an r-process powered transient, the explanation for the blue emission from high velocity material over the first few days is not as clear. Here we show that the power-law evolution of the luminosity, temperature, and photospheric radius during these early times can be explained by cooling of shock heated material around the neutron star merger. This heating is likely from the interaction of the gamma-ray burst jet with merger debris, so-called cocoon emission. We summarize the properties of this emission and provide formulae that can be used to study future detections of shock cooling from merging neutron stars. This argues that optical transient surveys should search for such early, blue light if they wish to find off-axis gamma-ray bursts and double neutron star gravitational wave events as soon as possible after the merger.

We describe and employ a stacking procedure to investigate abundances derived from the low S/N spectra obtained in the Carnegie RR Lyrae Survey (CARRS; Kollmeier et al. 2013). We find iron metallicities that extend from [Fe/H] ~ -2.5 to values at least as large as [Fe/H] ~ -0.5 in the 274-spectrum CARRS RRc data set. We consider RRc sample contamination by high amplitude solar metallicity delta Scuti stars (HADS) at periods less than 0.3 days, where photometric discrimination between RRc and delta Scuti stars has proven to be problematic. We offer a spectroscopic discriminant, the well-marked overabundance of heavy elements, principally [Ba/H], that is a common, if not universal, characteristic of HADS of all periods and axial rotations. No bona fide RRc stars known to us have verified heavy-element overabundances. Three out of 34 stars in our sample with [Fe/H] > -0.7 exhibit anomalously strong features of Sr, Y, Zr, Ba, and many rare earths. However, carbon is not enhanced in these three stars, and we conclude that their elevated n-capture abundances have not been generated in interior neutron-capture nucleosynthesis. Contamination by HADS appears to be unimportant, and metal-rich RRc stars occur in approximately the same proportion in the Galactic field as do metal-rich RRab stars. An apparent dearth of metal-rich RRc is probably a statistical fluke. Finally we show that RRc stars have a similar inverse period-metallicity relationship as has been found for RRab stars.

Stellar tidal streams provide an opportunity to study the motion and structure of the disrupting galaxy as well as the gravitational potential of its host. Streams around the Milky Way are especially promising as phase space positions of individual stars will be measured by ongoing or upcoming surveys. Nevertheless, it remains a challenge to accurately assess distances to stars farther than 10 kpc from the Sun, where we have the poorest knowledge of the Galaxy's mass distribution. To address this we present observations of 32 candidate RR Lyrae stars in the Orphan tidal stream taken as part of the Spitzer Merger History and Shape of the Galactic Halo (SMHASH) program. The extremely tight correlation between the periods, luminosities, and metallicities of RR Lyrae variable stars in the Spitzer IRAC $\mathrm{3.6 \mu m}$ band allows the determination of precise distances to individual stars; the median statistical distance uncertainty to each RR Lyrae star is $2.5\%$. By fitting orbits in an example potential we obtain an upper limit on the mass of the Milky Way interior to 60 kpc of $\mathrm{5.6_{-1.1}^{+1.2}\times 10^{11}\ M_\odot}$, bringing estimates based on the Orphan Stream in line with those using other tracers. The SMHASH data also resolve the stream in line--of--sight depth, allowing a new perspective on the internal structure of the disrupted dwarf galaxy. Comparing with N--body models we find that the progenitor had an initial dark halo mass of approximately $\mathrm{3.2 \times 10^{9}\ M_\odot}$, placing the Orphan Stream's progenitor amongst the classical dwarf spheroidals.

SDSS-V will be an all-sky, multi-epoch spectroscopic survey of over six million objects. It is designed to decode the history of the Milky Way, trace the emergence of the chemical elements, reveal the inner workings of stars, and investigate the origin of planets. It will also create an integral-field spectroscopic map of the gas in the Galaxy and the Local Group that is 1,000x larger than the current state of the art and at high enough spatial resolution to reveal the self-regulation mechanisms of galactic ecosystems. SDSS-V will pioneer systematic, spectroscopic monitoring across the whole sky, revealing changes on timescales from 20 minutes to 20 years. The survey will thus track the flickers, flares, and radical transformations of the most luminous persistent objects in the universe: massive black holes growing at the centers of galaxies. The scope and flexibility of SDSS-V will be unique among extant and future spectroscopic surveys: it is all-sky, with matched survey infrastructures in both hemispheres; it provides near-IR and optical multi-object fiber spectroscopy that is rapidly reconfigurable to serve high target densities, targets of opportunity, and time-domain monitoring; and it provides optical, ultra-wide-field integral field spectroscopy. SDSS-V, with its programs anticipated to start in 2020, will be well-timed to multiply the scientific output from major space missions (e.g., TESS, Gaia, eROSITA) and ground-based projects. SDSS-V builds on the 25-year heritage of SDSS's advances in data analysis, collaboration infrastructure, and product deliverables. The project is now refining its science scope, optimizing the survey strategies, and developing new hardware that builds on the SDSS-IV infrastructure. We present here an overview of the current state of these developments as we seek to build our worldwide consortium of institutional and individual members.

On 2017 August 17, Swope Supernova Survey 2017a (SSS17a) was discovered as the optical counterpart of the binary neutron star gravitational wave event GW170817. We report time-series spectroscopy of SSS17a from 11.75 hours until 8.5 days after merger. Over the first hour of observations the ejecta rapidly expanded and cooled. Applying blackbody fits to the spectra, we measure the photosphere cooling from $11,000^{+3400}_{-900}$ K to $9300^{+300}_{-300}$ K, and determine a photospheric velocity of roughly 30% of the speed of light. The spectra of SSS17a begin displaying broad features after 1.46 days, and evolve qualitatively over each subsequent day, with distinct blue (early-time) and red (late-time) components. The late-time component is consistent with theoretical models of r-process-enriched neutron star ejecta, whereas the blue component requires high velocity, lanthanide-free material.

On 2017 August 17, gravitational waves were detected from a binary neutron star merger, GW170817, along with a coincident short gamma-ray burst, GRB170817A. An optical transient source, Swope Supernova Survey 17a (SSS17a), was subsequently identified as the counterpart of this event. We present ultraviolet, optical and infrared light curves of SSS17a extending from 10.9 hours to 18 days post-merger. We constrain the radioactively-powered transient resulting from the ejection of neutron-rich material. The fast rise of the light curves, subsequent decay, and rapid color evolution are consistent with multiple ejecta components of differing lanthanide abundance. The late-time light curve indicates that SSS17a produced at least ~0.05 solar masses of heavy elements, demonstrating that neutron star mergers play a role in r-process nucleosynthesis in the Universe.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) began observations in July 2014. It pursues three core programs: APOGEE-2, MaNGA, and eBOSS. In addition, eBOSS contains two major subprograms: TDSS and SPIDERS. This paper describes the first data release from SDSS-IV, Data Release 13 (DR13), which contains new data, reanalysis of existing data sets and, like all SDSS data releases, is inclusive of previously released data. DR13 makes publicly available 1390 spatially resolved integral field unit observations of nearby galaxies from MaNGA, the first data released from this survey. It includes new observations from eBOSS, completing SEQUELS. In addition to targeting galaxies and quasars, SEQUELS also targeted variability-selected objects from TDSS and X-ray selected objects from SPIDERS. DR13 includes new reductions of the SDSS-III BOSS data, improving the spectrophotometric calibration and redshift classification. DR13 releases new reductions of the APOGEE-1 data from SDSS-III, with abundances of elements not previously included and improved stellar parameters for dwarf stars and cooler stars. For the SDSS imaging data, DR13 provides new, more robust and precise photometric calibrations. Several value-added catalogs are being released in tandem with DR13, in particular target catalogs relevant for eBOSS, TDSS, and SPIDERS, and an updated red-clump catalog for APOGEE. This paper describes the location and format of the data now publicly available, as well as providing references to the important technical papers that describe the targeting, observing, and data reduction. The SDSS website, http://www.sdss.org, provides links to the data, tutorials and examples of data access, and extensive documentation of the reduction and analysis procedures. DR13 is the first of a scheduled set that will contain new data and analyses from the planned ~6-year operations of SDSS-IV.

APOGEE-2 is a high-resolution, near-infrared spectroscopic survey observing roughly 300,000 stars across the entire sky. It is the successor to APOGEE and is part of the Sloan Digital Sky Survey IV (SDSS-IV). APOGEE-2 is expanding upon APOGEE's goals of addressing critical questions of stellar astrophysics, stellar populations, and Galactic chemodynamical evolution using (1) an enhanced set of target types and (2) a second spectrograph at Las Campanas Observatory in Chile. APOGEE-2 is targeting red giant branch (RGB) and red clump (RC) stars, RR Lyrae, low-mass dwarf stars, young stellar objects, and numerous other Milky Way and Local Group sources across the entire sky from both hemispheres. In this paper, we describe the APOGEE-2 observational design, target selection catalogs and algorithms, and the targeting-related documentation included in the SDSS data releases.

Zahn (1975) first put forward and calculated in detail the torque experienced by stars in a close binary systems due to dynamical tides. His widely used formula for stars with radiative envelopes and convective cores is expressed in terms of the stellar radius, even though the torque is actually being applied to the convective core at the core radius. This results in a large prefactor, which is very sensitive to the global properties of the star, that multiplies the torque. This large factor is compensated by a very small multiplicative factor, $E_{2}$. Although this is mathematically accurate, depending on the application this can lead to significant errors. The problem is even more severe, since the calculation of $E_{2}$ itself is non-trivial, and different authors have obtained inconsistent values of $E_{2}$. Moreover, many codes (e.g. BSE, StarTrack, MESA) interpolate (and sometimes extrapolate) a fit of $E_{2}$ values to the stellar mass, often in regimes where this is not sound practice. We express the torque in an alternate form, cast in terms of parameters at the envelope-core boundary and a dimensionless coefficient, $\beta_{2}$. Previous attempts to express the torque in such a form are either missing an important factor, which depends on the density profile of the star, or are not easy to implement. We show that $\beta_{2}$ is almost independent of the properties of the star and its value is approximately unity. Our formula for the torque is simple to implement and avoids the difficulties associated with the classic expression.

We present a new method for calculating loops in cosmological perturbation theory. This method is based on approximating a $\Lambda$CDM-like cosmology as a finite sum of complex power-law universes. The decomposition is naturally achieved using an FFTLog algorithm. For power-law cosmologies, all loop integrals are formally equivalent to loop integrals of massless quantum field theory. These integrals have analytic solutions in terms of generalized hypergeometric functions. We provide explicit formulae for the one-loop and the two-loop power spectrum and the one-loop bispectrum. A chief advantage of our approach is that the difficult part of the calculation is cosmology independent, need be done only once, and can be recycled for any relevant predictions. Evaluation of standard loop diagrams then boils down to a simple matrix multiplication. We demonstrate the promise of this method for applications to higher multiplicity/loop correlation functions.

We examine the statistics of the low-redshift Ly-alpha forest in an adaptive mesh refinement hydrodynamic cosmological simulation of sufficient volume to include distinct large-scale environments. We compare our HI column density distribution of absorbers both with recent work and between two highly-refined regions of our simulation: a large-scale overdensity and a large-scale underdensity (on scales of approximately 20 Mpc). We recover the average results presented in Kollmeier et al. (2014) using different simulation methods. We further break down these results as a function of environment to examine the detailed dependence of absorber statistics on large-scale density. We find that the slope of the HI column density distribution in the 10$^{12.5}$ $\le$ N$_{HI}$/cm$^{-2}$ $\le$ 10$^{14.5}$ range depends on environment such that the slope becomes steeper for higher environmental density, and this difference reflects distinct physical conditions of the intergalactic medium on these scales. We track this difference to the different temperature structures of filaments in varying environments. Specifically, filaments in the overdensity are hotter and, correspondingly, are composed of gas with lower HI fractions than those in underdense environments. Our results highlight that in order to understand the physics driving the HI CDD, we need not only improved accounting of the sources of ionizing UV photons, but also of the physical conditions of the IGM and how this may vary as a function of large-scale environment.

We describe the Sloan Digital Sky Survey IV (SDSS-IV), a project encompassing three major spectroscopic programs. The Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) is observing hundreds of thousands of Milky Way stars at high resolution and high signal-to-noise ratio in the near-infrared. The Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey is obtaining spatially-resolved spectroscopy for thousands of nearby galaxies (median redshift of z = 0.03). The extended Baryon Oscillation Spectroscopic Survey (eBOSS) is mapping the galaxy, quasar, and neutral gas distributions between redshifts z = 0.6 and 3.5 to constrain cosmology using baryon acoustic oscillations, redshift space distortions, and the shape of the power spectrum. Within eBOSS, we are conducting two major subprograms: the SPectroscopic IDentification of eROSITA Sources (SPIDERS), investigating X-ray AGN and galaxies in X-ray clusters, and the Time Domain Spectroscopic Survey (TDSS), obtaining spectra of variable sources. All programs use the 2.5-meter Sloan Foundation Telescope at Apache Point Observatory; observations there began in Summer 2014. APOGEE-2 also operates a second near-infrared spectrograph at the 2.5-meter du Pont Telescope at Las Campanas Observatory, with observations beginning in early 2017. Observations at both facilities are scheduled to continue through 2020. In keeping with previous SDSS policy, SDSS-IV provides regularly scheduled public data releases; the first one, Data Release 13, was made available in July 2016.

We present a simple model for low-mass planet formation and subsequent evolution within "transition" discs. We demonstrate quantitatively that the predicted and observed structure of such discs are prime birthsites of planets. Planet formation is likely to proceed through pebble accretion, should a planetary embryo ($M\gtrsim 10^{-4}\,$M$_\oplus$) form. Efficient pebble accretion is likely to be unavoidable in transition disc dust traps, as the size of the dust particles required for pebble accretion are those which are most efficiently trapped in the transition disc dust trap. Rapid pebble accretion within the dust trap gives rise, not only to low-mass planets, but to a large accretion luminosity. This accretion luminosity is sufficient to heat the disc outside the gravitational influence of the planet and makes the disc locally baroclinic, and a source of vorticity. Using numerical simulations we demonstrate that this source of vorticity can lead to the growth of a single large scale vortex in $\sim 100$ orbits, which is capable of trapping particles. Finally, we suggest an evolutionary cycle: planet formation proceeds through pebble accretion, followed by vortex formation and particle trapping in the vortex quenching the planetary accretion and thus removing the vorticity source. After the vortex is destroyed the process can begin anew. This means transition discs should present with large scale vortices for a significant fraction of their lifetimes and remnant planets at large ~10AU separation should be a common outcome of this cycle.

We present a multi-wavelength compilation of new and previously-published photometry for 55 Galactic field RR Lyrae variables. Individual studies, spanning a time baseline of up to 30 years, are self-consistently phased to produce light curves in 10 photometric bands covering the wavelength range from 0.4 to 4.5 microns. Data smoothing via the GLOESS technique is described and applied to generate high-fidelity light curves, from which mean magnitudes, amplitudes, rise-times, and times of minimum and maximum light are derived. 60,000 observations were acquired using the new robotic Three-hundred MilliMeter Telescope (TMMT), which was first deployed at the Carnegie Observatories in Pasadena, CA, and is now permanently installed and operating at Las Campanas Observatory in Chile. We provide a full description of the TMMT hardware, software, and data reduction pipeline. Archival photometry contributed approximately 31,000 observations. Photometric data are given in the standard Johnson UBV, Kron-Cousins RI, 2MASS JHK, and Spitzer [3.6] & [4.5] bandpasses.

We explore the consequences of dynamical evolution of field binaries composed of a primary black hole (BH) and a Wolf-Rayet (WR) star in the context of gravitational wave (GW) source progenitors. We argue, from general considerations, that the spin of the WR-descendent BH will be maximal in a significant number of cases due to dynamical effects. In other cases, the spin should reflect the natal spin of the primary BH which are currently theoretically unconstrained. We argue that the three currently published LIGO systems (GW150914, GW151226, LVT151012) suggest that this spin is small. The resultant effective spin distribution of gravitational wave sources should thus be bi-model if this classic GW progenitor channel is indeed dominant. While this is consistent with the LIGO detections thus far, it is in contrast to the three best-measured high-mass x-ray binary (HMXB) systems. A comparison of the spin distribution of HMXBs and GW sources should ultimately reveal whether or not these systems arise from similar astrophysical channels.

We present the Magellan/FIRE detection of highly-ionized CIV 1550 and OIII] 1666 in a deep infrared spectrum of the z=6.11 gravitationally lensed low-mass galaxy RXC J2248.7-4431-ID3, which has previously-known Lyman-alpha. No corresponding emission is detected at the expected location of HeII 1640. The upper limit on HeII paired with detection of OIII] and CIV constrains possible ionization scenarios. Production of CIV and OIII] requires ionizing photons of 2.5-3.5 Ryd, but once in that state their multiplet emission is powered by collisional excitation at lower energies (~0.5 Ryd). As a pure recombination line, HeII emission is powered by 4 Ryd ionizing photons. The data therefore require a spectrum with significant power at 3.5 Ryd but a rapid drop toward 4.0 Ryd. This hard spectrum with a steep drop is characteristic of low-metallicity stellar populations, and less consistent with soft AGN excitation, which features more 4 Ryd photons and hence higher HeII flux. The conclusions based on ratios of metal line detections to Helium non-detection are strengthened if the gas metallicity is low. RXJ2248-ID3 adds to the growing handful of reionization-era galaxies with UV emission line ratios distinct from the general z=2-3 population, in a way that suggests hard ionizing spectra that do not necessarily originate in AGN.

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.

The Gaia team has applied a renormalization to their internally-derived parallax errors $\sigma_{\rm int}(\pi)$ $$ \sigma_{tgas}(\pi) = \sqrt{[A\sigma_{int}(\pi)]^2 + \sigma_0^2}; \ \ \ \ (A,\sigma_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,\sigma_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: $$ \sigma_{true}(\pi) = \sqrt{(0.79\sigma_{tgas}(\pi))^2 - (0.10\ \rm mas)^2}. $$

The Gaia-Tycho release, scheduled for 14 September, is forecast to yield parallax errors of $\sigma(\pi)\sim 300\,\mu$as for about 2 million Tycho stars. We show analytically that the actual performance should be $$ \sigma(\pi) = {\rm max}(\sigma_{1991}/96,20\,\mu{\rm as}) $$ where $\sigma_{1991}$ is the positional error from the Hipparcos mission. For typical Tycho stars, $\sigma_{1991}\sim 30\,$mas, so this reproduces the usual claims. However, for the 100,000 star Hipparcos subset of this sample, $\sigma_{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.

Like Hipparcos, Gaia is designed to give absolute parallaxes, independent of any astrophysical reference system. And indeed, Gaia's internal zero-point error for parallaxes is likely to be smaller than any individual parallax error. Nevertheless, due in part to mechanical issues of unknown origin, there are many astrophysical questions for which the parallax zero-point error $\sigma(\pi_0)$ will be the fundamentally limiting constraint. These include the distance to the Large Magellanic Cloud and the Galactic Center. We show that by using the photometric parallax estimates for RR Lyrae stars (RRL) within 8kpc, via the ultra-precise infrared period-luminosity relation, one can independently determine a hyper-precise value for $\pi_{0}$. Despite their paucity relative to bright quasars, we show that RRL are competitive due to their order-of-magnitude improved parallax precision for each individual object relative to bright quasars. We show that this method is mathematically robust and well-approximated by analytic formulae over a wide range of relevant distances.

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.

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.

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$\beta$ reverberation mapped active galactic nuclei (AGN). We focus on the MgII, CIV, and CIII] broad emission lines and compare them to both H$\alpha$ and H$\beta$. 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$\alpha$ and H$\beta$, 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.

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.

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.

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.

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$\AA$. 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$\AA$ 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.

A robust prediction of $\Lambda{\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 $\Lambda{\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 $\Lambda{\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 $\Lambda{\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.

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.

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 $\chi^{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.

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.

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.

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.

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.

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.

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.

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%

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]

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.

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)

Using a sample of 92 UV continuum-selected, spectroscopically identified galaxies with <z> = 2.65, all of which have been imaged in the Ly-a line with extremely deep narrow-band imaging, we examine galaxy Ly-a emission profiles to very faint surface brightness limits. The galaxies are representative of spectroscopic samples of LBGs at similar redshifts in terms of apparent magnitude, UV luminosity, inferred extinction, and star formation rate, and were selected without regard to Ly-a emission properties. We use extremely deep stacks of UV continuum and Ly-a emission line images to show that all sub-samples exhibit diffuse Ly-a emission to radii of at least 10" (80 physical kpc), including galaxies whose spectra exhibit Ly-a in net absorption. The intensity scaling, but not the surface brightness distribution, is strongly correlated with the emission observed in the central ~1". The characteristic scale length for Ly-a line emission exceeds that of the UV continuum light for the same galaxies by factors of at least 5-10, regardless of the spectral morphology of Ly-a. Including the extended Ly-a halos increases the total Ly-a flux [and rest equivalent width W_0(Lya)] by an average factor of 5. We argue that most, if not all, of the observed Ly-a emission in the diffuse halos originates in the galaxy H II regions and is scattered in our direction by H I gas in the galaxy's circum-galactic medium (CGM). We show that whether or not a galaxy is classified as a giant "Lyman-a Blob" (LAB) depends sensitively on the Ly-a surface brightness threshold reached by an observation. Accounting for diffuse Ly-a halos, all LBGs would be LABs if surveys were routinely sensitive to 10 times lower surface brightness thresholds; also, essentially all LBGs would qualify as LAEs (W_0(Lya) > 20 A).

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.

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.

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\'e & 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).

We present new limits on the ejection of metal-rich old-population hypervelocity stars from the Galactic center (GC) as probed by the SEGUE-2 survey. Our limits are a factor of 3-10 more stringent than previously reported, depending on stellar type. Compared to the known population of B-star ejectees, there can be no more than 30 times more metal-rich old-population F/G stars ejected from the GC. Because B stars comprise a tiny fraction of a normal stellar population, this places significant limits on a combination of the GC mass function and the ejection mechanism for hypervelocity stars. In the presence of a normal GC mass function, our results require an ejection mechanism that is about 5.5 times more efficient at ejecting B-stars compared to low-mass F/G stars.

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]

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.

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.

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.

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.

We present a method for obtaining accurate black hole (BH) mass estimates from the MgII emission line in active galactic nuclei (AGNs). Employing the large database of AGN measurements from the Sloan Digital Sky Survey (SDSS) presented by Shen et al., we find that AGNs in the redshift range 0.3-0.9, for which a given object can have both H-beta and MgII line widths measured, display a modest but correctable discrepancy in MgII-based masses that correlates with the Eddington ratio. We use the SDSS database to estimate the probability distribution of the true (i.e., H-beta-based) mass given a measured MgII line width. These probability distributions are then applied to the SDSS measurements from Shen et al. across the entire MgII-accessible redshift range (0.3-2.2). We find that accounting for this residual correlation generally increases the dispersion of Eddington ratios by a small factor (~0.09 dex for the redshift and luminosity bins we consider). We continue to find that the intrinsic distribution of Eddington ratios for luminous AGNs is extremely narrow, 0.3-0.4 dex, as demonstrated by Kollmeier et al. Using the method we describe, MgII emission lines can be used with confidence to obtain BH mass estimates.

We report on the discovery of a bright Lyman alpha blob associated with the z=3 quasar SDSSJ124020.91+145535.6 which is also coincident with strong damped Lyman alpha absorption from a foreground galaxy (a so-called proximate damped Lyman alpha system; PDLA). The one dimensional spectrum acquired by the Sloan Digital Sky Survey (SDSS) shows a broad Lyman alpha emission line with a FWHM ~ 500 km/s and a luminosity of L_{Lya} = 3.9e43 erg/s superposed on the trough of the PDLA. Mechanisms for powering this large Lyman alpha luminosity are discussed. We argue against emission from HII regions in the PDLA galaxy since this requires an excessive star-formation rate ~ 500 Msun/yr and would correspond to the largest Lyman alpha luminosity ever measured from a damped Lyman alpha system or starburst galaxy. We use a Monte Carlo radiative transfer simulation to investigate the possibility that the line emission is fluorescent recombination radiation from the PDLA galaxy powered by the ionizing flux of the quasar, but find that the predicted Lyman alpha flux is several orders of magnitude lower than observed. We conclude that the Lyman alpha emission is not associated with the PDLA galaxy at all, but instead is intrinsic to the quasar's host and similar to the extended Lyman alpha 'fuzz' which is detected around many AGN. PDLAs are natural coronagraphs that block their background quasar at Lyman alpha, and we discuss how systems similar to SDSSJ124020.91+145535.6 might be used to image the neutral hydrogen in the PDLA galaxy in silhouette against the screen of extended Lyman alpha emission from the background quasar.

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.

The host galaxies of the five local, z<0.25, long-duration gamma-ray bursts, each of which had a well-documented associated supernova, are all faint and metal-poor compared to the population of local star-forming galaxies. We quantify this statement by using a previous analysis of star-forming galaxies (z<0.2) from the Sloan Digital Sky Survey to estimate the fraction of local star formation as a function of host galaxy oxygen abundance. We find that only a small fraction (<25%) of current star formation occurs in galaxies with oxygen abundance 12+log(O/H)<8.6, i.e., about half that of the Milky Way. However, all five low-z GRB hosts have oxygen abundance below this limit, in three cases very significantly so. If GRBs traced local star formation independent of metallicity, the probability of obtaining such low abundances for all five hosts would be P~0.1%. We conclude that GRBs trace only low-metallicity star formation, and that the Milky Way has been too metal rich to host long GRBs for at least the last several billion years. This result has implications for the potential role of GRBs in mass extinctions, for searches for recent burst remnants in the Milky Way and other large galaxies, for non-detections of late radio emission from local core-collapse supernovae, and for the production of cosmic rays in the local Universe. We also find that the isotropic energy release of these five GRBs, E_iso, steeply decreases with increasing host oxygen abundance, suggesting an upper metallicity limit for ``cosmological'' GRBs at ~0.15 Z_solar. (Abridged)

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.

We have detected Lyman-alpha emission from a damped Lyman-alpha system (DLA) that lies near the bright quasar HS1549+1919. The DLA has the same redshift as HS1549+1919 and was discovered in the spectrum of a faint QSO that lies 49" away (380 proper kpc). The emission line's luminosity, double-peaked profile, and small spatial separation from the DLA suggest that it may be fluorescent Lyman-alpha emission from gas that is absorbing the nearby QSO's radiation. If this is the case, our observations show that the DLA has a size of at least 1.5" and that the QSO's luminosity one million years ago was similar to its luminosity today. A survey for similar systems within 1' of bright QSOs would put interesting limits on the mean quasar lifetime.

We present metallicity estimates for seven open clusters based on spectrophotometric indices from moderate-resolution spectroscopy. Observations of field giants of known metallicity provide a correlation between the spectroscopic indices and the metallicity of open cluster giants. We use \chi^2 analysis to fit the relation of spectrophotometric indices to metallicity in field giants. The resulting function allows an estimate of the target-cluster giants' metallicities with an error in the method of \pm0.08 dex. We derive the following metallicities for the seven open clusters: NGC 1245, [m/H]=-0.14\pm0.04; NGC 2099, [m/H]=+0.05\pm0.05; NGC 2324, [m/H]=-0.06\pm0.04; NGC 2539, [m/H]=-0.04\pm0.03; NGC 2682 (M67), [m/H]=-0.05\pm0.02; NGC 6705, [m/H]=+0.14\pm0.08; NGC 6819, [m/H]=-0.07\pm0.12. These metallicity estimates will be useful in planning future extra-solar planet transit searches since planets may form more readily in metal-rich environments.

We present predictions of Ly-alpha forest-galaxy correlations at z=3 from Eulerian simulations that include the effects of galactic winds, driven primarily by supernova explosions. Galactic winds produce expanding bubbles of shock-heated gas within ~0.5 comoving Mpc/h of luminous galaxies in the simulation, which have space density similar to that of observed Lyman break galaxies (LBGs). However, most of the low-density intergalactic gas that determines the observed properties of the Ly-alpha forest is unaffected by winds. The impact of winds on the Ly-alpha optical depth near galaxies is less dramatic than their impact on gas temperature because winds heat only a small fraction of the gas present in the turnaround regions surrounding galaxies. Hence, Ly-alpha absorption from gas outside the wind bubbles is spread out over the same velocity range occupied by the wind-heated gas. In general, Ly-alpha absorption is expected to be stronger than average near galaxies because of the high gas density. Winds result in a modest reduction of this expected increase of Ly-alpha absorption. Our predictions can be compared to future observations to detect the wind effects and infer their strength, although with the caveat that the results are still dependent on the correspondence of simulated galaxies and observed LBGs. We find that wind effects in our simulations are not strong enough to reproduce the high Ly-alpha transmission within 0.5 Mpc/h comoving of galaxies that has been suggested by recent observations; powerful galactic explosions or ejecta with hyper-escape velocities would be required, but these are unlikely to be produced by ordinary star formation and supernovae alone

We calibrate the photographic photometry of the revised New Luyten Two-Tenths catalog (rNLTT) by matching 3448 rNLTT stars to the Sloan Digital Sky Survey (SDSS). The correction is linear in magnitude and goes from zero at V ~ 14 to 0.32 mag at V=19, in the sense that rNLTT was previously too bright. The correction implies that the underlying USNO-A2.0 photometry, on which rNLTT photometry is based, is non-linear. The new calibration somewhat improves the appearance of the (V,V-J) reduced proper motion diagram in the sense of better separation between disk and halo stars. We repeat Gould's analysis of 5000 halo stars in rNLTT. The most important change is to move the peak of the halo luminosity function about 0.5 mag dimmer, from M_V=10.5 to M_V=11, putting it into good agreement with the parallax-based determination of Dahn et al.

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.

We present a catalog of 345,000 stars with proper motions \mu > 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 \sigma_\mu=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.

We use hydrodynamic simulations to predict correlations between Lya forest absorption and galaxies at redshift z~3. The probability distribution function (PDF) of Lya flux decrements shifts systematically towards higher values in the vicinity of galaxies, reflecting the overdense environments in which these galaxies reside. The predicted signal remains strong in spectra smoothed over 50-200 km/s, allowing tests with moderate resolution quasar spectra. The strong bias of high redshift galaxies towards high density regions imprints a clear signature on the flux PDF, but the predictions are not sensitive to galaxy baryon mass or star formation rate, and they are similar for galaxies and for dark matter halos. The dependence of the flux PDF on galaxy proximity is sensitive to redshift determination errors, with rms errors of 150-300 km/s substantially weakening the predicted trends. On larger scales, the mean galaxy overdensity in a cube of 5 or 10 Mpc/h (comoving) is strongly correlated with the mean Lya flux decrement on a line of sight through the cube center. The slope of the correlation is ~3 times steeper for galaxies than for dark matter as a result of galaxy bias. The predicted large scale correlation is in qualitative agreement with recently reported observational results. However, observations also show a drop in absorption in the immediate vicinity of galaxies, which our models do not predict even if we allow the galaxies or AGNs within them to be ionizing sources. This decreased absorption could be a signature of galaxy feedback on the surrounding IGM, perhaps via galactic winds. Peculiar velocities often allow gas at comoving distances ~1.5 Mpc/h to produce saturated absorption at the galaxy redshift, so any feedback mechanism must suppress neutral hydrogen out to these radii to match the data. (Abridged)

We review recent developments in the theory of the Lyman-alpha forest and their implications for the role of the forest as a test of cosmological models. Simulations predict a relatively tight correlation between the local Lya optical depth and the local gas or dark matter density. Statistical properties of the transmitted flux can constrain the amplitude and shape of the matter power spectrum at high redshift, test the assumption of Gaussian initial conditions, and probe the evolution of dark energy by measuring the Hubble parameter H(z). Simulations predict increased Lya absorption in the vicinity of galaxies, but observations show a Lya deficit within Delta_r ~ 0.5 Mpc/h (comoving). We investigate idealized models of "winds" and find that they must eliminate neutral hydrogen out to comoving radii ~1.5 Mpc/h to marginally explain the data. Winds of this magnitude suppress the flux power spectrum by \~0.1 dex but have little effect on the distribution function or threshold crossing frequency. In light of the stringent demands on winds, we consider the alternative possibility that extended Lya emission from target galaxies replaces absorbed flux, but we conclude that this explanation is unlikely. Taking full advantage of the data coming from large telescopes and from the Sloan Digital Sky Survey will require more complete understanding of the galaxy proximity effect, careful attention to continuum determination, and more accurate numerical predictions, with the goal of reaching 5-10% precision on key cosmological quantities.

Hydrodynamic cosmological simulations predict that the average opacity of the Ly-alpha forest should increase in the neighborhood of galaxies because galaxies form in dense environments. Recent observations (Adelberger et al. 2002) confirm this expectation at large scales, but they show a decrease of absorption at comoving separations Delta_r <~ 1 Mpc/h. We show that this discrepancy is statistically significant, especially for the innermost data point at Delta_r <= 0.5 Mpc/h, even though this data point rests on three galaxy-quasar pairs. Galaxy redshift errors of the expected magnitude are insufficient to resolve the conflict. Peculiar velocities allow gas at comoving distances >~ 1 Mpc/h to produce saturated absorption at the galaxy redshift, putting stringent requirements on any ``feedback'' solution. Local photoionization is insufficient, even if we allow for recurrent AGN activity that keeps the neutral hydrogen fraction below its equilibrium value. A simple ``wind'' model that eliminates all neutral hydrogen in spheres around the observed galaxies can marginally explain the data, but only if the winds extend to comoving radii ~1.5 Mpc/h.

We present complementary data on 5 intermediate redshift (0.44 < z < 0.66) MgII absorbing galaxies, combining high spatial resolution imaging from HST, high--resolution QSO spectroscopy from Keck/HIRES, and galaxy kinematics from intermediate resolution spectroscopy using Keck/LRIS. These data allow a direct comparison of the kinematics of gas at large galactocentric impact parameters with the galaxy kinematics obtained from the faint galaxy spectroscopy. All 5 galaxies appear to be relatively normal spirals, with measured rotation curves yielding circular velocities in the range 100 < v_c < 260 km/s. We find that in 4 of the 5 cases examined, the velocities of all of the Mg II absorption components lie entirely to one side of the galaxy systemic redshift. These observations are consistent with rotation being dominant for the absorbing gas kinematics; however, the total range of velocities observed is inconsistent with simple disk rotation in every case. Simple kinematic models that simultaneously explain both the systemic offset of the absorbing material relative to the galaxy redshifts, and the total velocity width spanned by the absorption, require either extremely thick rotating gas layers, rotation velocities that vary with height above the extrapolation of the galactic plane, or both. In any case, our small sample suggests that rotating ``halo'' gas is a common feature of intermediate redshift spiral galaxies, and that the kinematic signature of rotation dominates over radial infall or outflow even for gas well away from the galactic plane. We discuss possible explanations for this behavior, and compare our observations to possible local analogs.