Source author record

Ying Zu

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

23works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

23 published item(s)

preprint2025arXiv

Introduction to the Chinese Space Station Survey Telescope (CSST)

The Chinese Space Station Survey Telescope (CSST) is an upcoming Stage-IV sky survey telescope, distinguished by its large field of view (FoV), high image quality, and multi-band observation capabilities. It can simultaneously conduct precise measurements of the Universe by performing multi-color photometric imaging and slitless spectroscopic surveys. The CSST is equipped with five scientific instruments, i.e. Multi-band Imaging and Slitless Spectroscopy Survey Camera (SC), Multi-Channel Imager (MCI), Integral Field Spectrograph (IFS), Cool Planet Imaging Coronagraph (CPI-C), and THz Spectrometer (TS). Using these instruments, CSST is expected to make significant contributions and discoveries across various astronomical fields, including cosmology, galaxies and active galactic nuclei (AGN), the Milky Way and nearby galaxies, stars, exoplanets, Solar System objects, astrometry, and transients and variable sources. This review aims to provide a comprehensive overview of the CSST instruments, observational capabilities, data products, and scientific potential.

preprint2024arXiv

AGN STORM 2: V. Anomalous Behavior of the CIV Light Curve in Mrk 817

An intensive reverberation mapping campaign on the Seyfert 1 galaxy Mrk817 using the Cosmic Origins Spectrograph (COS) on the Hubble Space Telescope (HST) revealed significant variations in the response of the broad UV emission lines to fluctuations in the continuum emission. The response of the prominent UV emission lines changes over a $\sim$60-day duration, resulting in distinctly different time lags in the various segments of the light curve over the 14 months observing campaign. One-dimensional echo-mapping models fit these variations if a slowly varying background is included for each emission line. These variations are more evident in the CIV light curve, which is the line least affected by intrinsic absorption in Mrk817 and least blended with neighboring emission lines. We identify five temporal windows with distinct emission line response, and measure their corresponding time delays, which range from 2 to 13 days. These temporal windows are plausibly linked to changes in the UV and X-ray obscuration occurring during these same intervals. The shortest time lags occur during periods with diminishing obscuration, whereas the longest lags occur during periods with rising obscuration. We propose that the obscuring outflow shields the ultraviolet broad lines from the ionizing continuum. The resulting change in the spectral energy distribution of the ionizing continuum, as seen by clouds at a range of distances from the nucleus, is responsible for the changes in the line response.

preprint2022arXiv

\textsc{The Three Hundred} project: The \textsc{Gizmo-Simba} run

We introduce \textsc{Gizmo-Simba}, a new suite of galaxy cluster simulations within \textsc{The Three Hundred} project. \textsc{The Three Hundred} consists of zoom re-simulations of 324 clusters with $M_{200}\gtrsim 10^{14.8}M_\odot$ drawn from the MultiDark-Planck $N$-body simulation, run using several hydrodynamic and semi-analytic codes. The \textsc{Gizmo-Simba} suite adds a state-of-the-art galaxy formation model based on the highly successful {\sc Simba} simulation, mildly re-calibrated to match $z=0$ cluster stellar properties. Comparing to \textsc{The Three Hundred} zooms run with \textsc{Gadget-X}, we find intrinsic differences in the evolution of the stellar and gas mass fractions, BCG ages, and galaxy colour-magnitude diagrams, with \textsc{Gizmo-Simba} generally providing a good match to available data at $z \approx 0$. \textsc{Gizmo-Simba}'s unique black hole growth and feedback model yields agreement with the observed BH scaling relations at the intermediate-mass range and predicts a slightly different slope at high masses where few observations currently lie. \textsc{Gizmo-Simba} provides a new and novel platform to elucidate the co-evolution of galaxies, gas, and black holes within the densest cosmic environments.

preprint2022arXiv

An Extended Halo-based Group/Cluster finder: application to the DESI legacy imaging surveys DR8

We extend the halo-based group finder developed by \citet[][]{Yang2005a} to use data {\it simultaneously} with either photometric or spectroscopic redshifts. A mock galaxy redshift survey constructed from a high-resolution N-body simulation is used to evaluate the performance of this extended group finder. For galaxies with magnitude ${\rm z\le 21}$ and redshift $0<z\le 1.0$ in the DESI legacy imaging surveys (the Legacy Surveys), our group finder successfully identifies more than 60\% of the members in about $90\%$ of halos with mass $\ga 10^{12.5}\msunh$. Detected groups with mass $\ga 10^{12.0}\msunh$ have a purity (the fraction of true groups) greater than 90\%. The halo mass assigned to each group has an uncertainty of about 0.2 dex at the high mass end $\ga 10^{13.5}\msunh$ and 0.40 dex at the low mass end. Groups with more than 10 members have a redshift accuracy of $\sim 0.008$. We apply this group finder to the Legacy Surveys DR8 and find 5.2 Million groups with at least 3 members. About 387,000 of these groups have at least 10 members. The resulting catalog containing 3D coordinates, richness, halo masses, and total group luminosities, is made publicly available.

preprint2022arXiv

Elucidating Galaxy Assembly Bias in SDSS

We investigate the level of galaxy assembly bias in the Sloan Digital Sky Survey (SDSS) main galaxy sample using ELUCID, a state-of-the-art constrained simulation that accurately reconstructed the initial density perturbations within the SDSS volume. On top of the ELUCID haloes, we develop an extended HOD model that includes the assembly bias of central and satellite galaxies, parameterized as $\mathcal{Q}_\mathrm{cen}$ and $\mathcal{Q}_\mathrm{sat}$, respectively, to predict a suite of one- and two-point observables. In particular, our fiducial constraint employs the probability distribution of the galaxy number counts measured on $8\,\mathrm{Mpc}\,h^{-1}$ scales $N_8^g$ and the projected cross-correlation functions of quintiles of galaxies selected by $N_8^g$ with our entire galaxy sample. We perform extensive tests of the efficacy of our method by fitting the same observables to mock data using both constrained and non-constrained simulations. We discover that in many cases the level of cosmic variance between the two simulations can produce biased constraints that lead to an erroneous detection of galaxy assembly bias if the non-constrained simulation is used. When applying our method to the SDSS data, the ELUCID reconstruction effectively removes an otherwise strong degeneracy between cosmic variance and galaxy assembly bias in SDSS, enabling us to derive an accurate and stringent constraint on the latter. Our fiducial ELUCID constraint, for galaxies above a stellar mass threshold $M_*{=}10^{10.2}\,h^{-2}\,M_\odot$, is $\mathcal{Q}_\mathrm{cen}{=}{-}0.09\pm{0.05}$ and $\mathcal{Q}_\mathrm{sat}{=}0.09\pm{0.10}$, indicating no evidence for a significant~($>2σ$) galaxy assembly bias in the local Universe probed by SDSS. Finally, our method provides a promising path to the robust modelling of the galaxy-halo connection within future surveys like DESI and PFS.

preprint2022arXiv

Strong Conformity and Assembly Bias: Towards a Physical Understanding of the Galaxy-Halo Connection in SDSS Clusters

Understanding the physical connection between cluster galaxies and massive haloes is key to mitigating systematic uncertainties in next-generation cluster cosmology. We develop a novel method to infer the level of conformity between the stellar mass of the brightest central galaxies~(BCGs) $M_*^{BCG}$ and the satellite richness $λ$, defined as their correlation coefficient $ρ_{cc}$ at fixed halo mass, using the abundance and weak lensing of SDSS clusters as functions of $M_*^{BCG}$ and $λ$. We detect a halo mass-dependent conformity as $ρ_{cc}{=}0.60{+}0.08\ln(M_h/3{\times}10^{14}M_{\odot}/h)$. The strong conformity successfully resolves the "halo mass equality" conundrum discovered in Zu et al. 2021 -- when split by $M_*^{BCG}$ at fixed $λ$, the low and high-$M_*^{BCG}$ clusters have the same average halo mass despite having a $0.34$ dex discrepancy in average $M_*^{BCG}$. On top of the best-fitting conformity model, we develop a cluster assembly bias~(AB) prescription calibrated against the CosmicGrowth simulation, and build a conformity+AB model for the cluster weak lensing measurements. Our model predicts that with a ${\sim}20\%$ lower halo concentration $c$, the low-$M_*^{BCG}$ clusters are ${\sim}10\%$ more biased than the high-$M_*^{BCG}$ systems, in excellent agreement with the observations. We also show that the observed conformity and assembly bias are unlikely due to projection effects. Finally, we build a toy model to argue that while the early-time BCG-halo co-evolution drives the $M_*^{BCG}$-$c$ correlation, the late-time dry merger-induced BCG growth naturally produces the $M_*^{BCG}$-$λ$ conformity despite the well-known anti-correlation between $λ$ and $c$. Our method paves the path towards simultaneously constraining cosmology and cluster formation with future cluster surveys.

preprint2022arXiv

The Sphere of Influence of the Bright Central Galaxies in the Diffuse Light of SDSS Clusters

The bright central galaxies (BCGs) dominate the inner portion of the diffuse cluster light, but it is still unclear where the intracluster light (ICL) takes over. To investigate the BCG-ICL transition, we stack the images of ${\sim}3000$ clusters between $0.2{<}z{<}0.3$ in the SDSS $gri$ bands, and measure their BCG+ICL stellar surface mass profile $Σ_{*}^{\texttt{B+I}}$ down to $3{\times}10^4\,M_{\odot}/\mathrm{kpc}^{2}$ at $R{\simeq}1\,\mathrm{Mpc}$ (${\sim}32$ mag/arcsec$^2$ in the $r$-band). We develop a physically-motivated method to decompose $Σ_{*}^{\texttt{B+I}}$ into three components, including an inner de Vaucouleurs' profile, an outer ICL that follows the dark matter distribution measured from weak lensing, and an intriguing transitional component between 70 and 200 kpc. To investigate the origin of this transition, we split the clusters into two subsamples by their BCG stellar mass $M_*^{\mathrm{BCG}}$ (mass enclosed roughly within 50 kpc) while making sure they have the same distribution of satellite richness. The $Σ_{*}^{\texttt{B+I}}$ profiles of the two subsamples differ by more than a factor of two at $R{<}50\,\mathrm{kpc}$, consistent with their 0.34 dex difference in $M_*^{\mathrm{BCG}}$, whereas on scales beyond 400 kpc the two profiles converge to the same amplitudes, suggesting a satellite-stripping origin of the outer ICL. Remarkably, however, the discrepancy between the two $Σ_{*}^{\texttt{B+I}}$ profiles persist at above $50\%$ level on all scales below 200 kpc, thereby revealing the BCG sphere of influence with radius $R_{\mathrm{SOI}}{\simeq}$ 200 kpc. Finally, we speculate that the surprisingly large sphere of influence of the BCG is tied to the elevated escape velocity profile within $r_s$, the characteristic radius of the dark matter haloes.

preprint2020arXiv

HI gas content of SDSS galaxies revealed by ALFALFA: implications for the mass-metallicity relation and the environmental dependence of HI in the local Universe

The neutral hydrogen~(HI) gas is an important barometer of recent star formation and metal enrichment activities in galaxies. I develop a novel statistical method for predicting the HI-to-stellar mass ratio $f_{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. I eliminate the impact of the Malmquist bias against HI-deficient systems on the $f_{gas}$ predictor by properly accounting for the HI detection probability of each galaxy in the analysis. The best-fitting $f_{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_{gas}$ from the expected value, I 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, I show that the metallicity dependence on HI is likely more fundamental than that on specific star formation rate. In addition, I find 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 $ρ_{cc}{=}-0.18$. This weak anti-correlation also successfully explains the observed dependence of HI clustering on $f_{gas}$. My method provides a useful framework for learning HI gas evolution from the synergy between future HI and optical galaxy surveys.

preprint2016arXiv

Application of Stochastic Modeling to Analysis of Photometric Reverberation Mapping Data

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.

preprint2016arXiv

Strong bimodality in the host halo mass of central galaxies from galaxy-galaxy lensing

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σ$ 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.

preprint2016arXiv

The DESI Experiment Part I: Science,Targeting, and Survey Design

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-$α$ 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.

preprint2016arXiv

The DESI Experiment Part II: Instrument Design

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= λ/Δλ$ 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.

preprint2015arXiv

Mapping stellar content to dark matter halos using galaxy clustering and galaxy-galaxy lensing in the SDSS DR7

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.

preprint2015arXiv

Mapping stellar content to dark matter halos. II. Halo mass is the main driver of galaxy quenching

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.

preprint2015arXiv

The Sloan Digital Sky Survey Reverberation Mapping Project: First Broad-line Hbeta and MgII Lags at z>~0.3 from six-Month Spectroscopy

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.

preprint2014arXiv

The Sloan Digital Sky Survey Reverberation Mapping Project: Technical Overview

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.

preprint2013arXiv

Galaxy Infall Kinematics as a Test of Modified Gravity

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_π), which shows ~ 1-2 Mpc/h increase in the characteristic line-of-sight distance r_π^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.

preprint2013arXiv

Is Quasar Optical Variability a Damped Random Walk?

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.

preprint2013arXiv

The Redshift-Space Cluster-Galaxy Cross-Correlation Function: I. Modeling Galaxy Infall onto Millennium Simulation Clusters and SDSS Groups

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 ξ_{cg}(r_p,r_π). 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 ξ_{cg}, and we show that measurements of ξ_{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 ξ_{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.

preprint2012arXiv

Cosmological Constraints from the Large Scale Weak Lensing of SDSS MaxBCG Clusters

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.

preprint2011arXiv

An Alternative Approach To Measuring Reverberation Lags in Active Galactic Nuclei

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.

preprint2011arXiv

Cosmological Constraints from Galaxy Clustering and the Mass-to-Number Ratio of Galaxy Clusters

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.

preprint2011arXiv

Intergalactic Dust Extinction in Hydrodynamic Cosmological Simulations

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.