Branimir Sesar

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

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

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.

LIGHT CURVE TEMPLATES AND GALACTIC DISTRIBUTION OF RR LYRAE STARS FROM SLOAN DIGITAL SKY SURVEY STRIPE 82

We present an improved analysis of halo substructure traced by RR Lyrae stars in the Sloan Digital Sky Survey (SDSS) stripe 82 region. With the addition of SDSS-II data, a revised selection method based on new ugriz light curve templates results in a sample of 483 RR Lyrae stars that is essentially free of contamination. The main result from our first study persists: the spatial distribution of halo stars at galactocentric distances 5-100 kpc is highly inhomogeneous. At least 20% of halo stars within 30 kpc from the Galactic center can be statistically associated with substructure. We present strong direct evidence, based on both RR Lyrae stars and main-sequence stars, that the halo stellar number density profile significantly steepens beyond a Galactocentric distance of ~30 kpc, and a larger fraction of the stars are associated with substructure. By using a novel method that simultaneously combines data for RR Lyrae and main-sequence stars, and using photometric metallicity estimates for main-sequence stars derived from deep co-added u-band data, we measure the metallicity of the Sagittarius dSph tidal stream (trailing arm) toward R.A. ~2h-3h and decl. ~ 0? to be 0.3 dex higher ([Fe/H] = ?1.2) than that of surrounding halo field stars. Together with a similar result for another major halo substructure, the Monoceros stream, these results support theoretical predictions that an early forming, smooth inner halo, is metal-poor compared to high surface brightness material that have been accreted onto a later-forming outer halo. The mean metallicity of stars in the outer halo that are not associated with detectable clumps may still be more metal-poor than the bulk of inner-halo stars, as has been argued from other data sets.

A DESCRIPTION OF QUASAR VARIABILITY MEASURED USING REPEATED SDSS AND POSS IMAGING

We provide a quantitative description and statistical interpretation of the optical continuum variability of quasars. The Sloan Digital Sky Survey (SDSS) has obtained repeated imaging in five UV-to-IR photometric bands for 33,881 spectroscopically confirmed quasars. About 10,000 quasars have an average of 60 observations in each band obtained over a decade along Stripe 82 (S82), whereas the remaining ~25,000 have 2-3 observations due to scan overlaps. The observed time lags span the range from a day to almost 10 years, and constrain quasar variability at rest-frame time lags of up to 4 years, and at rest-frame wavelengths from 1000 A to 6000 A. We publicly release a user-friendly catalog of quasars from the SDSS Data Release 7 that have been observed at least twice in SDSS or once in both SDSS and the Palomar Observatory Sky Survey, and we use it to analyze the ensemble properties of quasar variability. Based on a damped random walk (DRW) model defined by a characteristic timescale and an asymptotic variability amplitude that scale with the luminosity, black hole mass, and rest wavelength for individual quasars calibrated in S82, we can fully explain the ensemble variability statistics of the non-S82 quasars such as the exponential distribution of large magnitude changes. All available data are consistent with the DRW model as a viable description of the optical continuum variability of quasars on timescales of ~5-2000 days in the rest frame. We use these models to predict the incidence of quasar contamination in transient surveys such as those from the Palomar Transient Factory and Large Synoptic Survey Telescope.

The palomar transient factory photometric calibration

The Palomar Transient Factory (PTF) provides multiple epoch imaging for a large fraction of the celestial sphere. Here, we describe the photometric calibration of the PTF data products that allows the PTF magnitudes to be related to other magnitude systems. The calibration process utilizes Sloan Digital Sky Survey (SDSS) r ∼ 16 mag point-source objects as photometric standards. During photometric conditions, this allows us to solve for the extinction coefficients and color terms and to estimate the camera illumination correction. This also enables the calibration of fields that are outside the SDSS footprint. We test the precision and repeatability of the PTF photometric calibration. Given that PTF is observing in a single filter each night, we define a PTF calibrated magnitude system for the R band and g band. We show that, in this system, ≈59% (47%) of the photometrically calibrated PTF R-band (g-band) data achieve a photometric precision of 0.02–0.04 mag and have color terms and extinction coefficients that are close to their average values. Given the objects’ color, the PTF magnitude system can be converted to other systems. Moreover, a night-by-night comparison of the calibrated magnitudes of individual stars observed on multiple nights shows that they are consistent to a level of ≈0.02 mag. Most of the data that were taken under nonphotometric conditions can be calibrated relative to other epochs of the same sky footprint obtained during photometric conditions. We provide a concise guide describing how to use the PTF photometric-calibration data products, as well as the transformations between the PTF magnitude system and the SDSS and Johnson-Cousins systems.

Sagittarius II, Draco II and Laevens 3: Three New Milky Way Satellites Discovered in the Pan-STARRS 1 3π Survey

Author(s): Laevens, BPM; Martin, NF; Bernard, EJ; Schlafly, EF; Sesar, B; Rix, HW; Bell, EF; Ferguson, AMN; Slater, CT; Sweeney, WE; Wyse, RFG; Huxor, AP; Burgett, WS; Chambers, KC; Draper, PW; Hodapp, KA; Kaiser, N; Magnier, EA; Metcalfe, N; Tonry, JL; Wainscoat, RJ; Waters, C | Abstract: