We obtained thorough photometric observations of two binary near-Earth asteroids (66391) Moshup = 1999 KW4 and (88710) 2001 SL9 from June 2000 to June 2019 and derived physical and dynamical properties of the binary systems. We found that the data for 1999 KW4 are inconsistent with a constant orbital period and we obtained unique solution with a quadratic drift of the mean anomaly of the satellite of -0.65 +/- 0.16 deg/yr2 (all quoted uncertainties are 3sigma). This means that the semimajor axis of the mutual orbit of the components of this binary system increases in time with a mean rate of 1.2 +/- 0.3 cm/yr. The data for 2001 SL9 are also inconsistent with a constant orbital period and we obtained two solutions for the quadratic drift of the mean anomaly: 2.8 +/- 0.2 and 5.2 +/- 0.2 deg/yr2, implying that the semimajor axis of the mutual orbit of the components decreases in time with a mean rate of -2.8 +/- 0.2 or -5.1 +/- 0.2 cm/yr for the two solutions, respectively. The expanding orbit of 1999 KW4 may be explained by mutual tides interplaying with binary YORP (BYORP) effect (McMahon and Scheeres, 2010). However, a modeling of the BYORP drift using radar-derived shapes of the binary components predicted a much higher value of the orbital drift than the observed one. It suggests that either the radar-derived shape model of the secondary is inadequate for computing the BYORP effect, or the present theory of BYORP overestimates it. It is possible that the BYORP coefficient has instead an opposite sign than predicted; in that case, the system may be moving into an equilibrium between the BYORP and the tides. In the case of 2001 SL9, the BYORP effect is the only known physical mechanism that can cause the inward drift of its mutual orbit.

Prompted by peculiar spectroscopic variability observed in SDSS/APOGEE $H$-band spectra, we monitored the Be star HD 55606 using optical spectroscopy and found that it is an exotic double-lined spectroscopic binary (SB2) consisting of a Be star and a hot, compact companion that is probably an OB subdwarf (sdOB) star. Motion of the sdOB star is traced by its impact on the strong He~I lines, observed as radial velocity ($V_{r}$) variable, double-peaked emission profiles with narrow central absorption cores. Weak He II 4686 {\AA} absorption associated with the companion star is detected in most spectra. Use of the emission peaks of low-ionization emission lines to trace the Be star $V_{r}$ and the He I lines to trace the companion star $V_{r}$ yields a circular orbital solution with a 93.8-day period and masses of $M_{\rm Be}=6.2$ $M_{\rm \odot}$ and $M_{\rm sdOB}=0.9$ $M_{\rm \odot}$ in the case of $i=80^{\circ}$. HD 55606 exhibits a variety of phase-locked variability, including the development of shell lines twice per orbit. The shell phases coincide with variation in the double emission peak separations, and both forms of variability are likely caused by a two-armed spiral density perturbation in the Be disk. The intensity ratios of the double emission peaks are also phase-locked, possibly indicating heating by the sdOB star of the side of the Be disk facing it. HD 55606 is a new member of the growing sample of Be+sdOB binaries, in which the Be star's rapid rotation and ability to form a disk can be attributed to past mass transfer.

We demonstrate a path to hitherto unachievable differential photometric precisions from the ground, both in the optical and near-infrared (NIR), using custom-fabricated beam-shaping diffusers produced using specialized nanofabrication techniques. Such diffusers mold the focal plane image of a star into a broad and stable top-hat shape, minimizing photometric errors due to non-uniform pixel response, atmospheric seeing effects, imperfect guiding, and telescope-induced variable aberrations seen in defocusing. This PSF reshaping significantly increases the achievable dynamic range of our observations, increasing our observing efficiency and thus better averages over scintillation. Diffusers work in both collimated and converging beams. We present diffuser-assisted optical observations demonstrating $62^{+26}_{-16}$ppm precision in 30 minute bins on a nearby bright star 16-Cygni A (V=5.95) using the ARC 3.5m telescope---within a factor of $\sim$2 of Kepler's photometric precision on the same star. We also show a transit of WASP-85-Ab (V=11.2) and TRES-3b (V=12.4), where the residuals bin down to $180^{+66}_{-41}$ppm in 30 minute bins for WASP-85-Ab---a factor of $\sim$4 of the precision achieved by the K2 mission on this target---and to 101ppm for TRES-3b. In the NIR, where diffusers may provide even more significant improvements over the current state of the art, our preliminary tests have demonstrated $137^{+64}_{-36}$ppm precision for a $K_S =10.8$ star on the 200" Hale Telescope. These photometric precisions match or surpass the expected photometric precisions of TESS for the same magnitude range. This technology is inexpensive, scalable, easily adaptable, and can have an important and immediate impact on the observations of transits and secondary eclipses of exoplanets.

Lunar laser ranging provides a number of leading experimental tests of gravitation -- important in our quest to unify General Relativity and the Standard Model of physics. The Apache Point Observatory Lunar Laser-ranging Operation (APOLLO) has for years achieved median range precision at the ~2 mm level. Yet residuals in model-measurement comparisons are an order-of-magnitude larger, raising the question of whether the ranging data are not nearly as accurate as they are precise, or if the models are incomplete or ill-conditioned. This paper describes a new absolute calibration system (ACS) intended both as a tool for exposing and eliminating sources of systematic error, and also as a means to directly calibrate ranging data in-situ. The system consists of a high-repetition-rate (80 MHz) laser emitting short (< 10 ps) pulses that are locked to a cesium clock. In essence, the ACS delivers photons to the APOLLO detector at exquisitely well-defined time intervals as a "truth" input against which APOLLO's timing performance may be judged and corrected. Preliminary analysis indicates no inaccuracies in APOLLO data beyond the ~3 mm level, suggesting that historical APOLLO data are of high quality and motivating continued work on model capabilities. The ACS provides the means to deliver APOLLO data both accurate and precise below the 2 mm level.

The 163 comets observed during the WISE/NEOWISE prime mission represent the largest infrared survey to date of comets, providing constraints on dust, nucleus sizes, and CO+CO2 production. We present detailed analyses of the WISE/NEOWISE comet discoveries, and discuss observations of the active comets showing 4.6 $\mu$m band excess. We find a possible relation between dust and CO+CO2 production, as well as possible differences in the sizes of long and short period comet nuclei.

Laser ranging measurements during the total lunar eclipse on 2010 December 21 verify previously suspected thermal lensing in the retroreflectors left on the lunar surface by the Apollo astronauts. Signal levels during the eclipse far exceeded those historically seen at full moon, and varied over an order of magnitude as the eclipse progressed. These variations can be understood via a straightforward thermal scenario involving solar absorption by a ~50% covering of dust that has accumulated on the front surfaces of the reflectors. The same mechanism can explain the long-term degradation of signal from the reflectors as well as the acute signal deficit observed near full moon.

The Z Cam stars IW And and V513 Cas are unusual in having outbursts following their standstills in contrast to the usual Z Cam behavior of quiescence following standstills. In order to gain further understanding of these little-studied systems, we obtained spectra correlated with photometry from the AAVSO throughout a 3-4 month interval in 2011. In addition, time-resolved spectra were obtained in 2012 that provided orbital periods of 3.7 hrs for IW And and 5.2 hrs for V513 Cas. The photometry of V513 Cas revealed a regular pattern of standstills and outbursts with little time at quiescence, while IW And underwent many excursions from quiescence to outburst to short standstills. The spectra of IW And are similar to normal dwarf novae, with strong Balmer emission at quiescence and absorption at outburst. In contrast, V513 Cas shows a much flatter/redder spectrum near outburst with strong HeII emission and prominent emission cores in the Balmer lines. Part of this continuum difference may be due to reddening effects. While our attempts to model the outburst and standstill states of IW And indicate a mass accretion rate near 3E-9 solar masses per year, we could find no obvious reason why these systems behave differently following standstill compared to normal Z Cam stars.

Here we present observations of 7 large Kuiper Belt Objects. From these observations, we extract a point source catalog with $\sim0.01"$ precision, and astrometry of our target Kuiper Belt Objects with $0.04-0.08"$ precision within that catalog. We have developed a new technique to predict the future occurrence of stellar occultations by Kuiper Belt Objects. The technique makes use of a maximum likelihood approach which determines the best-fit adjustment to cataloged orbital elements of an object. Using simulations of a theoretical object, we discuss the merits and weaknesses of this technique compared to the commonly adopted ephemeris offset approach. We demonstrate that both methods suffer from separate weaknesses, and thus, together provide a fair assessment of the true uncertainty in a particular prediction. We present occultation predictions made by both methods for the 7 tracked objects, with dates as late as 2015. Finally, we discuss observations of three separate close passages of Quaoar to field stars, which reveal the accuracy of the element adjustment approach, and which also demonstrate the necessity of considering the uncertainty in stellar position when assessing potential occultations.

We describe the Pan-STARRS Moving Object Processing System (MOPS), a modern software package that produces automatic asteroid discoveries and identifications from catalogs of transient detections from next-generation astronomical survey telescopes. MOPS achieves > 99.5% efficiency in producing orbits from a synthetic but realistic population of asteroids whose measurements were simulated for a Pan-STARRS4-class telescope. Additionally, using a non-physical grid population, we demonstrate that MOPS can detect populations of currently unknown objects such as interstellar asteroids. MOPS has been adapted successfully to the prototype Pan-STARRS1 telescope despite differences in expected false detection rates, fill-factor loss and relatively sparse observing cadence compared to a hypothetical Pan-STARRS4 telescope and survey. MOPS remains >99.5% efficient at detecting objects on a single night but drops to 80% efficiency at producing orbits for objects detected on multiple nights. This loss is primarily due to configurable MOPS processing limits that are not yet tuned for the Pan-STARRS1 mission. The core MOPS software package is the product of more than 15 person-years of software development and incorporates countless additional years of effort in third-party software to perform lower-level functions such as spatial searching or orbit determination. We describe the high-level design of MOPS and essential subcomponents, the suitability of MOPS for other survey programs, and suggest a road map for future MOPS development.

Enhancements to the science data processing pipeline of NASA's Wide-field Infrared Explorer (WISE) mission, collectively known as NEOWISE, resulted in the detection of $>$158,000 minor planets in four infrared wavelengths during the fully cryogenic portion of the mission. Following the depletion of its cryogen, NASA's Planetary Science Directorate funded a four month extension to complete the survey of the inner edge of the Main Asteroid Belt and to detect and discover near-Earth objects (NEOs). This extended survey phase, known as the NEOWISE Post-Cryogenic Survey, resulted in the detection of $\sim$6500 large Main Belt asteroids and 88 NEOs in its 3.4 and 4.6 $\mu$m channels. During the Post-Cryogenic Survey, NEOWISE discovered and detected a number of asteroids co-orbital with the Earth and Mars, including the first known Earth Trojan. We present preliminary thermal fits for these and other NEOs detected during the 3-Band Cryogenic and Post-Cryogenic Surveys.

We present the preliminary results of an analysis of the sub-populations within the near-Earth asteroids, including the Atens, Apollos, Amors, and those that are considered potentially hazardous using data from the Wide-field Infrared Survey Explorer (WISE). In order to extrapolate the sample of objects detected by WISE to the greater population, we determined the survey biases for asteroids detected by the project's automated moving object processing system (known as NEOWISE) as a function of diameter, visible albedo, and orbital elements. Using this technique, we are able to place constraints on the number of potentially hazardous asteroids (PHAs) larger than 100 m and find that there are $\sim4700\pm1450$ such objects. As expected, the Atens, Apollos, and Amors are revealed by WISE to have somewhat different albedo distributions, with the Atens being brighter than the Amors. The cumulative size distributions of the various near-Earth object (NEO) subgroups vary slightly between 100 m and 1 km. A comparison of the observed orbital elements of the various sub-populations of the NEOs with the current best model is shown.

We present a reanalysis of the relationship between asteroid albedo and polarization properties using the albedos derived from the Wide-field Infrared Survey Explorer. We find that the function that best describes this relation is a three-dimensional linear fit in the space of log(albedo)-log(polarization slope)-log(minimum polarization). When projected to two dimensions the parameters of the fit are consistent with those found in previous work. We also define p* as the quantity of maximal polarization variation when compared with albedo and present the best fitting albedo-p* relation. Some asteroid taxonomic types stand out in this three-dimensional space, notably the E, B, and M Tholen types, while others cluster in clumps coincident with the S- and C-complex bodies. We note that both low albedo and small (D<30 km) asteroids are under-represented in the polarimetric sample, and we encourage future polarimetric surveys to focus on these bodies.

We report results based on mid-infrared photometry of 5 active main belt objects (AMBOs) detected by the Wide-field Infrared Survey Explorer (WISE) spacecraft. Four of these bodies, P/2010 R2 (La Sagra), 133P/Elst-Pizarro, (596) Scheila, and 176P/LINEAR, showed no signs of activity at the time of the observations, allowing the WISE detections to place firm constraints on their diameters and albedos. Geometric albedos were in the range of a few percent, and on the order of other measured comet nuclei. P/2010 A2 was observed on April 2-3, 2010, three months after its peak activity. Photometry of the coma at 12 and 22 {\mu}m combined with ground-based visible-wavelength measurements provides constraints on the dust particle mass distribution (PMD), dlogn/dlogm, yielding power-law slope values of {\alpha} = -0.5 +/- 0.1. This PMD is considerably more shallow than that found for other comets, in particular inbound particle fluence during the Stardust encounter of comet 81P/Wild 2. It is similar to the PMD seen for 9P/Tempel 1 in the immediate aftermath of the Deep Impact experiment. Upper limits for CO2 & CO production are also provided for each AMBO and compared with revised production numbers for WISE observations of 103P/Hartley 2.

We have combined the NEOWISE and Sloan Digital Sky Survey data to study the albedos of 24,353 asteroids with candidate taxonomic classifications derived using Sloan photometry. We find a wide range of moderate to high albedos for candidate S-type asteroids that are analogous to the S-complex defined by previous spectrophotometrically-based taxonomic systems. The candidate C-type asteroids, while generally very dark, have a tail of higher albedos that overlaps the S types. The albedo distribution for asteroids with a photometrically derived Q classification is extremely similar to those of the S types. Asteroids with similar colors to (4) Vesta have higher albedos than the S types, and most have orbital elements similar to known Vesta family members. Finally, we show that the relative reflectance at 3.4 and 4.6 $\mu$m is higher for D-type asteroids and suggest that their red visible and near-infrared spectral slope extends out to these wavelengths. Understanding the relationship between size, albedo, and taxonomic classification is complicated by the fact that the objects with classifications were selected from the visible/near-infrared Sloan Moving Object Catalog, which is biased against fainter asteroids, including those with lower albedos.

We present the preliminary analysis of 1023 known asteroids in the Hilda region of the Solar System observed by the NEOWISE component of the Wide-field Infrared Survey Explorer (WISE). The sizes of the Hildas observed range from $\sim 3 - 200$km. We find no size - albedo dependency as reported by other projects. The albedos of our sample are low, with a weighted mean value $p_V = 0.055\pm0.018$, for all sizes sampled by the NEOWISE survey. We observed a significant fraction of the objects in the two known collisional families in the Hilda population. It is found that the Hilda collisional family is brighter, with weighted mean albedo of $p_V = 0.061\pm0.011$, than the general population and dominated by D-type asteroids, while the Schubart collisional family is darker, with weighted mean albedo of ($p_V = 0.039\pm0.013$). Using the reflected sunlight in the two shortest WISE bandpasses we are able to derive a method for taxonomic classification of $\sim 10%$ of the Hildas detected in the NEOWISE survey. For the Hildas with diameter larger than 30km there are $67^{+7}_{-15}%$ D-type asteroids and $26^{+17}_{-5}%$ C-/P-type asteroids (with the majority of these being P-types).

We present the preliminary analysis of over 1739 known and 349 candidate Jovian Trojans observed by the NEOWISE component of the Wide-field Infrared Survey Explorer (WISE). With this survey the available diameters, albedos and beaming parameters for the Jovian Trojans have been increased by more than an order of magnitude compared to previous surveys. We find that the Jovian Trojan population is very homogenous for sizes larger than $\sim10$km (close to the detection limit of WISE for these objects). The observed sample consists almost exclusively of low albedo objects, having a mean albedo value of $0.07\pm0.03$. The beaming parameter was also derived for a large fraction of the observed sample, and it is also very homogenous with an observed mean value of $0.88\pm0.13$. Preliminary debiasing of the survey shows our observed sample is consistent with the leading cloud containing more objects than the trailing cloud. We estimate the fraction to be N(leading)/N(trailing) $\sim 1.4 \pm 0.2$, lower than the $1.6 \pm 0.1$ value derived by others.

With the NEOWISE portion of the \emph{Wide-field Infrared Survey Explorer} (WISE) project, we have carried out a highly uniform survey of the near-Earth object (NEO) population at thermal infrared wavelengths ranging from 3 to 22 $\mu$m, allowing us to refine estimates of their numbers, sizes, and albedos. The NEOWISE survey detected NEOs the same way whether they were previously known or not, subject to the availability of ground-based follow-up observations, resulting in the discovery of more than 130 new NEOs. The survey's uniformity in sensitivity, observing cadence, and image quality have permitted extrapolation of the 428 near-Earth asteroids (NEAs) detected by NEOWISE during the fully cryogenic portion of the WISE mission to the larger population. We find that there are 981$\pm$19 NEAs larger than 1 km and 20,500$\pm$3000 NEAs larger than 100 m. We show that the Spaceguard goal of detecting 90% of all 1 km NEAs has been met, and that the cumulative size distribution is best represented by a broken power law with a slope of 1.32$\pm$0.14 below 1.5 km. This power law slope produces $\sim13,200\pm$1,900 NEAs with $D>$140 m. Although previous studies predict another break in the cumulative size distribution below $D\sim$50-100 m, resulting in an increase in the number of NEOs in this size range and smaller, we did not detect enough objects to comment on this increase. The overall number for the NEA population between 100-1000 m are lower than previous estimates. The numbers of near-Earth comets will be the subject of future work.

The NEOWISE dataset offers the opportunity to study the variations in albedo for asteroid classification schemes based on visible and near-infrared observations for a large sample of minor planets. We have determined the albedos for nearly 1900 asteroids classified by the Tholen, Bus and Bus-DeMeo taxonomic classification schemes. We find that the S-complex spans a broad range of bright albedos, partially overlapping the low albedo C-complex at small sizes. As expected, the X-complex covers a wide range of albedos. The multi-wavelength infrared coverage provided by NEOWISE allows determination of the reflectivity at 3.4 and 4.6 $\mu$m relative to the visible albedo. The direct computation of the reflectivity at 3.4 and 4.6 $\mu$m enables a new means of comparing the various taxonomic classes. Although C, B, D and T asteroids all have similarly low visible albedos, the D and T types can be distinguished from the C and B types by examining their relative reflectance at 3.4 and 4.6 $\mu$m. All of the albedo distributions are strongly affected by selection biases against small, low albedo objects, as all objects selected for taxonomic classification were chosen according to their visible light brightness. Due to these strong selection biases, we are unable to determine whether or not there are correlations between size, albedo and space weathering. We argue that the current set of classified asteroids makes any such correlations difficult to verify. A sample of taxonomically classified asteroids drawn without significant albedo bias is needed in order to perform such an analysis.

We present initial results from the Wide-field Infrared Survey Explorer (WISE), a four-band all-sky thermal infrared survey that produces data well suited to measuring the physical properties of asteroids, and the NEOWISE enhancement to the WISE mission allowing for detailed study of Solar system objects. Using a NEATM thermal model fitting routine we compute diameters for over 100,000 Main Belt asteroids from their IR thermal flux, with errors better than 10%. We then incorporate literature values of visible measurements (in the form of the H absolute magnitude) to determine albedos. Using these data we investigate the albedo and diameter distributions of the Main Belt. As observed previously, we find a change in the average albedo when comparing the inner, middle, and outer portions of the Main Belt. We also confirm that the albedo distribution of each region is strongly bimodal. We observe groupings of objects with similar albedos in regions of the Main Belt associated with dynamical breakup families. Asteroid families typically show a characteristic albedo for all members, but there are notable exceptions to this. This paper is the first look at the Main Belt asteroids in the WISE data, and only represents the preliminary, observed raw size and albedo distributions for the populations considered. These distributions are subject to survey biases inherent to the NEOWISE dataset and cannot yet be interpreted as describing the true populations; the debiased size and albedo distributions will be the subject of the next paper in this series.

We report results based on mid-infrared photometry of comet 103P/Hartley 2 taken during May 4-13, 2010 (when the comet was at a heliocentric distance of 2.3 AU, and an observer distance of 2.0 AU) by the Wide-field Infrared Survey Explorer (Wright et al. 2010). Photometry of the coma at 22 microns and data from the University of Hawaii 2.2-m telescope obtained on May 22, 2010 provide constraints on the dust particle size distribution, dlogn/dlogm, yielding power-law slope values of alpha = -0.97 +/- 0.10, steeper than that found for the inbound particle fluence during the Stardust encounter of comet 81P/Wild 2 (Green et al. 2004). The extracted nucleus signal at 12 microns is consistent with a body of average spherical radius of 0.6 +/- 0.2 km (one standard deviation), assuming a beaming parameter of 1.2. The 4.6 micron-band signal in excess of dust and nucleus reflected and thermal contributions may be attributed to carbon monoxide or carbon dioxide emission lines and provides limits and estimates of species production. Derived carbon dioxide coma production rates are 3.5(+/- 0.9) \times 10^24 molecules per second. Analyses of the trail signal present in the stacked image with an effective exposure time of 158.4 seconds yields optical-depth values near 9 x 10^-10 at a delta mean anomaly of 0.2 deg trailing the comet nucleus, in both 12 and 22 {\mu}m bands. A minimum chi-squared analysis of the dust trail position yields a beta-parameter value of 1.0 x10^-4, consistent with a derived mean trail-grain diameter of 1.1/{\rho} cm for grains of {\rho} g/cm^3 density. This leads to a total detected trail mass of at least 4 x 10^10 {\rho} kg.

With thermal infrared observations detected by the NEOWISE project, we have measured diameters for 1742 minor planets that were also observed by the Infrared Astronomical Satellite (IRAS). We have compared the diameters and albedo derived by applying a spherical thermal model to the objects detected by NEOWISE and find that they are in generally good agreement with the IRAS values. We have shown that diameters computed from NEOWISE data are often less systematically biased than those found with IRAS. This demonstrates that the NEOWISE dataset can provide accurate physical parameters for the >157,000 minor planets that were detected by NEOWISE.

With the Wide-field Infrared Survey Explorer (WISE; Wright et al. 2010), we have observed over 157,000 minor planets (Mainzer et al. 2011). Included in these are a number of near-Earth objects, Main Belt Asteroids, and irregular satellites which have well-measured physical properties via radar, occultation and in situ imaging. We have used these objects to validate models of thermal models using the WISE measurements, as well as the color corrections derived in Wright et al. (2010) for the four WISE bandpasses as a function of effective temperature. We have used 50 objects with diameters measured by radar, occultation or in situ imaging to characterize the systematic errors implicit in using the WISE data with a faceted spherical NEATM model to compute diameters and albedos. By using the previously measured diameters and H magnitudes with a spherical NEATM model, we compute the predicted fluxes after applying the color corrections given in Wright et al. (2010) in the WISE bands and compare them to the measured magnitudes. We find minimum systematic flux errors of 5-10%, yielding minimum relative diameter and albedo errors of ~10% and ~20%, respectively. Visible albedos for the objects are computed and compared to the albedos at 3.4 and 4.6 microns, which contain a mix of reflected sunlight and thermal emission for most asteroids. We derive a linear relationship between subsolar temperature and effective temperature, which allows the color corrections given in Wright et al. (2010) to be used for asteroids by computing only subsolar temperature instead of a faceted thermal model. The thermal models derived in this paper are not intended to supplant previous measurements made using radar or spacecraft imaging; rather, we have used them to characterize the errors that should be expected when computing diameters and albedos of WISE asteroids using a spherical NEATM model.

We present orbital parameters for six double-lined spectroscopic binaries (iota Pegasi, omega Draconis, 12 Bootis, V1143 Cygni, beta Aurigae, and Mizar A) and two double-lined triple star systems (kappa Pegasi and eta Virginis). The orbital fits are based upon high-precision radial velocity observations made with a dispersed Fourier Transform Spectrograph, or dFTS, a new instrument which combines interferometric and dispersive elements. For some of the double-lined binaries with known inclination angles, the quality of our RV data permits us to determine the masses M_1 and M_2 of the stellar components with relative errors as small as 0.2%.

The \emph{Wide-field Infrared Survey Explorer} has surveyed the entire sky at four infrared wavelengths with greatly improved sensitivity and spatial resolution compared to its predecessors, the \emph{Infrared Astronomical Satellite} and the \emph{Cosmic Background Explorer}. NASA's Planetary Science Division has funded an enhancement to the \WISE\ data processing system called "NEOWISE" that allows detection and archiving of moving objects found in the \WISE\ data. NEOWISE has mined the \WISE\ images for a wide array of small bodies in our Solar System, including Near-Earth Objects (NEOs), Main Belt asteroids, comets, Trojans, and Centaurs. By the end of survey operations in February 2011, NEOWISE identified over 157,000 asteroids, including more than 500 NEOs and $\sim$120 comets. The NEOWISE dataset will enable a panoply of new scientific investigations.

The all sky surveys done by the Palomar Observatory Schmidt, the European Southern Observatory Schmidt, and the United Kingdom Schmidt, the InfraRed Astronomical Satellite and the 2 Micron All Sky Survey have proven to be extremely useful tools for astronomy with value that lasts for decades. The Wide-field Infrared Survey Explorer is mapping the whole sky following its launch on 14 December 2009. WISE began surveying the sky on 14 Jan 2010 and completed its first full coverage of the sky on July 17. The survey will continue to cover the sky a second time until the cryogen is exhausted (anticipated in November 2010). WISE is achieving 5 sigma point source sensitivities better than 0.08, 0.11, 1 and 6 mJy in unconfused regions on the ecliptic in bands centered at wavelengths of 3.4, 4.6, 12 and 22 microns. Sensitivity improves toward the ecliptic poles due to denser coverage and lower zodiacal background. The angular resolution is 6.1, 6.4, 6.5 and 12.0 arc-seconds at 3.4, 4.6, 12 and 22 microns, and the astrometric precision for high SNR sources is better than 0.15 arc-seconds.

In 1970, the Soviet Lunokhod 1 rover delivered a French-built laser reflector to the Moon. Although a few range measurements were made within three months of its landing, these measurements---and any that may have followed---are unpublished and unavailable. The Lunokhod 1 reflector was, therefore, effectively lost until March of 2010 when images from the Lunar Reconnaissance Orbiter (LRO) provided a positive identification of the rover and determined its coordinates with uncertainties of about 100 m. This allowed the Apache Point Observatory Lunar Laser-ranging Operation (APOLLO) to quickly acquire a laser signal. The reflector appears to be in excellent condition, delivering a signal roughly four times stronger than its twin reflector on the Lunokhod 2 rover. The Lunokhod 1 reflector is especially valuable for science because it is closer to the Moon's limb than any of the other reflectors and, unlike the Lunokhod 2 reflector, we find that it is usable during the lunar day. We report the selenographic position of the reflector to few-centimeter accuracy, comment on the health of the reflector, and illustrate the value of this reflector for achieving science goals.

Forty years ago, Apollo astronauts placed the first of several retroreflector arrays on the lunar surface. Their continued usefulness for laser-ranging might suggest that the lunar environment does not damage optical devices. However, new laser ranging data reveal that the efficiency of the three Apollo reflector arrays is now diminished by a factor of ten at all lunar phases and by an additional factor of ten when the lunar phase is near full moon. These deficits did not exist in the earliest years of lunar ranging, indicating that the lunar environment damages optical equipment on the timescale of decades. Dust or abrasion on the front faces of the corner-cube prisms may be responsible, reducing their reflectivity and degrading their thermal performance when exposed to face-on sunlight at full moon. These mechanisms can be tested using laboratory simulations and must be understood before designing equipment destined for the moon.

We have designed and constructed a second-generation version of the Dispersed Fourier Transform Spectrograph, or dFTS. This instrument combines a spectral interferometer with a dispersive spectrograph to provide high-accuracy, high-resolution optical spectra of stellar targets. The new version, dFTS2, is based upon the design of our prototype, with several modifications to improve the system throughput and performance. We deployed dFTS2 to the Steward Observatory 2.3-meter Bok Telescope from June 2007 to June 2008, and undertook an observing program on spectroscopic binary stars, with the goal of constraining the velocity amplitude K of the binary orbits with 0.1% accuracy, a significant improvement over most of the orbits reported in the literature. We present results for radial velocity reference stars and orbit solutions for single-lined spectroscopic binaries.

We present measurements of the Hubble diagram for 103 Type Ia supernovae (SNe) with redshifts 0.04 < z < 0.42, discovered during the first season (Fall 2005) of the Sloan Digital Sky Survey-II (SDSS-II) Supernova Survey. These data fill in the redshift "desert" between low- and high-redshift SN Ia surveys. We combine the SDSS-II measurements with new distance estimates for published SN data from the ESSENCE survey, the Supernova Legacy Survey, the Hubble Space Telescope, and a compilation of nearby SN Ia measurements. Combining the SN Hubble diagram with measurements of Baryon Acoustic Oscillations from the SDSS Luminous Red Galaxy sample and with CMB temperature anisotropy measurements from WMAP, we estimate the cosmological parameters w and Omega_M, assuming a spatially flat cosmological model (FwCDM) with constant dark energy equation of state parameter, w. For the FwCDM model and the combined sample of 288 SNe Ia, we find w = -0.76 +- 0.07(stat) +- 0.11(syst), Omega_M = 0.306 +- 0.019(stat) +- 0.023(syst) using MLCS2k2 and w = -0.96 +- 0.06(stat) +- 0.12(syst), Omega_M = 0.265 +- 0.016(stat) +- 0.025(syst) using the SALT-II fitter. We trace the discrepancy between these results to a difference in the rest-frame UV model combined with a different luminosity correction from color variations; these differences mostly affect the distance estimates for the SNLS and HST supernovae. We present detailed discussions of systematic errors for both light-curve methods and find that they both show data-model discrepancies in rest-frame $U$-band. For the SALT-II approach, we also see strong evidence for redshift-dependence of the color-luminosity parameter (beta). Restricting the analysis to the 136 SNe Ia in the Nearby+SDSS-II samples, we find much better agreement between the two analysis methods but with larger uncertainties.

We present a measurement of the rate of type Ia supernovae (SNe Ia) from the first of three seasons of data from the SDSS-II Supernova Survey. For this measurement, we include 17 SNe Ia at redshift $z\le0.12$. Assuming a flat cosmology with $\Omega_m = 0.3=1-\Omega_\Lambda$, we find a volumetric SN Ia rate of $[2.93^{+0.17}_{-0.04}({\rm systematic})^{+0.90}_{-0.71}({\rm statistical})] \times 10^{-5} {\rm SNe} {\rm Mpc}^{-3} h_{70}^3 {\rm year}^{-1}$, at a volume-weighted mean redshift of 0.09. This result is consistent with previous measurements of the SN Ia rate in a similar redshift range. The systematic errors are well controlled, resulting in the most precise measurement of the SN Ia rate in this redshift range. We use a maximum likelihood method to fit SN rate models to the SDSS-II Supernova Survey data in combination with other rate measurements, thereby constraining models for the redshift-evolution of the SN Ia rate. Fitting the combined data to a simple power-law evolution of the volumetric SN Ia rate, $r_V \propto (1+z)^{\beta}$, we obtain a value of $\beta = 1.5 \pm 0.6$, i.e. the SN Ia rate is determined to be an increasing function of redshift at the $\sim 2.5 \sigma$ level. Fitting the results to a model in which the volumetric SN rate, $r_V=A\rho(t)+B\dot \rho(t)$, where $\rho(t)$ is the stellar mass density and $\dot \rho(t)$ is the star formation rate, we find $A = (2.8 \pm 1.2) \times 10^{-14} \mathrm{SNe} \mathrm{M}_{\sun}^{-1} \mathrm{year}^{-1}$, $B = (9.3^{+3.4}_{-3.1})\times 10^{-4} \mathrm{SNe} \mathrm{M}_{\sun}^{-1}$.

This paper presents spectroscopy of supernovae discovered in the first season of the Sloan Digital Sky Survey-II Supernova Survey. This program searches for and measures multi-band light curves of supernovae in the redshift range z = 0.05 - 0.4, complementing existing surveys at lower and higher redshifts. Our goal is to better characterize the supernova population, with a particular focus on SNe Ia, improving their utility as cosmological distance indicators and as probes of dark energy. Our supernova spectroscopy program features rapid-response observations using telescopes of a range of apertures, and provides confirmation of the supernova and host-galaxy types as well as precise redshifts. We describe here the target identification and prioritization, data reduction, redshift measurement, and classification of 129 SNe Ia, 16 spectroscopically probable SNe Ia, 7 SNe Ib/c, and 11 SNe II from the first season. We also describe our efforts to measure and remove the substantial host galaxy contamination existing in the majority of our SN spectra.

We use the high magnification event seen in the 1999 OGLE campaign light curve of image C of the quadruply imaged gravitational lens Q2237+0305 to study the structure of the quasar engine. We have obtained g'- and r'-band photometry at the Apache Point Observatory 3.5m telescope where we find that the event has a smaller amplitude in the r'-band than in the g'- and OGLE V-bands. By comparing the light curves with microlensing simulations we obtain constraints on the sizes of the quasar regions contributing to the g'- and r'-band flux. Assuming that most of the surface mass density in the central kiloparsec of the lensing galaxy is due to stars and by modeling the source with a Gaussian profile, we obtain for the Gaussian width 1.20 x 10^15 sqrt(M/0.1M_sun)cm < sigma_g' < 7.96 x 10^15 sqrt(M/0.1Msun) cm, where M is the mean microlensing mass, and a ratio sigma_r'/sigma_g'=1.25^{+0.45}_{-0.15}. With the limits on the velocity of the lensing galaxy from Gil-Merino et al. (2005) as our only prior, we obtain 0.60 x 10^15 sqrt(M/0.1Msun) cm < sigma_g' < 1.57 x 10^15 sqrt(M/0.1Msun) cm and a ratio sigma_r'/sigma_g'=1.45^{+0.90}_{-0.25} (all values at 68 percent confidence). Additionally, from our microlensing simulations we find that, during the chromatic microlensing event observed, the continuum emitting region of the quasar crossed a caustic at >72 percent confidence.

The Sloan Digital Sky Survey-II Supernova Survey has identified a large number of new transient sources in a 300 sq. deg. region along the celestial equator during its first two seasons of a three-season campaign. Multi-band (ugriz) light curves were measured for most of the sources, which include solar system objects, Galactic variable stars, active galactic nuclei, supernovae (SNe), and other astronomical transients. The imaging survey is augmented by an extensive spectroscopic follow-up program to identify SNe, measure their redshifts, and study the physical conditions of the explosions and their environment through spectroscopic diagnostics. During the survey, light curves are rapidly evaluated to provide an initial photometric type of the SNe, and a selected sample of sources are targeted for spectroscopic observations. In the first two seasons, 476 sources were selected for spectroscopic observations, of which 403 were identified as SNe. For the Type Ia SNe, the main driver for the Survey, our photometric typing and targeting efficiency is 90%. Only 6% of the photometric SN Ia candidates were spectroscopically classified as non-SN Ia instead, and the remaining 4% resulted in low signal-to-noise, unclassified spectra. This paper describes the search algorithm and the software, and the real-time processing of the SDSS imaging data. We also present the details of the supernova candidate selection procedures and strategies for follow-up spectroscopic and imaging observations of the discovered sources.

The Sloan Digital Sky Survey-II (SDSS-II) has embarked on a multi-year project to identify and measure light curves for intermediate-redshift (0.05 < z < 0.35) Type Ia supernovae (SNe Ia) using repeated five-band (ugriz) imaging over an area of 300 sq. deg. The survey region is a stripe 2.5 degrees wide centered on the celestial equator in the Southern Galactic Cap that has been imaged numerous times in earlier years, enabling construction of a deep reference image for discovery of new objects. Supernova imaging observations are being acquired between 1 September and 30 November of 2005-7. During the first two seasons, each region was imaged on average every five nights. Spectroscopic follow-up observations to determine supernova type and redshift are carried out on a large number of telescopes. In its first two three-month seasons, the survey has discovered and measured light curves for 327 spectroscopically confirmed SNe Ia, 30 probable SNe Ia, 14 confirmed SNe Ib/c, 32 confirmed SNe II, plus a large number of photometrically identified SNe Ia, 94 of which have host-galaxy spectra taken so far. This paper provides an overview of the project and briefly describes the observations completed during the first two seasons of operation.

We present extensive ugrizYHJK photometry and optical spectroscopy of SN 2005gj obtained by the SDSS-II and CSP Supernova Projects, which give excellent coverage during the first 150 days after the time of explosion. These data show that SN 2005gj is the second clear case, after SN 2002ic, of a thermonuclear explosion in a dense circumstellar environment. Both the presence of singly and doubly ionized iron-peak elements (FeIII and weak SII, SiII) near maximum light as well as the spectral evolution show that SN 2002ic-like events are Type Ia explosions. Independent evidence comes from the exponential decay in luminosity of SN 2005gj, pointing to an exponential density distribution of the ejecta. The interaction of the supernova ejecta with the dense circumstellar medium is stronger than in SN 2002ic: (1) the supernova lines are weaker; (2) the Balmer emission lines are more luminous; and (3) the bolometric luminosity is higher close to maximum light. The velocity evolution of the Halpha components suggest that the CSM around SN 2005gj is clumpy and it has a flatter density distribution compared with the steady wind solution, in agreement with SN 2002ic. An early X-ray observation with Chandra gives an upper-limit on the mass loss rate from the companion of < 2x10^{-4} Msun/yr.

We have designed and constructed a ``dispersed Fourier Transform Spectrometer'' (dFTS), consisting of a conventional FTS followed by a grating spectrometer. By combining these two devices, we negate a substantial fraction of the sensitivity disadvantage of a conventional FTS for high resolution, broadband, optical spectroscopy, while preserving many of the advantages inherent to interferometric spectrometers. In addition, we have implemented a simple and inexpensive laser metrology system, which enables very precise calibration of the interferometer wavelength scale. The fusion of interferometric and dispersive technologies with a laser metrology system yields an instrument well-suited to stellar spectroscopy, velocimetry, and extrasolar planet detection, which is competitive with existing high-resolution, high accuracy stellar spectrometers. In this paper, we describe the design of our prototype dFTS, explain the algorithm we use to efficiently reconstruct a broadband spectrum from a sequence of narrowband interferograms, and present initial observations and resulting velocimetry of stellar targets.

We present 426 epochs of optical monitoring data spanning 1000 days from December 2003 to June 2006 for the gravitationally lensed quasar SDSS J1004+4112. The time delay between the A and B images is 38.4+/-2.0 days in the expected sense that B leads A and the overall time ordering is C-B-A-D-E. The measured delay invalidates all published models. The models failed because they neglected the perturbations from cluster member galaxies. Models including the galaxies can fit the data well, but strong conclusions about the cluster mass distribution should await the measurement of the longer, and less substructure sensitive, delays of the C and D images. For these images, a CB delay of 681+/-15 days is plausible but requires confirmation, while CB and AD delays of >560 days and > 800 days are required. We clearly detect microlensing of the A/B images, with the delay-corrected flux ratios changing from B-A=0.44+/-0.01 mag in the first season to 0.29+/-0.01 mag in the second season and 0.32+/-0.01 mag in the third season.

We present Follow-Up Network for Gamma-Ray Bursts (FUN GRB) Collaboration observations of the optical afterglow of GRB 021211 made between 143 seconds and 102 days after the burst. Our unique data set includes the earliest filtered detections and color information for an afterglow in the pre-Swift era. We find that the afterglow is best described by (1) propagation through a wind-swept medium, (2) a cooling break that is blueward of the observed optical frequencies, and (3) a hard electron energy distribution. However, superimposed on this "standard model" behavior we find significant chromatic variations within the first few hours after the burst. We consider possible reasons for these variations, including the possibility that they are due to a dust echo. Finally, we constrain physical parameters that describe the afterglow and surrounding medium for a variety of scenarios and find that GRB 021211's afterglow is faint for a combination of 3-4 reasons: (1) a low fraction of energy in relativistic electrons, (2) a low density for the wind-swept medium, implying either a low mass-loss rate and/or a high wind velocity for the progenitor, (3) a wide opening/viewing angle for the jet, and possibly (4) moderate source frame extinction. The jet appears to be significantly far from equipartition and magnetically dominated. More extreme versions of this might explain the darkness of many afterglows in the Swift era.

We present 39 nights of optical photometry, 34 nights of infrared photometry, and 4 nights of optical spectroscopy of the Type Ia SN 1999ac. This supernova was discovered two weeks before maximum light, and observations were begun shortly thereafter. At early times its spectra resembled the unusual SN 1999aa and were characterized by very high velocities in the Ca II H and K lines, but very low velocities in the Si II 6355 A line. The optical photometry showed a slow rise to peak brightness but, quite peculiarly, was followed by a more rapid decline from maximum. Thus, the B- and V-band light curves cannot be characterized by a single stretch factor. We argue that the best measure of the nature of this object is not the decline rate parameter Delta m_15 (B). The B-V colors were unusual from 30 to 90 days after maximum light in that they evolved to bluer values at a much slower rate than normal Type Ia supernovae. The spectra and bolometric light curve indicate that this event was similar to the spectroscopically peculiar slow decliner SN 1999aa.

We present previously unpublished BVRI photometry of the Type Ia supernovae 1999cc and 2000cf along with revised photometry of SN 1999cl. We confirm that SN 1999cl is reddened by highly non-standard dust, with R_V = 1.55 +/- 0.08. Excepting two quasar-lensing galaxies whose low values of R_V are controversial, this is the only known object with a published value of R_V less than 2.0. SNe 1999cl and 2000cf have near-infrared absolute magnitudes at maximum in good agreement with other Type Ia SNe of mid-range decline rates.

In preparation for the Supernova Survey of the Sloan Digital Sky Survey (SDSS) II, a proposed 3-year extension to the SDSS, we have conducted an early engineering and science run during the fall of 2004, which consisted of approximately 20 scheduled nights of repeated imaging of half of the southern equatorial stripe. Transient supernova-like events were detected in near real-time and photometric measurements were made in the five SDSS filter bandpasses with a cadence of ~2 days. Candidate type Ia supernovae (SNe) were pre-selected based on their colors, light curve shape, and the properties of the host galaxy. Follow-up spectroscopic observations were performed with the Astrophysical Research Consortium 3.5m telescope and the 9.2m Hobby-Eberly Telescope to confirm their types and measure the redshifts. The 2004 campaign resulted in 22 spectroscopically confirmed SNe, which includes 16 type Ia, 5 type II, and 1 type Ib/c. These SN Ia will help fill in the sparsely sampled redshift interval of z = 0.05 - 0.35, the so-called 'redshift desert', in the Hubble diagram. Detailed investigation of the spectral properties of these moderate-redshift SNe Ia will also provide a bridge between local SNe and high-redshift objects, and will help us understand the systematics for future cosmological applications that require high photometric precision. Finally, the large survey volume also provides the opportunity to select unusual supernovae for spectroscopic study that are poorly sampled in other surveys. We report on some of the early results from this program and discuss potential future applications.

The redshifts of faint radio galaxies identified with giant radio source candidates selected from the sample of Machalski et al. (2001) have been measured. Given the redshift, the projected linear size and radio luminosity are then determined. The above, supplemented with the axial ratio of the sources (evaluated from the radio maps) allows to constrain their jet power and the dynamical age using the analytical model of Kaiser et al. (1997) but modified by allowing the axial ratio of the source's cocoon to evolve in time.

We present optical and infrared photometry of the unusual Type Ia supernova 2000cx. With the data of Li et al. (2001) and Jha (2002), this comprises the largest dataset ever assembled for a Type Ia SN, more than 600 points in UBVRIJHK. We confirm the finding of Li et al. regarding the unusually blue B-V colors as SN 2000cx entered the nebular phase. Its I-band secondary hump was extremely weak given its B-band decline rate. The V minus near infrared colors likewise do not match loci based on other slowly declining Type Ia SNe, though V-K is the least ``abnormal''. In several ways SN 2000cx resembles other slow decliners, given its B-band decline rate (Delta m_15(B) = 0.93), the appearance of Fe III lines and weakness of Si II in its pre-maximum spectrum, the V-K colors and post-maximum V-H colors. If the distance modulus derived from Surface Brightness Fluctuations of the host galaxy is correct, we find that the rate of light increase prior to maximum, the characteristics of the bolometric light curve, and the implied absolute magnitude at maximum are all consistent with a sub-luminous object with Delta m_15(B) ~ 1.6-1.7 having a higher than normal kinetic energy.

The discovery of an optical counterpart to GRB010222 (detected by BeppoSAX; Piro 2001) was announced 4.4 hrs after the burst by Henden (2001a). The Sloan Digital Sky Survey's 0.5m photometric telescope (PT) and 2.5m survey telescope were used to observe the afterglow of GRB010222 starting 4.8 hours after the GRB. The 0.5m PT observed the afterglow in five, 300 sec g' band exposures over the course of half an hour, measuring a temporal decay rate in this short period of F_nu \propto t^{-1.0+/-0.5}. The 2.5m camera imaged the counterpart nearly simultaneously in five filters (u' g' r' i' z'), with r' = 18.74+/-0.02 at 12:10 UT. These multicolor observations, corrected for reddening and the afterglow's temporal decay, are well fit by the power-law F_nu \propto nu^{-0.90+/-0.03} with the exception of the u' band UV flux which is 20% below this slope. We examine possible interpretations of this spectral shape, including source extinction in a star forming region.

From an analysis of SNe 1972E, 1980N, 1981B, 1981D, 1983R, 1998bu, 1999cl, and 1999cp we find that the intrinsic V-K colors of Type Ia SNe with multi-color light curve shape (MLCS) parameter -0.4 < Delta < +0.2 suggest a uniform color curve. V-K colors become bluer linearly with time from roughly one week before B-band maximum until one week after maximum, after which they redden linearly until four weeks after maximum. V-H colors exhibit very similar color evolution. V-J colors exhibit slightly more complex evolution, with greater scatter. The existence of V minus near infrared color relations allows the construction of near infrared light curve templates that are an improvement on those of Elias et al. (1985). We provide optical BVRI and infrared JHK photometry of the Type Ia supernovae 1999aa, 1999cl, and 1999cp. SN 1999aa is an overluminous "slow decliner" (with Delta = -0.47 mag). SN 1999cp is a moderately bright SN unreddened in its host. SN 1999cl is extremely reddened in its host. The V minus near infrared colors of SN 1999cl yield A_V = 2.01 +/- 0.11 mag. This leads to a distance for its host galaxy (M 88) in agreement with other distance measurements for members of the Virgo cluster.