Nicholas A. Bond

The Milky Way Tomography with SDSS. I. Stellar Number Density Distribution

Using the photometric parallax method we estimate the distances to ~48 million stars detected by the Sloan Digital Sky Survey (SDSS) and map their three-dimensional number density distribution in the Galaxy. The currently available data sample the distance range from 100 pc to 20 kpc and cover 6500 deg2 of sky, mostly at high Galactic latitudes (|b| > 25). These stellar number density maps allow an investigation of the Galactic structure with no a priori assumptions about the functional form of its components. The data show strong evidence for a Galaxy consisting of an oblate halo, a disk component, and a number of localized overdensities. The number density distribution of stars as traced by M dwarfs in the solar neighborhood (D < 2 kpc) is well fit by two exponential disks (the thin and thick disk) with scale heights and lengths, bias corrected for an assumed 35% binary fraction, of H1 = 300 pc and L1 = 2600 pc, and H2 = 900 pc and L2 = 3600 pc, and local thick-to-thin disk density normalization ρthick(R☉)/ρthin(R☉) = 12% . We use the stars near main-sequence turnoff to measure the shape of the Galactic halo. We find a strong preference for oblate halo models, with best-fit axis ratio c/a = 0.64, ρH ∝ r−2.8 power-law profile, and the local halo-to-thin disk normalization of 0.5%. Based on a series of Monte Carlo simulations, we estimate the errors of derived model parameters not to be larger than ~20% for the disk scales and ~10% for the density normalization, with largest contributions to error coming from the uncertainty in calibration of the photometric parallax relation and poorly constrained binary fraction. While generally consistent with the above model, the measured density distribution shows a number of statistically significant localized deviations. In addition to known features, such as the Monoceros stream, we detect two overdensities in the thick disk region at cylindrical galactocentric radii and heights (R,Z) ~ (6.5,1.5) kpc and (R,Z) ~ (9.5,0.8) kpc and a remarkable density enhancement in the halo covering over 1000 deg2 of sky toward the constellation of Virgo, at distances of ~6-20 kpc. Compared to counts in a region symmetric with respect to the l = 0° line and with the same Galactic latitude, the Virgo overdensity is responsible for a factor of 2 number density excess and may be a nearby tidal stream or a low-surface brightness dwarf galaxy merging with the Milky Way. The u − g color distribution of stars associated with it implies metallicity lower than that of thick disk stars and consistent with the halo metallicity distribution. After removal of the resolved overdensities, the remaining data are consistent with a smooth density distribution; we detect no evidence of further unresolved clumpy substructure at scales ranging from ~50 pc in the disk to ~1-2 kpc in the halo.

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

Sloan Digital Sky Survey Standard Star Catalog for Stripe 82: The Dawn of Industrial 1% Optical Photometry

We describe a standard star catalog constructed using multiple SDSS photometric observations (at least four per band, with a median of 10) in the ugriz system. The catalog includes 1.01 million nonvariable unresolved objects from the equatorial stripe 82 (|δJ2000.0| < 1.266°) in the right ascension range 20h34m-4h00m and with the corresponding r-band (approximately Johnson V-band) magnitudes in the range 14-22. The distributions of measurements for individual sources demonstrate that the photometric pipeline correctly estimates random photometric errors, which are below 0.01 mag for stars brighter than 19.5, 20.5, 20.5, 20, and 18.5 in ugriz, respectively (about twice as good as for individual SDSS runs). Several independent tests of the internal consistency suggest that the spatial variation of photometric zero points is not larger than ~0.01 mag (rms). In addition to being the largest available data set with optical photometry internally consistent at the ~1% level, this catalog provides a practical definition of the SDSS photometric system. Using this catalog, we show that photometric zero points for SDSS observing runs can be calibrated within a nominal uncertainty of 2% even for data obtained through 1 mag thick clouds, and we demonstrate the existence of He and H white dwarf sequences using photometric data alone. Based on the properties of this catalog, we conclude that upcoming large-scale optical surveys such as the Large Synoptic Survey Telescope will be capable of delivering robust 1% photometry for billions of sources.

Exploring the Variable Sky with the Sloan Digital Sky Survey

We quantify the variability of faint unresolved optical sources using a catalog based on multiple SDSS imaging observations. The catalog covers SDSS stripe 82, which lies along the celestial equator in the southern Galactic hemisphere (22h24m < ?J2000.0 < 04h08m, -1.27? < ?J2000.0 < +1.27?, ~290 deg2), and contains 34 million photometric observations in the SDSS ugriz system for 748,084 unresolved sources at high Galactic latitudes (b < -20?) that were observed at least four times in each of the ugri bands (with a median of 10 observations obtained over ~6 yr). In each photometric bandpass we compute various low-order light-curve statistics, such as rms scatter, ?2 per degree of freedom, skewness, and minimum and maximum magnitude, and use them to select and study variable sources. We find that 2% of unresolved optical sources brighter than g = 20.5 appear variable at the 0.05 mag level (rms) simultaneously in the g and r bands (at high Galactic latitudes). The majority (2 out of 3) of these variable sources are low-redshift (<2) quasars, although they represent only 2% of all sources in the adopted flux-limited sample. We find that at least 90% of quasars are variable at the 0.03 mag level (rms) and confirm that variability is as good a method for finding low-redshift quasars as the UV excess color selection (at high Galactic latitudes). We analyze the distribution of light-curve skewness for quasars and find that it is centered on zero. We find that about one-fourth of the variable stars are RR Lyrae stars, and that only 0.5% of stars from the main stellar locus are variable at the 0.05 mag level. The distribution of light-curve skewness in the g - r versus u - g color-color diagram on the main stellar locus is found to be bimodal (with one mode consistent with Algol-like behavior). Using over 600 RR Lyrae stars, we demonstrate rich halo substructure out to distances of 100 kpc. We extrapolate these results to the expected performance by the Large Synoptic Survey Telescope and estimate that it will obtain well-sampled, 2% accurate, multicolor light curves for ~2 million low-redshift quasars and discover at least 50 million variable stars.

High-Redshift Superwinds as the Source of the Strongest Mg II Absorbers: A Feasibility Analysis

We present High Resolution Echelle Spectrometer/Keck profiles of four extremely strong (Wr > 1.8 A) Mg II absorbers at 1 1.8 A evolve away from z = 2 to the present. We propose that a substantial fraction of these very strong absorbers are due to superwinds and that their evolution is related to the redshift evolution of star-forming galaxies. Based on the observed redshift number density of Wr > 1.8 A Mg II absorbers at 1 < z < 2, we explore whether it is realistic that superwinds from starbursting galaxies could give rise to these absorbers. Finally, we do an analysis of the superwind connection to damped Lyα absorbers (DLAs). DLAs and superwinds evolve differently and usually have different kinematic structure, indicating that superwinds probably do not give rise to the majority of DLAs.

THE EVOLUTION OF Lyα-EMITTING GALAXIES BETWEEN z = 2.1 AND z = 3.1

We describe the results of a new, wide-field survey for z=3.1 Ly-alpha emission-line galaxies (LAEs) in the Extended Chandra Deep Field South (ECDF-S). By using a nearly top-hat 5010 Angstrom filter and complementary broadband photometry from the MUSYC survey, we identify a complete sample of 141 objects with monochromatic fluxes brighter than 2.4E-17 ergs/cm^2/s and observers-frame equivalent widths greater than ~ 80 Angstroms (i.e., 20 Angstroms in the rest-frame of Ly-alpha). The bright-end of this dataset is dominated by x-ray sources and foreground objects with GALEX detections, but when these interlopers are removed, we are still left with a sample of 130 LAE candidates, 39 of which have spectroscopic confirmations. This sample overlaps the set of objects found in an earlier ECDF-S survey, but due to our filters redder bandpass, it also includes 68 previously uncataloged sources. We confirm earlier measurements of the z=3.1 LAE emission-line luminosity function, and show that an apparent anti-correlation between equivalent width and continuum brightness is likely due to the effect of correlated errors in our heteroskedastic dataset. Finally, we compare the properties of z=3.1 LAEs to LAEs found at z=2.1. We show that in the ~1 Gyr after z~3, the LAE luminosity function evolved significantly, with L* fading by ~0.4 mag, the number density of sources with L > 1.5E42 ergs/s declining by ~50%, and the equivalent width scale-length contracting from 70^{+7}_{-5} Angstroms to 50^{+9}_{-6} Angstroms. When combined with literature results, our observations demonstrate that over the redshift range z~0 to z~4, LAEs contain less than ~10% of the star-formation rate density of the universe.

On the Origin of Intrinsic Narrow Absorption Lines in z ≲ 1 QSOs

We present an exhaustive statistical analysis of the associated (Δvabs -0.5], and mediocre C IV FWHM (6000 km s-1) do not have detectable associated NALs, down to Wr(C ) = 0.35 A. We also find that broad absorption line (BAL) QSOs have an enhanced probability of hosting detectable NAL gas. In addition, we find that the velocity distribution of associated NALs is peaked around the emission redshifts rather than the systemic redshifts of the QSOs. Finally, we find only one strong NAL [Wr(C ) 1.5 A] in our low-redshift sample. A comparison with previous higher redshift surveys reveals evolution in the number of strong NAL systems with redshift. We interpret these results in the context of an accretion disk model. We propose that NAL gas hugs the streamlines of the faster, denser, low-latitude wind, which has been associated with BALs. In the framework of this scenario, we can explain the observational clues as resulting from differences in orientation and wind properties, the latter presumably associated with the QSO radio properties.

Ly alpha-EMITTING GALAXIES AT z=2.1 : STELLAR MASSES, DUST, AND STAR FORMATION HISTORIES FROM SPECTRAL ENERGY DISTRIBUTION FITTING

We study the physical properties of 216 z similar or equal to 2.1 Ly alpha-emitting galaxies (LAEs) discovered in an ultra-deep narrow-MUSYC image of the ECDF-S. We fit their stacked spectral energ ...

Crawling the cosmic network: identifying and quantifying filamentary structure

We present the Smoothed Hessian Major Axis Filament Finder (SHMAFF), an algorithm that uses the eigenvectors of the Hessian matrix of the smoothed galaxy distribution to identify individual filamentary structures. Filaments are traced along the Hessian eigenvector corresponding to the largest eigenvalue and are stopped when the axis orientation changes more rapidly than a preset threshold. In both N-body simulations and the Sloan Digital Sky Survey (SDSS) main galaxy redshift survey data, the resulting filament length distributions are approximately exponential. In the SDSS galaxy distribution, using smoothing lengths of 10 and 15 h -1 Mpc, we find filament lengths per unit volume of 1.9 x 10 -3 and 7.6 × 10- 4 h 2 Mpc -2 , respectively. The filament width distributions, which are much more sensitive to non-linear growth, are also consistent between the real and mock galaxy distributions using a standard cosmology. In SDSS, we find mean filament widths of 5.5 and 8.4 h -1 Mpc on 10 and 15 h -1 Mpc smoothing scales, with standard deviations of 1.1 and 1.4 h -1 Mpc, respectively. Finally, the spatial distribution of filamentary structure in simulations is very similar between z = 3 and z = 0 on smoothing scales as large as 15 h -1 Mpc, suggesting that the outline of filamentary structure is already in place at high redshift.

UVUDF: ULTRAVIOLET THROUGH NEAR-INFRARED CATALOG AND PHOTOMETRIC REDSHIFTS OF GALAXIES IN THE HUBBLE ULTRA DEEP FIELD

We present photometry and derived redshifts from up to eleven bandpasses for 9927 galaxies in the Hubble Ultra Deep field (UDF), covering an observed wavelength range from the near-ultraviolet (NUV) to the near-infrared (NIR) with Hubble Space Telescope observations. Our Wide Field Camera 3 (WFC3)/UV F225W, F275W, and F336W image mosaics from the ultra-violet UDF (UVUDF) imaging campaign are newly calibrated to correct for charge transfer inefficiency, and use new dark calibrations to minimize background gradients and pattern noise. Our NIR WFC3/IR image mosaics combine the imaging from the UDF09 and UDF12 campaigns with CANDELS data to provide NIR coverage for the entire UDF field of view. We use aperture-matched point-spread function corrected photometry to measure photometric redshifts in the UDF, sampling both the Lyman break and Balmer break of galaxies at z~0.8-3.4, and one of the breaks over the rest of the redshift range. Our comparison of these results with a compilation of robust spectroscopic redshifts shows an improvement in the galaxy photometric redshifts by a factor of two in scatter and a factor three in outlier fraction over previous UDF catalogs. The inclusion of the new NUV data is responsible for a factor of two decrease in the outlier fraction compared to redshifts determined from only the optical and NIR data, and improves the scatter at z 2. The panchromatic coverage of the UDF from the NUV through the NIR yields robust photometric redshifts of the UDF, with the lowest outlier fraction available.