We report velocity-delay maps for prominent broad emission lines, Ly_alpha, CIV, HeII and H_beta, in the spectrum of NGC5548. The emission-line responses inhabit the interior of a virial envelope. The velocity-delay maps reveal stratified ionization structure. The HeII response inside 5-10 light-days has a broad single-peaked velocity profile. The Ly_alpha, CIV, and H_beta responses peak inside 10 light-days, extend outside 20 light-days, and exhibit a velocity profile with two peaks separated by 5000 km/s in the 10 to 20 light-day delay range. The velocity-delay maps show that the M-shaped lag vs velocity structure found in previous cross-correlation analysis is the signature of a Keplerian disk with a well-defined outer edge at R=20 light-days. The outer wings of the M arise from the virial envelope, and the U-shaped interior of the M is the lower half of an ellipse in the velocity-delay plane. The far-side response is weaker than that from the near side, so that we see clearly the lower half, but only faintly the upper half, of the velocity--delay ellipse. The delay tau=(R/c)(1-sin(i))=5 light-days at line center is from the near edge of the inclined ring, giving the inclination i=45 deg. A black hole mass of M=7x10^7 Msun is consistent with the velocity-delay structure. A barber-pole pattern with stripes moving from red to blue across the CIV and possibly Ly_alpha line profiles suggests the presence of azimuthal structure rotating around the far side of the broad-line region and may be the signature of precession or orbital motion of structures in the inner disk. Further HST observations of NGC 5548 over a multi-year timespan but with a cadence of perhaps 10 days rather than 1 day could help to clarify the nature of this new AGN phenomenon.

We broadly explore the effects of systematic errors on reverberation mapping lag uncertainty estimates from {\tt JAVELIN} and the interpolated cross-correlation function (ICCF) method. We focus on simulated lightcurves from random realizations of the lightcurves of five intensively monitored AGNs. Both methods generally work well even in the presence of systematic errors, although {\tt JAVELIN} generally provides better error estimates. Poorly estimated lightcurve uncertainties have less effect on the ICCF method because, unlike {\tt JAVELIN}, it does not explicitly assume Gaussian statistics. Neither method is sensitive to changes in the stochastic process driving the continuum or the transfer function relating the line lightcurve to the continuum. The only systematic error we considered that causes significant problems is if the line lightcurve is not a smoothed and shifted version of the continuum lightcurve but instead contains some additional sources of variability.

Next-generation optical imaging surveys will revolutionise the observations of weak gravitational lensing by galaxy clusters and provide stringent constraints on growth of structure and cosmic acceleration. In these experiments, accurate modelling of covariance matrices of cluster weak lensing plays the key role in obtaining robust measurements of the mean mass of clusters and cosmological parameters. We use a combination of analytical calculations and high-resolution N-body simulations to derive accurate covariance matrices that span from the virial regime to linear scales of the cluster-matter cross-correlation. We validate this calculation using a public ray-tracing lensing simulation and provide a software package for calculating covariance matrices for a wide range of cluster and source sample choices. We discuss the relative importance of shape noise and density fluctuations, the impact of radial bin size, and the impact of off-diagonal elements. For a weak lensing source density 10 per square arcmin, shape noise typically dominates the variance on comoving scales less than 5 Mpc/h. However, for 60 per square arcmin, potentially achievable with future weak lensing experiments, density fluctuations typically dominate the variance at scales greater than 1 Mpc/h and remain comparable to shape noise on smaller scales.

We present the status of the Dark Energy Spectroscopic Instrument (DESI) and its plans and opportunities for the coming decade. DESI construction and its initial five years of operations are an approved experiment of the US Department of Energy and is summarized here as context for the Astro2020 panel. Beyond 2025, DESI will require new funding to continue operations. We expect that DESI will remain one of the world's best facilities for wide-field spectroscopy throughout the decade. More about the DESI instrument and survey can be found at https://www.desi.lbl.gov.

We model the ultraviolet spectra of the Seyfert 1 galaxy NGC~5548 obtained with the Hubble Space Telescope during the 6-month reverberation-mapping campaign in 2014. Our model of the emission from NGC 5548 corrects for overlying absorption and deblends the individual emission lines. Using the modeled spectra, we measure the response to continuum variations for the deblended and absorption-corrected individual broad emission lines, the velocity-dependent profiles of Ly$\alpha$ and C IV, and the narrow and broad intrinsic absorption features. We find that the time lags for the corrected emission lines are comparable to those for the original data. The velocity-binned lag profiles of Ly$\alpha$ and C IV have a double-peaked structure indicative of a truncated Keplerian disk. The narrow absorption lines show delayed response to continuum variations corresponding to recombination in gas with a density of $\sim 10^5~\rm cm^{-3}$. The high-ionization narrow absorption lines decorrelate from continuum variations during the same period as the broad emission lines. Analyzing the response of these absorption lines during this period shows that the ionizing flux is diminished in strength relative to the far-ultraviolet continuum. The broad absorption lines associated with the X-ray obscurer decrease in strength during this same time interval. The appearance of X-ray obscuration in $\sim\,2012$ corresponds with an increase in the luminosity of NGC 5548 following an extended low state. We suggest that the obscurer is a disk wind triggered by the brightening of NGC 5548 following the decrease in size of the broad-line region during the preceding low-luminosity state.

We develop a novel method to explore the galaxy-halo connection using the galaxy imaging surveys by modeling the projected two-point correlation function measured from the galaxies with reasonable photometric redshift measurements. By assuming a Gaussian form of the photometric redshift errors, we are able to simultaneously constrain the halo occupation distribution (HOD) models and the effective photometric redshift uncertainties. Tests with mock galaxy catalogs demonstrate that this method can successfully recover (within $\sim 1\sigma$) the intrinsic large-scale galaxy bias, as well as the HOD models and the effective photometric redshift uncertainty. This method also works well even for galaxy samples with 10 per cent catastrophic photometric redshift errors.

We formulate a model of the conditional colour-magnitude distribution (CCMD) to describe the distribution of galaxy luminosity and colour as a function of halo mass. It consists of two populations of different colour distributions, dubbed pseudo-blue and pseudo-red, respectively, with each further separated into central and satellite galaxies. We define a global parameterization of these four colour-magnitude distributions and their dependence on halo mass, and we infer parameter values by simultaneously fitting the space densities and auto-correlation functions of 79 galaxy samples from the Sloan Digital Sky Survey defined by fine bins in the colour-magnitude diagram (CMD). The model deprojects the overall galaxy CMD, revealing its tomograph along the halo mass direction. The bimodality of the colour distribution is driven by central galaxies at most luminosities, though at low luminosities it is driven by the difference between blue centrals and red satellites. For central galaxies, the two pseudo-colour components are distinct and orthogonal to each other in the CCMD: at fixed halo mass, pseudo-blue galaxies have a narrow luminosity range and broad colour range, while pseudo-red galaxies have a narrow colour range and broad luminosity range. For pseudo-blue centrals, luminosity correlates tightly with halo mass, while for pseudo-red galaxies colour correlates more tightly (redder galaxies in more massive haloes). The satellite fraction is higher for redder and for fainter galaxies, with colour a stronger indicator than luminosity. We discuss the implications of the results and further applications of the CCMD model.

The neutral hydrogen (HI) gas is an important barometer of recent star formation and metal enrichment activities in galaxies. We develop a novel statistical method for predicting the HI-to-stellar mass ratio $f_{\mathrm{gas}}$ of galaxies from their stellar mass and optical colour, and apply it to a volume-limited galaxy sample jointly observed by the Sloan Digital Sky Survey and the Arecibo Legacy Fast ALFA survey. We eliminate the impact of the Malmquist bias against HI-deficient systems on the $f_{\mathrm{gas}}$ predictor by properly accounting for the HI detection probability of each galaxy in the analysis. The best-fitting $f_{\mathrm{gas}}$ predictor, with an estimated scatter of $0.272$ dex, provides excellent description to the observed HI mass function. After defining an HI excess parameter as the deviation of the observed $f_{\mathrm{gas}}$ from the expected value, we confirm that there exists a strong secondary dependence of the mass-metallicity relation on HI excess. By further examining the 2D metallicity distribution on the specific star formation rate vs. HI excess plane, we demonstrate that the metallicity dependence on HI is more fundamental than that on specific star formation rate. In addition, we show that the environmental dependence of HI in the local Universe can be effectively described by the cross-correlation coefficient between HI excess and the red galaxy overdensity $\rho_{cc}{=}-0.18$. This weak anti-correlation also successfully explains the observed dependence of HI clustering on $f_{\mathrm{gas}}$. Our method provides a useful framework for learning HI gas evolution from the synergy between future HI and optical galaxy surveys.

We present accretion disk size measurements for 15 luminous quasars at $0.7 \leq z \leq 1.9$ derived from $griz$ light curves from the Dark Energy Survey. We measure the disk sizes with continuum reverberation mapping using two methods, both of which are derived from the expectation that accretion disks have a radial temperature gradient and the continuum emission at a given radius is well-described by a single blackbody. In the first method we measure the relative lags between the multiband light curves, which provides the relative time lag between shorter and longer wavelength variations. From this, we are only able to constrain upper limits on disk sizes, as many are consistent with no lag the 2$\sigma$ level. The second method fits the model parameters for the canonical thin disk directly rather than solving for the individual time lags between the light curves. Our measurements demonstrate good agreement with the sizes predicted by this model for accretion rates between 0.3-1 times the Eddington rate. Given our large uncertainties, our measurements are also consistent with disk size measurements from gravitational microlensing studies of strongly lensed quasars, as well as other photometric reverberation mapping results, that find disk sizes that are a factor of a few ($\sim$3) larger than predictions.

We present the first results from a reverberation-mapping campaign undertaken during the first half of 2012, with additional data on one AGN (NGC 3227) from a 2014 campaign. Our main goals are (1) to determine the black hole masses from continuum-Hbeta reverberation signatures, and (2) to look for velocity-dependent time delays that might be indicators of the gross kinematics of the broad-line region. We successfully measure Hbeta time delays and black hole masses for five AGNs, four of which have previous reverberation mass measurements. The values measured here are in agreement with earlier estimates, though there is some intrinsic scatter beyond the formal measurement errors. We observe velocity dependent Hbeta lags in each case, and find that the patterns have changed in the intervening five years for three AGNs that were also observed in 2007.

Cosmic acceleration is the most surprising cosmological discovery in many decades. Testing and distinguishing among possible explanations requires cosmological measurements of extremely high precision probing the full history of cosmic expansion and structure growth and, ideally, compare and contrast matter and relativistic tracers of the gravity potential. This program is one of the defining objectives of the Wide-Field Infrared Survey Telescope (WFIRST), as set forth in the New Worlds, New Horizons report (NWNH) in 2010. The WFIRST mission has the ability to improve these measurements by 1-2 orders of magnitude compared to the current state of the art, while simultaneously extending their redshift grasp, greatly improving control of systematic effects, and taking a unified approach to multiple probes that provide complementary physical information and cross-checks of cosmological results. We describe in this annual report the activities of the Science Investigation Team (SIT) "Cosmology with the High Latitude Survey (HLS)" during the year 2017. This team was selected by NASA in December 2015 in order to address the stringent challenges of the WFIRST dark energy (DE) program through the Project's formulation phase. This SIT has elected to jointly address Galaxy Redshift Survey, Weak Lensing and Cluster Growth and thus fully embrace the fact that the imaging and spectroscopic elements of the HLS will be realized as an integrated observing program, and they jointly impose requirements on performance and operations. WFIRST is designed to be able to deliver a definitive result on the origin of cosmic acceleration. It is not optimized for Figure of Merit sensitivity but for control of systematic uncertainties and for having multiple techniques each with multiple cross-checks. Our SIT work focuses on understanding the potential systematics in the WFIRST DE measurements.

Galaxies exhibit different clustering and quenching properties in clusters, filaments, and the field, but it is still uncertain whether such differences are imprints of the tidal environment on galaxy formation, or if they reflect the variation of the underlying halo mass function across the cosmic web. We measure the dependence of galaxy clustering and quenching on the cosmic web in the Sloan Digital Sky Survey, characterized by the combination of spherical overdensity $\delta_8$ and tidal anisotropy $\alpha_5$ centred on each galaxy. We find that galaxy clustering is a strong function of either $\delta_8$ or $\alpha_5$, and the large-scale galaxy bias shows complex and rich behaviour on the $\delta_8$ vs. $\alpha_5$ plane. Using the mean galaxy colour as a proxy for the average quenched level of galaxies, we find that galaxy quenching is primarily a function of $\delta_8$, with some subtle yet non-trivial dependence on $\alpha_5$ at fixed $\delta_8$. The quenched galaxies generally show stronger small-scale clustering than the star-forming ones at fixed $\delta_8$ or $\alpha_5$, while the characteristic scale at which the amplitude of clustering becomes comparable for both galaxy populations varies with $\delta_8$ and $\alpha_5$. We compare these observed cosmic web dependences of galaxy clustering and quenching with a mock galaxy catalogue constructed from the iHOD model, which places quenched and star-forming galaxies inside dark matter haloes based on the stellar-to-halo mass relation and the halo quenching model --- the $\delta_8$ and $\alpha_5$ dependences of iHOD galaxies are thus solely derived from the cosmic web modulation of the halo mass function. The main observed trends are accounted for extremely well by the iHOD model. Thus any additional direct effect of the large-scale~(${>}5\,h^{-1}{\rm Mpc}$) tidal field on galaxy formation must be extremely weak. [Abridged]

The use of high-quality simulated sky catalogs is essential for the success of cosmological surveys. The catalogs have diverse applications, such as investigating signatures of fundamental physics in cosmological observables, understanding the effect of systematic uncertainties on measured signals and testing mitigation strategies for reducing these uncertainties, aiding analysis pipeline development and testing, and survey strategy optimization. The list of applications is growing with improvements in the quality of the catalogs and the details that they can provide. Given the importance of simulated catalogs, it is critical to provide rigorous validation protocols that enable both catalog providers and users to assess the quality of the catalogs in a straightforward and comprehensive way. For this purpose, we have developed the DESCQA framework for the Large Synoptic Survey Telescope Dark Energy Science Collaboration as well as for the broader community. The goal of DESCQA is to enable the inspection, validation, and comparison of an inhomogeneous set of synthetic catalogs via the provision of a common interface within an automated framework. In this paper, we present the design concept and first implementation of DESCQA. In order to establish and demonstrate its full functionality we use a set of interim catalogs and validation tests. We highlight several important aspects, both technical and scientific, that require thoughtful consideration when designing a validation framework, including validation metrics and how these metrics impose requirements on the synthetic sky catalogs.

We study the dependence of galaxy clustering on atomic gas mass using a sample of $\sim$16,000 galaxies with redshift in the range of $0.0025<z<0.05$ and HI mass of $M_{\rm HI}>10^8M_{\odot}$, drawn from the 70% complete sample of the Arecibo Legacy Fast ALFA survey. We construct subsamples of galaxies with $M_{\rm HI}$ above different thresholds, and make volume-limited clustering measurements in terms of three statistics: the projected two-point correlation function, the projected cross-correlation function with respect to a reference sample selected from the Sloan Digital Sky Survey, and the redshift-space monopole moment. In contrast to previous studies, which found no/weak HI-mass dependence, we find both the clustering amplitude on scales above a few Mpc and the bias factors to increase significantly with increasing HI mass for subsamples with HI mass thresholds above $10^9M_{\odot}$. For HI mass thresholds below $10^9M_{\odot}$, while the measurements have large uncertainties caused by the limited survey volume and sample size, the inferred galaxy bias factors are systematically lower than the minimum halo bias factor from mass-selected halo samples. The simple halo model, in which galaxy content is only determined by halo mass, has difficulties in interpreting the clustering measurements of the HI-selected samples. We extend the simple model by including the halo formation time as an additional parameter. A model that puts HI-rich galaxies into halos that formed late can reproduce the clustering measurements reasonably well. We present the implications of our best-fitting model on the correlation of HI mass with halo mass and formation time, as well as the halo occupation distributions and HI mass functions for central and satellite galaxies. These results are compared with the predictions from semi-analytic galaxy formation models and hydrodynamic galaxy formation simulations.

During the Space Telescope and Optical Reverberation Mapping Project (STORM) observations of NGC 5548, the continuum and emission-line variability became de-correlated during the second half of the 6-month long observing campaign. Here we present Swift and Chandra X-ray spectra of NGC 5548 obtained as a part of the campaign. The Swift spectra show that excess flux (relative to a power-law continuum) in the soft X-ray band appears before the start of the anomalous emission-line behavior, peaks during the period of the anomaly, and then declines. This is a model-independent result suggesting that the soft excess is related to the anomaly. We divide the Swift data into on- and off-anomaly spectra to characterize the soft excess via spectral fitting. The cause of the spectral differences is likely due to a change in the intrinsic spectrum rather than being due to variable obscuration or partial covering. The Chandra spectra have lower signal-to-noise ratios, but are consistent with Swift data. Our preferred model of the soft excess is emission from an optically thick, warm Comptonizing corona, the effective optical depth of which increases during the anomaly. This model simultaneously explains all the three observations: the UV emission line flux decrease, the soft-excess increase, and the emission line anomaly.

Recently, several studies have discovered a strong discrepancy between the large-scale clustering biases of two subsamples of galaxy clusters at the same halo mass, split by their average projected membership distances $R_{\mathrm{mem}}$. The level of this discrepancy significantly exceeds the maximum halo assembly bias signal predicted by LCDM. In this study, we explore whether some of the clustering bias differences could be caused by biases in $R_{\mathrm{mem}}$ due to projection effects from other systems along the line-of-sight. We thoroughly investigate the halo assembly bias of the photometrically-detected redMaPPer clusters in SDSS, by defining a new variant of the average membership distance estimator $\tilde{R}_{\mathrm{mem}}$ that is more robust against projection effects in the cluster membership identification. Using the angular mark correlation functions of clusters, we show that the large-scale bias differences when splitting by $R_{\mathrm{mem}}$ can be largely attributed to such projection effects. After splitting by $\tilde{R}_{\mathrm{mem}}$, the anomalously large signal is reduced, giving a ratio of $1.02\pm0.14$ between the two clustering biases as measured from weak lensing. Using a realistic mock cluster catalog, we predict that the bias ratio between two $\tilde{R}_{\mathrm{mem}}$-split subsamples should be $<1.10$, which is at least 60% weaker than the maximum halo assembly bias signal (1.24) when split by halo concentration. Therefore, our results demonstrate that the level of halo assembly bias exhibited by redMaPPer clusters in SDSS is consistent with the LCDM prediction. With a ten-fold increase in cluster numbers, deeper ongoing surveys will enable a more robust detection of halo assembly bias. Our findings also have important implications for how projection effects and their impact on cluster cosmology can be quantified in photometric cluster catalogs.

Recent studies suggest that the quenching properties of galaxies are correlated over several mega-parsecs. The large-scale "galactic conformity" phenomenon around central galaxies has been regarded as a potential signature of "galaxy assembly bias" or "pre-heating", both of which interpret conformity as a result of direct environmental effects acting on galaxy formation. Building on the iHOD halo quenching framework developed in Zu & Mandelbaum (2015, 2016), we discover that our fiducial halo mass quenching model, without any galaxy assembly bias, can successfully explain the overall environmental dependence and the conformity of galaxy colours in SDSS, as measured by the mark correlation functions of galaxy colours and the red galaxy fractions around isolated primaries, respectively. Our fiducial iHOD halo quenching mock also correctly predicts the differences in the spatial clustering and galaxy-galaxy lensing signals between the more vs. less red galaxy subsamples, split by the red-sequence ridge-line at fixed stellar mass. Meanwhile, models that tie galaxy colours fully or partially to halo assembly bias have difficulties in matching all these observables simultaneously. Therefore, we demonstrate that the observed environmental dependence of galaxy colours can be naturally explained by the combination of 1) halo quenching and 2) the variation of halo mass function with environment --- an indirect environmental effect mediated by two separate physical processes.

We present the results of an optical spectroscopic monitoring program targeting NGC 5548 as part of a larger multi-wavelength reverberation mapping campaign. The campaign spanned six months and achieved an almost daily cadence with observations from five ground-based telescopes. The H$\beta$ and He II $\lambda$4686 broad emission-line light curves lag that of the 5100 $\AA$ optical continuum by $4.17^{+0.36}_{-0.36}$ days and $0.79^{+0.35}_{-0.34}$ days, respectively. The H$\beta$ lag relative to the 1158 $\AA$ ultraviolet continuum light curve measured by the Hubble Space Telescope is roughly $\sim$50% longer than that measured against the optical continuum, and the lag difference is consistent with the observed lag between the optical and ultraviolet continua. This suggests that the characteristic radius of the broad-line region is $\sim$50% larger than the value inferred from optical data alone. We also measured velocity-resolved emission-line lags for H$\beta$ and found a complex velocity-lag structure with shorter lags in the line wings, indicative of a broad-line region dominated by Keplerian motion. The responses of both the H$\beta$ and He II $\lambda$4686 emission lines to the driving continuum changed significantly halfway through the campaign, a phenomenon also observed for C IV, Ly $\alpha$, He II(+O III]), and Si IV(+O IV]) during the same monitoring period. Finally, given the optical luminosity of NGC 5548 during our campaign, the measured H$\beta$ lag is a factor of five shorter than the expected value implied by the $R_\mathrm{BLR} - L_\mathrm{AGN}$ relation based on the past behavior of NGC 5548.

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$\alpha$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

We conduct a multiwavelength continuum variability study of the Seyfert 1 galaxy NGC 5548 to investigate the temperature structure of its accretion disk. The 19 overlapping continuum light curves (1158 to 9157 angstroms) combine simultaneous HST , Swift , and ground-based observations over a 180 day period from 2014 January to July. Light-curve variability is interpreted as the reverberation response of the accretion disk to irradiation by a central time-varying point source. Our model yields the disk inclination, i, temperature T1 at 1 light day from the black hole, and a temperature-radius slope, alpha. We also infer the driving light curve and find that it correlates poorly with both the hard and soft X-ray light curves, suggesting that the X-rays alone may not drive the ultraviolet and optical variability over the observing period. We also decompose the light curves into bright, faint, and mean accretion-disk spectra. These spectra lie below that expected for a standard blackbody accretion disk accreting at L/LEdd = 0.1

During an intensive Hubble Space Telescope (HST) Cosmic Origins Spectrograph (COS) UV monitoring campaign of the Seyfert~1 galaxy NGC 5548 performed from 2014 February to July, the normally highly correlated far-UV continuum and broad emission-line variations decorrelated for ~60 to 70 days, starting ~75 days after the first HST/COS observation. Following this anomalous state, the flux and variability of the broad emission lines returned to a more normal state. This transient behavior, characterised by significant deficits in flux and equivalent width of the strong broad UV emission lines, is the first of its kind to be unambiguously identified in an active galactic nucleus reverberation mapping campaign. The largest corresponding emission-line flux deficits occurred for the high-ionization collisionally excited lines, C IV and Si IV(+O IV]), and also He II(+O III]), while the anomaly in Ly-alpha was substantially smaller. This pattern of behavior indicates a depletion in the flux of photons with E_{\rm ph} > 54 eV, relative to those near 13.6 eV. We suggest two plausible mechanisms for the observed behavior: (i) temporary obscuration of the ionizing continuum incident upon BLR clouds by a moving veil of material lying between the inner accretion disk and inner BLR, perhaps resulting from an episodic ejection of material from the disk, or (ii) a temporary change in the intrinsic ionizing continuum spectral energy distribution resulting in a deficit of ionizing photons with energies > 54 eV, possibly due to a transient restructuring of the Comptonizing atmosphere above the disk. Current evidence appears to favor the latter explanation.

We use both simulated and real quasar light curves to explore modeling photometric reverberation-mapping (RM) data as a stochastic process. We do this using modifications to our previously developed RM method based on modeling quasar variability as a damped random walk. We consider the feasibility of one and two-band photometric RM and compare the results with those fr om spectroscopic RM. We find that our method for two-band photometric RM can be competitive with spectroscopic RM only for strong (large equivalent width) lines like Halpha and Hbeta, and that the one-band method is also feasible, but requires very high-precision photometry. We fail to robustly detect Halpha lags in single-band quasar light curves from OGLE-III and IV despite the outstanding cadence and time span of the data, on account of photometric uncertainties in the range 0.02-0.04 mag. Simulations suggest that success could be achieved if the photometric uncertainties were of order 0.01 mag. Single-band RM for all lines and two-band RM for lower equivalent width lines are likely only feasible for statistical estimate s of mean lags for large samples of AGN of similar properties (e.g., luminosity) rather than for individual quasars. Our approach is directly applicable to the time-domain programs within ongoing and future wide-field imaging surveys, and could provide robust lag measurements for an unprecedented number of systems.

We present ground-based optical photometric monitoring data for NGC 5548, part of an extended multi-wavelength reverberation mapping campaign. The light curves have nearly daily cadence from 2014 January to July in nine filters (\emph{BVRI} and \emph{ugriz}). Combined with ultraviolet data from the \emph{Hubble Space Telescope} and \emph{Swift}, we confirm significant time delays between the continuum bands as a function of wavelength, extending the wavelength coverage from 1158\,\AA\ to the $z$ band ($\sim\!9160$\,\AA). We find that the lags at wavelengths longer than the {\it V} band are equal to or greater than the lags of high-ionization-state emission lines (such as He\,{\sc ii}\,$\lambda 1640$ and $\lambda 4686$), suggesting that the continuum-emitting source is of a physical size comparable to the inner broad-line region (BLR). The trend of lag with wavelength is broadly consistent with the prediction for continuum reprocessing by an accretion disk with $\tau \propto \lambda^{4/3}$. However, the lags also imply a disk radius that is 3 times larger than the prediction from standard thin-disk theory, assuming that the bolometric luminosity is 10\% of the Eddington luminosity ($L = 0.1L_{\rm Edd}$). Using optical spectra from the Large Binocular Telescope, we estimate the bias of the interband continuum lags due to BLR emission observed in the filters. We find that the bias for filters with high levels of BLR contamination ($\sim\! 20\%$) can be important for the shortest continuum lags, and likely has a significant impact on the {\it u} and {\it U} bands owing to Balmer continuum emission.

We use galaxy-galaxy lensing to study the dark matter halos surrounding a sample of Locally Brightest Galaxies (LBGs) selected from the Sloan Digital Sky Survey. We measure mean halo mass as a function of the stellar mass and colour of the central galaxy. Mock catalogues constructed from semi-analytic galaxy formation simulations demonstrate that most LBGs are the central objects of their halos, greatly reducing interpretation uncertainties due to satellite contributions to the lensing signal. Over the full stellar mass range, $10.3 < \log [M_*/M_\odot] < 11.6$, we find that passive central galaxies have halos that are at least twice as massive as those of star-forming objects of the same stellar mass. The significance of this effect exceeds $3\sigma$ for $\log [M_*/M_\odot] > 10.7$. Tests using the mock catalogues and on the data themselves clarify the effects of LBG selection and show that it cannot artificially induce a systematic dependence of halo mass on LBG colour. The bimodality in halo mass at fixed stellar mass is reproduced by the astrophysical model underlying our mock catalogue, but the sign of the effect is inconsistent with recent, nearly parameter-free age-matching models. The sign and magnitude of the effect can, however, be reproduced by halo occupation distribution models with a simple (few-parameter) prescription for type-dependence.

Reverberation mapping (RM) measurements of broad-line region (BLR) lags in z>0.3 quasars are important for directly measuring black hole masses in these distant objects, but so far there have been limited attempts and success given the practical difficulties of RM in this regime. Here we report preliminary results of 15 BLR lag measurements from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project, a dedicated RM program with multi-object spectroscopy designed for RM over a wide redshift range. The lags are based on the 2014 spectroscopic light curves alone (32 epochs over 6 months) and focus on the Hbeta and MgII broad lines in the 100 lowest-redshift (z<0.8) quasars included in SDSS-RM; they represent a small subset of the lags that SDSS-RM (including 849 quasars to z~4.5) is expected to deliver. The reported preliminary lag measurements are for intermediate-luminosity quasars at 0.3<~z<0.8, including 9 Hbeta lags and 6 MgII lags, for the first time extending RM results to this redshift-luminosity regime and providing direct quasar black hole mass estimates over ~ half of cosmic time. The MgII lags also increase the number of known MgII lags by several-fold, and start to explore the utility of MgII for RM at high redshift. The location of these new lags at higher redshifts on the observed BLR size-luminosity relationship is statistically consistent with previous Hbeta results at z<0.3. However, an independent constraint on the relationship slope at z>0.3 is not yet possible due to the limitations in our current sample. Our results demonstrate the general feasibility and potential of multi-object RM for z>0.3 quasars.

We develop a simple yet comprehensive method to distinguish the underlying drivers of galaxy quenching, using the clustering and galaxy-galaxy lensing of red and blue galaxies in SDSS. Building on the iHOD framework developed by Zu & Mandelbaum (2015a), we consider two quenching scenarios: 1) a "halo" quenching model in which halo mass is the sole driver for turning off star formation in both centrals and satellites; and 2) a "hybrid" quenching model in which the quenched fraction of galaxies depends on their stellar mass while the satellite quenching has an extra dependence on halo mass. The two best-fit models describe the red galaxy clustering and lensing equally well, but halo quenching provides significantly better fits to the blue galaxies above $10^{11} M_\odot/h^2$. The halo quenching model also correctly predicts the average halo mass of the red and blue centrals, showing excellent agreement with the direct weak lensing measurements of locally brightest galaxies. Models in which quenching is not tied to halo mass, including an age-matching model in which galaxy colour depends on halo age at fixed $M_*$, fail to reproduce the observed halo mass for massive blue centrals. We find similar critical halo masses responsible for the quenching of centrals and satellites (~$1.5\times10^{12} M\odot/h^2$), hinting at a uniform quenching mechanism for both, e.g., the virial shock-heating of infalling gas. The success of the iHOD halo quenching model provides strong evidence that the physical mechanism that quenches star formation in galaxies is tied principally to the masses of their dark matter halos rather than the properties of their stellar components.

The mapping between the distributions of the observed galaxy stellar mass and the underlying dark matter halos provides the crucial link from theories of large-scale structure formation to interpreting the complex phenomena of galaxy formation and evolution. We develop a novel statistical method, based on the Halo Occupation Distribution model (HOD), to solve for this mapping by jointly fitting the galaxy clustering and the galaxy-galaxy lensing measured from the Sloan Digital Sky Survey (SDSS). The method, called the iHOD model, extracts maximum information from the survey by including ~80% more galaxies than the traditional HOD methods, and takes into account the incompleteness of the stellar mass samples in a statistically consistent manner. The derived stellar-to-halo mass relation not only explains the clustering and lensing of SDSS galaxies over almost four decades in stellar mass, but also successfully predicts the stellar mass functions observed in SDSS. Due to its capability of modelling significantly more galaxies, the iHOD is able to break the degeneracy between the logarithmic scatter in the stellar mass at fixed halo mass and the slope of the stellar-to-halo mass relation at high mass end, without the need to assume a strong prior on the scatter and/or use the stellar mass function as an input. We detect a decline of the scatter with halo mass, from 0.22 dex at below 10^{12} Msun to 0.18 dex at 10^{14} Msun. The model also enables stringent constraints on the satellite stellar mass functions at fixed halo mass, predicting a departure from the Schechter functional form in high mass halos. The iHOD model can be easily applied to existing and future spectroscopic datasets, greatly improving the statistical constraint on the stellar-to-halo mass relation compared to the traditional HOD methods within the same survey.

We describe the first results from a six-month long reverberation-mapping experiment in the ultraviolet based on 170 observations of the Seyfert 1 galaxy NGC 5548 with the Cosmic Origins Spectrograph on the Hubble Space Telescope. Significant correlated variability is found in the continuum and broad emission lines, with amplitudes ranging from ~30% to a factor of two in the emission lines and a factor of three in the continuum. The variations of all the strong emission lines lag behind those of the continuum, with He II 1640 lagging behind the continuum by ~2.5 days and Lyman alpha 1215, C IV 1550, and Si IV 1400 lagging by ~5-6 days. The relationship between the continuum and emission lines is complex. In particular, during the second half of the campaign, all emission-line lags increased by a factor of 1.3-2 and differences appear in the detailed structure of the continuum and emission-line light curves. Velocity-resolved cross-correlation analysis shows coherent structure in lag versus line-of-sight velocity for the emission lines; the high-velocity wings of C IV respond to continuum variations more rapidly than the line core, probably indicating higher velocity BLR clouds at smaller distances from the central engine. The velocity-dependent response of Lyman alpha, however, is more complex and will require further analysis.

Recent intensive Swift monitoring of the Seyfert 1 galaxy NGC 5548 yielded 282 usable epochs over 125 days across six UV/optical bands and the X-rays. This is the densest extended AGN UV/optical continuum sampling ever obtained, with a mean sampling rate <0.5 day. Approximately daily HST UV sampling was also obtained. The UV/optical light curves show strong correlations (r_max = 0.57 - 0.90) and the clearest measurement to date of interband lags. These lags are well-fit by a \tau propto \lambda^4/3 wavelength dependence, with a normalization that indicates an unexpectedly large disk radius of 0.35 +/- 0.05 lt-day at 1367 A, assuming a simple face-on model. The U-band shows a marginally larger lag than expected from the fit and surrounding bands, which could be due to Balmer continuum emission from the broad-line region as suggested by Korista and Goad. The UV/X-ray correlation is weaker (r_max < 0.45) and less consistent over time. This indicates that while Swift is beginning to measure UV/optical lags in general agreement with accretion disk theory (although the derived size is larger than predicted), the relationship with X-ray variability is less well understood. Combining this accretion disk size estimate with those from quasar microlensing studies suggests that AGN disk sizes scale approximately linearly with central black hole mass over a wide range of masses.

A large reverberation mapping study of the Seyfert 1 galaxy NGC 7469 has yielded emission-line lags for Hbeta 4861 and He II 4686 and a central black hole mass measurement of about 10 million solar masses, consistent with previous measurements. A very low level of variability during the monitoring campaign precluded meeting our original goal of recovering velocity-delay maps from the data, but with the new Hbeta measurement, NGC 7469 is no longer an outlier in the relationship between the size of the Hbeta-emitting broad-line region and the AGN luminosity. It was necessary to detrend the continuum and Hbeta and He II 4686 line light curves and those from archival UV data for different time-series analysis methods to yield consistent results.

The Sloan Digital Sky Survey Reverberation Mapping project (SDSS-RM) is a dedicated multi-object RM experiment that has spectroscopically monitored a sample of 849 broad-line quasars in a single 7 deg$^2$ field with the SDSS-III BOSS spectrograph. The RM quasar sample is flux-limited to i_psf=21.7 mag, and covers a redshift range of 0.1<z<4.5. Optical spectroscopy was performed during 2014 Jan-Jul dark/grey time, with an average cadence of ~4 days, totaling more than 30 epochs. Supporting photometric monitoring in the g and i bands was conducted at multiple facilities including the CFHT and the Steward Observatory Bok telescopes in 2014, with a cadence of ~2 days and covering all lunar phases. The RM field (RA, DEC=14:14:49.00, +53:05:00.0) lies within the CFHT-LS W3 field, and coincides with the Pan-STARRS 1 (PS1) Medium Deep Field MD07, with three prior years of multi-band PS1 light curves. The SDSS-RM 6-month baseline program aims to detect time lags between the quasar continuum and broad line region (BLR) variability on timescales of up to several months (in the observed frame) for ~10% of the sample, and to anchor the time baseline for continued monitoring in the future to detect lags on longer timescales and at higher redshift. SDSS-RM is the first major program to systematically explore the potential of RM for broad-line quasars at z>0.3, and will investigate the prospects of RM with all major broad lines covered in optical spectroscopy. SDSS-RM will provide guidance on future multi-object RM campaigns on larger scales, and is aiming to deliver more than tens of BLR lag detections for a homogeneous sample of quasars. We describe the motivation, design and implementation of this program, and outline the science impact expected from the resulting data for RM and general quasar science.

After the All-Sky Automated Survey for SuperNovae (ASAS-SN) discovered a significant brightening of the inner region of NGC 2617, we began a ~70 day photometric and spectroscopic monitoring campaign from the X-ray through near-infrared (NIR) wavelengths. We report that NGC 2617 went through a dramatic outburst, during which its X-ray flux increased by over an order of magnitude followed by an increase of its optical/ultraviolet (UV) continuum flux by almost an order of magnitude. NGC 2617, classified as a Seyfert 1.8 galaxy in 2003, is now a Seyfert 1 due to the appearance of broad optical emission lines and a continuum blue bump. Such "changing look Active Galactic Nuclei (AGN)" are rare and provide us with important insights about AGN physics. Based on the Hbeta line width and the radius-luminosity relation, we estimate the mass of central black hole to be (4 +/- 1) x 10^7 M_sun. When we cross-correlate the light curves, we find that the disk emission lags the X-rays, with the lag becoming longer as we move from the UV (2-3 days) to the NIR (6-9 days). Also, the NIR is more heavily temporally smoothed than the UV. This can largely be explained by a simple model of a thermally emitting thin disk around a black hole of the estimated mass that is illuminated by the observed, variable X-ray fluxes.

Infrared modifications of General Relativity (GR) can be revealed by comparing the mass of galaxy clusters estimated from weak lensing to that from infall kinematics. We measure the 2D galaxy velocity distribution in the cluster infall region by applying the galaxy infall kinematics (GIK) model developed by Zu and Weinberg (2013) to two suites of f(R) and Galileon modified gravity simulations. Despite having distinct screening mechanisms, namely, the Chameleon and the Vainshtein effects, the f(R) and Galileon clusters exhibit very similar deviations in their GIK profiles from GR, with ~ 100-200 k/s enhancement in the characteristic infall velocity at r=5 Mpc/h and 50-100 km/s broadening in the radial and tangential velocity dispersions across the entire infall region, for clusters with mass ~ 10^{14} Msol/h at z=0.25. These deviations are detectable via the GIK reconstruction of the redshift--space cluster-galaxy cross-correlation function, xi_cg^s(r_p,r_\pi), which shows ~ 1-2 Mpc/h increase in the characteristic line-of-sight distance r_\pi^c at r_p<6 Mpc/h from GR predictions. With overlapping deep imaging and large redshift surveys in the future, we expect that the GIK modelling of xi_cg^s, in combination with the stacked weak lensing measurements, will provide powerful diagnostics of modified gravity theories and the origin of cosmic acceleration.

We present new stellar velocity dispersion measurements for four luminous quasars with the NIFS instrument and the ALTAIR laser guide star adaptive optics system on the Gemini North 8-m telescope. Stellar velocity dispersion measurements and measurements of the supermassive black hole masses in luminous quasars are necessary to investigate the coevolution of black holes and galaxies, trace the details of accretion, and probe the nature of feedback. We find that higher-luminosity quasars with higher-mass black holes are not offset with respect to the MBH-sigma relation exhibited by lower-luminosity AGNs with lower-mass black holes, nor do we see correlations with galaxy morphology. As part of this analysis, we have recalculated the virial products for the entire sample of reverberation-mapped AGNs and used these data to redetermine the mean virial factor hfi that places the reverberation data on the quiescent M_BH-sigma relation. With our updated measurements and new additions to the AGN sample, we obtain <f> = 4.31 +/- 1.05, which is slightly lower than, but consistent with, most previous determinations.

The Astro2010 Decadal Survey recommended a Wide Field Infrared Survey Telescope (WFIRST) as its top priority for a new large space mission. The report of the WFIRST-AFTA Science Definition Team (SDT) presents a Design Reference Mission for WFIRST that employs one of the 2.4-m, Hubble-quality mirror assemblies recently made available to NASA. The 2.4-m primary mirror enables a mission with greater sensitivity and higher angular resolution than the smaller aperture designs previously considered for WFIRST, increasing both the science return of the primary surveys and the capabilities of WFIRST as a Guest Observer facility. The option of adding an on-axis, coronagraphic instrument would enable imaging and spectroscopic studies of planets around nearby stars. This short article, produced as a companion to the SDT report, summarizes the key points of the WFIRST-2.4 DRM. It highlights the remarkable opportunity that the 2.4-m telescope affords for advances in many fields of astrophysics and cosmology, including dark energy, the demographics and characterization of exoplanets, the evolution of galaxies and quasars, and the stellar populations of the Milky Way and its neighbors.

The Astro2010 Decadal Survey recommended a Wide Field Infrared Survey Telescope (WFIRST) as its top priority for a new large space mission. As conceived by the decadal survey, WFIRST would carry out a dark energy science program, a microlensing program to determine the demographics of exoplanets, and a general observing program utilizing its ultra wide field. In October 2012, NASA chartered a Science Definition Team (SDT) to produce, in collaboration with the WFIRST Project Office at GSFC and the Program Office at JPL, a Design Reference Mission (DRM) for an implementation of WFIRST using one of the 2.4-m, Hubble-quality mirror assemblies recently made available to NASA. This DRM builds on the work of the earlier WFIRST SDT, reported by Green et al. (2012). The 2.4-m primary mirror enables a mission with greater sensitivity and higher angular resolution than the 1.3-m and 1.1-m designs considered previously, increasing both the science return of the primary surveys and the capabilities of WFIRST as a Guest Observer facility. The option of adding an on-axis, coronagraphic instrument would enable imaging and spectroscopic studies of planets around nearby stars. This document presents the final report of the SDT.

The large scale infall of galaxies around massive clusters provides a potentially powerful diagnostic of structure growth, dark energy, and cosmological deviations from General Relativity. We develop and test a method to recover galaxy infall kinematics (GIK) from measurements of the redshift-space cluster-galaxy cross-correlation function \xi_{cg}(r_p,r_\pi). Using galaxy and halo samples from the Millennium simulation, we calibrate an analytic model of the galaxy kinematic profiles comprised of a virialized component with an isotropic Gaussian velocity distribution and an infall component described by a skewed 2D t-distribution with a characteristic infall velocity v_r and separate radial and tangential dispersions. We show that convolving the real-space cross-correlation function with this velocity distribution accurately predicts the redshift-space \xi_{cg}, and we show that measurements of \xi_{cg} can be inverted to recover the four distinct elements of the GIK profiles. These in turn provide diagnostics of cluster mass profiles, and we expect the characteristic infall velocity v_r(r) in particular to be insensitive to galaxy formation physics that can affect velocity dispersions within halos. As a proof of concept we measure \xi_{cg} for rich galaxy groups in the Sloan Digital Sky Survey and recover GIK profiles for groups in two bins of central galaxy stellar mass. The higher mass bin has a v_r(r) curve very similar to that of 10^{14} Msun halos in the Millennium simulation, and the recovered kinematics follow the expected trends with mass. GIK modeling of cluster-galaxy cross-correlations can be a valuable complement to stacked weak lensing analyses, allowing novel tests of modified gravity theories that seek to explain cosmic acceleration.

The damped random walk (DRW) model is increasingly used to model the variability in quasar optical light curves, but it is still uncertain whether the DRW model provides an adequate description of quasar optical variability across all time scales. Using a sample of OGLE quasar light curves, we consider four modifications to the DRW model by introducing additional parameters into the covariance function to search for deviations from the DRW model on both short and long time scales. We find good agreement with the DRW model on time scales that are well sampled by the data (from a month to a few years), possibly with some intrinsic scatter in the additional parameters, but this conclusion depends on the statistical test employed and is sensitive to whether the estimates of the photometric errors are correct to within ~10%. On very short time scales (below a few months), we see some evidence of the existence of a cutoff below which the correlation is stronger than the DRW model, echoing the recent finding of Mushotzky et al. (2011) using quasar light curves from Kepler. On very long time scales (> a few years), the light curves do not constrain models well, but are consistent with the DRW model.

We present velocity-resolved reverberation results for five active galactic nuclei. We recovered velocity-delay maps using the maximum-entropy method for four objects: Mrk 335, Mrk 1501, 3C120, and PG2130+099. For the fifth, Mrk 6, we were only able to measure mean time delays in different velocity bins of the Hbeta emission line. The four velocity-delay maps show unique dynamical signatures for each object. For 3C120, the Balmer lines show kinematic signatures consistent with both an inclined disk and infalling gas, but the HeII 4686 emission line is suggestive only of inflow. The Balmer lines in Mrk 335, Mrk 1501, and PG 2130+099 show signs of infalling gas, but the HeII emission in Mrk 335 is consistent with an inclined disk. We also see tentative evidence of combined virial motion and infalling gas from the velocity-binned analysis of Mrk 6. The maps for 3C120 and Mrk 335 are two of the most clearly defined velocity-delay maps to date. These maps constitute a large increase in the number of objects for which we have resolved velocity-delay maps and provide evidence supporting the reliability of reverberation-based black hole mass measurements.

We derive constraints on the matter density \Om and the amplitude of matter clustering \sig8 from measurements of large scale weak lensing (projected separation R=5-30\hmpc) by clusters in the Sloan Digital Sky Survey MaxBCG catalog. The weak lensing signal is proportional to the product of \Om and the cluster-mass correlation function \xicm. With the relation between optical richness and cluster mass constrained by the observed cluster number counts, the predicted lensing signal increases with increasing \Om or \sig8, with mild additional dependence on the assumed scatter between richness and mass. The dependence of the signal on scale and richness partly breaks the degeneracies among these parameters. We incorporate external priors on the richness-mass scatter from comparisons to X-ray data and on the shape of the matter power spectrum from galaxy clustering, and we test our adopted model for \xicm against N-body simulations. Using a Bayesian approach with minimal restrictive priors, we find \sig8(\Om/0.325)^{0.501}=0.828 +/- 0.049, with marginalized constraints of \Om=0.325_{-0.067}^{+0.086} and \sig8=0.828_{-0.097}^{+0.111}, consistent with constraints from other MaxBCG studies that use weak lensing measurements on small scales (R<=2\hmpc). The (\Om,\sig8) constraint is consistent with and orthogonal to the one inferred from WMAP CMB data, reflecting agreement with the structure growth predicted by GR for an LCDM cosmological model. A joint constraint assuming LCDM yields \Om=0.298 +/- 0.020 and \sig8=0.831 +/- 0.020. Our cosmological parameter errors are dominated by the statistical uncertainties of the large scale weak lensing measurements, which should shrink sharply with current and future imaging surveys.

We present the results from a detailed analysis of photometric and spectrophotometric data on five Seyfert 1 galaxies observed as a part of a recent reverberation mapping program. The data were collected at several observatories over a 140-day span beginning in 2010 August and ending in 2011 January. We obtained high sampling-rate light curves for Mrk 335, Mrk 1501, 3C120, Mrk 6, and PG2130+099, from which we have measured the time lag between variations in the 5100 Angstrom continuum and the H-beta broad emission line. We then used these measurements to calculate the mass of the supermassive black hole at the center of each of these galaxies. Our new measurements substantially improve previous measurements of MBH and the size of the broad line-emitting region for four sources and add a measurement for one new object. Our new measurements are consistent with photoionization physics regulating the location of the broad line region in active galactic nuclei.

We present the first results from a detailed analysis of photometric and spectrophotometric data on the narrow-line Seyfert 1 galaxy Mrk 335, collected over a 120-day span in the fall of 2010. From these data we measure the lag in the He II 4686 broad emission line relative to the optical continuum to be 2.7 \pm 0.6 days and the lag in the H\beta 4861 broad emission line to be 13.9 \pm 0.9 days. Combined with the line width, the He II lag yields a black hole mass, MBH = (2.6 \pm 0.8)\times 10^7 Msun. This measurement is consistent with measurements made using the H\beta 4861 line, suggesting that the He II emission originates in the same region as H\beta, but at a much smaller radius. This constitutes the first robust lag measurement for a high-ionization line in a narrow-line Seyfert 1 galaxy.

Motivated by recent progress in the statistical modeling of quasar variability, we develop a new approach to measuring emission-line reverberation lags to estimate the size of broad-line regions (BLRs) in active galactic nuclei. Assuming that all emission-line light curves are scaled, smoothed, and displaced versions of the continuum, this alternative approach fits the light curves directly using a damped random walk model and aligns them to recover the time lag and its statistical confidence limits. We introduce the mathematical formalism of this approach and demonstrate its ability to cope with some of the problems for traditional methods, such as irregular sampling, correlated errors, and seasonal gaps. We redetermine the lags for 87 emission lines in 31 quasars and reassess the BLR size--luminosity relationship using 60 H-beta lags. We confirm the general results from the traditional cross-correlation methods, with a few exceptions. Our method, however, also supports a broad range of extensions. In particular, it can simultaneously fit multiple lines and continuum light curves which improves the lag estimate for the lines and provides estimates of the error correlations between them. Determining these correlations is of particular importance for interpreting emission-line velocity--delay maps. We can also include parameters for luminosity-dependent lags or line responses. We use this to detect the scaling of the BLR size with continuum luminosity in NGC 5548.

We place constraints on the average density (Omega_m) and clustering amplitude (sigma_8) of matter using a combination of two measurements from the Sloan Digital Sky Survey: the galaxy two-point correlation function, w_p, and the mass-to-galaxy-number ratio within galaxy clusters, M/N, analogous to cluster M/L ratios. Our w_p measurements are obtained from DR7 while the sample of clusters is the maxBCG sample, with cluster masses derived from weak gravitational lensing. We construct non-linear galaxy bias models using the Halo Occupation Distribution (HOD) to fit both w_p and M/N for different cosmological parameters. HOD models that match the same two-point clustering predict different numbers of galaxies in massive halos when Omega_m or sigma_8 is varied, thereby breaking the degeneracy between cosmology and bias. We demonstrate that this technique yields constraints that are consistent and competitive with current results from cluster abundance studies, even though this technique does not use abundance information. Using w_p and M/N alone, we find Omega_m^0.5*sigma_8=0.465+/-0.026, with individual constraints of Omega_m=0.29+/-0.03 and sigma_8=0.85+/-0.06. Combined with current CMB data, these constraints are Omega_m=0.290+/-0.016 and sigma_8=0.826+/-0.020. All errors are 1-sigma. The systematic uncertainties that the M/N technique are most sensitive to are the amplitude of the bias function of dark matter halos and the possibility of redshift evolution between the SDSS Main sample and the maxBCG sample. Our derived constraints are insensitive to the current level of uncertainties in the halo mass function and in the mass-richness relation of clusters and its scatter, making the M/N technique complementary to cluster abundances as a method for constraining cosmology with future galaxy surveys.

Recently Menard et al. detected a subtle but systematic change in the mean color of quasars as a function of their projected separation from foreground galaxies, extending to comoving separations of ~10Mpc/h, which they interpret as a signature of reddening by intergalactic dust. We present theoretical models of this remarkable observation, using SPH cosmological simulations of a (50Mpc/h)^3 volume. Our primary model uses a simulation with galactic winds and assumes that dust traces the intergalactic metals. The predicted galaxy-dust correlation function is similar in form to the galaxy-mass correlation function, and reproducing the MSFR data requires a dust-to-metal mass ratio of 0.24, about half the value in the Galactic ISM. Roughly half of the reddening arises in dust that is more than 100Kpc/h from the nearest massive galaxy. We also examine a simulation with no galactic winds, which predicts a much smaller fraction of intergalactic metals (3% vs. 35%) and therefore requires an unphysical dust-to-metal ratio of 2.18 to reproduce the MSFR data. In both models, the signal is dominated by sightlines with E(g-i)=0.001-0.1. The no-wind simulation can be reconciled with the data if we also allow reddening to arise in galaxies up to several x 10^10 Msun. The wind model predicts a mean visual extinction of A_V ~0.0133 mag out to z=0.5, with a sightline-to-sightline dispersion similar to the mean, which could be significant for future supernova cosmology studies. Reproducing the MSFR results in these simulations requires that a large fraction of ISM dust survive its expulsion from galaxies and its residence in the intergalactic medium. Future observational studies that provide higher precision and measure the dependence on galaxy type and environment will allow detailed tests for models of enriched galactic outflows and the survival of IG dust.

Galaxy formation inside dark matter halos, as well as the halo formation itself, can be affected by large-scale environments. Evaluating the imprints of environmental effects on galaxy clustering is crucial for precise cosmological constraints with data from galaxy redshift surveys. We investigate such an environmental impact on both real-space and redshift-space galaxy clustering statistics using a semi-analytic model derived from the Millennium Simulation. We compare clustering statistics from original SAM galaxy samples and shuffled ones with environmental influence on galaxy properties eliminated. Among the luminosity-threshold samples examined, the one with the lowest threshold luminosity (~0.2L_*) is affected by environmental effects the most, which has a ~10% decrease in the real-space two-point correlation function (2PCF) after shuffling. By decomposing the 2PCF into five different components based on the source of pairs, we show that the change in the 2PCF can be explained by the age and richness dependence of halo clustering. The 2PCFs in redshift space are found to change in a similar manner after shuffling. If the environmental effects are neglected, halo occupation distribution modeling of the real-space and redshift-space clustering may have a less than 6.5% systematic uncertainty in constraining beta from the most affected SAM sample and have substantially smaller uncertainties from the other, more luminous samples. We argue that the effect could be even smaller in reality. In the Appendix, we present a method to decompose the 2PCF, which can be applied to measure the two-point auto-correlation functions of galaxy sub-samples in a volume-limited galaxy sample and their two-point cross-correlation functions in a single run utilizing only one random catalog.