We present visible and near-infrared photometric and spectroscopic observations of interstellar object 2I/Borisov taken from 2019 September 10 to 2019 November 29 using the GROWTH, the APO ARC 3.5 m and the NASA/IRTF 3.0 m combined with post and pre-discovery observations of 2I obtained by ZTF from 2019 March 17 to 2019 May 5. Comparison with imaging of distant Solar System comets shows an object very similar to mildly active Solar System comets with an out-gassing rate of $\sim$10$^{27}$ mol/sec. The photometry, taken in filters spanning the visible and NIR range shows a gradual brightening trend of $\sim0.03$ mags/day since 2019 September 10 UTC for a reddish object becoming neutral in the NIR. The lightcurve from recent and pre-discovery data reveals a brightness trend suggesting the recent onset of significant H$_2$O sublimation with the comet being active with super volatiles such as CO at heliocentric distances $>$6 au consistent with its extended morphology. Using the advanced capability to significantly reduce the scattered light from the coma enabled by high-resolution NIR images from Keck adaptive optics taken on 2019 October 04, we estimate a diameter of 2I's nucleus of $\lesssim$1.4 km. We use the size estimates of 1I/'Oumuamua and 2I/Borisov to roughly estimate the slope of the ISO size-distribution resulting in a slope of $\sim$3.4$\pm$1.2, similar to Solar System comets and bodies produced from collisional equilibrium.

The identification of the electromagnetic counterpart of the gravitational wave event, GW170817, and discovery of neutrinos and gamma-rays from TXS 0506+056 heralded the new era of multi-messenger astrophysics. As the number of multi-messenger events rapidly grow over the next decade, the cyberinfrastructure requirements to handle the increase in data rates, data volume, need for event follow up, and analysis across the different messengers will also explosively grow. The cyberinfrastructure requirements to enhance multi-messenger astrophysics will both be a major challenge and opportunity for astronomers, physicists, computer scientists and cyberinfrastructure specialists. Here we outline some of these requirements and argue for a distributed cyberinfrastructure institute for multi-messenger astrophysics to meet these challenges.

We have obtained repeated images of 6 fields towards the Galactic bulge in 5 passbands (u, g, r, i, z) with the DECam imager on the Blanco 4m telescope at CTIO. From over 1.6 billion individual photometric measurements in the field centered on Baade's window, we have detected 4877 putative variable stars. 474 of these have been confirmed as fundamental mode RR Lyrae stars, whose colors at minimum light yield line-of-sight reddening determinations as well as a reddening law towards the Galactic Bulge which differs significantly from the standard R_V = 3.1 formulation. Assuming that the stellar mix is invariant over the 3 square-degree field, we are able to derive a line-of-sight reddening map with sub-arcminute resolution, enabling us to obtain de-reddened and extinction corrected color-magnitude diagrams (CMD's) of this bulge field using up to 2.5 million well-measured stars. The corrected CMD's show unprecedented detail and expose sparsely populated sequences: e.g., delineation of the very wide red giant branch, structure within the red giant clump, the full extent of the horizontal branch, and a surprising bright feature which is likely due to stars with ages younger than 1 Gyr. We use the RR Lyrae stars to trace the spatial structure of the ancient stars, and find an exponential decline in density with Galactocentric distance. We discuss ways in which our data products can be used to explore the age and metallicity properties of the bulge, and how our larger list of all variables is useful for learning to interpret future LSST alerts.

The Zwicky Transient Facility (ZTF) survey generates real-time alerts for optical transients, variables, and moving objects discovered in its wide-field survey. We describe the ZTF alert stream distribution and processing (filtering) system. The system uses existing open-source technologies developed in industry: Kafka, a real-time streaming platform, and Avro, a binary serialization format. The technologies used in this system provide a number of advantages for the ZTF use case, including (1) built-in replication, scalability, and stream rewind for the distribution mechanism; (2) structured messages with strictly enforced schemas and dynamic typing for fast parsing; and (3) a Python-based stream processing interface that is similar to batch for a familiar and user-friendly plug-in filter system, all in a modular, primarily containerized system. The production deployment has successfully supported streaming up to 1.2 million alerts or roughly 70 GB of data per night, with each alert available to a consumer within about 10 s of alert candidate production. Data transfer rates of about 80,000 alerts/minute have been observed. In this paper, we discuss this alert distribution and processing system, the design motivations for the technology choices for the framework, performance in production, and how this system may be generally suitable for other alert stream use cases, including the upcoming Large Synoptic Survey Telescope.

The Zwicky Transient Facility (ZTF), a public-private enterprise, is a new time domain survey employing a dedicated camera on the Palomar 48-inch Schmidt telescope with a 47 deg$^2$ field of view and 8 second readout time. It is well positioned in the development of time domain astronomy, offering operations at 10% of the scale and style of the Large Synoptic Survey Telescope (LSST) with a single 1-m class survey telescope. The public surveys will cover the observable northern sky every three nights in g and r filters and the visible Galactic plane every night in g and r. Alerts generated by these surveys are sent in real time to brokers. A consortium of universities which provided funding ("partnership") are undertaking several boutique surveys. The combination of these surveys producing one million alerts per night allows for exploration of transient and variable astrophysical phenomena brighter than r $\sim$ 20.5 on timescales of minutes to years. We describe the primary science objectives driving ZTF including the physics of supernovae and relativistic explosions, multi-messenger astrophysics, supernova cosmology, active galactic nuclei and tidal disruption events, stellar variability, and Solar System objects.

The Zwicky Transient Facility (ZTF) is a new optical time-domain survey that uses the Palomar 48-inch Schmidt telescope. A custom-built wide-field camera provides a 47 deg$^2$ field of view and 8 second readout time, yielding more than an order of magnitude improvement in survey speed relative to its predecessor survey, the Palomar Transient Factory (PTF). We describe the design and implementation of the camera and observing system. The ZTF data system at the Infrared Processing and Analysis Center provides near-real-time reduction to identify moving and varying objects. We outline the analysis pipelines, data products, and associated archive. Finally, we present on-sky performance analysis and first scientific results from commissioning and the early survey. ZTF's public alert stream will serve as a useful precursor for that of the Large Synoptic Survey Telescope.

The ESA Euclid mission is a space telescope that will survey ~15,000 square degrees of the sky, primarily to study the distant universe (constraining cosmological parameters through the lensing of galaxies). It is also expected to observe ~150,000 Solar System Objects (SSOs), primarily in poorly understood high inclination populations, as it will mostly avoid +/-15 degrees from the ecliptic plane. With a launch date of 2022 and a 6 year survey, Euclid and LSST will operate at the same time, and have complementary capabilities. We propose a LSST mini-survey to coordinate quasi-simultaneous observations between these two powerful observatories, when possible, with the primary aim of greatly improving the orbits of SSOs discovered by these facilities. As Euclid will operate from a halo orbit around the Sun-Earth L2 Lagrangian point, there will be significant parallax between observations from Earth and Euclid (0.01 AU). This means that simultaneous observations will give an independent distance measurement to SSOs, giving additional constraints on orbits compared to single Euclid visits.

Making an inventory of the Solar System is one of the four fundamental science requirements for the Large Synoptic Survey Telescope (LSST). The current baseline footprint for LSST's main Wide-Fast-Deep (WFD) Survey observes the sky below 0$^\circ$ declination, which includes only half of the ecliptic plane. Critically, key Solar System populations are asymmetrically distributed on the sky: they will be entirely missed, or only partially mapped, if only the WFD occurs. We propose a Northern Ecliptic Spur (NES) mini survey, observing the northern sky up to +10$^\circ$ ecliptic latitude, to maximize Solar System science with LSST. The mini survey comprises a total area of $\sim$5800 deg$^2$/604 fields, with 255 observations/field over the decade, split between g,r, and z bands. Our proposed survey will 1) obtain a census of main-belt comets; 2) probe Neptune's past migration history, by exploring the resonant structure of the Kuiper belt and the Neptune Trojan population; 3) explore the origin of Inner Oort cloud objects and place significant constraints on the existence of a hypothesized planet beyond Neptune; and 4) enable precise predictions of KBO stellar occultations. These high-ranked science goals of the Solar System Science Collaboration are only achievable with this proposed northern survey.

We propose a LSST Solar System near-Sun Survey, to be implemented during twilight hours, that extends the seasonal reach of LSST to its maximum as fresh sky is uncovered at about 50 square degrees per night (1500 sq. deg. per lunation) in the morning eastern sky, and surveyable sky is lost at the same rate to the western evening sky due to the Earth's synodic motion. By establishing near-horizon fence post picket lines to the far west and far east we address Solar System science use cases (including Near Earth Objects, Interior Earth Objects, Potentially Hazardous Asteroids, Earth Trojans, near-Sun asteroids, sun-grazing comets, and dormant comets) as well as provide the first look and last look that LSST will have at the transient and variable objects within each survey field. This proposed near-Sun Survey will also maximize the overlap with the field of regard of the proposed NEOCam spacecraft that will be stationed at the Earth's L1 Lagrange point and survey near quadrature with the Sun. This will allow LSST to incidently follow-up NEOCam targets and vice-versa (as well as targets from missions such as Euclid), and will roughly correspond to the Earth's L4 and L5 regions.

We report on $g$, $r$ and $i$ band observations of the Interstellar Object 'Oumuamua (1I) taken on 2017 October 29 from 04:28 to 08:40 UTC by the Apache Point Observatory (APO) 3.5m telescope's ARCTIC camera. We find that 1I's colors are $g-r=0.41\pm0.24$ and $r-i=0.23\pm0.25$, consistent with the visible spectra of Masiero (2017), Ye et al. (2017) and Fitzsimmons et al. (2017), and most comparable to the population of Solar System C/D asteroids, Trojans, or comets. We find no evidence of any cometary activity at a heliocentric distance of 1.46 au, approximately 1.5 months after 1I's closest approach distance to the Sun. Significant brightness variability was seen in the $r$ observations, with the object becoming notably brighter towards the end of the run. By combining our APO photometric time series data with the Discovery Channel Telescope (DCT) data of Knight et al. (2017), taken 20 h later on 2017 October 30, we construct an almost complete light curve with a most probable lightcurve period of $P \simeq 4~{\rm h}$. Our results imply a double peaked rotation period of 8.1 $\pm$ 0.02 h, with a peak-to-peak amplitude of 1.5 - 2.1 mags. Assuming that 1I's shape can be approximated by an ellipsoid, the amplitude constraint implies that 1I has an axial ratio of 3.5 to 10.3, which is strikingly elongated. Assuming that 1I is rotating above its critical break up limit, our results are compatible with 1I having having modest cohesive strength and may have obtained its elongated shape during a tidal disruption event before being ejected from its home system. Astrometry useful for constraining 1I's orbit was also obtained and published in Weaver et al. (2017).

Few observational constraints exist for the tidal synchronization rate of late-type stars, despite its fundamental role in binary evolution. We visually inspected the light curves of 2278 eclipsing binaries (EBs) from the Kepler Eclipsing Binary Catalog to identify those with starspot modulations, as well as other types of out-of-eclipse variability. We report rotation periods for 816 EBs with starspot modulations, and find that 79% of EBs with orbital periods less than ten days are synchronized. However, a population of short period EBs exists with rotation periods typically 13% slower than synchronous, which we attribute to the differential rotation of high latitude starspots. At 10 days, there is a transition from predominantly circular, synchronized EBs to predominantly eccentric, pseudosynchronized EBs. This transition period is in good agreement with the predicted and observed circularization period for Milky Way field binaries. At orbital periods greater than about 30 days, the amount of tidal synchronization decreases. We also report 12 previously unidentified candidate $\delta$ Scuti and $\gamma$ Doradus pulsators, as well as a candidate RS CVn system with an evolved primary that exhibits starspot occultations. For short period contact binaries, we observe a period-color relation, and compare it to previous studies. As a whole, these results represent the largest homogeneous study of tidal synchronization of late-type stars.

The recently discovered fast radio bursts (FRBs), presumably of extra-galactic origin, have the potential to become a powerful probe of the intergalactic medium (IGM). We point out a few such potential applications. We provide expressions for the dispersion measure and rotation measure as a function of redshift, and we discuss the sensitivity of these measures to the HeII reionization and the IGM magnetic field. Finally we calculate the microlensing effect from an isolate, extragalctic stellar-mass compact object on the FRB spectrum. The time delays between the two lensing images will induce constructive and destructive interference, leaving a specific imprint on the spectra of FRBs. With a high all-sky rate, a large statistical sample of FRBs is expected to make these applications feasible.

We present a recalibration of the Sloan Digital Sky Survey (SDSS) photometry with new flat fields and zero points derived from Pan-STARRS1 (PS1). Using PSF photometry of 60 million stars with $16 < r < 20$, we derive a model of amplifier gain and flat-field corrections with per-run RMS residuals of 3 millimagnitudes (mmag) in $griz$ bands and 15 mmag in $u$ band. The new photometric zero points are adjusted to leave the median in the Galactic North unchanged for compatibility with previous SDSS work. We also identify transient non-photometric periods in SDSS ("contrails") based on photometric deviations co-temporal in SDSS bands. The recalibrated stellar PSF photometry of SDSS and PS1 has an RMS difference of {9,7,7,8} mmag in $griz$, respectively, when averaged over $15'$ regions.

The Large Synoptic Survey Telescope (LSST) will be a ground-based, optical, all-sky, rapid cadence survey project with tremendous potential for discovering and characterizing asteroids. With LSST's large 6.5m diameter primary mirror, a wide 9.6 square degree field of view 3.2 Gigapixel camera, and rapid observational cadence, LSST will discover more than 5 million asteroids over its ten year survey lifetime. With a single visit limiting magnitude of 24.5 in r-band, LSST will be able to detect asteroids in the Main Belt down to sub-kilometer sizes. The current strawman for the LSST survey strategy is to obtain two visits (each visit being a pair of back-to-back 15s exposures) per field, separated by about 30 minutes, covering the entire visible sky every 3-4 days throughout the observing season, for ten years. The catalogs generated by LSST will increase the known number of small bodies in the Solar System by a factor of 10-100 times, among all populations. The median number of observations for Main Belt asteroids will be on the order of 200-300, with Near Earth Objects receiving a median of 90 observations. These observations will be spread among ugrizy bandpasses, providing photometric colors and allowing sparse lightcurve inversion to determine rotation periods, spin axes, and shape information. These catalogs will be created using automated detection software, the LSST Moving Object Processing System (MOPS), that will take advantage of the carefully characterized LSST optical system, cosmetically clean camera, and recent improvements in difference imaging. Tests with the prototype MOPS software indicate that linking detections (and thus discovery) will be possible at LSST depths with our working model for the survey strategy, but evaluation of MOPS and improvements in the survey strategy will continue. All data products and software created by LSST will be publicly available.

We present robust constraints from the Sloan Digital Sky Survey (SDSS) on the shape and distribution of the dark matter halo within the Milky Way (MW). Using the number density distribution and kinematics of SDSS halo stars, we probe the dark matter distribution to heliocentric distances exceeding 10 kpc and galactocentric distances exceeding 20 kpc. Our analysis utilizes Jeans equations to generate two-dimensional acceleration maps throughout the volume; this approach is thoroughly tested on a cosmologically derived N-body+SPH simulation of a MW-like galaxy. We show that the known accelerations (gradients of the gravitational potential) can be successfully recovered in such a realistic system. Leveraging the baryonic gravitational potential derived by Bovy & Rix (2013), we show that the gravitational potential implied by the SDSS observations cannot be explained, assuming Newtonian gravity, by visible matter alone: the gravitational force experienced by stars at galactocentric distances of 20 kpc is as much as three times stronger than what can be attributed to purely visible matter. We also show that the SDSS data provide a strong constraint on the shape of the dark matter halo potential. Within galactocentric distances of 20 kpc, the dark matter halo potential is well described as an oblate halo with axis ratio qDM=0.7+/-0.1; this corresponds to an axis ratio qDM=0.4+/-0.1 for the dark matter density distribution. Because of our precise two-dimensional measurements of the acceleration of the halo stars, we can reject several MOND models as an explanation of the observed behavior.

Millisecond duration bright radio pulses at 1.4-GHz with high dispersion measures (DM) were reported by Lorimer et al., Keane et al., and Thornton et al. Their all-sky rate is $\approx 10^4$/day above $\sim$1 Jy. Related events are "Perytons" -- similar pulsed, dispersed sources, but most certainly local. Suggested models of fast radio bursts (FRBs) can originate in the Earth's atmosphere, in stellar coronae, in other galaxies, and even at cosmological distances. Using physically motivated assumptions combined with observed properties, we explore these models. In our analysis, we focus on the Lorimer event: a 30 Jy, 5-ms duration burst with DM$=$ 375 cm$^{-3}$ pc, exhibiting a steep frequency-dependent pulse width (the {\it Sparker}). To be complete, we drop the assumption that high DMs are produced by plasma propagation and assume that the source produces pulses with frequency-dependent arrival time ("chirped signals"). Within this framework we explore a scenario in which Perytons, the {\it Sparker}, and the FRBs are all atmospheric phenomenon occurring at different heights. This model is {\it ad hoc} in that we cannot explain why Perytons at higher altitudes show greater DMs or exhibit narrower pulses. Nonetheless, we argue the {\it Sparker} may be a Peryton. We end with two remarks. First, the detection of a single FRB by an interferometer with a kilometer (or longer) baseline will prove that FRBs are of extra-terrestrial origin. Second, we urge astronomers to pursue observations and understanding of Perytons since they form (at least) a formidable foreground for the FRBs.

We present a map of the dust reddening to 4.5 kpc derived from Pan-STARRS1 stellar photometry. The map covers almost the entire sky north of declination -30 degrees at a resolution of 7' to 14', and is based on the estimated distances and reddenings to more than 500 million stars. The technique is designed to map dust in the Galactic plane, where many other techniques are stymied by the presence of multiple dust clouds at different distances along each line of sight. This reddening-based dust map agrees closely with the Schlegel, Finkbeiner, and Davis (SFD; 1998) far-infrared emission-based dust map away from the Galactic plane, and the most prominent differences between the two maps stem from known limitations of SFD in the plane. We also compare the map with Planck, finding likewise good agreement in general at high latitudes. The use of optical data from Pan-STARRS1 yields reddening uncertainty as low as 25 mmag E(B-V).

A community meeting on the topic of "Radio Astronomy in the LSST Era" was hosted by the National Radio Astronomy Observatory in Charlottesville, VA (2013 May 6--8). The focus of the workshop was on time domain radio astronomy and sky surveys. For the time domain, the extent to which radio and visible wavelength observations are required to understand several classes of transients was stressed, but there are also classes of radio transients for which no visible wavelength counterpart is yet known, providing an opportunity for discovery. From the LSST perspective, the LSST is expected to generate as many as 1 million alerts nightly, which will require even more selective specification and identification of the classes and characteristics of transients that can warrant follow up, at radio or any wavelength. The LSST will also conduct a deep survey of the sky, producing a catalog expected to contain over 38 billion objects in it. Deep radio wavelength sky surveys will also be conducted on a comparable time scale, and radio and visible wavelength observations are part of the multi-wavelength approach needed to classify and understand these objects. Radio wavelengths are valuable because they are unaffected by dust obscuration and, for galaxies, contain contributions both from star formation and from active galactic nuclei. The workshop touched on several other topics, on which there was consensus including the placement of other LSST "Deep Drilling Fields," inter-operability of software tools, and the challenge of filtering and exploiting the LSST data stream. There were also topics for which there was insufficient time for full discussion or for which no consensus was reached, which included the procedures for following up on LSST observations and the nature for future support of researchers desiring to use LSST data products.

Studies of stellar populations, understood to mean collections of stars with common spatial, kinematic, chemical, and/or age distributions, have been reinvigorated during the last decade by the advent of large-area sky surveys such as SDSS, 2MASS, RAVE, and others. We review recent analyses of these data that, together with theoretical and modeling advances, are revolutionizing our understanding of the nature of the Milky Way, and galaxy formation and evolution in general. The formation of galaxies like the Milky Way was long thought to be a steady process leading to a smooth distribution of stars. However, the abundance of substructure in the multi-dimensional space of various observables, such as position, kinematics, and metallicity, is by now proven beyond doubt, and demonstrates the importance of mergers in the growth of galaxies. Unlike smooth models that involve simple components, the new data reviewed here clearly show many irregular structures, such as the Sagittarius dwarf tidal stream and the Virgo and Pisces overdensities in the halo, and the Monoceros stream closer to the Galactic plane. These recent developments have made it clear that the Milky Way is a complex and dynamical structure, one that is still being shaped by the merging of neighboring smaller galaxies. We also briefly discuss the next generation of wide-field sky surveys, such as SkyMapper, Pan-STARRS, Gaia and LSST, which will improve measurement precision manyfold, and comprise billions of individual stars. The ultimate goal, development of a coherent and detailed story of the assembly and evolutionary history of the Milky Way and other large spirals like it, now appears well within reach.

The Hyades cluster is an ideal target to study the dynamical evolution of a star cluster over the entire mass range due to its intermediate age and proximity to the Sun. We wanted to extend the Hyades mass function towards lower masses down to 0.1 Msol and to use the full three-dimensional spatial information to characterize the dynamical evolution of the cluster. We performed a kinematic and photometric selection using the PPMXL and Pan-STARRS1 sky surveys, to search for cluster members up to 30 pc from the cluster centre. We determined our detection efficiency and field star contamination rate to derive the cluster luminosity and mass functions down to masses of 0.1 Msol. The thorough astrometric and photometric constraints minimized the contamination. A minimum spanning tree algorithm was used to quantify the mass segregation. We discovered 43 new Hyades member candidates with velocity perpendicular to the Hyades motion up to 2 km/s. They have mass estimates between 0.43 and 0.09 Msol, for a total mass of 10 Msol. This doubles the number of Hyades candidates with masses smaller than 0.15 Msol. We provide an additional list of 11 possible candidates with velocity perpendicular to the Hyades motion up to 4 km/s. The cluster is significantly mass segregated. The extension of the mass function towards lower masses provided an even clearer signature than estimated in the past. We also identified as likely Hyades member an L0 dwarf previously assumed to be a field dwarf. Finally we question the membership of a number of previously published candidates, including a L2.5-type dwarf.

As of 2012 Jan 21, the Pan-STARRS 1 3\pi Survey has observed the 3/4 of the sky visible from Hawaii with a minimum of 2 and mean of 7.6 observations in 5 filters, g_p1, r_p1, i_p1, z_p1, y_p1. Now at the end of the second year of the mission, we are in a position to make an initial public release of a portion of this unprecedented dataset. This article describes the PS1 Photometric Ladder, Release 12.01 This is the first of a series of data releases to be generated as the survey coverage increases and the data analysis improves. The Photometric Ladder has rungs every hour in RA and at 4 intervals in declination. We will release updates with increased area coverage (more rungs) from the latest dataset until the PS1 survey and the final re-reduction are completed. The currently released catalog presents photometry of \approx 1000 objects per square degree in the rungs of the ladder. Saturation occurs at g_p1, r_p1, i_p1 \approx 13.5; z_p1 \approx 13.0; and y_p1 \approx 12.0. Photometry is provided for stars down to g_p1, r_p1, i_p1 \approx 19.1 in the AB system. This data release depends on the rigid `Ubercal' photometric calibration using only the photometric nights, with systematic uncertainties of (8.0, 7.0, 9.0, 10.7, 12.4) millimags in (g_p1, r_p1, i_p1, z_p1, y_p1). Areas covered only with lower quality nights are also included, and have been tied to the Ubercal solution via relative photometry; photometric accuracy of the non-photometric regions is lower and should be used with caution.

We use data from the Pan-STARRS1 survey to present a panoramic view of the Sagittarius tidal stream in the southern Galactic hemisphere. As a result of the extensive sky coverage of Pan-STARRS1, the southern stream is visible along more than 60 degrees of its orbit, nearly double the length seen by the SDSS. The recently discovered southern bifurcation of the stream is also apparent, with the fainter branch of the stream visible over at least 30 degrees. Using a combination of fitting both the main sequence turn-off and the red clump, we measure the distance to both arms of the stream in the south. We find that the distances to the bright arm of the stream agree very well with the N-body models of Law & Majewski (2010). We also find that the faint arm lies ~5 kpc closer to the Sun than the bright arm, similar to the behavior seen in the northern hemisphere.

We present Swarm-NG, a C++ library for the efficient direct integration of many n-body systems using highly-parallel Graphics Processing Unit (GPU), such as NVIDIA's Tesla T10 and M2070 GPUs. While previous studies have demonstrated the benefit of GPUs for n-body simulations with thousands to millions of bodies, Swarm-NG focuses on many few-body systems, e.g., thousands of systems with 3...15 bodies each, as is typical for the study of planetary systems. Swarm-NG parallelizes the simulation, including both the numerical integration of the equations of motion and the evaluation of forces using NVIDIA's "Compute Unified Device Architecture" (CUDA) on the GPU. Swarm-NG includes optimized implementations of 4th order time-symmetrized Hermite integration and mixed variable symplectic integration, as well as several sample codes for other algorithms to illustrate how non-CUDA-savvy users may themselves introduce customized integrators into the Swarm-NG framework. To optimize performance, we analyze the effect of GPU-specific parameters on performance under double precision. Applications of Swarm-NG include studying the late stages of planet formation, testing the stability of planetary systems and evaluating the goodness-of-fit between many planetary system models and observations of extrasolar planet host stars (e.g., radial velocity, astrometry, transit timing). While Swarm-NG focuses on the parallel integration of many planetary systems,the underlying integrators could be applied to a wide variety of problems that require repeatedly integrating a set of ordinary differential equations many times using different initial conditions and/or parameter values.

We search for evidence of dark matter in the Milky Way by utilizing the stellar number density distribution and kinematics measured by the Sloan Digital Sky Survey (SDSS) to heliocentric distances exceeding ~10 kpc. We employ the cylindrically symmetric form of Jeans equations and focus on the morphology of the resulting acceleration maps, rather than the normalization of the total mass as done in previous, mostly local, studies. Jeans equations are first applied to a mock catalog based on a cosmologically derived N-body + SPH simulation, and the known acceleration (gradient of gravitational potential) is successfully recovered. The same simulation is also used to quantify the impact of dark matter on the total acceleration. We use Galfast, a code designed to quantitatively reproduce SDSS measurements and selection effects, to generate a synthetic stellar catalog. We apply Jeans equations to this catalog and produce two-dimensional maps of stellar acceleration. These maps reveal that in a Newtonian framework, the implied gravitational potential cannot be explained by visible matter alone. The acceleration experienced by stars at galactocentric distances of ~20 kpc is three times larger than what can be explained by purely visible matter. The application of an analytic method for estimating the dark matter halo axis ratio to SDSS data implies an oblate halo with q_DM = 0.47 +/- 0.14 within the same distance range. These techniques can be used to map the dark matter halo to much larger distances from the Galactic center using upcoming deep optical surveys, such as LSST.

We present a precise photometric calibration of the first 1.5 years of science imaging from the Pan-STARRS1 survey (PS1), an ongoing optical survey of the entire sky north of declination -30 degrees in five bands. Building on the techniques employed by Padmanabhan et al. (2008) in the Sloan Digital Sky Survey (SDSS), we use repeat PS1 observations of stars to perform the relative calibration of PS1 in each of its five bands, solving simultaneously for the system throughput, the atmospheric transparency, and the large-scale detector flat field. Both internal consistency tests and comparison against the SDSS indicate that we achieve relative precision of <10 mmag in g, r, and i_P1, and ~10 mmag in z and y_P1. The spatial structure of the differences with the SDSS indicates that errors in both the PS1 and SDSS photometric calibration contribute similarly to the differences. The analysis suggests that both the PS1 system and the Haleakala site will enable <1% photometry over much of the sky.

We report the discovery of two new halo velocity groups (Cancer groups A and B) traced by 8 distant RR Lyrae stars and observed by the Palomar Transient Factory (PTF) survey at R.A.~129 deg, Dec~20 deg (l~205 deg, b~32 deg). Located at 92 kpc from the Galactic center (86 kpc from the Sun), these are some of the most distant substructures in the Galactic halo known to date. Follow-up spectroscopic observations with the Palomar Observatory 5.1-m Hale telescope and W. M. Keck Observatory 10-m Keck I telescope indicate that the two groups are moving away from the Galaxy at v_{gsr} = 78.0+-5.6 km/s (Cancer group A) and v_{gsr} = 16.3+-7.1 km/s (Cancer group B). The groups have velocity dispersions of \sigma_{v_{gsr}}=12.4+-5.0 km/s and \sigma_{v_{gsr}}=14.9+-6.2 km/s, and are spatially extended (about several kpc) making it very unlikely that they are bound systems, and are more likely to be debris of tidally disrupted dwarf galaxies or globular clusters. Both groups are metal-poor (median metallicities of [Fe/H] = -1.6 dex and [Fe/H] =-2.1 dex), and have a somewhat uncertain (due to small sample size) metallicity dispersion of ~0.4 dex, suggesting dwarf galaxies as progenitors. Two additional RR Lyrae stars with velocities consistent with those of the Cancer groups have been observed ~25 deg east, suggesting possible extension of the groups in that direction.

We present low-resolution MMT Hectospec spectroscopy of 594 candidate Monoceros stream member stars. Based on strong color-magnitude diagram overdensities, we targeted three fields within the stream's footprint, with 178 deg < l < 203 deg and -25 deg < b < 25 deg. By comparing the measured iron abundances with those expected from smooth Galactic components alone, we measure, for the first time, the spectroscopic metallicity distribution function for Monoceros. We find the stream to be chemically distinct from both the thick disk and halo, with [Fe/H] = -1, and do not detect a trend in the stream's metallicity with Galactic longitude. Passing from b = +25 deg to b = -25 deg the median Monoceros metallicity trends upward by 0.1 dex, though uncertainties in modeling sample contamination by the disk and halo make this a marginal detection. In each field, we find Monoceros to have an intrinsic [Fe/H] dispersion of 0.10-0.22 dex. From the CaII K line, we measure [Ca/Fe] for a subsample of metal poor program stars with -1.1 < [Fe/H] < -0.5. In two of three fields, we find calcium deficiencies qualitatively similar to previously reported [Ti/Fe] underabundances in Monoceros and the Sagittarius tidal stream. Further, using 90 spectra of thick disk stars in the Monoceros pointings with b ~ +/-25 deg, we detect a 0.22 dex north/south metallicity asymmetry coincident with known stellar density asymmetry at R_GC ~ 12 kpc and |Z| ~ 1.7 kpc. Our median Monoceros [Fe/H] = -1.0 and its relatively low dispersion naturally fit the expectation for an appropriately luminous M_V ~ -13 dwarf galaxy progenitor.

We use the Eighth Data Release of Sloan Digital Sky Survey (SDSS DR8) catalog with its additional sky coverage of the southern Galactic hemisphere, to measure the extent and study the nature of the Virgo Overdensity (VOD; Juric et al. 2008). The data show that the VOD extends over no less than 2000 deg^2, with its true extent likely closer to 3000 deg^2. We test whether the VOD can be attributed to a tilt in the stellar halo ellipsoid with respect to the plane of the Galactic disk and find that the observed symmetry of the north-south Galactic hemisphere star counts excludes this possibility. We argue that the Virgo Overdensity, in spite of its wide area and cloud-like appearance, is still best explained by a minor merger. Its appearance and position is qualitatively similar to a near perigalacticon merger event and, assuming that the VOD and the Virgo Stellar Stream share the same progenitor, consistent with the VSS orbit determined by Casetti-Dinescu et al. (2009).

We present the ensemble properties of 31 comets (27 resolved and 4 unresolved) observed by the Sloan Digital Sky Survey (SDSS). This sample of comets represents about 1 comet per 10 million SDSS photometric objects. Five-band (u,g,r,i,z) photometry is used to determine the comets' colors, sizes, surface brightness profiles, and rates of dust production in terms of the Af{\rho} formalism. We find that the cumulative luminosity function for the Jupiter Family Comets in our sample is well fit by a power law of the form N(< H) \propto 10(0.49\pm0.05)H for H < 18, with evidence of a much shallower fit N(< H) \propto 10(0.19\pm0.03)H for the faint (14.5 < H < 18) comets. The resolved comets show an extremely narrow distribution of colors (0.57 \pm 0.05 in g - r for example), which are statistically indistinguishable from that of the Jupiter Trojans. Further, there is no evidence of correlation between color and physical, dynamical, or observational parameters for the observed comets.

We use Canada-France-Hawaii Telescope Legacy Survey data for 170 deg^2, recalibrated and transformed to the Sloan Digital Sky Survey ugri photometric system, to study the distribution of near-turnoff main-sequence stars in the Galactic halo along four lines of sight to heliocentric distances of ~35 kpc. We find that the halo stellar number density profile becomes steeper at Galactocentric distances greater than R_{gal}~28 kpc, with the power law index changing from n_{inner}=-2.62+-0.04 to n_{outer}=-3.8+-0.1. In particular, we test a series of single power law models and find them to be strongly disfavored by the data. The parameters for the best-fit Einasto profile are n=2.2+-0.2 and R_e=22.2+-0.4 kpc. We measure the oblateness of the halo to be q=c/a=0.70+-0.01 and detect no evidence of it changing across the range of probed distances. The Sagittarius stream is detected in the l=173 and b=-62 direction as an overdensity of [Fe/H] ~ -1.5 dex stars at R_{gal}~32 kpc, providing a new constraint for the Sagittarius stream and dark matter halo models. We also detect the Monoceros stream as an overdensity of [Fe/H] > -1.5 dex stars in the l=232 and b=26 direction at R_{gal}<25 kpc. In the two sightlines where we do not detect significant substructure, the median metallicity is found to be independent of distance within systematic uncertainties ([Fe/H] ~ -1.5+-0.1 dex).

The advent of deep, wide, accurate, digital photometric surveys exemplified by the Sloan Digital Sky Survey (SDSS) has had a profound impact on studies of the Milky Way. In the past decade, we have transitioned from a scarcity to an (over)abundance of precise, well calibrated, observations of stars over a large fraction of the Galaxy. The avalanche of data will continue throughout this decade, culminating with Gaia and LSST. This new reality will necessitate changes in methodology, habits, and expectations both on the side of the large survey projects as well as the astrophysics community at large. We argue, based on the experience with SDSS, that surveys should release data as early and often as possible incorporating incremental improvements in each subsequent release, as opposed to holding off for a single, big, final release. The scientific community will need to reciprocate by performing analyses and (re-analyses) appropriate to the current fidelity of the released data, understanding that these are continually evolving and improving products.

We study Milky Way kinematics using a sample of 18.8 million main-sequence stars with r<20 and proper-motion measurements derived from SDSS and POSS astrometry, including ~170,000 stars with radial-velocity measurements from the SDSS spectroscopic survey. Distances to stars are determined using a photometric parallax relation, covering a distance range from ~100 pc to 10 kpc over a quarter of the sky at high Galactic latitudes (|b|>20 degrees). We find that in the region defined by 1 kpc <Z< 5 kpc and 3 kpc <R< 13 kpc, the rotational velocity for disk stars smoothly decreases, and all three components of the velocity dispersion increase, with distance from the Galactic plane. In contrast, the velocity ellipsoid for halo stars is aligned with a spherical coordinate system and appears to be spatially invariant within the probed volume. The velocity distribution of nearby ($Z<1$ kpc) K/M stars is complex, and cannot be described by a standard Schwarzschild ellipsoid. For stars in a distance-limited subsample of stars (<100 pc), we detect a multimodal velocity distribution consistent with that seen by HIPPARCOS. This strong non-Gaussianity significantly affects the measurements of the velocity ellipsoid tilt and vertex deviation when using the Schwarzschild approximation. We develop and test a simple descriptive model for the overall kinematic behavior that captures these features over most of the probed volume, and can be used to search for substructure in kinematic and metallicity space. We use this model to predict further improvements in kinematic mapping of the Galaxy expected from Gaia and LSST.

A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy.

We present an improved analysis of halo substructure traced by RR Lyrae stars in the 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) towards R.A.2h-3h and Dec~0 deg 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.

In order to attract and retain excellent researchers and diverse individuals in astrophysics, we recommend action be taken in several key areas impacting young scientists: (1) Maintain balance between large collaborations and individual projects through distribution of funding; encourage public releases of observational and simulation data for use by a broader community. (2) Improve the involvement of women, particularly at leading institutions. (3) Address the critical shortage of child care options and design reasonable profession-wide parental leave policies. (4) Streamline the job application and hiring process. We summarize our reasons for bringing these areas to the attention of the committee, and we suggest several practical steps that can be taken to address them.

Using SDSS Data Release 6, we construct two independent samples of candidate stellar wide binaries selected as i) pairs of unresolved sources with angular separation in the range 3'' - 16'', ii) common proper motion pairs with 5'' - 30'' angular separation, and make them publicly available. These samples are dominated by disk stars, and we use them to constrain the shape of the main-sequence photometric parallax relation M_r(r-i) and to study the properties of wide binary systems. We estimate M_r(r-i) by searching for a relation that minimizes the difference between distance moduli of primary and secondary components of wide binary candidates. We model M_r(r-i) by a fourth degree polynomial and determine the coefficients using Markov Chain Monte Carlo fitting, independently for each sample. Aided by the derived photometric parallax relation, we construct a series of high-quality catalogs of candidate main-sequence binary stars. Using these catalogs, we study the distribution of semi-major axes of wide binaries, a, in the 2,000 < a < 47,000 AU range. We find the observations to be well described by the Opik distribution, f(a)~1/a, for a<a_{break}, where a_{break} increases roughly linearly with the height Z above the Galactic plane (a_{break}~12,300 Z[kpc]^0.7 AU). The number of wide binary systems with 100 AU < a < a_{break}, as a fraction of the total number of stars, decreases from 0.9% at Z=0.5 kpc to 0.5% at Z=3 kpc. The probability for a star to be in a wide binary system is independent of its color. Given this color, the companions of red components seem to be drawn randomly from the stellar luminosity function, while blue components have a larger blue-to-red companion ratio than expected from luminosity function.

Asteroid families, traditionally defined as clusters of objects in orbital parameter space, often have distinctive optical colors. We show that the separation of family members from background interlopers can be improved with the aid of SDSS colors as a qualifier for family membership. Based on an ~88,000 object subset of the Sloan Digital Sky Survey Moving Object Catalog 4 with available proper orbital elements, we define 37 statistically robust asteroid families with at least 100 members using a simple Gaussian distribution model in both orbital and color space. The interloper rejection rate based on colors is typically ~10% for a given orbital family definition, with four families that can be reliably isolated only with the aid of colors. About 50% of all objects in this data set belong to families, and this fraction varies from about 35% for objects brighter than an H magnitude of 13 and rises to 60% for objects fainter than this. The fraction of C-type objects in families decreases with increasing H magnitude for H > 13, while the fraction of S-type objects above this limit remains effectively constant. This suggests that S-type objects require a shorter timescale for equilibrating the background and family size distributions via collisional processing. The size distributions for 15 families display a well-defined change of slope and can be modeled as a "broken" double power-law. Such "broken" size distributions are twice as likely for S-type familes than for C-type families, and are dominated by dynamically old families. The remaining families with size distributions that can be modeled as a single power law are dominated by young families. When size distribution requires a double power-law model, the two slopes are correlated and are steeper for S-type families.

We explore the possibility that the observed eccentricity distribution of extrasolar planets arose through planet-planet interactions, after the initial stage of planet formation was complete. Our results are based on ~3250 numerical integrations of ensembles of randomly constructed planetary systems, each lasting 100 Myr. We find that for a remarkably wide range of initial conditions the eccentricity distributions of dynamically active planetary systems relax towards a common final equilibrium distribution, well described by the fitting formula dn ~ e exp[-1/2 (e/0.3)^2] de. This distribution agrees well with the observed eccentricity distribution for e > 0.2, but predicts too few planets at lower eccentricities, even when we exclude planets subject to tidal circularization. These findings suggest that a period of large-scale dynamical instability has occurred in a significant fraction of newly formed planetary systems, lasting 1--2 orders of magnitude longer than the ~1 Myr interval in which gas-giant planets are assembled. This mechanism predicts no (or weak) correlations between semimajor axis, eccentricity, inclination, and mass in dynamically relaxed planetary systems. An additional observational consequence of dynamical relaxation is a significant population of planets (>10%) that are highly inclined (>25deg) with respect to the initial symmetry plane of the protoplanetary disk; this population may be detectable in transiting planets through the Rossiter-McLaughlin effect.

Using effective 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 for estimating these parameters from the SDSS u-g and g-r colors. We apply this method to SDSS photometric data for about 2 million F/G stars and measure 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 two components also possess the kinematics expected for disk and halo stars. The metallicity of the halo component is spatially invariant, while the median disk metallicity smoothly decreases with distance from the Galactic plane from -0.6 at 500 pc to -0.8 beyond several kpc. The absence of a correlation between metallicity and kinematics for disk stars is in a conflict with the traditional decomposition in terms of thin and thick disks. We detect coherent substructures in the kinematics--metallicity space, such as the Monoceros stream, which rotates faster than the LSR, and has a median metallicity of [Fe/H]=-0.96, 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 the LSST will obtain metallicity measurements accurate to 0.2 dex or better, with proper motion measurements accurate to ~0.2 mas/yr, for about 200 million F/G dwarf stars within a distance limit of ~100 kpc (g<23.5). [abridged]

Abridged: We estimate the distances to ~48 million stars detected by the Sloan Digital Sky Survey and map their 3D number density distribution in 100 < D < 20 kpc range over 6,500 deg^2 of sky. The data show strong evidence for a Galaxy consisting of an oblate halo, a disk component, and a number of localized overdensities with exponential disk parameters (bias-corrected for an assumed 35% binary fraction) H_1 = 300 pc, L_1 = 2600 pc, H_2 = 900 pc, L_2 = 3600 pc, and local density normalization of 12%. We find the halo to be oblate, with best-fit axis ratio c/a = 0.64, r^{-2.8} profile, and the local halo-to-thin disk normalization of 0.5%. We estimate the errors of derived model parameters to be no larger than ~20% (disk scales) and ~10% (thick disk normalization). While generally consistent with the above model, the density distribution shows a number of statistically significant localized deviations. We detect two overdensities in the thick disk region at (R, Z) ~ (6.5, 1.5)kpc and (R, Z) ~ (9.5, 0.8) kpc, and a remarkable density enhancement in the halo covering >1000deg^2 of sky towards the constellation of Virgo, at distances of ~6-20 kpc. Compared to a region symmetric with respect to the l=0 line, 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. 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 ~50pc in the disk, to ~1 - 2 kpc in the halo.

We quantify the variability of faint unresolved optical sources using a catalog based on multiple SDSS imaging observations. The catalog covers SDSS Stripe 82, and contains 58 million photometric observations in the SDSS ugriz system for 1.4 million unresolved sources. In each photometric bandpass we compute various low-order lightcurve statistics 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. The majority (2/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 is the UV excess color selection (at high Galactic latitudes). We analyze the distribution of lightcurve skewness for quasars and find that is centered on zero. We find that about 1/4 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 lightcurve skewness in the g-r vs. 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 six hundred RR Lyrae stars, we demonstrate rich halo substructure out to distances of 100 kpc. We extrapolate these results to expected performance by the Large Synoptic Survey Telescope and estimate that it will obtain well-sampled 2% accurate, multi-color lightcurves for ~2 million low-redshift quasars, and will discover at least 50 million variable stars.

We analyze 1187 observations of about 860 unique candidate Jovian Trojan asteroids listed in the 3rd release of Sloan Digital Sky Survey (SDSS) Moving Object Catalog. The sample is complete at the faint end to r=21.2 mag (apparent brightness) and H=13.8 (absolute brightness, approximately corresponding to 10 km diameter). A subset of 297 detections of previously known Trojans were used to design and optimize a selection method based on observed angular velocity that resulted in the remaining objects. Using a sample of objects with known orbits, we estimate that the candidate sample contamination is about 3%. The well-controlled selection effects, the sample size, depth and accurate five-band UV-IR photometry enabled several new findings and the placement of older results on a firmer statistical footing. We find that there are significantly more asteroids in the leading swarm (L4) than in the trailing swarm (L5): N(L4)/N(L5)=1.6$\pm$0.1, independently of limiting object's size. The overall counts normalization suggests that there are about as many Jovians Trojans as there are main-belt asteroids down to the same size limit, in agreement with earlier estimates. We find that Trojan asteroids have a remarkably narrow color distribution (root-mean-scatter of only $\sim$0.05 mag) that is significantly different from the color distribution of the main-belt asteroids. The color of Trojan asteroids is correlated with their orbital inclination, in a similar way for both swarms, but appears uncorrelated with the object's size. We extrapolate the results presented here and estimate that Large Synoptic Survey Telescope will determine orbits, accurate colors and measure light curves in six photometric bandpasses for about 100,000 Jovian Trojan asteroids.

The recently discovered Virgo stellar over-density, which expands over \~1000deg^2 perpendicularly to the Galactic disk plane (7< Z <15 kpc, R~7 kpc), is the largest clump of tidal debris ever detected in the outer halo and is likely related with the accretion of a nearby dwarf galaxy by the Milky Way. We carry out N-body simulations of the Sagittarius stream to show that this giant stellar over-density is a confirmation of theoretical model predictions for the leading tail of the Sagittarius stream to cross the Milky Way plane in the Solar neighborhood. Radial velocity measurements are needed to confirm this association and to further constrain the shape of the Milky Way dark matter halo through a new generation of theoretical models. If the identification of Virgo over-density and the Sagittarius leading arm is correct, we predict highly negative radial velocities for the stars of Virgo over-density. The detection of this new portion of the Sagittarius tidal stream would represent an excellent target for the on-going and future kinematic surveys and for dark matter direct detection experiments in the proximity of the Sun.

In addition to optical photometry of unprecedented quality, the Sloan Digital Sky Survey (SDSS) is also producing a massive spectroscopic database. We discuss determination of stellar parameters, such as effective temperature, gravity and metallicity from SDSS spectra, describe correlations between kinematics and metallicity, and study their variation as a function of the position in the Galaxy. We show that stellar parameter estimates by Beers et al. show a good correlation with the position of a star in the g-r vs. u-g color-color diagram, thereby demonstrating their robustness as well as a potential for photometric parameter estimation methods. Using Beers et al. parameters, we find that the metallicity distribution of the Milky Way stars at a few kpc from the galactic plane is bimodal with a local minimum at [Z/Zo]~ -1.3. The median metallicity for the low-metallicity [Z/Zo]<-1.3 subsample is nearly independent of Galactic cylindrical coordinates R and z, while it decreases with z for the high-metallicity [Z/Zo]> -1.3 sample. We also find that the low-metallicity sample has ~2.5 times larger velocity dispersion and that it does not rotate (at the ~10 km/s level), while the rotational velocity of the high-metallicity sample decreases smoothly with the height above the galactic plane.

(Abridged) The Large Synoptic Survey Telescope (LSST) is currently by far the most ambitious proposed ground-based optical survey. Solar System mapping is one of the four key scientific design drivers, with emphasis on efficient Near-Earth Object (NEO) and Potentially Hazardous Asteroid (PHA) detection, orbit determination, and characterization. In a continuous observing campaign of pairs of 15 second exposures of its 3,200 megapixel camera, LSST will cover the entire available sky every three nights in two photometric bands to a depth of V=25 per visit (two exposures), with exquisitely accurate astrometry and photometry. Over the proposed survey lifetime of 10 years, each sky location would be visited about 1000 times. The baseline design satisfies strong constraints on the cadence of observations mandated by PHAs such as closely spaced pairs of observations to link different detections and short exposures to avoid trailing losses. Equally important, due to frequent repeat visits LSST will effectively provide its own follow-up to derive orbits for detected moving objects. Detailed modeling of LSST operations, incorporating real historical weather and seeing data from LSST site at Cerro Pachon, shows that LSST using its baseline design cadence could find 90% of the PHAs with diameters larger than 250 m, and 75% of those greater than 140 m within ten years. However, by optimizing sky coverage, the ongoing simulations suggest that the LSST system, with its first light in 2013, can reach the Congressional mandate of cataloging 90% of PHAs larger than 140m by 2020.

We discuss several results made possible by accurate SDSS astrometric measurements in a large sky area, with emphasis on asteroids and stellar proper motions obtained by comparing POSS and SDSS. SDSS has observed over 200,000 moving objects in five photometric bands, corresponding to about two orders of magnitude increase over previous multi--color surveys. These data were used to extend the measurement of asteroid size distribution to a smaller size limit, to demonstrate that asteroid dynamical families, defined as clusters in orbital parameter space, also strongly segregate in color space, and to discover a correlation between asteroid age and colors. A preliminary analysis of SDSS-POSS proper motions for about 1 million M dwarf stars demonstrates that, in the 0.1-1 kpc distance range, the rotational velocity and its dispersion for disk stars increase with the distance from the Galactic plane.

We compare Stetson's photometric standards with measurements listed in a standard star catalog constructed using repeated SDSS imaging observations. The SDSS catalog includes over 700,000 candidate standard stars from the equatorial stripe 82 (|Dec|<1.266 deg) in the RA range 20h 34' to 4h 00', and with the $r$ band magnitudes in the range 14--21. The distributions of measurements for individual sources demonstrate that the SDSS photometric pipeline correctly estimates random photometric errors, which are below 0.01 mag for stars brighter than (19.5, 20.5, 20.5, 20, 18.5) in ugriz, respectively (about twice as good as for individual SDSS runs). We derive mean photometric transformations between the SDSS gri and the BVRI system using 1165 Stetson stars found in the equatorial stripe 82, and then study the spatial variation of the difference in zeropoints between the two catalogs. Using third order polynomials to describe the color terms, we find that photometric measurements for main-sequence stars can be transformed between the two systems with systematic errors smaller than a few millimagnitudes. The spatial variation of photometric zeropoints in the two catalogs typically does not exceed 0.01 magnitude. Consequently, the SDSS Standard Star Catalog for Stripe 82 can be used to calibrate new data in both the SDSS ugriz and the BVRI systems with a similar accuracy.

We discuss the panchromatic properties of 99,088 galaxies selected from the SDSS Data Release 1 spectroscopic sample (a flux-limited sample for 1360 deg^2). These galaxies are positionally matched to sources detected by ROSAT, GALEX, 2MASS, IRAS, GB6, FIRST, NVSS and WENSS. We find strong correlations between the detection fraction at other wavelengths and optical properties such as flux, colors, and emission-line strengths. Using GALEX, SDSS, and 2MASS data, we construct the UV-IR broad-band spectral energy distributions for various types of galaxies, and find that they form a nearly one-parameter family. For example, based on SDSS u- and r-band data, supplemented with redshift, the K-band 2MASS magnitudes can be "predicted" with an rms scatter of only 0.2 mag. When a dust content estimate determined from SDSS data by Kauffmann et al. (2003) is also utilized, this scatter decreases to 0.1 mag. We demonstrate that this dust content is indeed higher for galaxies detected by IRAS and that it can be used to "predict" measured IRAS 60 micron flux density within a factor of two using only SDSS data. We also show that the position of a galaxy in the emission-line-based Baldwin-Phillips-Terlevich diagram is correlated with the optical light concentration index and u-r color determined from the SDSS broad-band imaging data, and discuss changes in the morphology of this diagram induced by requiring detections at other wavelengths. We study the IR-radio correlation and find evidence that its slope may be different for AGN and star-forming galaxies and related to the H_alpha/H_beta line strength ratio.

We analyze the properties of quasar variability using repeated SDSS imaging data in five UV-to-far red photometric bands, accurate to 0.02 mag, for 13,000 spectroscopically confirmed quasars. The observed time lags span the range from 3 hours to over 3 years, and constrain the quasar variability for rest-frame time lags of up to two years, and at rest-frame wavelengths from 1000 Ang. to 6000 Ang. We demonstrate that 66,000 SDSS measurements of magnitude differences can be described within the measurement noise by a simple function of only three free parameters. The addition of POSS data constrains the long-term behavior of quasar variability and provides evidence for a turn-over in the structure function. This turn-over indicates that the characteristic time scale for optical variability of quasars is of the order 1 year.

We positionally match the 6 cm GB6, 20 cm FIRST and NVSS, and 92 cm WENSS radio catalogs and find 16,500 matches in ~3,000 deg2 of sky. Using this unified radio database, we construct radio "color-magnitude-morphology" diagrams and find that they display a clear structure, rather than a random scatter. We propose a simple, yet powerful, method for morphological classification of radio sources based on FIRST and NVSS measurements. For a subset of matched sources, we find optical identifications using the SDSS Data Release 1 catalogs, and separate them into quasars and galaxies. Compact radio sources with flat radio spectra are dominated by quasars, while compact sources with steep spectra, and resolved radio sources, contain substantial numbers of both quasars and galaxies.

We report a detection of statistically significant color variations for a sample of 7,531 multiply observed asteroids that are listed in the Sloan Digital Sky Survey Moving Object Catalog. Using 5-band photometric observations accurate to ~0.02 mag, we detect color variations in the range 0.06-0.11 mag (rms). These variations appear uncorrelated with asteroids physical characteristics such as diameter (in the probed 1-10 km range), taxonomic class, and family membership. Despite such a lack of correlation, which implies a random nature for the cause of color variability, a suite of tests suggest that the detected variations are not instrumental effects. In particular, the observed color variations are incompatible with photometric errors, and, for objects observed at least four times, the color change in the first pair of observations is correlated with the color change in the second pair. These facts strongly suggest that the observed effect is real, and also indicate that for some asteroids color variations are larger than for others. The detected color variations can be explained as due to inhomogeneous albedo distribution over an asteroid surface. Although relatively small, these variations suggest that fairly large patches with different color than their surroundings exist on a significant fraction of asteroids. This conclusion is in agreement with spatially resolved color images of several large asteroids obtained by NEAR spacecraft and HST.

The potential of the Sloan Digital Sky Survey for wide-field variability studies is illustrated using multi-epoch observations for 3,000,000 point sources observed in 700 deg2 of sky, with time spans ranging from 3 hours to 3 years. These repeated observations of the same sources demonstrate that SDSS delivers ~0.02 mag photometry with well behaved and understood errors. We show that quasars dominate optically faint (r > 18) point sources that are variable on time scales longer than a few months, while for shorter time scales, and at bright magnitudes, most variable sources are stars.

We announce the first public release of the SDSS Moving Object Catalog, with SDSS observations for 58,117 asteroids. The catalog lists astrometric and photometric data for moving objects observed prior to Dec 15, 2001, and also includes orbital elements for 10,592 previously known objects. We analyze the correlation between the orbital parameters and optical colors for the known objects, and confirm that asteroid dynamical families, defined as clusters in orbital parameter space, also strongly segregate in color space. Their distinctive optical colors indicate that the variations in chemical composition within a family are much smaller than the compositional differences between families, and strongly support earlier suggestions that asteroids belonging to a particular family have a common origin.

We discuss optical colors of 10,592 asteroids with known orbits selected from a sample of 58,000 moving objects observed by the Sloan Digital Sky Survey (SDSS). This is more than ten times larger sample that includes both orbital parameters and multi-band photometric measurements than previously available. We confirm that asteroid dynamical families, defined as clusters in orbital parameter space, also strongly segregate in color space. In particular, we demonstrate that the three major asteroid families (Eos, Koronis, and Themis), together with the Vesta family, represent four main asteroid color types. Their distinctive optical colors indicate that the variations in chemical composition within a family are much smaller than the compositional differences between families, and strongly support earlier suggestions that asteroids belonging to a particular family have a common origin. We estimate that over 90% of asteroids belong to families.

We positionally correlate asteroids from existing catalogs with a sample of $\about$18,000 asteroids detected by the Sloan Digital Sky Survey (SDSS, Ivezi\'{c} {\em et al.} 2001). We find 2641 unique matches, which represent the largest sample of asteroids with both accurate multi-color photometry and known orbital parameters. The matched objects are predominantly bright, and demonstrate that the SDSS photometric pipeline recovers \about90% of the known asteroids in the observed region. For the recovered asteroids we find a large offset (\about 0.4 mag) between Johnson V magnitudes derived from SDSS photometry and the predicted catalog-based visual magnitudes. This offset varies with the asteroid color from 0.34 mag for blue asteroids to 0.44 mag for red asteroids, and is probably caused by the use of unfiltered CCD observations in the majority of recent asteroid surveys. This systematic photometric error leads to an overestimate of the number of asteroids brighter than a given absolute magnitude limit by a factor of \about 1.7. The distribution of the matched asteroids in orbital parameter space indicates strong color segregation. We confirm that some families are dominated by a single asteroid type (e.g. the Koronis family by red asteroids and the Themis family by blue asteroids), while others appear to be a mixture of blue and red objects (e.g. the Nysa/Polana family). Asteroids with the bluest \iz colors, which can be associated with the Vesta family, show particularly striking localization in orbital parameter space.