Showing 1–50 of 89 results for author: Jedicke, R

NEOMOD 3: The Debiased Size Distribution of Near Earth Objects

Authors: David Nesvorny, David Vokrouhlicky, Frank Shelly, Rogerio Deienno, William F. Bottke, Carson Fuls, Robert Jedicke, Shantanu Naidu, Steven R. Chesley, Paul W. Chodas, Davide Farnocchia, Marco Delbo

Abstract: Our previous model (NEOMOD2) for the orbital and absolute magnitude distribution of Near Earth Objects (NEOs) was calibrated on the Catalina Sky Survey observations between 2013 and 2022. Here we extend NEOMOD2 to include visible albedo information from the Wide-Field Infrared Survey Explorer. The debiased albedo distribution of NEOs can be approximated by the sum of two Rayleigh distributions wit… ▽ More Our previous model (NEOMOD2) for the orbital and absolute magnitude distribution of Near Earth Objects (NEOs) was calibrated on the Catalina Sky Survey observations between 2013 and 2022. Here we extend NEOMOD2 to include visible albedo information from the Wide-Field Infrared Survey Explorer. The debiased albedo distribution of NEOs can be approximated by the sum of two Rayleigh distributions with the scale parameters p_V,dark=0.03 and p_V,bright=0.17. We find evidence for smaller NEOs having (on average) higher albedos than larger NEOs; this is likely a consequence of the size-dependent sampling of different main belt sources. These inferences and the absolute magnitude distribution from NEOMOD2 are used to construct the debiased size distribution of NEOs. We estimate 830+/-60 NEOs with diameters D>1 km and 20,000+/-2,000 NEOs with D>140 m. The new model, NEOMOD3, is available via the NEOMOD Simulator -- an easy-to-operate code that can be used to generate user-defined samples (orbits, sizes and albedos) from the model. △ Less

Submitted 29 April, 2024; originally announced April 2024.

Comments: Icarus, in press

Numerical behavior of the Keplerian Integral methods for initial orbit determination

Authors: Óscar Rodríguez, Giovanni F. Gronchi, Giulio Baù, Robert Jedicke

Abstract: We investigate the behaviour of two recent methods for the computation of preliminary orbits. These methods are based on the conservation laws of Kepler's problem, and enable the linkage of very short arcs of optical observations even when they are separated in time by a few years. Our analysis is performed using both synthetic and real data of 822 main belt asteroids. The differences between comp… ▽ More We investigate the behaviour of two recent methods for the computation of preliminary orbits. These methods are based on the conservation laws of Kepler's problem, and enable the linkage of very short arcs of optical observations even when they are separated in time by a few years. Our analysis is performed using both synthetic and real data of 822 main belt asteroids. The differences between computed and true orbital elements have been analysed for the true linkages, as well as the occurrence of alternative solutions. Some metrics have been introduced to quantify the results, with the aim of discarding as many of the false linkages as possible and keeping the vast majority of true ones. These numerical experiments provide thresholds for the metrics which take advantage of the knowledge of the \emph{ground truth}: the values of these thresholds can be used in normal operation mode, when we do not know the correct values of the orbital elements and whether the linkages are true or false. △ Less

Submitted 2 January, 2024; originally announced January 2024.

Linking tracklets over the years in large datasets

Authors: Óscar Rodríguez, Giovanni F. Gronchi, Giulio Baù, Robert Jedicke

Abstract: We present a new procedure to identify observations of known objects in large data sets of unlinked detections. It begins with a Keplerian integrals method that allows us to link two tracklets, computing preliminary orbits, even when the tracklets are separated in time by a few years. In the second step, we represent the results in a `graph' where the tracklets are the nodes and the preliminary or… ▽ More We present a new procedure to identify observations of known objects in large data sets of unlinked detections. It begins with a Keplerian integrals method that allows us to link two tracklets, computing preliminary orbits, even when the tracklets are separated in time by a few years. In the second step, we represent the results in a `graph' where the tracklets are the nodes and the preliminary orbits are the edges. Then, acceptable `3-cycles' are identified and a least squares orbit is computed for each of them. Finally, we construct sequences of $n \geq 4$ tracklets by searching through the orbits of nearby 3-cycles and attempting to attribute the remaining tracklets. We calculate the technique's efficiency at identifying unknown objects using real detections that attempt to mimic key parameters of the Minor Planet Center's Isolated Tracklet File (ITF) and then apply the procedure to the ITF to identify tens of thousands of new objects. △ Less

Submitted 2 January, 2024; originally announced January 2024.

NEOMOD 2: An Updated Model of Near-Earth Objects from a Decade of Catalina Sky Survey Observations

Authors: David Nesvorny, David Vokrouhlicky, Frank Shelly, Rogerio Deienno, William F. Bottke, Eric Christensen, Robert Jedicke, Shantanu Naidu, Steven R. Chesley, Paul W. Chodas, Davide Farnocchia, Mikael Granvik

Abstract: Catalina Sky Survey (CSS) is a major survey of Near-Earth Objects (NEOs). In a recent work, we used CSS observations from 2005-2012 to develop a new population model of NEOs (NEOMOD). CSS's G96 telescope was upgraded in 2016 and detected over 10,000 unique NEOs since then. Here we characterize the NEO detection efficiency of G96 and use G96's NEO detections from 2013-2022 to update NEOMOD. This re… ▽ More Catalina Sky Survey (CSS) is a major survey of Near-Earth Objects (NEOs). In a recent work, we used CSS observations from 2005-2012 to develop a new population model of NEOs (NEOMOD). CSS's G96 telescope was upgraded in 2016 and detected over 10,000 unique NEOs since then. Here we characterize the NEO detection efficiency of G96 and use G96's NEO detections from 2013-2022 to update NEOMOD. This resolves previous model inconsistencies related to the population of large NEOs. We estimate there are 936+/-29 NEOs with absolute magnitude H<17.75 (diameter D>1 km for the reference albedo p_V=0.14). The slope of the NEO size distribution for H=25-28 is found to be relatively shallow (cumulative index 2.6) and the number of H<28 NEOs (D>9 m) is determined to be (1.20+/-0.04)x10^7. Small NEOs have a different orbital distribution and higher impact probabilities than large NEOs. We estimate 0.034+/-0.002 impacts of H<28 NEOs on the Earth per year, which is near the low end of the impact flux range inferred from atmospheric bolide observations. Relative to a model where all NEOs are delivered directly from the main belt, the population of small NEOs detected by G96 shows an excess of low-eccentricity orbits with a=1--1.6 au that appears to increase with H. We suggest that the population of very small NEOs is boosted by tidal disruption of large NEOs during close encounters to the terrestrial planets. When the effect of tidal disruption is (approximately) accounted for in the model, we estimate 0.06+/-0.01 impacts of H<28 NEOs on the Earth per year, which is more in line with the bolide data. △ Less

Submitted 14 December, 2023; originally announced December 2023.

Comments: Icarus, in press. arXiv admin note: text overlap with arXiv:2306.09521

Rotation periods and colours of 10-m scale near-Earth asteroids from CFHT target of opportunity streak photometry

Authors: B. T. Bolin, M. Ghosal, R. Jedicke

Abstract: The rotational properties of $\sim$10~m-scale asteroids are poorly understood with only a few measurements. Additionally, collisions or thermal recoil can spin their rotations to periods less than a few seconds obfuscating their study due to the observational cadence imposed by the long read-out times of charge-coupled device imagers. We present a method to measure the rotation periods of 10~m-sca… ▽ More The rotational properties of $\sim$10~m-scale asteroids are poorly understood with only a few measurements. Additionally, collisions or thermal recoil can spin their rotations to periods less than a few seconds obfuscating their study due to the observational cadence imposed by the long read-out times of charge-coupled device imagers. We present a method to measure the rotation periods of 10~m-scale asteroids using the target of opportunity capability of the Canada France Hawaii Telescope and its MegaCam imager by intentionally streaking their detections in single exposures when they are at their brightest. Periodic changes in brightness as small as $\sim$0.05 mag along the streak can be measured as short as a few seconds. Additionally, the streak photometry is taken in multiple g, r, and i filter exposures enabling the measurement of asteroid colours. The streak photometry method was tested on CFHT observations of three 10~m-scale asteroids, 2016 GE$_1$, 2016 CG$_{18}$, and 2016 EV$_{84}$. Our 3 targets are among the smallest known asteroids with measured rotation periods/colours having some of the shortest known rotation periods. We compare our rotation period and taxonomic results with independent data from the literature and discuss applications of the method to future small asteroid observations. △ Less

Submitted 13 October, 2023; originally announced October 2023.

Comments: Revised version, MNRAS:L, 13 pages, 10 figures, 3 tables

Interstellar Comets from Post-Main Sequence Systems as Tracers of Extrasolar Oort Clouds

Authors: W. Garrett Levine, Aster G. Taylor, Darryl Z. Seligman, Devin J. Hoover, Robert Jedicke, Jennifer B. Bergner, Gregory P. Laughlin

Abstract: Interstellar small bodies are unique probes into the histories of exoplanetary systems. One hypothesized class of interlopers are "Jurads," exo-comets released into the Milky Way during the post-main sequence as the thermally-pulsing asymptotic giant branch (AGB) host stars lose mass. In this study, we assess the prospects for the Legacy Survey of Space and Time (LSST) to detect a Jurad and examin… ▽ More Interstellar small bodies are unique probes into the histories of exoplanetary systems. One hypothesized class of interlopers are "Jurads," exo-comets released into the Milky Way during the post-main sequence as the thermally-pulsing asymptotic giant branch (AGB) host stars lose mass. In this study, we assess the prospects for the Legacy Survey of Space and Time (LSST) to detect a Jurad and examine whether such an interloper would be observationally distinguishable from exo-comets ejected during the (pre-)main sequence. Using analytic and numerical methods, we estimate the fraction of exo-Oort Cloud objects that are released from 1-8 solar mass stars during post-main sequence evolution. We quantify the extent to which small bodies are altered by the increased luminosity and stellar outflows during the AGB, finding that some Jurads may lack hypervolatiles and that stellar winds could deposit dust that covers the entire exo-comet surface. Next, we construct models of the interstellar small body reservoir for various size-frequency distribution slopes, characteristic sizes, and the total mass sequestered in the minor planets of exo-Oort Clouds. Even with the LSST's increased search volume compared to contemporary surveys, we find that detecting a Jurad is unlikely but not infeasible given the current understanding of (exo)planet formation. △ Less

Submitted 21 June, 2023; originally announced June 2023.

Comments: 28 pages, 13 figures; accepted to PSJ

NEOMOD: A New Orbital Distribution Model for Near Earth Objects

Authors: David Nesvorny, Rogerio Deienno, William F. Bottke, Robert Jedicke, Shantanu Naidu, Steven R. Chesley, Paul W. Chodas, Mikael Granvik, David Vokrouhlicky, Miroslav Broz, Alessandro Morbidelli, Eric Christensen, Bryce T. Bolin

Abstract: Near Earth Objects (NEOs) are a transient population of small bodies with orbits near or in the terrestrial planet region. They represent a mid-stage in the dynamical cycle of asteroids and comets, which starts with their removal from the respective source regions -- the main belt and trans-Neptunian scattered disk -- and ends as bodies impact planets, disintegrate near the Sun, or are ejected fro… ▽ More Near Earth Objects (NEOs) are a transient population of small bodies with orbits near or in the terrestrial planet region. They represent a mid-stage in the dynamical cycle of asteroids and comets, which starts with their removal from the respective source regions -- the main belt and trans-Neptunian scattered disk -- and ends as bodies impact planets, disintegrate near the Sun, or are ejected from the Solar System. Here we develop a new orbital model of NEOs by numerically integrating asteroid orbits from main belt sources and calibrating the results on observations of the Catalina Sky Survey. The results imply a size-dependent sampling of the main belt with the $ν_6$ and 3:1 resonances producing $\simeq 30$\% of NEOs with absolute magnitudes $H = 15$ and $\simeq 80$\% of NEOs with $H = 25$. Hence, the large and small NEOs have different orbital distributions. The inferred flux of $H<18$ bodies into the 3:1 resonance can be sustained only if the main-belt asteroids near the resonance drift toward the resonance at the maximal Yarkovsky rate ($\simeq 2 \times 10^{-4}$ au Myr$^{-1}$ for diameter $D=1$ km and semimajor axis $a=2.5$~au). This implies obliquities $θ\simeq 0^\circ$ for $a<2.5$~au and $θ\simeq 180^\circ$ for $a>2.5$~au, both in the immediate neighborhood of the resonance (the same applies to other resonances as well). We confirm the size-dependent disruption of asteroids near the Sun found in previous studies. An interested researcher can use the publicly available NEOMOD Simulator to generate user-defined samples of NEOs from our model. △ Less

Submitted 15 June, 2023; originally announced June 2023.

Comments: AJ

Mineralogical Characterization and Phase Angle Study of Two Binary Near-Earth Asteroids, Potential Targets for NASA's Janus Mission

Authors: Lucille Le Corre, Juan A. Sanchez, Vishnu Reddy, Adam Battle, David Cantillo, Benjamin Sharkey, Robert Jedicke, Daniel Scheeres

Abstract: Ground-based characterization of spacecraft targets prior to mission operations is critical to properly plan and execute measurements. Understanding surface properties, like mineralogical composition and phase curves (expected brightness at different viewing geometries) informs data acquisition during the flybys. Binary near-Earth asteroids (NEA) (35107) 1991 VH and (175706) 1996 FG3 were selected… ▽ More Ground-based characterization of spacecraft targets prior to mission operations is critical to properly plan and execute measurements. Understanding surface properties, like mineralogical composition and phase curves (expected brightness at different viewing geometries) informs data acquisition during the flybys. Binary near-Earth asteroids (NEA) (35107) 1991 VH and (175706) 1996 FG3 were selected as potential targets of the National Aeronautics and Space Administration's (NASA) dual spacecraft Janus mission. We observed 1991 VH using the 3-m NASA Infrared Telescope Facility (IRTF) on Mauna Kea, Hawaii, on July 26, 2008. 1996 FG3 was observed with the IRTF for seven nights during the spring of 2022. Compositional analysis of 1991 VH revealed that this NEA is classified as an Sq-type in the Bus-DeMeo taxonomy classification, with a composition consistent with LL ordinary chondrites. Using thermal modeling, we computed the thermally corrected spectra for 1996 FG3 and the corresponding best fit albedo of about 2-3% for the best spectra averaged for each night. Our spectral analysis indicates that this NEA is a Ch-type. The best possible meteorite analogs for 1996 FG3, based on curve matching, are two carbonaceous chondrites, Y-86789 and Murchison. No rotational variation was detected in the spectra of 1996 FG3, which means there may not be any heterogeneities on the surface of the primary. However, a clear phase reddening effect was observed in our data, confirming findings from previous ground-based studies. △ Less

Submitted 26 April, 2023; originally announced April 2023.

Comments: 18 pages, 8 figures, 1 table, accepted for publication in the Planetary Science Journal

Journal ref: Planet. Sci. J. 4 91 (2023)

Shape-Driven Selection Effects for Aspherical Near-Earth Objects in Systematic Surveys

Authors: W. Garrett Levine, Robert Jedicke

Abstract: The apparent magnitude of elongated small bodies is time-dependent over their rotation phase. Therefore, previously undiscovered aspherical minor planets may experience a shape-driven selection effect in systematic surveys versus their spherical counterparts. In this study, we conduct injection-recovery exercises of synthetic asteroid lightcurves using a simple model to quantify the effect of vary… ▽ More The apparent magnitude of elongated small bodies is time-dependent over their rotation phase. Therefore, previously undiscovered aspherical minor planets may experience a shape-driven selection effect in systematic surveys versus their spherical counterparts. In this study, we conduct injection-recovery exercises of synthetic asteroid lightcurves using a simple model to quantify the effect of varying axial ratio on detection efficiencies. We find that high-amplitude lightcurves are confronted with adverse selection effects for survey cadences and discovery thresholds for constructing tracklets that are representative of modern and proposed future NEO searches. Furthermore, we illustrate the possible hazards of drawing population-level inferences on an underlying reservoir of elongated small bodies. If physical size and characteristic axial ratios are correlated, then size-frequency distributions may require revision at small diameters. In particular, this effect could alter the estimated populations of near-Earth objects. We conclude by discussing the applicability of our results to various other classes of solar system minor planets and interstellar interlopers, as well as discuss future work that may further interrogate this detection bias. △ Less

Submitted 6 March, 2023; v1 submitted 2 March, 2023; originally announced March 2023.

Comments: 21 pages, 10 figure, 1 table. Accepted for publication in Icarus. (small typo fixed in v2)

R2-D2: Roman and Rubin -- From Data to Discovery

Authors: Suvi Gezari, Misty Bentz, Kishalay De, K. Decker French, Aaron Meisner, Michelle Ntampaka, Robert Jedicke, Ekta Patel, Daniel Perley, Robyn Sanderson, Christian Aganze, Igor Andreoni, Eric F. Bell, Edo Berger, Ian Dell'Antonio, Ryan Foley, Henry Hsieh, Mansi Kasliwal, Joel Kastner, Charles D. Kilpatrick, J. Davy Kirkpatrick, Casey Lam, Karen Meech, Dante Minniti, Ethan O. Nadler , et al. (6 additional authors not shown)

Abstract: The NASA Nancy Grace Roman Space Telescope (Roman) and the Vera C. Rubin Observatory Legacy Survey of Space and Time (Rubin), will transform our view of the wide-field sky, with similar sensitivities, but complementary in wavelength, spatial resolution, and time domain coverage. Here we present findings from the AURA Roman+Rubin Synergy Working group, charged by the STScI and NOIRLab Directors to… ▽ More The NASA Nancy Grace Roman Space Telescope (Roman) and the Vera C. Rubin Observatory Legacy Survey of Space and Time (Rubin), will transform our view of the wide-field sky, with similar sensitivities, but complementary in wavelength, spatial resolution, and time domain coverage. Here we present findings from the AURA Roman+Rubin Synergy Working group, charged by the STScI and NOIRLab Directors to identify frontier science questions in General Astrophysics, beyond the well-covered areas of Dark Energy and Cosmology, that can be uniquely addressed with Roman and Rubin synergies in observing strategy, data products and archiving, joint analysis, and community engagement. This analysis was conducted with input from the community in the form of brief (1-2 paragraph) "science pitches" (see Appendix), and testimony from "outside experts" (included as co-authors). We identify a rich and broad landscape of potential discoveries catalyzed by the combination of exceptional quality and quantity of Roman and Rubin data, and summarize implementation requirements that would facilitate this bounty of additional science with coordination of survey fields, joint coverage of the Galactic plane, bulge, and ecliptic, expansion of General Investigator and Target of Opportunity observing modes, co-location of Roman and Rubin data, and timely distribution of data, transient alerts, catalogs, value-added joint analysis products, and simulations to the broad astronomical community. △ Less

Submitted 24 February, 2022; originally announced February 2022.

Comments: 29 pages, 12 figures, Table of Implementation Recommendations, Appendix of Community Science Pitches, AURA-commissioned whitepaper submitted to the Director of STScI (Ken Sembach) and the Director of NOIRLab (Pat McCarthy)

A Sublime Opportunity: The Dynamics of Transitioning Cometary Bodies and the Feasibility of $\textit{In Situ}$ Observations of The Evolution of Their Activity

Authors: Darryl Z. Seligman, Kaitlin M. Kratter, W. Garrett Levine, Robert Jedicke

Abstract: The compositional and morphological evolution of minor bodies in the Solar System is primarily driven by the evolution of their heliocentric distances, as the level of incident solar radiation regulates cometary activity. We investigate the dynamical transfer of Centaurs into the inner Solar System, facilitated by mean motion resonances with Jupiter and Saturn. The recently discovered object, P/20… ▽ More The compositional and morphological evolution of minor bodies in the Solar System is primarily driven by the evolution of their heliocentric distances, as the level of incident solar radiation regulates cometary activity. We investigate the dynamical transfer of Centaurs into the inner Solar System, facilitated by mean motion resonances with Jupiter and Saturn. The recently discovered object, P/2019 LD2, will transition from the Centaur region to the inner Solar System in 2063. In order to contextualize LD2, we perform N-body simulations of a population of Centaurs and JFCs. Objects between Jupiter and Saturn with Tisserand parameter $T_J\sim$3 are transferred onto orbits with perihelia $q<4$au within the next 1000 years with notably high efficiency. Our simulations show that there may be additional LD2-like objects transitioning into the inner Solar System in the near-term future, all of which have low $Δ$V with respect to Jupiter. We calculate the distribution of orbital elements resulting from a single Jovian encounter and show that objects with initial perihelia close to Jupiter are efficiently scattered to $q<4$au. Moreover, approximately $55\%$ of the transitioning objects in our simulated population experience at least 1 Jovian encounter prior to reaching $q<4$au. We demonstrate that a spacecraft stationed near Jupiter would be well-positioned to rendezvous, orbit match, and accompany LD2 into the inner Solar System, providing an opportunity to observe the onset of intense activity in a pristine comet $\textit{in situ}$. Finally, we discuss the prospect of identifying additional targets for similar measurements with forthcoming observational facilities. △ Less

Submitted 6 October, 2021; originally announced October 2021.

Comments: 23 pages, 17 figures, accepted for publication at PSJ

Generalization of a method by Mossotti for initial orbit determination

Authors: Giovanni F. Gronchi, Giulio Baù, Óscar Rodríguez, Robert Jedicke, Joachim Moeyens

Abstract: Here we revisit an initial orbit determination method introduced by O. F. Mossotti employing four geocentric sky-plane observations and a linear equation to compute the angular momentum of the observed body. We then extend the method to topocentric observations, yielding a quadratic equation for the angular momentum. The performance of the two versions are compared through numerical tests with syn… ▽ More Here we revisit an initial orbit determination method introduced by O. F. Mossotti employing four geocentric sky-plane observations and a linear equation to compute the angular momentum of the observed body. We then extend the method to topocentric observations, yielding a quadratic equation for the angular momentum. The performance of the two versions are compared through numerical tests with synthetic asteroid data using different time intervals between consecutive observations and different astrometric errors. We also show a comparison test with Gauss's method using simulated observations with the expected cadence of the VRO-LSST telescope. △ Less

Submitted 1 September, 2021; v1 submitted 1 April, 2021; originally announced April 2021.

Comments: 22 pages, 9 figures

MSC Class: 70-XX

Establishing Earth's Minimoon Population through Characterization of Asteroid 2020 CD$_3$

Authors: Grigori Fedorets, Marco Micheli, Robert Jedicke, Shantanu P. Naidu, Davide Farnocchia, Mikael Granvik, Nicholas Moskovitz, Megan E. Schwamb, Robert Weryk, Kacper Wierzchoś, Eric Christensen, Theodore Pruyne, William F. Bottke, Quanzhi Ye, Richard Wainscoat, Maxime Devogèle, Laura E. Buchanan, Anlaug Amanda Djupvik, Daniel M. Faes, Dora Föhring, Joel Roediger, Tom Seccull, Adam B. Smith

Abstract: We report on our detailed characterization of Earth's second known temporary natural satellite, or minimoon, asteroid 2020 CD3. An artificial origin can be ruled out based on its area-to-mass ratio and broadband photometry, which suggest that it is a silicate asteroid belonging to the S or V complex in asteroid taxonomy. The discovery of 2020 CD3 allows for the first time a comparison between know… ▽ More We report on our detailed characterization of Earth's second known temporary natural satellite, or minimoon, asteroid 2020 CD3. An artificial origin can be ruled out based on its area-to-mass ratio and broadband photometry, which suggest that it is a silicate asteroid belonging to the S or V complex in asteroid taxonomy. The discovery of 2020 CD3 allows for the first time a comparison between known minimoons and theoretical models of their expected physical and dynamical properties. The estimated diameter of 1.2+0.4-0.2 m and geocentric capture approximately a decade after the first known minimoon, 2006 RH120, are in agreement with theoretical predictions. The capture duration of 2020 CD3 of at least 2.7 yr is unexpectedly long compared to the simulation average, but it is in agreement with simulated minimoons that have close lunar encounters, providing additional support for the orbital models. 2020 CD3's atypical rotation period, significantly longer than theoretical predictions, suggests that our understanding of meter-scale asteroids needs revision. More discoveries and a detailed characterization of the population can be expected with the forthcoming Vera C. Rubin Observatory Legacy Survey of Space and Time. △ Less

Submitted 20 November, 2020; originally announced November 2020.

Comments: 22 pages, 6 Figures, 5 Tables, to appear in the Astronomical Journal

Debiased albedo distribution for Near Earth Objects

Authors: Alessandro Morbidelli, Marco Delbo, Mikael Granvik, William F. Bottke, Robert Jedicke, Bryce Bolin, Patrick Michel, David Vokrouhlický

Abstract: We extend the most recent orbital and absolute magnitude Near Earth Object (NEO) model (Granvik et al., 2018) to provide a statistical description of NEO geometric albedos. Our model is calibrated on NEOWISE albedo data for the NEO population and reproduces these data very well once a simple model for the NEOWISE observational biases is applied. The results are consistent with previous estimates.… ▽ More We extend the most recent orbital and absolute magnitude Near Earth Object (NEO) model (Granvik et al., 2018) to provide a statistical description of NEO geometric albedos. Our model is calibrated on NEOWISE albedo data for the NEO population and reproduces these data very well once a simple model for the NEOWISE observational biases is applied. The results are consistent with previous estimates. There are about 1,000 NEOs with diameter D>1km and the mean albedo to convert absolute magnitude into diameter is 0.147. We do not find any statistically significant evidence that the albedo distribution of NEOs depends on NEO size. Instead, we find evidence that the disruption of NEOs at small perihelion distances found in Granvik et al. (2016) occurs preferentially for dark NEOs. The interval between km-sized bodies striking the Earth should occur on average once every 750,000 years. Low and high albedo NEOs are well mixed in orbital space, but a trend remains with higher albedo objects being at smaller semimajor axes and lower albedo objects more likely found at larger semimajor axes. △ Less

Submitted 10 January, 2020; originally announced January 2020.

Comments: In press in Icarus

Discovering Earth's transient moons with the Large Synoptic Survey Telescope

Authors: Grigori Fedorets, Mikael Granvik, R. Lynne Jones, Mario Jurić, Robert Jedicke

Abstract: Earth's temporarily-captured orbiters (TCOs) are a sub-population of near-Earth objects (NEOs). TCOs can provide constraints for NEO population models in the 1--10-metre-diameter range, and they are outstanding targets for in situ exploration of asteroids due to a low requirement on $Δv$. So far there has only been a single serendipitous discovery of a TCO. Here we assess in detail the possibility… ▽ More Earth's temporarily-captured orbiters (TCOs) are a sub-population of near-Earth objects (NEOs). TCOs can provide constraints for NEO population models in the 1--10-metre-diameter range, and they are outstanding targets for in situ exploration of asteroids due to a low requirement on $Δv$. So far there has only been a single serendipitous discovery of a TCO. Here we assess in detail the possibility of their discovery with the upcoming Large Synoptic Survey Telescope (LSST), previously identified as the primary facility for such discoveries. We simulated observations of TCOs by combining a synthetic TCO population with an LSST survey simulation. We then assessed the detection rates, detection linking and orbit computation, and sources for confusion. Typical velocities of detectable TCOs will range from 1$^{\circ}$/day to 50$^{\circ}$/day, and typical apparent $V$ magnitudes from 21 to 23. Potentially-hazardous asteroids have observational characteristics similar to TCOs, but the two populations can be distinguished based on their orbits with LSST data alone. We predict that a TCO can be discovered once every year with the baseline moving-object processing system (MOPS). The rate can be increased to one TCO discovery every two months if tools complementary to the baseline MOPS are developed for the specific purpose of discovering these objects. △ Less

Submitted 5 November, 2019; originally announced November 2019.

Comments: 26 pages, 9 figures

Journal ref: Icarus 338, 113517 (2020)

The Orbit and Size-Frequency Distribution of Long Period Comets Observed by Pan-STARRS1

Authors: Benjamin Boe, Robert Jedicke, Karen J. Meech, Paul Wiegert, Robert J. Weryk, K. C. Chambers, L. Denneau, N. Kaiser, R. P. Kudritzki, E. A. Magnier, R. J. Wainscoat, C. Waters

Abstract: We introduce a new technique to estimate the comet nuclear size frequency distribution (SFD) that combines a cometary activity model with a survey simulation and apply it to 150 long period comets (LPC) detected by the Pan-STARRS1 near-Earth object survey. The debiased LPC size-frequency distribution is in agreement with previous estimates for large comets with nuclear diameter $>\sim 1$~km but we… ▽ More We introduce a new technique to estimate the comet nuclear size frequency distribution (SFD) that combines a cometary activity model with a survey simulation and apply it to 150 long period comets (LPC) detected by the Pan-STARRS1 near-Earth object survey. The debiased LPC size-frequency distribution is in agreement with previous estimates for large comets with nuclear diameter $>\sim 1$~km but we measure a significant drop in the SFD slope for small objects with diameters $<1$~km and approaching only $100$~m diameter. Large objects have a slope $α_{big} = 0.72 \pm 0.09 (stat.) \pm 0.15 (sys.)$ while small objects behave as $α_{small} = 0.07 \pm 0.03 (stat.) \pm 0.09 (sys.)$ where the SFD is $\propto 10^{αH_N}$ and $H_N$ represents the cometary nuclear absolute magnitude. The total number of LPCs that are $>1$~km diameter and have perihelia $q<10$~au is $0.46 \pm 0.15 \times 10^9$ while there are only $2.4 \pm 0.5 (stat.) \pm 2 (sys.) \times 10^9$ objects with diameters $>100$~m due to the shallow slope of the SFD for diameters $<1$~m. We estimate that the total number of `potentially active' objects with diameters $\ge 1$~km in the Oort cloud, objects that would be defined as LPCs if their perihelia evolved to $<10$~au, is $(1.5\pm1)\times10^{12}$ with a combined mass of $1.3\pm0.9 \, M_{Earth}$. The debiased LPC orbit distribution is broadly in agreement with expectations from contemporary dynamical models but there are discrepancies that could point towards a future ability to disentangle the relative importance of stellar perturbations and galactic tides in producing the LPC population. △ Less

Submitted 31 May, 2019; originally announced May 2019.

Comments: accepted for publication in Icarus

Earth's Minimoons: Opportunities for Science and Technology

Authors: Robert Jedicke, Bryce T. Bolin, William F. Bottke, Monique Chyba, Grigori Fedorets, Mikael Granvik, Lynne Jones, Hodei Urrutxua

Abstract: (modified from published version) Twelve years ago the Catalina Sky Survey discovered Earth's first known natural geocentric object other than the Moon, a few-meter diameter asteroid designated 2006 RH120. Despite significant improvements in ground-based telescope and detector technology in the past decade the asteroid surveys have not discovered another temporarily-captured orbiter (TCO; colloqui… ▽ More (modified from published version) Twelve years ago the Catalina Sky Survey discovered Earth's first known natural geocentric object other than the Moon, a few-meter diameter asteroid designated 2006 RH120. Despite significant improvements in ground-based telescope and detector technology in the past decade the asteroid surveys have not discovered another temporarily-captured orbiter (TCO; colloquially known as minimoons). Within a few years the Large Synoptic Survey Telescope (LSST) will either begin to regularly detect TCOs or force a re-analysis of the creation and dynamical evolution of small asteroids in the inner solar system. The first studies of the provenance, properties, and dynamics of Earth's minimoons suggested that there should be a steady state population with about one 1- to 2-meter diameter captured objects at any time. That model was then improved and extended to include the population of temporarily-captured flybys (TCFs), objects that fail to make an entire revolution around Earth while energetically bound to the Earth-Moon system. Several different techniques for discovering TCOs have been considered but their small diameters, proximity, and rapid motion make them challenging targets for existing ground-based optical, meteor, and radar surveys. We expect that if the TCO population is confirmed, and new objects are frequently discovered, they can provide new opportunities for 1) studying the dynamics of the Earth-Moon system, 2) testing models of the production and dynamical evolution of small asteroids from the asteroid belt, 3) rapid and frequent low delta-v missions to multiple minimoons, and 4) evaluating in-situ resource utilization techniques on asteroidal material. Here we review the past decade of minimoon studies in preparation for capitalizing on the scientific and commercial opportunities of TCOs in the first decade of LSST operations. △ Less

Submitted 31 May, 2019; originally announced May 2019.

Journal ref: Frontiers in Astronomy and Space Sciences (2018), 5, 13

The digest2 NEO Classification Code

Authors: Sonia Keys, Peter Vereš, Matthew J. Payne, Matthew J. Holman, Robert Jedicke, Gareth V. Williams, Tim Spahr, David J. Asher, Carl Hergenrother

Abstract: We describe the digest2 software package, a fast, short-arc orbit classifier for small Solar System bodies. The digest2 algorithm has been serving the community for more than 13 years. The code provides a score, D2, which represents a pseudo-probability that a tracklet belongs to a given Solar System orbit type. digest2 is primarily used as a classifier for Near-Earth Object (NEO) candidates, to i… ▽ More We describe the digest2 software package, a fast, short-arc orbit classifier for small Solar System bodies. The digest2 algorithm has been serving the community for more than 13 years. The code provides a score, D2, which represents a pseudo-probability that a tracklet belongs to a given Solar System orbit type. digest2 is primarily used as a classifier for Near-Earth Object (NEO) candidates, to identify those to be prioritized for follow-up observation. We describe the historical development of digest2 and demonstrate its use on real and synthetic data. We find that digest2 can accurately and precisely distinguish NEOs from non- NEOs. At the time of detection, 14% of NEO tracklets and 98.5% of non-NEOs tracklets have D2 below the critical value of D2 = 65. 94% of our simulated NEOs achieved the maximum D2 = 100 and 99.6% of NEOs achieved $D2 \ge 65$ at least once during the simulated 10-year timeframe. We demonstrate that D2 varies as a function of time, rate of motion, magnitude and sky-plane location, and show that NEOs tend to have lower D2 at low Solar elongations close to the ecliptic. We use our findings to recommend future development directions for the digest2 code. △ Less

Submitted 19 April, 2019; originally announced April 2019.

Comments: 19 figures, 7 tables

Journal ref: Publications of the Astronomical Society of the Pacific (2019), Vol 131, Number 1000

First Results from the rapid-response spectrophotometric characterization of Near-Earth Objects

Authors: Samuel Navarro-Meza, Michael Mommert, David Trilling, Nathaniel Butler, Mauricio Reyes-Ruiz, Barbara Pichardo, Tim Axelrod, Robert Jedicke, Nicholas Moskovitz

Abstract: As part of our multi-observatory, multi-filter campaign, we present \rmi color observations of 82 Near-Earth Objects (NEOs) obtained with the RATIR instrument on the 1.5m robotic telescope at the San Pedro Martir's National Observatory in Mexico. Our project is particularly focused on rapid response observations of small ($\lesssim 850$ m) NEOs. The rapid response and the use of spectrophotometry… ▽ More As part of our multi-observatory, multi-filter campaign, we present \rmi color observations of 82 Near-Earth Objects (NEOs) obtained with the RATIR instrument on the 1.5m robotic telescope at the San Pedro Martir's National Observatory in Mexico. Our project is particularly focused on rapid response observations of small ($\lesssim 850$ m) NEOs. The rapid response and the use of spectrophotometry allows us to constrain the taxonomic classification of NEOs with high efficiency. Here we present the methodology of our observations and our result, suggesting that the ratio of C-type to S-type asteroids in a size range of $\sim$30-850m is 1.1, which is in accordance with our previous results. We also find that 10$\%$ of all NEOs in our sample are neither C- nor S-type asteroids △ Less

Submitted 19 March, 2019; originally announced March 2019.

Comments: 31 pages, 4 tables, 10 figures

Extreme asteroids in the Pan-STARRS 1 Survey

Authors: Andrew McNeill, Alan Fitzsimmons, Robert Jedicke, Pedro Lacerda, Eva Lilly, Andrew Thompson, David E. Trilling, Ernst DeMooij, Matthew J. Hooton, Christopher A. Watson

Abstract: Using the first 18 months of the Pan-STARRS 1 survey we have identified 33 candidate high-amplitude objects for follow-up observations and carried out observations of 22 asteroids. 4 of the observed objects were found to have observed amplitude $A_{obs}\geq 1.0$ mag. We find that these high amplitude objects are most simply explained by single rubble pile objects with some density-dependent intern… ▽ More Using the first 18 months of the Pan-STARRS 1 survey we have identified 33 candidate high-amplitude objects for follow-up observations and carried out observations of 22 asteroids. 4 of the observed objects were found to have observed amplitude $A_{obs}\geq 1.0$ mag. We find that these high amplitude objects are most simply explained by single rubble pile objects with some density-dependent internal strength, allowing them to resist mass shedding even at their highly elongated shapes. 3 further objects although below the cut-off for 'high-amplitude' had a combination of elongation and rotation period which also may require internal cohesive strength, depending on the density of the body. We find that none of the 'high-amplitude asteroids' identified here require any unusual cohesive strengths to resist rotational fission. 3 asteroids were sufficiently observed to allow for shape and spin pole models to be determined through light curve inversion. 45864 was determined to have retrograde rotation with spin pole axes $λ=218\pm 10^{\circ}, β=-82\pm 5^{\circ}$ and asteroid 206167 was found to have best fit spin pole axes $λ= 57 \pm 5^{\circ}$, $β=-67 \pm 5^{\circ}$. An additional object not initially measured with $A_{obs}>1.0$ mag, 49257, was determined to have a shape model which does suggest a high-amplitude object. Its spin pole axes were best fit for values $λ=112\pm 6^{\circ}, β=6\pm 5^{\circ}$. In the course of this project to date no large super-fast rotators ($P_{rot} < 2.2$ h) have been identified. △ Less

Submitted 27 October, 2018; originally announced October 2018.

Comments: 31 pages; accepted by AJ

Debiased orbit and absolute-magnitude distributions for near-Earth objects

Authors: Mikael Granvik, Alessandro Morbidelli, Robert Jedicke, Bryce Bolin, William Bottke, Edward Beshore, David Vokrouhlicky, David Nesvorny, Patrick Michel

Abstract: The debiased absolute-magnitude and orbit distributions as well as source regions for near-Earth objects (NEOs) provide a fundamental frame of reference for studies of individual NEOs and more complex population-level questions. We present a new four-dimensional model of the NEO population that describes debiased steady-state distributions of semimajor axis, eccentricity, inclination, and absolute… ▽ More The debiased absolute-magnitude and orbit distributions as well as source regions for near-Earth objects (NEOs) provide a fundamental frame of reference for studies of individual NEOs and more complex population-level questions. We present a new four-dimensional model of the NEO population that describes debiased steady-state distributions of semimajor axis, eccentricity, inclination, and absolute magnitude $H$ in the range $17<H<25$. The modeling approach improves upon the methodology originally developed by Bottke et al. (2000; Science 288, 2190-2194) in that it is, for example, based on more realistic orbit distributions and uses source-specific absolute-magnitude distributions that allow for a power-law slope that varies with $H$. We divide the main asteroid belt into six different entrance routes or regions (ER) to the NEO region: the $ν_6$, 3:1J, 5:2J and 2:1J resonance complexes as well as Hungarias and Phocaeas. In addition we include the Jupiter-family comets as the primary cometary source of NEOs. We calibrate the model against NEO detections by Catalina Sky Surveys' stations 703 and G96 during 2005-2012, and utilize the complementary nature of these two systems to quantify the systematic uncertainties associated to the resulting model. The most important ERs are the $ν_6$ and 3:1J resonance complexes with JFCs contributing a few percent of NEOs on average. We predict that there are $962^{+52}_{-56}$ ($802^{+48}_{-42}\times10^3$) NEOs with $H<17.75$ ($H<25$) and these numbers are in agreement with the most recent estimates found in the literature (the uncertainty estimates only account for the random component). △ Less

Submitted 26 April, 2018; originally announced April 2018.

Comments: 33 pages, 28 figures, accepted for publication in Icarus

An Investigation of the Ranges of Validity of Asteroid Thermal Models for Near-Earth Asteroid Observations

Authors: Michael Mommert, Robert Jedicke, David E. Trilling

Abstract: The majority of known asteroid diameters are derived from thermal-infrared observations. Diameters are derived using asteroid thermal models that approximate their surface temperature distributions and compare the measured thermal-infrared flux with model-dependent predictions. The most commonly used thermal model is the Near-Earth Asteroid Thermal Model (NEATM), which is usually perceived as supe… ▽ More The majority of known asteroid diameters are derived from thermal-infrared observations. Diameters are derived using asteroid thermal models that approximate their surface temperature distributions and compare the measured thermal-infrared flux with model-dependent predictions. The most commonly used thermal model is the Near-Earth Asteroid Thermal Model (NEATM), which is usually perceived as superior to other models like the Fast-Rotating Model (FRM). We investigate the applicability of the NEATM and the FRM to thermal-infrared observations of Near-Earth Objects using synthetic asteroids with properties based on the real Near-Earth Asteroid (NEA) population. We find the NEATM to provide more accurate diameters and albedos than the FRM in most cases, with a few exceptions. The modeling results are barely affected by the physical properties of the objects, but we find a large impact of the solar phase angle on the modeling results. We conclude that the NEATM provides statistically more robust diameter estimates for NEAs observed at solar phase angles less than ~65°, while the FRM provides more robust diameter estimates for solar phase angles greater than ~65°. We estimate that <5% of all NEA diameters and albedos derived up to date are affected by systematic effects that are of the same order of magnitude as the typical thermal model uncertainties. We provide statistical correction functions for diameters and albedos derived using the NEATM and FRM as a function of solar phase angle. △ Less

Submitted 10 January, 2018; originally announced January 2018.

Comments: 25 pages, 8 figures, accepted for publication in AJ

Spectroscopy and thermal modelling of the first interstellar object 1I/2017 U1 'Oumuamua

Authors: Alan Fitzsimmons, Colin Snodgrass, Ben Rozitis, Bin Yang, Meabh Hyland, Tom Seccull, Michele T. Bannister, Wesley C. Fraser, Robert Jedicke, Pedro Lacerda

Abstract: During the formation and evolution of the Solar System, significant numbers of cometary and asteroidal bodies were ejected into interstellar space$^{1,2}$. It can be reasonably expected that the same happened for planetary systems other than our own. Detection of such Inter- stellar Objects (ISOs) would allow us to probe the planetesimal formation processes around other stars, possibly together wi… ▽ More During the formation and evolution of the Solar System, significant numbers of cometary and asteroidal bodies were ejected into interstellar space$^{1,2}$. It can be reasonably expected that the same happened for planetary systems other than our own. Detection of such Inter- stellar Objects (ISOs) would allow us to probe the planetesimal formation processes around other stars, possibly together with the effects of long-term exposure to the interstellar medium. 1I/2017 U1 'Oumuamua is the first known ISO, discovered by the Pan-STARRS1 telescope in October 2017$^3$.The discovery epoch photometry implies a highly elongated body with radii of $\sim 200 \times 20$ m when a comet-like geometric albedo of 0.04 is assumed. Here we report spectroscopic characterisation of 'Oumuamua, finding it to be variable with time but similar to organically rich surfaces found in the outer Solar System. The observable ISO population is expected to be dominated by comet-like bodies in agreement with our spectra, yet the reported inactivity implies a lack of surface ice. We show this is consistent with predictions of an insulating mantle produced by long-term cosmic ray exposure. An internal icy composition cannot therefore be ruled out by the lack of activity, even though 'Oumuamua passed within 0.25 au of the Sun. △ Less

Submitted 18 December, 2017; originally announced December 2017.

Comments: 4 figures, 1 table; Accepted for publication in Nature Astronomy, data are available at https://github.com/asteroidnerd/Oumuamua-spectroscopy

CO-Driven Activity in Comet C/2017 K2 (PANSTARRS)

Authors: Karen J. Meech, Jan T. Kleyna, Olivier Hainaut, Marco Micheli, James Bauer, Larry Denneau, Jacqueline V. Keane, Haynes Stephens, Robert Jedicke, Richard Wainscoat, Robert Weryk, Heather Flewelling, Eva Lilly, Eugene Magnier, Kenneth C. Chambers

Abstract: Comet C/2017 K2 (PANSTARRS) was discovered by the Pan-STARRS1 (PS1) Survey on 2017 May 21 at a distance 16.09 au from the Sun, the second most distant discovery of an active comet. Pre-discovery images in the PS1 archive back to 2014 and additional deep CFHT images between 2013 May 10-13 showed the comet to be active at 23.75 au. We derive an upper limit to the nucleus radius of $R_N$=80 km, assum… ▽ More Comet C/2017 K2 (PANSTARRS) was discovered by the Pan-STARRS1 (PS1) Survey on 2017 May 21 at a distance 16.09 au from the Sun, the second most distant discovery of an active comet. Pre-discovery images in the PS1 archive back to 2014 and additional deep CFHT images between 2013 May 10-13 showed the comet to be active at 23.75 au. We derive an upper limit to the nucleus radius of $R_N$=80 km, assuming a 4\% albedo. The spectral reflectivity of the comet surface is similar to "fresh" regions seen on comet 67P/Churyumov-Gerasimenko using the $Rosetta$ OSIRIS camera. Pre-discovery photometry combined with new data obtained with Megacam on the CFHT show that the activity is consistent with CO-ice sublimation and inconsistent with CO$_2$-ice sublimation. The ice sublimation models were run out to perihelion in 2022 at 1.8 au to predict the CO production rates, assuming that the outgassing area does not change. Assuming a canonical 4\% active surface area for water-ice sublimation, we present production rate ratios, $Q_{\rm CO}$/$Q_{\rm H2O}$, for a range of nucleus sizes. Comparing these results with other CO-rich comets we derive a lower limit to the nucleus radius of $\sim$14 km. We present predictions for $Q_{\rm CO}$ at a range of distances that will be useful for planning observations with JWST and large ground-based facilities. △ Less

Submitted 10 October, 2017; originally announced October 2017.

Comments: 5 pages, 4 figures

Distribution of shape elongations of main belt asteroids derived from Pan-STARRS1 photometry

Authors: H. Cibulková, H. Nortunen, J. Ďurech, M. Kaasalainen, P. Vereš, R. Jedicke, R. J. Wainscoat, M. Mommert, D. E. Trilling, E. Schunová-Lilly, E. A. Magnier, C. Waters, H. Flewelling

Abstract: Context. A lot of photometric data is produced by surveys such as Pan-STARRS, LONEOS, WISE or Catalina. These data are a rich source of information about the physical properties of asteroids. There are several possible approaches for utilizing these data. Lightcurve inversion is a typical method that works with individual asteroids. Our approach in this paper is statistical when we focused on larg… ▽ More Context. A lot of photometric data is produced by surveys such as Pan-STARRS, LONEOS, WISE or Catalina. These data are a rich source of information about the physical properties of asteroids. There are several possible approaches for utilizing these data. Lightcurve inversion is a typical method that works with individual asteroids. Our approach in this paper is statistical when we focused on large groups of asteroids like dynamical families and taxonomic classes, and the data were not sufficient for individual models. Aims. Our aim was to study the distributions of shape elongation $b/a$ and the spin axis latitude $β$ for various subpopulations of asteroids and to compare our results, based on Pan-STARRS1 survey, with statistics previously done using different photometric data (Lowell database, WISE data). Methods. We use the LEADER algorithm to compare the $b/a$ and $β$ distributions for different subpopulations of asteroids. The algorithm creates a cumulative distributive function (CDF) of observed brightness variations, and computes the $b/a$ and $β$ distributions using analytical basis functions that yield the observed CDF. A variant of LEADER is used to solve the joint distributions for synthetic populations to test the validity of the method. Results. When comparing distributions of shape elongation for groups of asteroids with different diameters $D$, we found that there are no differences for $D < 25$ km. We also constructed distributions for asteroids with different rotation periods and revealed that the fastest rotators with $P = 0 - 4$ h are more spheroidal than the population with $P = 4 - 8$ h. △ Less

Submitted 17 September, 2017; originally announced September 2017.

Comments: Accepted for publication in Astronomy&Astrophysics

Journal ref: A&A 611, A86 (2018)

A dwarf planet class object in the 21:5 resonance with Neptune

Authors: Matthew J. Holman, Matthew J. Payne, Wesley Fraser, Pedro Lacerda, Michele T. Bannister, Michael Lackner, Ying-Tung Chen, Hsing Wen Lin, Kenneth W. Smith, Rositako Kotanekova, David Young, K. Chambers, S. Chastel, L. Denneau, A. Fitzsimmons, H. Flewelling, Tommy Grav, M. Huber, Nick Induni, Rolf-Peter Kudritzki, Alex Krolewski, R. Jedicke, N. Kaiser, E. Lilly, E. Magnier , et al. (11 additional authors not shown)

Abstract: We report the discovery of a $H_r = 3.4\pm0.1$ dwarf planet candidate by the Pan-STARRS Outer Solar System Survey. 2010 JO$_{179}$ is red with $(g-r)=0.88 \pm 0.21$, roughly round, and slowly rotating, with a period of $30.6$ hr. Estimates of its albedo imply a diameter of 600--900~km. Observations sampling the span between 2005--2016 provide an exceptionally well-determined orbit for 2010 JO… ▽ More We report the discovery of a $H_r = 3.4\pm0.1$ dwarf planet candidate by the Pan-STARRS Outer Solar System Survey. 2010 JO$_{179}$ is red with $(g-r)=0.88 \pm 0.21$, roughly round, and slowly rotating, with a period of $30.6$ hr. Estimates of its albedo imply a diameter of 600--900~km. Observations sampling the span between 2005--2016 provide an exceptionally well-determined orbit for 2010 JO$_{179}$, with a semi-major axis of $78.307\pm0.009$ au, distant orbits known to this precision are rare. We find that 2010 JO$_{179}$ librates securely within the 21:5 mean-motion resonance with Neptune on hundred-megayear time scales, joining the small but growing set of known distant dwarf planets on metastable resonant orbits. These imply a substantial trans-Neptunian population that shifts between stability in high-order resonances, the detached population, and the eroding population of the scattering disk. △ Less

Submitted 15 September, 2017; originally announced September 2017.

Spin states of asteroids in the Eos collisional family

Authors: J. Hanuš, M. Delbo', V. Alí-Lagoa, B. Bolin, R. Jedicke, J. Ďurech, H. Cibulková, P. Pravec, P. Kušnirák, R. Behrend, F. Marchis, P. Antonini, L. Arnold, M. Audejean, M. Bachschmidt, L. Bernasconi, L. Brunetto, S. Casulli, R. Dymock, N. Esseiva, M. Esteban, O. Gerteis, H. de Groot, H. Gully, H. Hamanowa , et al. (12 additional authors not shown)

Abstract: Eos family was created during a catastrophic impact about 1.3 Gyr ago. Rotation states of individual family members contain information about the history of the whole population. We aim to increase the number of asteroid shape models and rotation states within the Eos collision family, as well as to revise previously published shape models from the literature. Such results can be used to constrain… ▽ More Eos family was created during a catastrophic impact about 1.3 Gyr ago. Rotation states of individual family members contain information about the history of the whole population. We aim to increase the number of asteroid shape models and rotation states within the Eos collision family, as well as to revise previously published shape models from the literature. Such results can be used to constrain theoretical collisional and evolution models of the family, or to estimate other physical parameters by a thermophysical modeling of the thermal infrared data. We use all available disk-integrated optical data (i.e., classical dense-in-time photometry obtained from public databases and through a large collaboration network as well as sparse-in-time individual measurements from a few sky surveys) as input for the convex inversion method, and derive 3D shape models of asteroids together with their rotation periods and orientations of rotation axes. We present updated shape models for 15 asteroids and new shape model determinations for 16 asteroids. Together with the already published models from the publicly available DAMIT database, we compiled a sample of 56 Eos family members with known shape models that we used in our analysis of physical properties within the family. Rotation states of asteroids smaller than ~20 km are heavily influenced by the YORP effect, whilst the large objects more or less retained their rotation state properties since the family creation. Moreover, we also present a shape model and bulk density of asteroid (423) Diotima, an interloper in the Eos family, based on the disk-resolved data obtained by the Near InfraRed Camera (Nirc2) mounted on the W.M. Keck II telescope. △ Less

Submitted 18 July, 2017; originally announced July 2017.

Comments: Accepted for publication in ICARUS Special Issue - Asteroids: Origin, Evolution & Characterization

The splitting of double-component active asteroid P/2016 J1 (PANSTARRS)

Authors: Fernando Moreno, Francisco Pozuelos, Bojan Novakovic, Javier Licandro, Antonio Cabrera-Lavers, Bryce Bolin, Robert Jedicke, Brett Gladman, Michele Bannister, Stephen Gwyn, Peter Veres, Kenneth Chambers, Serge Chastel, Larry Denneau, Heather Flewelling, Mark Huber, Eva Schunova-Lilly, Eugene Magnier, Richard Wainscoat, Christopher Waters, Robert Weryk, Davide Farnocchia, Marco Micheli

Abstract: We present deep imaging observations, orbital dynamics, and dust tail model analyses of the double-component asteroid P/2016 J1 (J1-A and J1-B). The observations were acquired at the Gran Telescopio Canarias (GTC) and the Canada-France-Hawaii Telescope (CFHT) from mid March to late July, 2016. A statistical analysis of backward-in-time integrations of the orbits of a large sample of clone objects… ▽ More We present deep imaging observations, orbital dynamics, and dust tail model analyses of the double-component asteroid P/2016 J1 (J1-A and J1-B). The observations were acquired at the Gran Telescopio Canarias (GTC) and the Canada-France-Hawaii Telescope (CFHT) from mid March to late July, 2016. A statistical analysis of backward-in-time integrations of the orbits of a large sample of clone objects of P/2016 J1-A and J1-B shows that the minimum separation between them occurred most likely $\sim$2300 days prior to the current perihelion passage, i.e., during the previous orbit near perihelion. This closest approach was probably linked to a fragmentation event of their parent body. Monte Carlo dust tail models show that those two components became active simultaneously $\sim$250 days before the current perihelion, with comparable maximum loss rates of $\sim$0.7 kg s$^{-1}$ and $\sim$0.5 kg s$^{-1}$, and total ejected masses of 8$\times$10$^{6}$ kg and 6$\times$10$^{6}$ kg for fragments J1-A and J1-B, respectively. In consequence, the fragmentation event and the present dust activity are unrelated. The simultaneous activation times of the two components and the fact that the activity lasted 6 to 9 months or longer, strongly indicate ice sublimation as the most likely mechanism involved in the dust emission process. △ Less

Submitted 13 February, 2017; originally announced February 2017.

Comments: Accepted by ApJ Letters, Feb. 17, 2017

An Observational Upper Limit on the Interstellar Number Density of Asteroids and Comets

Authors: Toni Engelhardt, Robert Jedicke, Peter Veres, Alan Fitzsimmons, Larry Denneau, Ed Beshore, Bonnie Meinke

Abstract: We derived 90% confidence limits (CL) on the interstellar number density ($ρ_{IS}^{CL}$) of interstellar objects (ISO; comets and asteroids) as a function of the slope of their size-frequency distribution and limiting absolute magnitude. To account for gravitational focusing, we first generated a quasi-realistic ISO population to ~750 au from the Sun and propagated it forward in time to generate a… ▽ More We derived 90% confidence limits (CL) on the interstellar number density ($ρ_{IS}^{CL}$) of interstellar objects (ISO; comets and asteroids) as a function of the slope of their size-frequency distribution and limiting absolute magnitude. To account for gravitational focusing, we first generated a quasi-realistic ISO population to ~750 au from the Sun and propagated it forward in time to generate a steady state population of ISOs with heliocentric distance <50 au. We then simulated the detection of the synthetic ISOs using pointing data for each image and average detection efficiencies for each of three contemporary solar system surveys --- PS1, the Mt. Lemmon Survey, and the Catalina Sky Survey. These simulations allowed us to determine the surveys' combined ISO detection efficiency under several different but realistic modes of identifying ISOs in the survey data. Some of the synthetic detected ISOs had eccentricities as small as 1.01 --- in the range of the largest eccentricities of several known comets. Our best CL of $ρ_{IS}^{CL} = 1.4 \times 10^{-4}$ au$^{-3}$ implies that the expectation that extra-solar systems form like our solar system, eject planetesimals in the same way, and then distribute them throughout the galaxy, is too simplistic, or that the SFD or behavior of ISOs as they pass through our solar system is far from expectations. △ Less

Submitted 7 February, 2017; originally announced February 2017.

Comments: accepted for publication in AJ

The Pan-STARRS1 Surveys

Authors: K. C. Chambers, E. A. Magnier, N. Metcalfe, H. A. Flewelling, M. E. Huber, C. Z. Waters, L. Denneau, P. W. Draper, D. Farrow, D. P. Finkbeiner, C. Holmberg, J. Koppenhoefer, P. A. Price, A. Rest, R. P. Saglia, E. F. Schlafly, S. J. Smartt, W. Sweeney, R. J. Wainscoat, W. S. Burgett, S. Chastel, T. Grav, J. N. Heasley, K. W. Hodapp, R. Jedicke , et al. (101 additional authors not shown)

Abstract: Pan-STARRS1 has carried out a set of distinct synoptic imaging sky surveys including the $3π$ Steradian Survey and the Medium Deep Survey in 5 bands ($grizy_{P1}$). The mean 5$σ$ point source limiting sensitivities in the stacked 3$π$ Steradian Survey in $grizy_{P1}$ are (23.3, 23.2, 23.1, 22.3, 21.4) respectively. The upper bound on the systematic uncertainty in the photometric calibration across… ▽ More Pan-STARRS1 has carried out a set of distinct synoptic imaging sky surveys including the $3π$ Steradian Survey and the Medium Deep Survey in 5 bands ($grizy_{P1}$). The mean 5$σ$ point source limiting sensitivities in the stacked 3$π$ Steradian Survey in $grizy_{P1}$ are (23.3, 23.2, 23.1, 22.3, 21.4) respectively. The upper bound on the systematic uncertainty in the photometric calibration across the sky is 7-12 millimag depending on the bandpass. The systematic uncertainty of the astrometric calibration using the Gaia frame comes from a comparison of the results with Gaia: the standard deviation of the mean and median residuals ($ Δra, Δdec $) are (2.3, 1.7) milliarcsec, and (3.1, 4.8) milliarcsec respectively. The Pan-STARRS system and the design of the PS1 surveys are described and an overview of the resulting image and catalog data products and their basic characteristics are described together with a summary of important results. The images, reduced data products, and derived data products from the Pan-STARRS1 surveys are available to the community from the Mikulski Archive for Space Telescopes (MAST) at STScI. △ Less

Submitted 28 January, 2019; v1 submitted 16 December, 2016; originally announced December 2016.

Comments: 38 pages, 29 figures, 12 tables

Pan-STARRS Pixel Processing: Detrending, Warping, Stacking

Authors: C. Z. Waters, E. A. Magnier, P. A. Price, K. C. Chambers, W. S. Burgett, P. Draper, H. A. Flewelling, K. W. Hodapp, M. E. Huber, R. Jedicke, N. Kaiser, R. -P. Kudritzki, R. H. Lupton, N. Metcalfe, A. Rest, W. E. Sweeney, J. L. Tonry, R. J. Wainscoat, W. M. Wood-Vasey, PS1 Builders

Abstract: The Pan-STARRS1 Science Consortium have carried out a set of imaging surveys using the 1.4 giga-pixel GPC1 camera on the PS1 telescope. As this camera is composed of many individual electronic readouts, and covers a very large field of view, great care was taken to ensure that the many instrumental effects were corrected to produce the most uniform detector response possible. We present the image… ▽ More The Pan-STARRS1 Science Consortium have carried out a set of imaging surveys using the 1.4 giga-pixel GPC1 camera on the PS1 telescope. As this camera is composed of many individual electronic readouts, and covers a very large field of view, great care was taken to ensure that the many instrumental effects were corrected to produce the most uniform detector response possible. We present the image detrending steps used as part of the processing of the data contained within the public release of the Pan-STARRS1 Data Release 1 (DR1). In addition to the single image processing, the methods used to transform the 375,573 individual exposures into a common sky-oriented grid are discussed, as well as those used to produce both the image stack and difference combination products. △ Less

Submitted 28 January, 2019; v1 submitted 15 December, 2016; originally announced December 2016.

Comments: 21 pages, 21 figures, Paper III of Pan-STARRS DR1 public data release documentation

Pan-STARRS Pixel Analysis : Source Detection and Characterization

Authors: Eugene A. Magnier, W. E. Sweeney, K. C. Chambers, H. A. Flewelling, M. E. Huber, P. A. Price, C. Z. Waters, L. Denneau, P. Draper, R. Jedicke, K. W. Hodapp, N. Kaiser, R. -P. Kudritzki, N. Metcalfe, C. W. Stubbs, R. J. Wainscoast

Abstract: Over 3 billion astronomical objects have been detected in the more than 22 million orthogonal transfer CCD images obtained as part of the Pan-STARRS1 $3π$ survey. Over 85 billion instances of those objects have been automatically detected and characterized by the Pan-STARRS Image Processing Pipeline photometry software, psphot. This fast, automatic, and reliable software was developed for the Pan-… ▽ More Over 3 billion astronomical objects have been detected in the more than 22 million orthogonal transfer CCD images obtained as part of the Pan-STARRS1 $3π$ survey. Over 85 billion instances of those objects have been automatically detected and characterized by the Pan-STARRS Image Processing Pipeline photometry software, psphot. This fast, automatic, and reliable software was developed for the Pan-STARRS project, but is easily adaptable to images from other telescopes. We describe the analysis of the astronomical objects by psphot in general as well as for the specific case of the 3rd processing version used for the first two public releases of the Pan-STARRS $3π$ survey data, DR1 & DR2. △ Less

Submitted 26 January, 2019; v1 submitted 15 December, 2016; originally announced December 2016.

Comments: Pan-STARRS Public Data Release 2 : Paper IV

The Pan-STARRS1 Database and Data Products

Authors: H. A. Flewelling, E. A. Magnier, K. C. Chambers, J. N. Heasley, C. Holmberg, M. E. Huber, W. Sweeney, C. Z. Waters, A. Calamida, S. Casertano, X. Chen, D. Farrow, G. Hasinger, R. Henderson, K. S. Long, N. Metcalfe, G. Narayan, M. A. Nieto-Santisteban, P. Norberg, A. Rest, R. P. Saglia, A. Szalay, A. R. Thakar, J. L. Tonry, J. Valenti , et al. (15 additional authors not shown)

Abstract: This paper describes the organization of the database and the catalog data products from the Pan-STARRS1 $3π$ Steradian Survey. The catalog data products are available in the form of an SQL-based relational database from MAST, the Mikulski Archive for Space Telescopes at STScI. The database is described in detail, including the construction of the database, the provenance of the data, the schema,… ▽ More This paper describes the organization of the database and the catalog data products from the Pan-STARRS1 $3π$ Steradian Survey. The catalog data products are available in the form of an SQL-based relational database from MAST, the Mikulski Archive for Space Telescopes at STScI. The database is described in detail, including the construction of the database, the provenance of the data, the schema, and how the database tables are related. Examples of queries for a range of science goals are included. The catalog data products are available in the form of an SQL-based relational database from MAST, the Mikulski Archive for Space Telescopes at STScI. △ Less

Submitted 29 January, 2019; v1 submitted 15 December, 2016; originally announced December 2016.

Comments: 63 pages, 59 tables, 16 figures

Pan-STARRS Photometric and Astrometric Calibration

Authors: Eugene. A. Magnier, Edward. F. Schlafly, Douglas P. Finkbeiner, J. L. Tonry, B. Goldman, S. Röser, E. Schilbach, K. C. Chambers, H. A. Flewelling, M. E. Huber, P. A. Price, W. E. Sweeney, C. Z. Waters, L. Denneau, P. Draper, K. W. Hodapp, R. Jedicke, R. -P. Kudritzki, N. Metcalfe, C. W. Stubbs, R. J. Wainscoast

Abstract: We present the details of the photometric and astrometric calibration of the Pan-STARRS1 $3π$ Survey. The photometric goals were to reduce the systematic effects introduced by the camera and detectors, and to place all of the observations onto a photometric system with consistent zero points over the entire area surveyed, the ~30,000 square degrees north of $δ$ = -30 degrees. The astrometric calib… ▽ More We present the details of the photometric and astrometric calibration of the Pan-STARRS1 $3π$ Survey. The photometric goals were to reduce the systematic effects introduced by the camera and detectors, and to place all of the observations onto a photometric system with consistent zero points over the entire area surveyed, the ~30,000 square degrees north of $δ$ = -30 degrees. The astrometric calibration compensates for similar systematic effects so that positions, proper motions, and parallaxes are reliable as well. The Pan-STARRS Data Release 2 (DR2) astrometry is tied to the Gaia DR1 release. △ Less

Submitted 26 January, 2019; v1 submitted 15 December, 2016; originally announced December 2016.

Comments: Pan-STARRS Public Data Release 2 : Paper V

The Pan-STARRS Data Processing System

Authors: Eugene A. Magnier, K. C. Chambers, H. A. Flewelling, J. C. Hoblitt, M. E. Huber, P. A. Price, W. E. Sweeney, C. Z. Waters, L. Denneau, P. Draper, K. W. Hodapp, R. Jedicke, N. Kaiser, R. -P. Kudritzki, N. Metcalfe, C. W. Stubbs, R. J. Wainscoast

Abstract: The Pan-STARRS Data Processing System is responsible for the steps needed to downloaded, archive, and process all images obtained by the Pan-STARRS telescopes, including real-time detection of transient sources such as supernovae and moving objects including potentially hazardous asteroids. With a nightly data volume of up to 4 terabytes and an archive of over 4 petabytes of raw imagery, Pan-STARR… ▽ More The Pan-STARRS Data Processing System is responsible for the steps needed to downloaded, archive, and process all images obtained by the Pan-STARRS telescopes, including real-time detection of transient sources such as supernovae and moving objects including potentially hazardous asteroids. With a nightly data volume of up to 4 terabytes and an archive of over 4 petabytes of raw imagery, Pan-STARRS is solidly in the realm of Big Data astronomy. The full data processing system consists of several subsystems covering the wide range of necessary capabilities. This article describes the Image Processing Pipeline and its connections to both the summit data systems and the outward-facing systems downstream. The latter include the Moving Object Processing System (MOPS) & the public database: the Published Science Products Subsystem (PSPS). △ Less

Submitted 26 January, 2019; v1 submitted 15 December, 2016; originally announced December 2016.

Comments: Pan-STARRS Public Data Release 2 : Paper II

The size-frequency distribution of H>13 NEOs and ARM targets detected by Pan-STARRS1

Authors: Eva Schunová Lilly, Robert Jedicke, Peter Vereš, Larry Denneau, Richard J. Wainscoat

Abstract: We determine the absolute magnitude (H) distribution (or size-frequency distribution, SFD; $N(H) \propto 10^{αH}$ where $α$ is the slope of the distribution) for near-Earth objects (NEO) with $13<H<30$ and Asteroid Retrieval Mission (ARM) targets with $27<H<31$ that were detected by the 1\st\ telescope of the Panoramic Survey Telescope and Rapid Response System - Pan-STARRS1 (e.g. Kaiser et al. 20… ▽ More We determine the absolute magnitude (H) distribution (or size-frequency distribution, SFD; $N(H) \propto 10^{αH}$ where $α$ is the slope of the distribution) for near-Earth objects (NEO) with $13<H<30$ and Asteroid Retrieval Mission (ARM) targets with $27<H<31$ that were detected by the 1\st\ telescope of the Panoramic Survey Telescope and Rapid Response System - Pan-STARRS1 (e.g. Kaiser et al. 2002, Kaiser 2004, Hodapp et al. 2004). The NEO and ARM target detection efficiencies were calculated using the Greenstreet et al. (2012) NEO orbit distribution. The debiased Pan-STARRS1 NEO absolute magnitude distribution is more complex than a single slope power law - it shows two transitions - at H$\sim$16 from steep to shallow slope, and in the $21<H<23$ interval from a shallow to steep slope, which is consistent with other recent works (e.g. Mainzer et al. 2011c, Brown et al. 2013, Harris and D`Abramo 2015). We fit $α= 0.48\pm0.02$ for NEOs with $13<H<16$, $α= 0.33\pm0.01$ for NEOs with $16<H<22$, and $α= 0.62\pm0.03$ for the smaller objects with $H>22$. There is also another change in slope from steep to shallow around H=27. The three ARM target candidates detected by Pan-STARRS1 in one year of surveying have a corrected SFD with slope $α= 0.40^{+0.33}_{-0.45}$. We also show that the window for follow up observations of small (H$\gtrsim$22) NEOs with the NASA IRTF telescope and Arecibo and Goldstone radars are extremely short - on order of days, and procedures for fast response must be implemented in order to measure physical characteristics of small Earth-approaching objects. CFHT's MegaCam and Pan-STARRS1 have longer observing windows and are capable of following-up more NEOs due to their deeper limiting magnitudes and wider fields of view. △ Less

Submitted 17 November, 2016; v1 submitted 9 November, 2016; originally announced November 2016.

Comments: 48 pages, 11 figures, 2 tables

The Pan-STARRS 1 Discoveries of five new Neptune Trojans

Authors: Hsing Wen Lin, Ying-Tung Chen, Matthew J. Holman, Wing-Huen Ip, M. J. Payne, P. Lacerda, W. C. Fraser, D. W. Gerdes, A. Bieryla, Z. -F. Sie, W. -P. Chen, W. S. Burgett, L. Denneau, R. Jedicke, N. Kaiser, E. A. Magnier, J. L. Tonry, R. J. Wainscoat, C. Waters

Abstract: In this work we report the detection of seven Neptune Trojans (NTs) in the Pan-STARRS 1 (PS1) survey. Five of these are new discoveries, consisting of four L4 Trojans and one L5 Trojan. Our orbital simulations show that the L5 Trojan stably librates for only several million years. This suggests that the L5 Trojan must be of recent capture origin. On the other hand, all four new L4 Trojans stably o… ▽ More In this work we report the detection of seven Neptune Trojans (NTs) in the Pan-STARRS 1 (PS1) survey. Five of these are new discoveries, consisting of four L4 Trojans and one L5 Trojan. Our orbital simulations show that the L5 Trojan stably librates for only several million years. This suggests that the L5 Trojan must be of recent capture origin. On the other hand, all four new L4 Trojans stably occupy the 1:1 resonance with Neptune for more than 1 Gyr. They can, therefore, be of primordial origin. Our survey simulation results show that the inclination width of the Neptune Trojan population should be between $7^{\circ}$ and $27^{\circ}$ at $>$ 95% confidence, and most likely $\sim 11^{\circ}$. In this paper, we describe the PS1 survey, the Outer Solar System pipeline, the confirming observations, and the orbital/physical properties of the new Neptune Trojans. △ Less

Submitted 15 September, 2016; originally announced September 2016.

Comments: 24 pages, 6 figures, AJ accepted

Discovery of A New Retrograde Trans-Neptunian Object: Hint of A Common Orbital Plane for Low Semi-Major Axis, High Inclination TNOs and Centaurs

Authors: Ying-Tung Chen, Hsing Wen Lin, Matthew J. Holman, Matthew J. Payne, Wesley C. Fraser, Pedro Lacerda, Wing-Huen Ip, Wen-Ping Chen, Rolf-Peter Kudritzki, Robert Jedicke, Richard J. Wainscoat, John L. Tonry, Eugene A. Magnier, Christopher Waters, Nick Kaiser, Shiang-Yu Wang, Matthew Lehner

Abstract: Although the majority of Centaurs are thought to have originated in the scattered disk, with the high-inclination members coming from the Oort cloud, the origin of the high inclination component of trans-Neptunian objects (TNOs) remains uncertain. We report the discovery of a retrograde TNO, which we nickname "Niku", detected by the Pan-STARRS 1 Outer Solar System Survey. Our numerical integration… ▽ More Although the majority of Centaurs are thought to have originated in the scattered disk, with the high-inclination members coming from the Oort cloud, the origin of the high inclination component of trans-Neptunian objects (TNOs) remains uncertain. We report the discovery of a retrograde TNO, which we nickname "Niku", detected by the Pan-STARRS 1 Outer Solar System Survey. Our numerical integrations show that the orbital dynamics of Niku are very similar to that of 2008 KV$_{42}$ (Drac), with a half-life of $\sim 500$ Myr. Comparing similar high inclination TNOs and Centaurs ($q > 10$ AU, $a < 100$ AU and $i > 60^\circ$), we find that these objects exhibit a surprising clustering of ascending node, and occupy a common orbital plane. This orbital configuration has high statistical significance: 3.8-$σ$. An unknown mechanism is required to explain the observed clustering. This discovery may provide a pathway to investigate a possible reservoir of high-inclination objects. △ Less

Submitted 5 August, 2016; originally announced August 2016.

Comments: 18 pages, 4 figures, 1 table, accepted for publication in ApJ Letters

Distant Solar System Objects identified in the Pan-STARRS1 survey

Authors: R. J. Weryk, E. Lilly, S. Chastel, L. Denneau, R. Jedicke, E. Magnier, R. J. Wainscoat, K. Chambers, H. Flewelling, M. E. Huber, C. Waters, the PS1 Builders

Abstract: We present a method to identify distant solar system objects in long-term wide-field asteroid survey data, and conduct a search for them in the Pan-STARRS1 (PS1) image data acquired from 2010 to mid-2015. We demonstrate that our method is able to find multi-opposition orbital links, and we present the resulting orbital distributions which consist of 154 Centaurs, 255 classical Trans-Neptunian Obje… ▽ More We present a method to identify distant solar system objects in long-term wide-field asteroid survey data, and conduct a search for them in the Pan-STARRS1 (PS1) image data acquired from 2010 to mid-2015. We demonstrate that our method is able to find multi-opposition orbital links, and we present the resulting orbital distributions which consist of 154 Centaurs, 255 classical Trans-Neptunian Objects (TNOs), 121 resonant TNOs, 89 Scattered Disc Objects (SDOs) and 10 comets. Our results show more than half of these are new discoveries, including a newly discovered 19th magnitude TNO. Our identified objects do not show clustering in their argument of perihelia, which if present, might support the existence of a large unknown planetary-sized object in the outer solar system. △ Less

Submitted 17 July, 2016; originally announced July 2016.

Comments: submitted to Icarus, 32 pages, 9 figures

Near-Earth asteroid (3200) Phaethon. Characterization of its orbit, spin state, and thermophysical parameters

Authors: J. Hanus, M. Delbo, D. Vokrouhlicky, P. Pravec, J. P. Emery, V. Ali-Lagoa, B. Bolin, M. Devogele, R. Dyvig, A. Galad, R. Jedicke, L. Kornos, P. Kusnirak, J. Licandro, V. Reddy, J-P. Rivet, J. Vilagi, B. D. Warner

Abstract: The near-Earth asteroid (3200) Phaethon is an intriguing object: its perihelion is at only 0.14 au and is associated with the Geminid meteor stream. We aim to use all available disk-integrated optical data to derive a reliable convex shape model of Phaethon. By interpreting the available space- and ground-based thermal infrared data and Spitzer spectra using a thermophysical model, we also aim to… ▽ More The near-Earth asteroid (3200) Phaethon is an intriguing object: its perihelion is at only 0.14 au and is associated with the Geminid meteor stream. We aim to use all available disk-integrated optical data to derive a reliable convex shape model of Phaethon. By interpreting the available space- and ground-based thermal infrared data and Spitzer spectra using a thermophysical model, we also aim to further constrain its size, thermal inertia, and visible geometric albedo. We applied the convex inversion method to the new optical data obtained by six instruments and to previous observations. The convex shape model was then used as input for the thermophysical modeling. We also studied the long-term stability of Phaethon's orbit and spin axis with a numerical orbital and rotation-state integrator. We present a new convex shape model and rotational state of Phaethon: a sidereal rotation period of 3.603958(2) h and ecliptic coordinates of the preferred pole orientation of (319$^{\circ}$, $-$39$^{\circ}$) with a 5$^{\circ}$ uncertainty. Moreover, we derive its size ($D$=5.1$\pm$0.2 km), thermal inertia ($Γ$=600$\pm$200 J m$^{-2}$ s$^{-1/2}$ K$^{-1}$), geometric visible albedo ($p_{\mathrm{V}}$=0.122$\pm$0.008), and estimate the macroscopic surface roughness. We also find that the Sun illumination at the perihelion passage during the past several thousand years is not connected to a specific area on the surface, which implies non-preferential heating. △ Less

Submitted 26 May, 2016; v1 submitted 17 May, 2016; originally announced May 2016.

Comments: Astronomy and Astrophysics. In press

Journal ref: A&A 592, A34 (2016)

Brightness variation distributions among main belt asteroids from sparse light curve sampling with Pan-STARRS 1

Authors: A. McNeill, A. Fitzsimmons, R. Jedicke, R. Wainscoat, L. Denneau, P. Veres, E. Magnier, K. C. Chambers, N. Kaiser, C. Waters

Abstract: The rotational state of asteroids is controlled by various physical mechanisms including collisions, internal damping and the Yarkovsky-O'Keefe-Radzievskii-Paddack (YORP) effect. We have analysed the changes in magnitude between consecutive detections of approximately 60,000 asteroids measured by the PanSTARRS 1 survey during its first 18 months of operations. We have attempted to explain the deri… ▽ More The rotational state of asteroids is controlled by various physical mechanisms including collisions, internal damping and the Yarkovsky-O'Keefe-Radzievskii-Paddack (YORP) effect. We have analysed the changes in magnitude between consecutive detections of approximately 60,000 asteroids measured by the PanSTARRS 1 survey during its first 18 months of operations. We have attempted to explain the derived brightness changes physically and through the application of a simple model. We have found a tendency toward smaller magnitude variations with decreasing diameter for objects of 1 < D < 8 km. Assuming the shape distribution of objects in this size range to be independent of size and composition our model suggests a population with average axial ratios 1 : 0.85 \pm 0.13 : 0.71 \pm 0.13, with larger objects more likely to have spin axes perpendicular to the orbital plane. △ Less

Submitted 15 April, 2016; originally announced April 2016.

Comments: 10 pages, 10 figures, accepted by MNRAS

First Results from the Rapid-Response Spectrophotometric Characterization of Near-Earth Objects using UKIRT

Authors: M. Mommert, D. E. Trilling, D. Borth, R. Jedicke, N. Butler, M. Reyes-Ruiz, B. Pichardo, E. Petersen, T. Axelrod, N. Moskovitz

Abstract: Using the Wide Field Camera for the United Kingdom Infrared Telescope, we measure the near-infrared colors of near-Earth objects (NEOs) in order to put constraints on their taxonomic classifications. The rapid-response character of our observations allows us to observe NEOs when they are close to the Earth and bright. Here we present near-infrared color measurements of 86 NEOs, most of which were… ▽ More Using the Wide Field Camera for the United Kingdom Infrared Telescope, we measure the near-infrared colors of near-Earth objects (NEOs) in order to put constraints on their taxonomic classifications. The rapid-response character of our observations allows us to observe NEOs when they are close to the Earth and bright. Here we present near-infrared color measurements of 86 NEOs, most of which were observed within a few days of their discovery, allowing us to characterize NEOs with diameters of only a few meters. Using machine-learning methods, we compare our measurements to existing asteroid spectral data and provide probabilistic taxonomic classifications for our targets. Our observations allow us to distinguish between S-complex, C/X-complex, D-type, and V-type asteroids. Our results suggest that the fraction of S-complex asteroids in the whole NEO population is lower than the fraction of ordinary chondrites in the meteorite fall statistics. Future data obtained with UKIRT will be used to investigate the significance of this discrepancy. △ Less

Submitted 18 February, 2016; originally announced February 2016.

Comments: 27 pages, 6 figures, accepted by AJ

Absolute magnitudes and slope parameters for 250,000 asteroids observed by Pan-STARRS PS1 - preliminary results

Authors: Peter Vereš, Robert Jedicke, Alan Fitzsimmons, Larry Denneau, Mikael Granvik, Bryce Bolin, Serge Chastel, Richard Wainscoat, William Burgett, Ken Chambers, Heather Flewelling, Nick Kaiser, Eugen Magnier, Jeff Morgan, Paul Price, John Tonry, Chris Waters

Abstract: We present the results of a Monte Carlo technique to calculate the absolute magnitudes (H) and slope parameters (G) of about 240000 asteroids observed by the Pan-STARRS1 telescope during the first 15 months of its 3-year all-sky survey mission. The system's exquisite photometry with photometric errors < 0.04 mag, and well-defined filter and photometric system, allowed us to derive accurate H and G… ▽ More We present the results of a Monte Carlo technique to calculate the absolute magnitudes (H) and slope parameters (G) of about 240000 asteroids observed by the Pan-STARRS1 telescope during the first 15 months of its 3-year all-sky survey mission. The system's exquisite photometry with photometric errors < 0.04 mag, and well-defined filter and photometric system, allowed us to derive accurate H and G even with a limited number of observations and restricted range in phase angles. Our Monte Carlo method simulates each asteroid's rotation period, amplitude and color to derive the most-likely H and G, but its major advantage is in estimating realistic statistical+systematic uncertainties and errors on each parameter. The method was tested by comparison with the well-established and accurate results for about 500 asteroids provided by Pravec et al. (2012) and then applied to determining H and G for the Pan-STARRS1 asteroids using both the Muinonen et al. (2010) and Bowell et al. (1989) phase functions. Our results confirm the bias in MPC photometry discovered by (Juric et al. 2002). △ Less

Submitted 25 August, 2015; v1 submitted 2 June, 2015; originally announced June 2015.

Comments: 13 figures, 7 tables

Surveys, Astrometric Follow-up & Population Statistics

Authors: Robert Jedicke, Mikael Granvik, Marco Micheli, Eileen Ryan, Timothy Spahr, Donald K. Yeomans

Abstract: Asteroid surveys are the backbone of asteroid science, and with this in mind we begin with a broad review of the impact of asteroid surveys on our field. We then provide a brief history of asteroid discoveries so as to place contemporary and future surveys in perspective. Surveys in the United States have discovered the vast majority of the asteroids and this dominance has been consolidated since… ▽ More Asteroid surveys are the backbone of asteroid science, and with this in mind we begin with a broad review of the impact of asteroid surveys on our field. We then provide a brief history of asteroid discoveries so as to place contemporary and future surveys in perspective. Surveys in the United States have discovered the vast majority of the asteroids and this dominance has been consolidated since the publication of Asteroids III. Our descriptions of the asteroid surveys that have been operational since that time are focussed upon those that have contributed the vast majority of asteroid observations and discoveries. We also provide some insight into upcoming next-generation surveys that are sure to alter our understanding of the small bodies in the inner solar system and provide evidence to untangle their complicated dynamical and physical histories. The Minor Planet Center, the nerve center of the asteroid discovery effort, has improved its operations significantly in the past decade so that it can manage the increasing discovery rate, and ensure that it is well-placed to handle the data rates expected in the next decade. We also consider the difficulties associated with astrometric follow-up of newly identified objects. It seems clear that both of these efforts must operate in new modes in order to keep pace with expected discovery rates of next-generation ground- and space-based surveys. △ Less

Submitted 14 March, 2015; originally announced March 2015.

Comments: Chapter to appear in the book ASTEROIDS IV, (University of Arizona Press) Space Science Series, edited by P. Michel, F. DeMeo and W. Bottke

Machine learning for transient discovery in Pan-STARRS1 difference imaging

Authors: D. E. Wright, S. J. Smartt, K. W. Smith, P. Miller, R. Kotak, A. Rest, W. S. Burgett, K. C. Chambers, H. Flewelling, K. W. Hodapp, M. Huber, R. Jedicke, N. Kaiser, N. Metcalfe, P. A. Price, J. L. Tonry, R. J. Wainscoat, C. Waters

Abstract: Efficient identification and follow-up of astronomical transients is hindered by the need for humans to manually select promising candidates from data streams that contain many false positives. These artefacts arise in the difference images that are produced by most major ground-based time domain surveys with large format CCD cameras. This dependence on humans to reject bogus detections is unsusta… ▽ More Efficient identification and follow-up of astronomical transients is hindered by the need for humans to manually select promising candidates from data streams that contain many false positives. These artefacts arise in the difference images that are produced by most major ground-based time domain surveys with large format CCD cameras. This dependence on humans to reject bogus detections is unsustainable for next generation all-sky surveys and significant effort is now being invested to solve the problem computationally. In this paper we explore a simple machine learning approach to real-bogus classification by constructing a training set from the image data of ~32000 real astrophysical transients and bogus detections from the Pan-STARRS1 Medium Deep Survey. We derive our feature representation from the pixel intensity values of a 20x20 pixel stamp around the centre of the candidates. This differs from previous work in that it works directly on the pixels rather than catalogued domain knowledge for feature design or selection. Three machine learning algorithms are trained (artificial neural networks, support vector machines and random forests) and their performances are tested on a held-out subset of 25% of the training data. We find the best results from the random forest classifier and demonstrate that by accepting a false positive rate of 1%, the classifier initially suggests a missed detection rate of around 10%. However we also find that a combination of bright star variability, nuclear transients and uncertainty in human labelling means that our best estimate of the missed detection rate is approximately 6%. △ Less

Submitted 3 March, 2015; v1 submitted 22 January, 2015; originally announced January 2015.

Comments: 17 pages, 16 figures, accepted for publication in MNRAS

The Main-Belt Comets: The Pan-STARRS1 Perspective

Authors: Henry H. Hsieh, Larry Denneau, Richard J. Wainscoat, Norbert Schorghofer, Bryce Bolin, Alan Fitzsimmons, Robert Jedicke, Jan Kleyna, Marco Micheli, Peter Veres, Nicholas Kaiser, Kenneth C. Chambers, William S. Burgett, Heather Flewelling, Klaus W. Hodapp, Eugene A. Magnier, Jeffrey S. Morgan, Paul A. Price, John L. Tonry, Christopher Waters

Abstract: We analyze 760475 observations of 333026 main-belt objects obtained by the Pan-STARRS1 (PS1) survey telescope between 2012 May 20 and 2013 November 9, a period during which PS1 discovered two main-belt comets, P/2012 T1 (PANSTARRS) and P/2013 R3 (Catalina-PANSTARRS). PS1 comet detection procedures currently consist of the comparison of the point spread functions (PSFs) of moving objects to those o… ▽ More We analyze 760475 observations of 333026 main-belt objects obtained by the Pan-STARRS1 (PS1) survey telescope between 2012 May 20 and 2013 November 9, a period during which PS1 discovered two main-belt comets, P/2012 T1 (PANSTARRS) and P/2013 R3 (Catalina-PANSTARRS). PS1 comet detection procedures currently consist of the comparison of the point spread functions (PSFs) of moving objects to those of reference stars, and the flagging of objects that show anomalously large radial PSF widths. Based on the number of missed discovery opportunities among comets discovered by other observers, we estimate an upper limit comet discovery efficiency rate of ~70% for PS1. Additional analyses that could improve comet discovery yields in future surveys include linear PSF analysis, modeling of trailed stellar PSFs for comparison to trailed moving object PSFs, searches for azimuthally localized activity, comparison of point-source-optimized photometry to extended-source-optimized photometry, searches for photometric excesses in objects with known absolute magnitudes, and crowd-sourcing. Analysis of PS1 survey statistics indicates an expected fraction of 59 MBCs per 10^6 outer main-belt asteroids, and a 95% confidence upper limit of 96 MBCs per 10^6 outer main-belt asteroids. We note that more sensitive future surveys could detect many more MBCs than estimated here. We find an excess of high eccentricities (0.1 < e < 0.3) among all known MBCs relative to the background asteroid population. Theoretical calculations show that, given these eccentricities, the sublimation rate for a typical MBC is orders of magnitude larger at perihelion than at aphelion, providing a plausible physical explanation for the observed behavior of MBCs peaking in observed activity strength near perihelion. △ Less

Submitted 19 October, 2014; originally announced October 2014.

Comments: 38 pages, 14 figures; accepted for publication in Icarus

Observational Constraints on the Catastrophic Disruption Rate of Small Main Belt Asteroids

Authors: Larry Denneau, Robert Jedicke, Alan Fitzsimmons, Henry Hsieh, Jan Kleyna, Mikael Granvik, Marco Micheli, T. Spahr, Peter Vereš, Richard Wainscoat, W. S. Burgett, K. C. Chambers, P. W. Draper, H. Flewelling, M. E. Huber, N. Kaiser, J. S. Morgan, J. L. Tonry

Abstract: We have calculated 90% confidence limits on the steady-state rate of catastrophic disruptions of main belt asteroids in terms of the absolute magnitude at which one catastrophic disruption occurs per year (HCL) as a function of the post-disruption increase in brightness (delta m) and subsequent brightness decay rate (tau). The confidence limits were calculated using the brightest unknown main belt… ▽ More We have calculated 90% confidence limits on the steady-state rate of catastrophic disruptions of main belt asteroids in terms of the absolute magnitude at which one catastrophic disruption occurs per year (HCL) as a function of the post-disruption increase in brightness (delta m) and subsequent brightness decay rate (tau). The confidence limits were calculated using the brightest unknown main belt asteroid (V = 18.5) detected with the Pan-STARRS1 (Pan-STARRS1) telescope. We measured the Pan-STARRS1's catastrophic disruption detection efficiency over a 453-day interval using the Pan-STARRS moving object processing system (MOPS) and a simple model for the catastrophic disruption event's photometric behavior in a small aperture centered on the catastrophic disruption event. Our simplistic catastrophic disruption model suggests that delta m = 20 mag and 0.01 mag d-1 < tau < 0.1 mag d-1 which would imply that H0 = 28 -- strongly inconsistent with H0,B2005 = 23.26 +/- 0.02 predicted by Bottke et al. (2005) using purely collisional models. We postulate that the solution to the discrepancy is that > 99% of main belt catastrophic disruptions in the size range to which this study was sensitive (100 m) are not impact-generated, but are instead due to fainter rotational breakups, of which the recent discoveries of disrupted asteroids P/2013 P5 and P/2013 R3 are probable examples. We estimate that current and upcoming asteroid surveys may discover up to 10 catastrophic disruptions/year brighter than V = 18.5. △ Less

Submitted 28 August, 2014; v1 submitted 28 August, 2014; originally announced August 2014.

Comments: 61 Pages, 10 Figures, 3 Tables

Search for the Return of Activity in Active Asteroid 176P/LINEAR

Authors: Henry H. Hsieh, Larry Denneau, Alan Fitzsimmons, Olivier R. Hainaut, Masateru Ishiguro, Robert Jedicke, Heather M. Kaluna, Jacqueline V. Keane, Jan Kleyna, Pedro Lacerda, Eric M. MacLennan, Karen J. Meech, Nick A. Moskovitz, Timm Riesen, Eva Schunova, Colin Snodgrass, Chadwick A. Trujillo, Laurie Urban, Peter Veres, Richard J. Wainscoat, Bin Yang

Abstract: We present the results of a search for the reactivation of active asteroid 176P/LINEAR during its 2011 perihelion passage using deep optical observations obtained before, during, and after that perihelion passage. Deep composite images of 176P constructed from data obtained between June 2011 and December 2011 show no visible signs of activity, while photometric measurements of the object during th… ▽ More We present the results of a search for the reactivation of active asteroid 176P/LINEAR during its 2011 perihelion passage using deep optical observations obtained before, during, and after that perihelion passage. Deep composite images of 176P constructed from data obtained between June 2011 and December 2011 show no visible signs of activity, while photometric measurements of the object during this period also show no significant brightness enhancements similar to that observed for 176P between November 2005 and December 2005 when it was previously observed to be active. An azimuthal search for dust emission likewise reveals no evidence for directed emission (i.e., a tail, as was previously observed for 176P), while a one-dimensional surface brightness profile analysis shows no indication of a spherically symmetric coma at any time in 2011. We conclude that 176P did not in fact exhibit activity in 2011, at least not on the level on which it exhibited activity in 2005, and suggest that this could be due to the devolatization or mantling of the active site responsible for its activity in 2005. △ Less

Submitted 20 August, 2014; originally announced August 2014.

Comments: 17 pages, 9 figures; published in AJ

Journal ref: Astronomical Journal (2014), 147, 89

Detecting Earth's Temporarily-Captured Natural Satellites - Minimoons

Authors: Bryce Bolin, Robert Jedicke, Mikael Granvik, Peter Brown, Ellen Howell, Michael C. Nolan, Peter Jenniskens, Monique Chyba, Geoff Patterson, Richard Wainscoat

Abstract: We present a study on the discoverability of temporarily captured orbiters (TCOs) by present day or near-term anticipated ground-based and space-based facilities. TCOs (Granvik et al. 2012) are potential targets for spacecraft ren- dezvous or human exploration (Chyba et al. 2014) and provide an opportunity to study the population of the smallest asteroids in the solar system. We find that present… ▽ More We present a study on the discoverability of temporarily captured orbiters (TCOs) by present day or near-term anticipated ground-based and space-based facilities. TCOs (Granvik et al. 2012) are potential targets for spacecraft ren- dezvous or human exploration (Chyba et al. 2014) and provide an opportunity to study the population of the smallest asteroids in the solar system. We find that present day ground-based optical surveys such as Pan-STARRS and ATLAS can discover the largest TCOs over years of operation. A targeted survey conducted with the Subaru telescope can discover TCOs in the 0.5 m to 1.0 m diameter size range in about 5 nights of observing. Furthermore, we discuss the application of space-based infrared surveys, such as NEOWISE, and ground-based meteor detection systems such as CAMS, CAMO and ASGARD in discovering TCOs. These systems can detect TCOs but at a uninteresting rate. Finally, we discuss the application of bi-static radar at Arecibo and Green Bank to discover TCOs. Our radar simulations are strongly dependent on the rotation rate distribution of the smallest asteroids but with an optimistic distribution we find that these systems have > 80% chance of detecting a > 10 cm diameter TCO in about 40 h of operation. △ Less

Submitted 16 June, 2014; v1 submitted 10 June, 2014; originally announced June 2014.

Comments: 44 Pages, 14 Figures, 1 Table

Properties and evolution of NEO families created by tidal disruption at Earth

Authors: Eva Schunová, Robert Jedicke, Kevin J. Walsh, Mikael Granvik, Richard J. Wainscoat, Nader Haghighipour

Abstract: We have calculated the coherence and detectable lifetimes of synthetic near-Earth object (NEO) families created by catastrophic disruption of a progenitor as it suffers a very close Earth approach. The closest or slowest approaches yield the most violent `s-class' disruption events. We found that the average slope of the absolute magnitude (H) distribution, $N(H)\propto10^{(0.55\pm0.04)\,H}$, for… ▽ More We have calculated the coherence and detectable lifetimes of synthetic near-Earth object (NEO) families created by catastrophic disruption of a progenitor as it suffers a very close Earth approach. The closest or slowest approaches yield the most violent `s-class' disruption events. We found that the average slope of the absolute magnitude (H) distribution, $N(H)\propto10^{(0.55\pm0.04)\,H}$, for the fragments in the s-class families is steeper than the slope of the NEO population \citep{mainzer2011} in the same size range. The families remain coherent as statistically significant clusters of orbits within the NEO population for an average of $\barτ_c = (14.7\pm0.6)\times10^3$ years after disruption. The s-class families are detectable with the techniques developed by \citet{fu2005} and \citet{Schunova2012} for an average duration ($\barτ_{det}$) ranging from about 2,000 to about 12,000 years for progenitors in the absolute magnitude ($H_p$) range from 20 to 13 corresponding to diameters in the range from about 0.5 to 10$\km$ respectively. The short detectability lifetime explains why zero NEO families have been discovered to-date. Nonetheless, every tidal disruption event of a progenitor with D$>0.5\km$ is capable of producing several million fragments in the $1\meter$ to $10\meter$ diameter range that can contribute to temporary local density enhancements of small NEOs in Earth's vicinity. We expect that there are about 1,200 objects in the steady state NEO population in this size range due to tidal disruption assuming that one $1\km$ diameter NEO tidally disrupts at Earth every 2,500 years. These objects may be suitable targets for asteroid retrieval missions due to their Earth-like orbits with corresponding low $v_{\infty}$. The fragments from the tidal disruptions at Earth have $\sim5\times$ the collision probability with Earth compared to the background NEO population. △ Less

Submitted 16 May, 2014; originally announced May 2014.

Comments: 60 pages, 16 figures

MSC Class: 85