The Cold Classical Kuiper Belt, a class of small bodies in undisturbed orbits beyond Neptune, are primitive objects preserving information about Solar System formation. The New Horizons spacecraft flew past one of these objects, the 36 km long contact binary (486958) Arrokoth (2014 MU69), in January 2019. Images from the flyby show that Arrokoth has no detectable rings, and no satellites (larger than 180 meters diameter) within a radius of 8000 km, and has a lightly-cratered smooth surface with complex geological features, unlike those on previously visited Solar System bodies. The density of impact craters indicates the surface dates from the formation of the Solar System. The two lobes of the contact binary have closely aligned poles and equators, constraining their accretion mechanism.

We report the observations of solar system objects during the 2015 campaign of the High cadence Transient Survey (HiTS). We found 5740 bodies (mostly Main Belt asteroids), 1203 of which were detected in different nights and in $g'$ and $r'$. Objects were linked in the barycenter system and their orbital parameters were computed assuming Keplerian motion. We identified 6 near Earth objects, 1738 Main Belt asteroids and 4 Trans-Neptunian objects. We did not find a $g'-r'$ color-size correlation for $14<H_{g'}<18$ ($1<D<10$ km) asteroids. We show asteroids' colors are disturbed by HiTS' 1.6 hour cadence and estimate that observations should be separated by at most 14 minutes to avoid confusion in future wide-field surveys like LSST. The size distribution for the Main Belt objects can be characterized as a simple power law with slope $\sim0.9$, steeper than in any other survey, while data from HiTS 2014's campaign is consistent with previous ones (slopes $\sim0.68$ at the bright end and $\sim0.34$ at the faint end). This difference is likely due to the ecliptic distribution of the Main Belt since 2015's campaign surveyed farther from the ecliptic than did 2014's and most previous surveys.

We propose an ecliptic Deep Drilling Field that will discover some 10,000~small and faint Kuiper Belt Objects (KBOs) --- primitive rocky/icy bodies that orbit at the outside of our Solar System and uniquely record the processes of planetary system formation and evolution. The primary goals are to measure the KBO size and shape distributions down to 25~km, a size that probes both the early and ongoing evolution of this population. These goals can be met with around 10~hours total of on-sky time (five separate fields that are observed for 2.1~hours each). Additional science will result from downstream observations that provide colors and orbit refinement, for a total time request of 40~hours over the ten year LSST main survey.

We report on the serendipitous observations of Solar System objects imaged during the High cadence Transient Survey (HiTS) 2014 observation campaign. Data from this high cadence, wide field survey was originally analyzed for finding variable static sources using Machine Learning to select the most-likely candidates. In this work we search for moving transients consistent with Solar System objects and derive their orbital parameters. We use a simple, custom detection algorithm to link trajectories and assume Keplerian motion to derive the asteroid's orbital parameters. We use known asteroids from the Minor Planet Center (MPC) database to assess the detection efficiency of the survey and our search algorithm. Trajectories have an average of nine detections spread over 2 days, and our fit yields typical errors of $\sigma_a\sim 0.07 ~{\rm AU}$, $\sigma_{\rm e} \sim 0.07 $ and $\sigma_i\sim 0.^{\circ}5~ {\rm deg}$ in semi-major axis, eccentricity, and inclination respectively for known asteroids in our sample. We extract 7,700 orbits from our trajectories, identifying 19 near Earth objects, 6,687 asteroids, 14 Centaurs, and 15 trans-Neptunian objects. This highlights the complementarity of supernova wide field surveys for Solar System research and the significance of machine learning to clean data of false detections. It is a good example of the data--driven science that LSST will deliver.

We identify sources with extremely hard X-ray spectra (i.e., with photon indices of Gamma<0.6 in the 13 sq. deg. NuSTAR serendipitous survey, to search for the most highly obscured AGNs detected at >10 keV. Eight extreme NuSTAR sources are identified, and we use the NuSTAR data in combination with lower energy X-ray observations (from Chandra, Swift XRT, and XMM-Newton) to characterize the broad-band (0.5-24 keV) X-ray spectra. We find that all of the extreme sources are highly obscured AGNs, including three robust Compton-thick (CT; N_H > 1.5e24 cm^-2) AGNs at low redshift (z<0.1), and a likely-CT AGN at higher redshift (z=0.16). Most of the extreme sources would not have been identified as highly obscured based on the low energy (<10 keV) X-ray coverage alone. The multiwavelength properties (e.g., optical spectra and X-ray/MIR luminosity ratios) provide further support for the eight sources being significantly obscured. Correcting for absorption, the intrinsic rest-frame 10-40 keV luminosities of the extreme sources cover a broad range, from ~ 5 x 10^42 to 10^45 erg s^-1. The estimated number counts of CT AGNs in the NuSTAR serendipitous survey are in broad agreement with model expectations based on previous X-ray surveys, except for the lowest redshifts (z<0.07) where we measure a high CT fraction of f_CT^obs = 30 (+16 -12) %. For the small sample of CT AGNs, we find a high fraction of galaxy major mergers (50 +/- 33%) compared to control samples of "normal" AGNs.

We analyzed data from the first year of a survey for Near Earth Objects (NEOs) that we are carrying out with the Dark Energy Camera (DECam) on the 4-meter Blanco telescope at the Cerro Tololo Inter-American Observatory. We implanted synthetic NEOs into the data stream to derive our nightly detection efficiency as a function of magnitude and rate of motion. Using these measured efficiencies and the Solar System absolute magnitudes derived by the Minor Planet Center for the 1377 measurements of 235 unique NEOs detected, we directly derive, for the first time from a single observational data set, the NEO size distribution from 1 km down to 10 meters. We find that there are 10^6.6 NEOs larger than 10 meters. This result implies a factor of ten fewer small NEOs than some previous results, though our derived size distribution is in good agreement with several other estimates.

We report the peculiar chemical abundance patterns of eleven atypical Milky Way (MW) field red giant stars observed by the Apache Point Observatory Galactic Evolution Experiment (APOGEE). These atypical giants exhibit strong Al and N enhancements accompanied by C and Mg depletions, strikingly similar to those observed in the so-called second-generation (SG) stars of globular clusters (GCs). Remarkably, we find low-Mg abundances ([Mg/Fe]$<$0.0) together with strong Al and N overabundances in the majority (5/7) of the metal-rich ([Fe/H]$\gtrsim - 1.0$) sample stars, which is at odds with actual observations of SG stars in Galactic CGs of similar metallicities. This chemical pattern is unique and unprecedented among MW stars, posing urgent questions about its origin. These atypical stars could be former SG stars of dissolved GCs formed with intrinsically lower abundances of Mg and enriched Al (subsequently self-polluted by massive AGB stars) or the result of exotic binary systems. We speculate that the stars Mg-deficiency as well as the orbital properties suggest that they could have an extragalactic origin. This discovery should guide future dedicated spectroscopic searches of atypical stellar chemical patterns in our Galaxy; a fundamental step forward to understand the Galactic formation and evolution.

Since the first discovery of a broad-lined type Ic supernova (SN) with a long-duration gamma-ray burst (GRB) in 1998, fewer than fifty gamma-ray burst supernovae (GRB-SNe) have been discovered. The intermediate-luminosity Swift GRB 161219B and its associated supernova SN 2016jca, which occurred at a redshift of z=0.1475, represents only the seventh GRB-SN to have been discovered within 1 Gpc, and hence provides an excellent opportunity to investigate the observational and physical properties of these very elusive and rare type of SN. As such, we present optical to near-infrared photometry and optical spectroscopy of GRB 161219B and SN 2016jca, spanning the first three months since its discovery. GRB 161219B exploded in the disk of an edge-on spiral galaxy at a projected distance of 3.4 kpc from the galactic centre. GRB 161219B itself is an outlier of the Amati relation, while SN 2016jca had a rest-frame, peak absolute V-band magnitude of M_V = -19.0, which it reached after 12.5 rest-frame days. We find that the bolometric properties of SN 2016jca are inconsistent with being powered solely by a magnetar central engine, as proposed by other authors, and demonstrate that it was likely powered exclusively by energy deposited by the radioactive decay of nickel and cobalt into their daughter products, which were nucleosynthesized when its progenitor underwent core collapse. We find that 0.22 solar masses of nickel is required to reproduce the peak luminosity of SN 2016jca, and we constrain an ejecta mass of 5.8 solar masses and a kinetic energy of ~5 x 10^52 erg. Finally, we report on a chromatic, pre-maximum bump in the g-band light curve, and discuss its possible origin. [Abridged]

We present the first full catalog and science results for the NuSTAR serendipitous survey. The catalog incorporates data taken during the first 40 months of NuSTAR operation, which provide ~20Ms of effective exposure time over 331 fields, with an areal coverage of 13 sq deg, and 497 sources detected in total over the 3-24 keV energy range. There are 276 sources with spectroscopic redshifts and classifications, largely resulting from our extensive campaign of ground-based spectroscopic followup. We characterize the overall sample in terms of the X-ray, optical, and infrared source properties. The sample is primarily comprised of active galactic nuclei (AGNs), detected over a large range in redshift from z = 0.002 - 3.4 (median of <z> = 0.56), but also includes 16 spectroscopically confirmed Galactic sources. There is a large range in X-ray flux, from log( f_3-24keV / erg s^-1 cm^-2 ) ~ -14 to -11, and in rest-frame 10-40 keV luminosity, from log( L_10-40keV / erg s^-1 ) ~ 39 to 46, with a median of 44.1. Approximately 79% of the NuSTAR sources have lower energy (<10 keV) X-ray counterparts from XMM-Newton, Chandra, and Swift/XRT. The mid-infrared (MIR) analysis, using WISE all-sky survey data, shows that MIR AGN color selections miss a large fraction of the NuSTAR-selected AGN population, from ~15% at the highest luminosities (Lx > 10^44 erg s^-1) to ~80% at the lowest luminosities (Lx < 10^43 erg s^-1). Our optical spectroscopic analysis finds that the observed fraction of optically obscured AGNs (i.e., the Type 2 fraction) is F_Type2 = 53(+14-15)%, for a well-defined subset of the 8-24 keV selected sample. This is higher, albeit at a low significance level, than the Type 2 fraction measured for redshift- and luminosity-matched AGNs selected by <10 keV X-ray missions.

Hi'iaka is the larger outer satellite of the dwarf planet Haumea. Using relative photometry from the Hubble Space Telescope and Magellan and a phase dispersion minimization analysis, we have identified the rotation period of Hi'iaka to be ~9.8 hrs (double-peaked). This is ~120 times faster than its orbital period, creating new questions about the formation of this system and possible tidal evolution. The rapid rotation suggests that Hi'iaka could have a significant obliquity and spin precession that could be visible in light curves within a few years. We then turn to an investigation of what we learn about the (presently unclear) formation of the Haumea system and family based on this unexpectedly rapid rotation rate. We explore the importance of the initial semi-major axis and rotation period in tidal evolution theory and find they strongly influence the time required to despin to synchronous rotation, relevant to understanding a wide variety of satellite and binary systems. We find that despinning tides do not necessarily lead to synchronous spin periods for Hi'iaka, even if it formed near the Roche limit. Therefore the short rotation period of Hi'iaka does not rule out significant tidal evolution. Hi'iaka's spin period is also consistent with formation near its current location and spin up due to Haumea-centric impactors.

We show, by comparing observations with theoretical models, that the observed Kuiper Belt size distribution is well matched by coagulation models, which start from an initial planetesimal population with radii of about 1km, and subsequent collisional evolution. We find that the observed size distribution for R > 30km has not been modified by collisional evolution over the age of the solar system, and that the size distribution below R ~ 30km has been modified by collisions and that its slope is well matched by collisional evolution models that use published strength laws. We investigate in detail the resulting size distribution of bodies ranging from 0.01km to 30km and find that its slope changes several times as a function of radius before approaching the expected value for an equilibrium collisional cascade of material strength dominated bodies for R < 0.1km. Compared to a single power law size distribution that would span the whole range from 0.01km to 30km, we find a strong deficit of bodies around R ~10km and a strong excess of bodies around 2km. This deficit and excess are caused by an excess mass in small planetesimals in the km size range that was left over from the runaway growth phase and that leaves a signature in the size distribution that is not erased after 4.5 Gyrs of collisional evolution. Observations of the small KBO size distribution can therefore test if large KBOs grew as a result of runaway growth and constrain the initial planetesimal sizes. We find that results from recent KBO occultation surveys and the observed KBO size distribution can be best matched by an initial planetesimal population that contained about equal mass per logarithmic mass bin in bodies ranging from 0.4km to 4km in radius. In addition, we find that we cannot match the observed KBO size distribution if most of the planetesimal mass was contained in bodies that were 10km in radius or larger.

We present the discovery of a long-term stable L5 (trailing) Neptune Trojan in data acquired to search for candidate Trans-Neptunian objects for the New Horizons spacecraft to fly by during an extended post-Pluto mission. This Neptune Trojan, 2011 HM102, has the highest inclination (29.4 degrees) of any known member of this population. It is intrinsically brighter than any single L5 Jupiter Trojan at H~8.18. We have determined its gri colors (a first for any L5 Neptune Trojan), which we find to be similar to the moderately red colors of the L4 Neptune Trojans, suggesting similar surface properties for members of both Trojan clouds. We also present colors derived from archival data for two L4 Neptune Trojans (2006 RJ103 and 2007 VL305), better refining the overall color distribution of the population. In this document we describe the discovery circumstances, our physical characterization of 2011 HM102, and this object's implications for the Neptune Trojan population overall. Finally, we discuss the prospects for detecting 2011 HM102 from the New Horizons spacecraft during their close approach in mid- to late-2013.

We present the faintest mid ecliptic latitude survey in the second part of HST archival search for outer Solar System bodies. We report the discovery of 28 new trans-Neptunian objects and 1 small centaur (R ~ 2km) in the band 5{\circ} - 20{\circ} off the ecliptic. The inclination distribution of these excited ob jects is consistent with the distribution derived from brighter ecliptic surveys. We suggest that the size and inclination distribution should be estimated consistently using suitable surveys with calibrated search algorithms and reliable orbital information.

We present photometry of 4 transits of the exoplanet WASP-4b, each with a precision of approximately 500 ppm and a time sampling of 40-60s. We have used the data to refine the estimates of the system parameters and ephemerides. During two of the transits we observed a short-lived, low-amplitude anomaly that we interpret as the occultation of a starspot by the planet. We also find evidence for a pair of similar anomalies in previously published photometry. The recurrence of these anomalies suggests that the stellar rotation axis is nearly aligned with the orbital axis, or else the star spot would not have remained on the transit chord. By analyzing the timings of the anomalies we find the sky-projected stellar obliquity to be -1_{-12}^{+14} degrees. This result is consistent with (and more constraining than) a recent observation of the Rossiter-McLaughlin effect. It suggests that the planet migration mechanism preserved the initially low obliquity, or else that tidal evolution has realigned the system. Future applications of this method using data from the Corot and Kepler missions will allow spin-orbit alignment to be probed for many other exoplanets.

Bernstein et al. (2004) found that the population of faint (R>26) trans-Neptunian objects (TNOs) known at that time was dominated by "Classical" objects, which have low inclinations (i<5 degrees) and distances 40--45 AU. Since those observations, the number of faint TNOs whose orbits are sufficiently well known to be classified as "Classical" or "Excited" has grown from seven to 39. We analyze the dynamical classifications of faint TNOs known today and find that this population is dominated by Excited objects. We discuss some implications of this result.

We introduce a novel search technique that can identify trans-neptunian objects in three to five exposures of a pointing within a single Hubble Space Telescope orbit. The process is fast enough to allow the discovery of candidates soon after the data are available. This allows sufficient time to schedule follow up observations with HST within a month. We report the discovery of 14 slow-moving objects found within 5\circ of the ecliptic in archival data taken with the Wide Field Channel of the Advanced Camera for Surveys. The luminosity function of these objects is consistent with previous ground-based and space-based results. We show evidence that the size distribution of both high and low inclination populations is similar for objects smaller than 100 km, as expected from collisional evolution models, while their size distribution differ for brighter objects. We suggest the two populations formed in different parts of the protoplanetary disk and after being dynamically mixed have collisionally evolved together. Among the objects discovered there is an equal mass binary with an angular separation ~ 0."53.

We present the results of an archival search for Trans-neptunian objects (TNOs) in an ecliptic field observed with Subaru in 2002. The depth of the search allowed us to find 20 new TNOs with magnitudes between R=24 and 27. We fit a double power law model to the data; the most likely values for the bright and faint power law exponents are alpha_1=0.73_{-0.09}^{+0.08} and alpha_2=0.20_{-0.14}^{+0.12}; the differential number density at R=23 is sigma_{23}=1.46_{-0.12}^{+0.14} and the break magnitude is R_{eq}=25.0_{-0.6}^{+0.8}. This is the most precise measurement of the break in the TNO luminosity function to date. The break in the size distribution corresponds to a diameter of D = 90+-30 km assuming a 4% albedo.

We present the results of a survey for trans-neptunian objects (TNOs) based on Subaru archival images, originally collected by Sheppard et al. (2005) as part of a search for irregular satellites of Uranus. The survey region covers 2.8 deg^2, centered on Uranus and observed near opposition on two adjacent nights. Our survey reaches half its maximum detection efficiency at R=25.69$\pm$0.01. The objects detected correspond to 82 TNOs, five Centaurs, and five irregular satellites. We model the cumulative number of TNOs brighter than a given apparent magnitude with both a single and double power law. The best fit single power law, with one object per square degree at magnitude $R_0 =22.6_{-0.4}^{+0.3}$ and a slope of $\alpha = 0.51_{-0.6}^{+0.5}$, is inconsistent with the results of similar searches with shallower limiting magnitudes. The best fit double power law, with a bright-end slope $\alpha_1 = 0.7_{-0.1}^{+0.2}$, a faint-end slope $\alpha_2=0.3_{-0.2}^{+0.2}$, a differential number density at R=23, $\sigma_{23} = 2.0_{-0.5}^{+0.5}$ and a magnitude break in the slope at $R_{eq} = 24.3_{-0.1}^{+0.8}$, is more likely than the single power law by a Bayes factor of ~26. This is the first survey with sufficient depth and areal coverage to identify the magnitude at which the break occurs without relying on the results of other surveys. ...

We present I and B photometry of five distinct transits of the exoplanet OGLE-TR-10b. By modeling the light curves, we find the planetary radius to be R_P = 1.06 +/- 0.08 R_Jup and the stellar radius to be R_S = 1.10 +/- 0.07 R_sun. The uncertainties are dominated by statistical errors in the photometry. Our estimate of the planetary radius is smaller than previous estimates that were based on lower-precision photometry, and hence the planet is not as anomalously large as was previously thought. We provide updated determinations of all the system parameters, including the transit ephemerides.

As part of our ongoing effort to measure exoplanet sizes and transit times with greater accuracy, we present I band observations of two transits of OGLE-TR-111b. The photometry has an accuracy of 0.15-0.20% and a cadence of 1-2 minutes. We derive a planetary radius of 1.067 +/- 0.054 Jupiter radii and a stellar radius of 0.831 +/- 0.031 solar radii. The uncertainties are dominated by errors in the photometry, rather than by systematic errors arising from uncertainties in the limb darkening function or the stellar mass. Both the stellar radius and the planetary radius are in agreement with theoretical expectations. The transit times are accurate to within 30 seconds, and allow us to refine the estimate of the mean orbital period: 4.0144479 +/- 0.0000041 days.

We present high-precision photometry of the hypervelocity star SDSS J090745.0+024507 (HVS), which has a Galactic rest-frame radial velocity of v=709 km/s, and so has likely been ejected from the supermassive black hole in the Galactic center. Our data were obtained on two nights using the MMT 6.5m telescope, and is supplemented by lower precision photometry obtained on four nights using the FLWO 1.2m telescope. The high-precision photometry indicates that the HVS is a short-period, low-amplitude variable, with period P=0.2-2 days and amplitude A = 2-10%. Together with the known effective temperature of T_eff ~ 10,500 K (spectral type B9), this variability implies that the HVS is a member of the class of slowly pulsating B-type main sequence stars, thus resolving the previously-reported two-fold degeneracy in the luminosity and distance of the star. The HVS has a heliocentric distance of 71 kpc, and an age of ~0.35 Gyr. The time of ejection from the center of the Galaxy is < 100 Myr, and thus the existence of the OS constitutes observational evidence of a population of young stars in the proximity of the central supermassive black hole ~0.1 Gyr ago. It is possible that the HVS was a member of a binary that was tidally disrupted by the central black hole; we discuss constraints on the properties of the companion's orbit.