Amy E. Kimball

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

MODELING THE TIME VARIABILITY OF SDSS STRIPE 82 QUASARS AS A DAMPED RANDOM WALK

We model the time variability of ~9000 spectroscopically confirmed quasars in SDSS Stripe 82 as a damped random walk (DRW). Using 2.7 million photometric measurements collected over 10 yr, we confirm the results of Kelly et al. and Kozlowski et al. that this model can explain quasar light curves at an impressive fidelity level (0.01-0.02 mag). The DRW model provides a simple, fast (O(N) for N data points), and powerful statistical description of quasar light curves by a characteristic timescale (τ) and an asymptotic rms variability on long timescales (SF∞). We searched for correlations between these two variability parameters and physical parameters such as luminosity and black hole mass, and rest-frame wavelength. Our analysis shows SF∞ to increase with decreasing luminosity and rest-frame wavelength as observed previously, and without a correlation with redshift. We find a correlation between SF∞ and black hole mass with a power-law index of 0.18 ± 0.03, independent of the anti-correlation with luminosity. We find that τ increases with increasing wavelength with a power-law index of 0.17, remains nearly constant with redshift and luminosity, and increases with increasing black hole mass with a power-law index of 0.21 ± 0.07. The amplitude of variability is anti-correlated with the Eddington ratio, which suggests a scenario where optical fluctuations are tied to variations in the accretion rate. However, we find an additional dependence on luminosity and/or black hole mass that cannot be explained by the trend with Eddington ratio. The radio-loudest quasars have systematically larger variability amplitudes by about 30%, when corrected for the other observed trends, while the distribution of their characteristic timescale is indistinguishable from that of the full sample. We do not detect any statistically robust differences in the characteristic timescale and variability amplitude between the full sample and the small subsample of quasars detected by ROSAT. Our results provide a simple quantitative framework for generating mock quasar light curves, such as currently used in LSST image simulations.

Stellar SEDs from 0.3 to 2.5 μm: Tracing the Stellar Locus and Searching for Color Outliers in the SDSS and 2MASS

The Sloan Digital Sky Survey (SDSS) and Two Micron All Sky Survey (2MASS) are rich resources for studying stellar astrophysics and the structure and formation history of the Galaxy. As new surveys and instruments adopt similar filter sets, it is increasingly important to understand the properties of the ugrizJHKs stellar locus, both to inform studies of ‘normal’ main sequence stars as well as for robust searches for point sources with unusual colors. Using a sample of � 600,000 point sources detected by SDSS and 2MASS, we tabulate the position and width of the ugrizJHKs stellar locus as a function of g i color, and provide accurate polynomial fits. We map the Morgan-Keenan spectral type sequence to the median stellar locus by using synthetic photometry of spectral standards and by analyzing 3000 SDSS stellar spectra with a custom spectral typing pipeline, described in full in an attached Appendix. Having characterized the properties of ‘normal’ main sequence stars, we develop an algorithm for identifying point sources whose colors differ significantly from those of normal stars. This algorithm calculates a point source’s minimum separation from the stellar locus in a seven-dimensional color space, and robustly identifies objects with unusual colors, as well as spurious SDSS/2MASS matches. Analysis of a final catalog of 2117 color outliers identifies 370 white-dwarf/M dwarf (WDMD) pairs, 93 QSOs, and 90 M giant/carbon star candidates, and demonstrates that WDMD pairs and QSOs can be distinguished on the basis of their J Ks and r z colors. We also identify a group of objects with correlated offsets in the u g vs. g r and g r vs. r i color-color spaces, but subsequent follow-up is required to reveal the nature of these objects. Future applications of this algorithm to a matched SDSS-UKIDSS catalog may well identify additional classes of objects with unusual colors by probing new areas of color-magnitude space. Subject headings: surveys — stars:late-type — stars:early-type — Galaxy:stellar content — infrared:stars

The Sloan Digital Sky Survey quasar catalog: ninth data release

We present the Data Release 9 Quasar (DR9Q) catalog from the Baryon Oscillation Spectroscopic Survey (BOSS) of the Sloan Digital Sky Survey III. The catalog includes all BOSS objects that were targeted as quasar candidates during the survey, are spectrocopically confirmed as quasars via visual inspection, have luminosities Mi[z = 2] 2.15 (61 931) is ~2.8 times larger than the number of z > 2.15 quasars previously known. Redshifts and FWHMs are provided for the strongest emission lines (C iv, C iii], Mg ii). The catalog identifies 7533 broad absorption line quasars and gives their characteristics. For each object the catalog presents five-band (u, g, r, i, z) CCD-based photometry with typical accuracy of 0.03 mag, and information on the morphology and selection method. The catalog also contains X-ray, ultraviolet, near-infrared, and radio emission properties of the quasars, when available, from other large-area surveys. The calibrated digital spectra cover the wavelength region 3600−10 500 A at a spectral resolution in the range 1300 < R < 2500; the spectra can be retrieved from the SDSS Catalog Archive Server. We also provide a supplemental list of an additional 949 quasars that have been identified, among galaxy targets of the BOSS or among quasar targets after DR9 was frozen.

EMU: Evolutionary Map of the Universe

EMU is a wide-field radio continuum survey planned for the new Australian Square Kilometre Array Pathfinder (ASKAP) telescope. The primary goal of EMU is to make a deep (rms ~10 μJy/beam) radio continuum survey of the entire Southern sky at 1.3 GHz, extending as far North as +30° declination, with a resolution of 10 arcsec. EMU is expected to detect and catalogue about 70 million galaxies, including typical star-forming galaxies up to z ~ 1, powerful starbursts to even greater redshifts, and active galactic nuclei to the edge of the visible Universe. It will undoubtedly discover new classes of object. This paper defines the science goals and parameters of the survey, and describes the development of techniques necessary to maximise the science return from EMU.

Measurement of Galaxy Cluster Sizes, Radial Profiles, and Luminosity Functions from SDSS Photometric Data

Imaging data from the Sloan Digital Sky Survey are used to measure the empirical size-richness relation for a large sample of galaxy clusters. Using population subtraction methods, we determine the radius at which the cluster galaxy number density is 200Ωm-1 times the mean galaxy density, without assuming a model for the radial distribution of galaxies in clusters. If these galaxies are unbiased on megaparsec scales, this galaxy density-based R200 reflects the characteristic radii of clusters. We measure the scaling of this characteristic radius with richness over an order of magnitude in cluster richness, from rich clusters to poor groups. We use this information to examine the radial profiles of galaxies in clusters as a function of cluster richness, finding that the concentration of the galaxy distribution decreases with richness and is systematically lower than the concentrations measured for dark matter profiles in N-body simulations. Using these scaled radii, we investigate the behavior of the cluster luminosity function and find that it is well matched by a Schechter function for galaxies brighter than Mr = -18 only after the central galaxy has been removed. We find that the luminosity function varies with richness and with distance from the cluster center, underscoring the importance of using an aperture that scales with cluster mass to compare physically equivalent regions of these different systems. We note that the lowest richness systems in our catalog have properties consistent with those expected of the earliest forming halos; our cluster-finding algorithm, in addition to reliably finding clusters, may be efficient at finding fossil groups.

Optically selected BL lacertae candidates from the sloan digital sky survey data release seven

We present a sample of 723 optically selected BL Lac candidates from the Sloan Digital Sky Survey Data Release 7 (SDSS DR7) spectroscopic database encompassing 8250?deg2 of sky; our sample constitutes one of the largest uniform BL Lac samples yet derived. Each BL Lac candidate has a high-quality SDSS spectrum from which we determine spectroscopic redshifts for ~60% of the objects. Redshift lower limits are estimated for the remaining objects utilizing the lack of host galaxy flux contamination in their optical spectra; we find that objects lacking spectroscopic redshifts are likely at systematically higher redshifts. Approximately 80% of our BL Lac candidates match to a radio source in FIRST/NVSS, and ~40% match to a ROSAT X-ray source. The homogeneous multiwavelength coverage allows subdivision of the sample into 637 radio-loud BL Lac candidates and 86 weak-featured radio-quiet objects. The radio-loud objects broadly support the standard paradigm unifying BL Lac objects with beamed radio galaxies. We propose that the majority of the radio-quiet objects may be lower-redshift (z < 2.2) analogs to high-redshift weak line quasars (i.e., active galactic nucleus with unusually anemic broad emission line regions). These would constitute the largest sample of such objects, being of similar size and complementary in redshift to the samples of high-redshift weak line quasars previously discovered by the SDSS. However, some fraction of the weak-featured radio-quiet objects may instead populate a rare and extreme radio-weak tail of the much larger radio-loud BL Lac population. Serendipitous discoveries of unusual white dwarfs, high-redshift weak line quasars, and broad absorption line quasars with extreme continuum dropoffs blueward of rest-frame 2800?? are also briefly described.

THE TWO-COMPONENT RADIO LUMINOSITY FUNCTION OF QUASI-STELLAR OBJECTS: STAR FORMATION AND ACTIVE GALACTIC NUCLEUS

Despite decades of study, it remains unclear whether there are distinct radio-loud and radio-quiet populations of quasi-stellar objects (QSOs). Early studies were limited by inhomogeneous QSO samples, inadequate sensitivity to probe the radio-quiet population, and degeneracy between redshift and luminosity for flux-density-limited samples. Our new 6 GHz EVLA observations allow us for the first time to obtain nearly complete (97%) radio detections in a volume-limited color-selected sample of 179 QSOs more luminous than M_i = -23 from the Sloan Digital Sky Survey (SDSS) Data Release Seven in the narrow redshift range 0.2 23.5, and (c) the total number of SDSS QSOs in our volume-limited sample. We show that the RLF can be explained as a superposition of two populations, dominated by AGNs at the bright end and star formation in the QSO host galaxies at the faint end.

Active Galactic Nucleus and Starburst Radio Emission from Optically Selected Quasi-stellar Objects

We used the 1.4 GHz NVSS to study radio sources in two color-selected QSO samples: a volumelimited sample of 1313 QSOs defined by Mi < −23 in the redshift range 0.2 < z < 0.45 and a magnitude-limited sample of 2471 QSOs with mr ≤ 18.5 and 1.8 < z < 2.5. About 10% were detected above the 2.4 mJy NVSS catalog limit and are powered primarily by AGNs. The space density of the low-redshift QSOs evolves as ρ ∝ (1 + z). In both redshift ranges the flux-density distributions and luminosity functions of QSOs stronger than 2.4 mJy are power laws, with no features to suggest more than one kind of radio source. Extrapolating the power laws to lower luminosities predicts the remaining QSOs should be extremely radio quiet, but they are not. Most were detected statistically on the NVSS images with median peak flux densities Sp(mJy beam ) ≈ 0.3 and 0.05 in the lowand high-redshift samples, corresponding to spectral luminosities log[L1.4GHz(W Hz )] ≈ 22.7 and 24.1, respectively. We suggest that the faint radio sources are powered by star formation at rates Ṁ ∼ 20M⊙ yr −1 in the moderate luminosity (median 〈Mi〉 ≈ −23.4) low-redshift QSOs and Ṁ ∼ 500M⊙ yr −1 in the very luminous (〈Mi〉 ≈ −27.5) high-redshift QSOs. Such luminous starbursts [〈log(LIR/L⊙)〉 ∼ 11.2 and 12.6, respectively] are consistent with “quasar mode” accretion in which cold gas flows fuel both AGN and starburst. Subject headings: galaxies: active—galaxies: starburst—quasars: general—radio continuum: galaxies

Spitzer Observations of a Gravitationally Lensed Quasar, QSO 2237+0305

The four-image gravitationally lensed quasar QSO 2237+0305 is microlensed by stars in the lens galaxy. The amplitude of microlensing variability can be used to infer the relative size of the quasar as a function of wavelength; this provides a test of quasar models. Toward this end, we present Spitzer Space Telescope Infrared Spectrograph and Infrared Array Camera (IRAC) observations of QSO 2237+0305, finding the following. (1) The infrared (IR) spectral energy distribution (SED) is similar to that of other bright radio-quiet quasars, contrary to an earlier claim. (2) A dusty torus model with a small opening angle fits the overall shape of the IR SED well, but the quantitative agreement is poor due to an offset in wavelength of the silicate feature. (3) The flux ratios of the four lensed images can be derived from the IRAC data despite being unresolved. We find that the near-IR fluxes are increasingly affected by microlensing toward shorter wavelengths. (4) The wavelength dependence of the IRAC flux ratios is consistent with the standard quasar model in which an accretion disk and a dusty torus both contribute near 1 μm in the rest frame. This is also consistent with recent IR spectropolarimetry of nearby quasars.