Scott F. Anderson

DETECTION OF THE BARYON ACOUSTIC PEAK IN THE LARGE-SCALE CORRELATION FUNCTION OF SDSS LUMINOUS RED GALAXIES

We present the large-scale correlation function measured from a spectroscopic sample of 46,748 luminous red galaxies from the Sloan Digital Sky Survey. The survey region covers 0.72h −3 Gpc 3 over 3816 square degrees and 0.16 < z < 0.47, making it the best sample yet for the study of large-scale structure. We find a well-detected peak in the correlation function at 100h −1 Mpc separation that is an excellent match to the predicted shape and location of the imprint of the recombination-epoch acoustic oscillations on the low-redshift clustering of matter. This detection demonstrates the linear growth of structure by gravitational instability between z ≈ 1000 and the present and confirms a firm prediction of the standard cosmological theory. The acoustic peak provides a standard ruler by which we can measure the ratio of the distances to z = 0.35 and z = 1089 to 4% fractional accuracy and the absolute distance to z = 0.35 to 5% accuracy. From the overall shape of the correlation function, we measure the matter density mh 2 to 8% and find agreement with the value from cosmic microwave background (CMB) anisotropies. Independent of the constraints provided by the CMB acoustic scale, we find m = 0.273 ±0.025+0.123(1+ w0)+0.137K. Including the CMB acoustic scale, we find that the spatial curvature is K = −0.010 ± 0.009 if the dark energy is a cosmological constant. More generally, our results provide a measurement of cosmological distance, and hence an argument for dark energy, based on a geometric method with the same simple physics as the microwave background anisotropies. The standard cosmological model convincingly passes these new and robust tests of its fundamental properties. Subject headings: cosmology: observations — large-scale structure of the universe — distance scale — cosmological parameters — cosmic microwave background — galaxies: elliptical and lenticular, cD

Cosmological constraints from the SDSS luminous red galaxies

We measure the large-scale real-space power spectrum P(k) using luminous red galaxies (LRGs) in the Sloan Digital Sky Survey (SDSS) and use this measurement to sharpen constraints on cosmological parameters from the Wilkinson Microwave Anisotropy Probe (WMAP). We employ a matrix-based power spectrum estimation method using Pseudo-Karhunen-Loeve eigenmodes, producing uncorrelated minimum-variance measurements in 20 k-bands of both the clustering power and its anisotropy due to redshift-space distortions, with narrow and well-behaved window functions in the range 0.01h/Mpc 0.1h/Mpc and associated nonlinear complications, yet agree well with more aggressive published analyses where nonlinear modeling is crucial.

Cosmological parameter analysis including SDSS lyα forest and galaxy bias : Constraints on the primordial spectrum of fluctuations, neutrino mass, and dark energy

We combine the constraints from the recent Lyα forest analysis of the Sloan Digital Sky Survey (SDSS) and the SDSS galaxy bias analysis with previous constraints from SDSS galaxy clustering, the latest supernovae, and 1st year WMAP cosmic microwave background anisotropies. We find significant improvements on all of the cosmological parameters compared to previous constraints, which highlights the importance of combining Lyα forest constraints with other probes. combining WMAP and the Lyα forest we find for the primordial slope ns = 0:98±0:02. We see no evidence of running, dn/=d lnk 0:003±0:010, a factor of 3 improvement over previous constraints. We also find no evidence of tensors, r < 0:36 (95% c.l.). Inflationary models predict the absence of running and many among them satisfy these constraints, particularly negative curvature models such as those based on spontaneous symmetry breaking. A positive correlation between tensors and primordial slope disfavors chaotic inflation-type models with steep slopes: while the V αo 2 model is within the 2-sigma contour, V αo4 is outside the 3- sigma contour. For the amplitude we find σ8 = 0:90 ± 0:03 from the Lyα forest and WMAP alone. We find no evidence of neutrino mass: for the case of 3 massive neutrino families with an inflationary prior, Σmv < 0:42 eV and the mass of lightest neutrino is m1 < 0:13 eV at 95% c.l. For the 3 massless +1 massive neutrino case we find mv < 0:79 eV for the massive neutrino, excluding at 95% c.l. all neutrino mass solutions compatible with the LSND results. We explore dark energy constraints in models with a fairly general time dependence of dark energy equation of state, finding Ωλ =0:72± 0:02, w(z = 0:3) = 0:98+0.10 -0.12,the latter changing to w(z = 0:3) = -0.92+0.09-0.10 if tensors are allowed. We find no evidence for variation of the equation of state with redshift, w(z = 1) = -1.03+0.21-0.28. These results rely on the current understanding of the Lyα forest and other probes, which need to be explored further both observationally and theoretically, but extensive tests reveal no evidence of inconsistency among different data sets used here.

Spectral Energy Distributions and Multiwavelength Selection of Type 1 Quasars

We present an analysis of the mid-infrared (MIR) and optical properties of type 1 (broad-line) quasars detected by the Spitzer Space Telescope. The MIR color-redshift relation is characterized to z ~ 3, with predictions to z = 7. We demonstrate how combining MIR and optical colors can yield even more efficient selection of active galactic nuclei (AGNs) than MIR or optical colors alone. Composite spectral energy distributions (SEDs) are constructed for 259 quasars with both Sloan Digital Sky Survey and Spitzer photometry, supplemented by near-IR, GALEX, VLA, and ROSAT data, where available. We discuss how the spectral diversity of quasars influences the determination of bolometric luminosities and accretion rates; assuming the mean SED can lead to errors as large as 50% for individual quasars when inferring a bolometric luminosity from an optical luminosity. Finally, we show that careful consideration of the shape of the mean quasar SED and its redshift dependence leads to a lower estimate of the fraction of reddened/obscured AGNs missed by optical surveys as compared to estimates derived from a single mean MIR to optical flux ratio.

A Survey of z > 5.8 Quasars in the Sloan Digital Sky Survey. I. Discovery of Three New Quasars and the Spatial Density of Luminous Quasars at z ∼ 6* **

We present the results from a survey of i-dropout objects selected from ~1550 deg2 of multicolor imaging data from the Sloan Digital Sky Survey to search for luminous quasars at z 5.8. Objects with i*-z* > 2.2 and z* 0.90. The ARC 3.5 m spectrum of SDSSp J103027.10+052455.0 shows that over a range of ~300 A immediately blueward of the Lyα emission, the average transmitted flux is only 0.003 ± 0.020 times that of the continuum level, consistent with zero flux over a ~300 A range of the Lyα forest region and suggesting a tentative detection of the complete Gunn-Peterson trough. The existence of strong metal lines in the quasar spectra suggests early metal enrichment in the quasar environment. The three new objects, together with the previously published z = 5.8 quasar SDSSp J104433.04-012502.2, form a complete color-selected flux-limited sample at z 5.8. We estimate the selection function of this sample, taking into account the estimated variations in the quasar spectral energy distribution, as well as observational photometric errors. We find that at z = 6, the comoving density of luminous quasars at M1450 < -26.8 (H0 = 50 km s-1 Mpc-1, Ω = 1) is 1.1 × 10-9 Mpc-3. This is a factor of ~2 lower than that at z ~ 5 and is consistent with an extrapolation of the observed quasar evolution at z < 5. Using the current sample, we discuss the constraint on the shape of the quasar luminosity function and the implications for the contribution of quasars to the ionizing background at z ~ 6. The luminous quasars discussed in the paper have central black hole masses of several times 109 M⊙ by the Eddington argument, with likely dark halo masses on the order of 1013 M⊙. Their observed space density provides a sensitive test of models of quasar and galaxy formation at high redshift.

Galaxy Clustering in Early Sloan Digital Sky Survey Redshift Data

We present the first measurements of clustering in the Sloan Digital Sky Survey (SDSS) galaxy redshift survey. Our sample consists of 29,300 galaxies with redshifts 5700 km s-1 ≤ cz ≤ 39,000 km s-1, distributed in several long but narrow (25-5°) segments, covering 690 deg2. For the full, flux-limited sample, the redshift-space correlation length is approximately 8 h-1 Mpc. The two-dimensional correlation function ξ(rp,π) shows clear signatures of both the small-scale, fingers-of-God distortion caused by velocity dispersions in collapsed objects and the large-scale compression caused by coherent flows, though the latter cannot be measured with high precision in the present sample. The inferred real-space correlation function is well described by a power law, ξ(r) = (r/6.1 ± 0.2 h-1 Mpc)-1.75±0.03, for 0.1 h-1 Mpc ≤ r ≤ 16 h-1 Mpc. The galaxy pairwise velocity dispersion is σ12 ≈ 600 ± 100 km s-1 for projected separations 0.15 h-1 Mpc ≤ rp ≤ 5 h-1 Mpc. When we divide the sample by color, the red galaxies exhibit a stronger and steeper real-space correlation function and a higher pairwise velocity dispersion than do the blue galaxies. The relative behavior of subsamples defined by high/low profile concentration or high/low surface brightness is qualitatively similar to that of the red/blue subsamples. Our most striking result is a clear measurement of scale-independent luminosity bias at r 10 h-1 Mpc: subsamples with absolute magnitude ranges centered on M* - 1.5, M*, and M* + 1.5 have real-space correlation functions that are parallel power laws of slope ≈-1.8 with correlation lengths of approximately 7.4, 6.3, and 4.7 h-1 Mpc, respectively.

A Survey of z > 5.7 Quasars in the Sloan Digital Sky Survey. II. Discovery of Three Additional Quasars at z > 6*

We present the discovery of three new quasars at z > 6 in ~ 1300 deg2 of Sloan Digital Sky Survey imaging data, J114816.64+525150.3 (z = 6.43), J104845.05+463718.3 (z = 6.23), and J163033.90+401209.6 (z = 6.05). The first two objects have weak Lyα emission lines; their redshifts are determined from the positions of the Lyman break. They are only accurate to ~0.05 and could be affected by the presence of broad absorption line systems. The last object has a Lyα strength more typical of lower redshift quasars. Based on a sample of six quasars at z > 5.7 that cover 2870 deg2 presented in this paper and in Paper I, we estimate the comoving density of luminous quasars at z ~ 6 and M1450 5.7 quasars and high-resolution, ground-based images (seeing ~04) of three additional z > 5.7 quasars show that none of them is gravitationally lensed. The luminosity distribution of the high-redshift quasar sample suggests the bright-end slope of the quasar luminosity function at z ~ 6 is shallower than Ψ ∝ L-3.5 (2 σ), consistent with the absence of strongly lensed objects.

The Sloan Digital Sky Survey Quasar Survey: Quasar Luminosity Function from Data Release 3

We determine the number counts and z = 0-5 luminosity function for a well-defined, homogeneous sample of quasars from the Sloan Digital Sky Survey (SDSS). We conservatively define the most uniform statistical sample possible, consisting of 15,343 quasars within an effective area of 1622 deg2 that was derived from a parent sample of 46,420 spectroscopically confirmed broad-line quasars in the 5282 deg2 of imaging data from SDSS Data Release 3. The sample extends from i = 15 to 19.1 at z 3 and to i = 20.2 for z 3. The number counts and luminosity function agree well with the results of the Two Degree Field QSO Redshift Survey (2QZ) at redshifts and luminosities at which the SDSS and 2QZ quasar samples overlap, but the SDSS data probe to much higher redshifts than does the 2QZ sample. The number density of luminous quasars peaks between redshifts 2 and 3, although uncertainties in the selection function in this range do not allow us to determine the peak redshift more precisely. Our best-fit model has a flatter bright-end slope at high redshift than at low redshift. For z < 2.4 the data are best fit by a redshift-independent slope of ? = -3.1 [?(L) ? L?]. Above z = 2.4 the slope flattens with redshift to ? -2.37 at z = 5. This slope change, which is significant at the 5 ? level, must be accounted for in models of the evolution of accretion onto supermassive black holes.

The Milky Way Tomography with SDSS. II. Stellar Metallicity

In addition to optical photometry of unprecedented quality, the Sloan Digital Sky Survey (SDSS) is producing a massive spectroscopic database which already contains over 280,000 stellar spectra. Using eectiv e temperature and metallicity derived from SDSS spectra for 60,000 F and G type main sequence stars (0:2 < g r < 0:6), we develop polynomial models, reminiscent of traditional methods based on the UBV photometry, for estimating these parameters from the SDSS u g and g r colors. These estimators reproduce SDSS spectroscopic parameters with a root-mean-square scatter of 100 K for eectiv e temperature, and 0.2 dex for metallicity (limited by photometric errors), which are similar to random and systematic uncertainties in spectroscopic determinations. We apply this method to a photometric catalog of coadded SDSS observations and study the photometric metallicity distribution of 200,000 F and G type stars observed in 300 deg 2 of high Galactic latitude sky. These deeper (g < 20:5) and photometrically precise ( 0.01 mag) coadded data enable an accurate measurement of the unbiased metallicity distribution for a complete volume-limited sample of stars at distances between 500 pc and 8 kpc. The metallicity distribution can be exquisitely modeled using two components with a spatially varying number ratio, that correspond to disk and halo. The best-t number ratio of the two components is consistent with that implied by the decomposition of stellar counts proles into exponential disk and power-law halo components by Juri c et al. (2008). The two components also possess the kinematics expected for disk and halo stars. The metallicity of the halo component can be modeled as a spatially invariant Gaussian distribution with a mean of [F e=H] = 1:46 and a standard deviation of 0.3 dex. The disk metallicity distribution is non-Gaussian, with a remarkably small scatter (rms 0.16 dex) and the median smoothly decreasing with distance from the plane from 0:6 at 500 pc to 0:8 beyond several kpc. Similarly, we nd using proper motion measurements that a nonGaussian rotational velocity distribution of disk stars shifts by 50 km/s as the distance from the plane increases from 500 pc to several kpc. Despite this similarity, the metallicity and rotational velocity distributions of disk stars are not correlated (Kendall’s = 0:017 0:018). This absence of a correlation between metallicity and kinematics for disk stars is in a conict with the traditional decomposition in terms of thin and thick disks, which predicts a strong correlation ( = 0:30 0:04) at 1 kpc from the mid-plane. Instead, the variation of the metallicity and rotational velocity distributions can be modeled using non-Gaussian functions that retain their shapes and only shift as the distance from the mid-plane increases. We also study the metallicity distribution using a shallower (g < 19:5) but much larger sample of close to three million stars in 8500 sq. deg. of sky included in SDSS Data Release 6. The large sky coverage enables the detection of coherent substructures in the kinematics{ metallicity space, such as the Monoceros stream, which rotates faster than the LSR, and has a median metallicity of [F e=H] = 0:95, with an rms scatter of only 0.15 dex. We extrapolate our results to the performance expected from the Large Synoptic Survey Telescope (LSST) and estimate that LSST will obtain metallicity measurements accurate to 0.2 dex or better, with proper motion measurements accurate to 0.2-0.5 mas/yr, for about 200 million F/G dwarf stars within a distance limit of 100 kpc (g < 23:5). Subject headings: methods: data analysis | stars: statistics | Galaxy: halo, kinematics and dynamics, stellar content, structure

The ensemble photometric variability of ∼25,000 quasars in the Sloan Digital Sky Survey

Using a sample of over 25,000 spectroscopically confirmed quasars from the Sloan Digital Sky Survey, we show how quasar variability in the rest-frame optical/UV regime depends on rest-frame time lag, luminosity, rest wavelength, redshift, the presence of radio and X-ray emission, and the presence of broad absorption line systems. Imaging photometry is compared with three-band spectrophotometry obtained at later epochs spanning time lags up to about 2 yr. The large sample size and wide range of parameter values allow the dependence of variability to be isolated as a function of many independent parameters. The time dependence of variability (the structure function) is well fitted by a single power law with an index γ = 0.246 ± 0.008, on timescales from days to years. There is an anticorrelation of variability amplitude with rest wavelength—e.g., quasars are about twice as variable at 1000 A as at 6000 A—and quasars are systematically bluer when brighter at all redshifts. There is a strong anticorrelation of variability with quasar luminosity—variability amplitude decreases by a factor of about 4 when luminosity increases by a factor of 100. There is also a significant positive correlation of variability amplitude with redshift, indicating evolution of the quasar population or the variability mechanism. We parameterize all of these relationships. Quasars with ROSAT All-Sky Survey X-ray detections are significantly more variable (at optical/UV wavelengths) than those without, and radio-loud quasars are marginally more variable than their radio-quiet counterparts. We find no significant difference in the variability of quasars with and without broad absorption line troughs. Currently, no models of quasar variability address more than a few of these relationships. Models involving multiple discrete events or gravitational microlensing are unlikely by themselves to account for the data. So-called accretion disk instability models are promising, but more quantitative predictions are needed.