Source author record

Larry Denneau

Larry Denneau appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

24works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

24 published item(s)

preprint2021arXiv

A Two-Stage Deep Learning Detection Classifier for the ATLAS Asteroid Survey

In this paper we present a two-step neural network model to separate detections of solar system objects from optical and electronic artifacts in data obtained with the "Asteroid Terrestrial-impact Last Alert System" (ATLAS), a near-Earth asteroid sky survey system [arXiv:1802.00879]. A convolutional neural network [arXiv:1807.10912] is used to classify small "postage-stamp" images of candidate detections of astronomical sources into eight classes, followed by a multi-layered perceptron that provides a probability that a temporal sequence of four candidate detections represents a real astronomical source. The goal of this work is to reduce the time delay between Near-Earth Object (NEO) detections and submission to the Minor Planet Center. Due to the rare and hazardous nature of NEOs [Harris and D'Abramo, 2015], a low false negative rate is a priority for the model. We show that the model reaches 99.6\% accuracy on real asteroids in ATLAS data with a 0.4\% false negative rate. Deployment of this model on ATLAS has reduced the amount of NEO candidates that astronomers must screen by 90%, thereby bringing ATLAS one step closer to full autonomy.

preprint2021arXiv

Characterizing the Manx Candidate A/2018 V3

Manx objects approach the inner solar system on long-period comet (LPC) orbits with the consequent high inbound velocities, but unlike comets, Manxes display very little to no activity even near perihelion. This suggests that they may have formed in circumstances different from typical LPCs; moreover, this lack of significant activity also renders them difficult to detect at large distances. Thus, analyzing their physical properties can help constrain models of solar system formation as well as sharpen detection methods for those classified as NEOs. Here, we focus on the Manx candidate A/2018 V3 as part of a larger effort to characterize Manxes as a whole. This particular object was observed to be inactive even at its perihelion at $q$ = 1.34 au in 2019 September. Its spectral reflectivity is consistent with typical organic-rich comet surfaces with colors of $g'-r'= 0.67\pm0.02$, $r'-i' = 0.26\pm0.02$, and $r'-z' = 0.45\pm0.02$, corresponding to a spectral reflectivity slope of $10.6\pm 0.9$ %/100nm. A least-squares fit of our constructed light curve to the observational data yields an average nucleus radius of $\approx$2 km assuming an albedo of 0.04. This is consistent with the value measured from NEOWISE. A surface brightness analysis for data taken 2020 July 13 indicated possible low activity ($\lesssim0.68$ g $\rm s^{-1}$), but not enough to lift optically significant amounts of dust. Finally, we discuss Manxes as a constraint on solar system dynamical models as well as their implications for planetary defense.

preprint2020arXiv

Beyond Gaia: Asteroseismic Distances of M giants using Ground-Based Transient Surveys

Evolved stars near the tip of the red giant branch (TRGB) show solar-like oscillations with periods spanning hours to months and amplitudes ranging from $\sim$1 mmag to $\sim$100 mmag. The systematic detection of the resulting photometric variations with ground-based telescopes would enable the application of asteroseismology to a much larger and more distant sample of stars than is currently accessible with space-based telescopes such as \textit{Kepler} or the ongoing Transiting Exoplanet Survey Satellite (\textit{TESS}) mission. We present an asteroseismic analysis of 493 M giants using data from two ground-based surveys: the Asteroid Terrestrial-impact Last Alert System (ATLAS) and the All-Sky Automated Survey for Supernovae (ASAS-SN). By comparing the extracted frequencies with constraints from \textit{Kepler}, the Sloan Digital Sky Survey Apache Point Observatory Galaxy Evolution Experiment (APOGEE), and Gaia we demonstrate that ground-based transient surveys allow accurate distance measurements to oscillating M giants with a precision of $\sim$15$\%$. Using stellar population synthesis models we predict that ATLAS and ASAS-SN can provide asteroseismic distances to $\sim$2$\times$10$^{6}$ galactic M giants out to typical distances of $20-50 \; \rm{kpc}$, vastly improving the reach of Gaia and providing critical constraints for Galactic archaeology and galactic dynamics.

preprint2020arXiv

The Outburst of the Young Star Gaia19bey

We report photometry and spectroscopy of the outburst of the young stellar object Gaia19bey. We have established the outburst light curve with archival Gaia G, ATLAS Orange, ZTF r-band and Pan-STARRS rizy-filter photometry, showing an outburst of approximately 4 years duration, longer than typical EXors but shorter than FUors. Its pre-outburst SED shows a flat far-infrared spectrum, confirming the early evolutionary state of Gaia19bey and its similarity to other deeply embedded young stars experiencing outbursts. A lower limit to the peak outburst luminosity is approximately 182 L_sun at an assumed distance of 1.4 kpc, the minimum plausible distance. Infrared and optical spectroscopy near maximum light showed an emission line spectrum, including HI lines, strong red CaII emission, other metal emission lines, infrared CO bandhead emission, and a strong infrared continuum. Towards the end of the outburst, the emission lines have all but disappeared and the spectrum has changed into an almost pure continuum spectrum. This indicates a cessation of magnetospheric accretion activity. The near-infrared colors have become redder as Gaia19bey has faded, indicating a cooling of the continuum component. Near the end of the outburst, the only remaining strong emission lines are forbidden shock-excited emission lines. Adaptive optics integral field spectroscopy shows the H_2 1--0 S(1) emission with the morphology of an outflow cavity and the extended emission in the [FeII] line at 1644 nm with the morphology of an edge-on disk. However, we do not detect any large-scale jet from Gaia19bey.

preprint2020arXiv

The Transient Jupiter Trojan-Like Orbit of P/2019 LD2 (ATLAS)

Comet P/2019 LD2 has orbital elements currently resembling those of a Jupiter Trojan, and therefore superficially appears to represent a unique opportunity to study the volatile content and active behavior of a member of this population for the first time. However, numerical integrations show that it was previously a Centaur before reaching its current Jupiter Trojan-like orbit in 2018 July, and is expected to return to being a Centaur in 2028 February, before eventually becoming a Jupiter-family comet in 2063 February. The case of P/2019 LD2 highlights the need for mechanisms to quickly and reliably dynamically classify small solar system bodies discovered in current and upcoming wide-field surveys.

preprint2019arXiv

The New EXor Outburst of ESO-H-alpha~99 observed by Gaia ATLAS and TESS

We report photometry and spectroscopy of the outburst of the young stellar object ESO-Halpha 99. The outburst was first noticed in Gaia alert Gaia18dvc and later by ATLAS. We have established the outburst light curve with archival ATLAS ``Orange'' filter photometry, Gaia data, new V-band photometry, and J, H, and K_s photometry from IRIS and UKIRT. The brightness has fluctuated several times near the light curve maximum. The TESS satellite observed ESO-Halpha 99 with high cadence during one of these minor minima and found brightness fluctuations on timescales of days and hours. Imaging with UKIRT shows the outline of an outflow cavity, and we find one knot of H_2 1-0 S(1) emission, now named MHO 1520, on the symmetry axis of this nebula, indicating recent collimated outflow activity from ESO-Halpha 99. Its pre-outburst SED shows a flat FIR spectrum, confirming its early evolutionary state and its similarity to other deeply embedded objects in the broader EXor class. The pre-outburst luminosity is 34 L_sun, a much higher luminosity than typical EXors, indicating that ESO-Halpha 99 may be a star of intermediate mass. Infrared and optical spectroscopy show a rich emission line spectrum, including HI lines, strong red CaII emission, as well as infrared CO bandhead emission, all characteristic EXors in the broadest sense. Comparison of the present spectra with an optical spectrum obtained in 1993, presumably in the quiescent state of the object, shows that during the present outburst the continuum component of the spectrum has increased notably more than the emission lines. The Halpha equivalent width during the outburst is down to one half of its 1993 level, and shock-excited emission lines are much less prominent.

preprint2016arXiv

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

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.

preprint2015arXiv

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

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).

preprint2015arXiv

Hypercalibration: A Pan-STARRS1-based recalibration of the Sloan Digital Sky Survey

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

preprint2015arXiv

Sublimation-Driven Activity in Main-Belt Comet 313P/Gibbs

We present an observational and dynamical study of newly discovered main-belt comet 313P/Gibbs. We find that the object is clearly active both in observations obtained in 2014 and in precovery observations obtained in 2003 by the Sloan Digital Sky Survey, strongly suggesting that its activity is sublimation-driven. This conclusion is supported by a photometric analysis showing an increase in the total brightness of the comet over the 2014 observing period, and dust modeling results showing that the dust emission persists over at least three months during both active periods, where we find start dates for emission no later than 2003 July 24+/-10 for the 2003 active period and 2014 July 28+/-10 for the 2014 active period. From serendipitous observations by the Subaru Telescope in 2004 when the object was apparently inactive, we estimate that the nucleus has an absolute R-band magnitude of H_R=17.1+/-0.3, corresponding to an effective nucleus radius of r_e~1.00+/-0.15 km. The object's faintness at that time means we cannot rule out the presence of activity, and so this computed radius should be considered an upper limit. We find that 313P's orbit is intrinsically chaotic, having a Lyapunov time of T_l=12000 yr and being located near two 3-body mean-motion resonances with Jupiter and Saturn, 11J-1S-5A and 10J+12S-7A, yet appears stable over >50 Myr in an apparent example of stable chaos. We furthermore find that 313P is the second main-belt comet, after P/2012 T1 (PANSTARRS), to belong to the ~155 Myr old Lixiaohua asteroid family.

preprint2014arXiv

A new distant Milky Way globular cluster in the Pan-STARRS1 3π survey

We present a new satellite in the outer halo of the Galaxy, the first Milky Way satellite found in the stacked photometric catalog of the Panoramic Survey Telescope and Rapid Response System 1 (Pan-STARRS1) Survey. From follow-up photometry obtained with WFI on the MPG/ESO 2.2m telescope, we argue that the object, located at a heliocentric distance of 145+/-17 kpc, is the most distant Milky Way globular cluster yet known. With a total magnitude of M_V=-4.3+/-0.2 and a half-light radius of 20+/-2 pc, it shares the properties of extended globular clusters found in the outer halo of our Galaxy and the Andromeda galaxy. The discovery of this distant cluster shows that the full spatial extent of the Milky Way globular cluster system has not yet been fully explored.

preprint2014arXiv

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

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.

preprint2014arXiv

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

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.

preprint2014arXiv

The Main-Belt Comets: The Pan-STARRS1 Perspective

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.

preprint2013arXiv

Main-Belt Comet P/2012 T1 (PANSTARRS)

We present initial results from observations and numerical analyses aimed at characterizing main-belt comet P/2012 T1 (PANSTARRS). Optical monitoring observations were made between October 2012 and February 2013 using the University of Hawaii 2.2 m telescope, the Keck I telescope, the Baade and Clay Magellan telescopes, Faulkes Telescope South, the Perkins Telescope at Lowell Observatory, and the Southern Astrophysical Research (SOAR) telescope. The object's intrinsic brightness approximately doubles from the time of its discovery in early October until mid-November and then decreases by ~60% between late December and early February, similar to photometric behavior exhibited by several other main-belt comets and unlike that exhibited by disrupted asteroid (596) Scheila. We also used Keck to conduct spectroscopic searches for CN emission as well as absorption at 0.7 microns that could indicate the presence of hydrated minerals, finding an upper limit CN production rate of QCN<1.5x10^23 mol/s, from which we infer a water production rate of QH2O<5x10^25 mol/s, and no evidence of the presence of hydrated minerals. Numerical simulations indicate that P/2012 T1 is largely dynamically stable for >100 Myr and is unlikely to be a recently implanted interloper from the outer solar system, while a search for potential asteroid family associations reveal that it is dynamically linked to the ~155 Myr-old Lixiaohua asteroid family.

preprint2013arXiv

NEOKepler: Discovering Near-Earth Objects Using the Kepler Spacecraft

We propose a new Kepler mission, called NEOKepler, that would survey near Earth's orbit to identify potentially hazardous objects (PHOs). To understand its surveying power, Kepler's large field of view produces an 'etendue' (A*Omega) that is 4.5 times larger than the best survey telescope currently in operation. In this paper, we investigate the feasibility of NEOKepler using a double "fence post" survey pattern that efficiently detects PHOs. In a simulated 12-month survey, we estimate that NEOKepler would detect ~150 new NEOs with absolute magnitudes of less than 21.5, ~50 of which would be new PHOs. This would increase the annual PHO discovery rate by at least 50% and improve upon our goal of discovering 90% of PHOs by the end of 2020. Due to its heliocentric orbit, Kepler would also be sensitive to objects inside Earth's orbit, discovering more objects in its first year than are currently known to exist. Understanding this undersampled sub-population of NEOs will reveal new insights into the actual PHO distribution by further constraining current NEO models. As an alternative science goal, NEOKepler could employ a different observing strategy to discover suitable targets for NASA's Asteroid Redirect Mission.

preprint2013arXiv

Outgassing Behavior of C/2012 S1 (ISON) From September 2011 to June 2013

We report photometric observations for comet C/2012 S1 (ISON) obtained during the time period immediately after discovery (r=6.28 AU) until it moved into solar conjunction in mid-2013 June using the UH2.2m, and Gemini North 8-m telescopes on Mauna Kea, the Lowell 1.8m in Flagstaff, the Calar Alto 1.2m telescope in Spain, the VYSOS-5 telescopes on Mauna Loa Hawaii and data from the CARA network. Additional pre-discovery data from the Pan STARRS1 survey extends the light curve back to 2011 September 30 (r=9.4 AU). The images showed a similar tail morphology due to small micron sized particles throughout 2013. Observations at sub-mm wavelengths using the JCMT on 15 nights between 2013 March 9 (r=4.52 AU) and June 16 (r=3.35 AU) were used to search for CO and HCN rotation lines. No gas was detected, with upper limits for CO ranging between (3.5-4.5)E27 molec/s. Combined with published water production rate estimates we have generated ice sublimation models consistent with the photometric light curve. The inbound light curve is likely controlled by sublimation of CO2. At these distances water is not a strong contributor to the outgassing. We also infer that there was a long slow outburst of activity beginning in late 2011 peaking in mid-2013 January (r~5 AU) at which point the activity decreased again through 2013 June. We suggest that this outburst was driven by CO injecting large water ice grains into the coma. Observations as the comet came out of solar conjunction seem to confirm our models.

preprint2013arXiv

The Pan-STARRS Moving Object Processing System

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

preprint2012arXiv

Discovery of Main-Belt Comet P/2006 VW139 by Pan-STARRS1

Main belt asteroid (300163) 2006 VW139 (later designated P/2006 VW139) was discovered to exhibit comet-like activity by the Pan-STARRS1 survey telescope using automated point-spread-function analyses performed by PS1's Moving Object Processing System. Deep follow-up observations show both a short (\sim 10") antisolar dust tail and a longer (\sim 60") dust trail aligned with the object's orbit plane, similar to the morphology observed for another main-belt comet, P/2010 R2 (La Sagra), and other well-established comets, implying the action of a long-lived, sublimation-driven emission event. Photometry showing the brightness of the near-nucleus coma remaining constant over \sim 30 days provides further evidence for this object's cometary nature, suggesting it is in fact a main-belt comet, and not a disrupted asteroid. A spectroscopic search for CN emission was unsuccessful, though we find an upper limit CN production rate of Q_CN < 1.3x10^24 mol/s, from which we infer a water production rate of Q_H2O < 10^26 mol/s. We also find an approximately linear optical spectral slope of 7.2%/1000A, similar to other cometary dust comae. Numerical simulations indicate that P/2006 VW139 is dynamically stable for > 100 Myr, while a search for a potential asteroid family around the object reveals a cluster of 24 asteroids within a cutoff distance of 68 m/s. At 70 m/s, this cluster merges with the Themis family, suggesting that it could be similar to the Beagle family to which another main-belt comet, 133P/Elst-Pizarro, belongs.

preprint2012arXiv

Identification of known objects in solar system surveys

The discovery of new objects in modern wide-field asteroid and comet surveys can be enhanced by first identifying observations belonging to known solar system objects. The assignation of new observations to a known object is an attribution problem that occurs when a least squares orbit already exists for the object but a separate fit is not possible to just the set of new observations. In this work we explore the strongly asymmetric attribution problem in which the existing least squares orbit is very well constrained and the new data are sparse. We describe an attribution algorithm that introduces new quality control metrics in the presence of strong biases in the astrometric residuals. The main biases arise from the stellar catalogs used in the reduction of asteroid observations and we show that a simple debiasing with measured regional catalog biases significantly improves the results. We tested the attribution algorithm using data from the PS1 survey that used the 2MASS star catalog for the astrometric reduction. We found small but statistically significant biases in the data of up to 0.1 arcsec that are relevant only when the observations reach the level of accuracy made possible by instruments like PS1. The false attribution rate was measured to be < 1/1,000 with a simple additional condition that can reduce it to zero while the attribution efficiency is consistent with 100%.

preprint2012arXiv

Improved Asteroid Astrometry and Photometry with Trail Fitting

Asteroid detections in astronomical images may appear as trails due to a combination of their apparent rate of motion and exposure duration. Nearby asteroids in particular typically have high apparent rates of motion and acceleration. Their recovery, especially on their discovery apparition, depends upon obtaining good astrometry from the trailed detections. We present an analytic function describing a trailed detection under the assumption of a Gaussian point spread function (PSF) and constant rate of motion. We have fit the function to both synthetic and real trailed asteroid detections from the Pan-STARRS1 survey telescope to obtain accurate astrometry and photometry. For short trails our trailing function yields the same astrometric and photometry accuracy as a functionally simpler 2-d Gaussian but the latter underestimates the length of the trail - a parameter that can be important for measuring the object's rate of motion and assessing its cometary activity. For trails longer than about 10 pixels (> 3xPSF) our trail fitting provides 3-times better astrometric accuracy and up to 2 magnitudes improvement in the photometry. The trail fitting algorithm can be implemented at the source detection level for all detections to provide trail length and position angle that can be used to reduce the false tracklet rate.

preprint2012arXiv

Observational and Dynamical Characterization of Main-Belt Comet P/2010 R2 (La Sagra)

We present observations of comet-like main-belt object P/2010 R2 (La Sagra) obtained by Pan-STARRS 1 and the Faulkes Telescope-North on Haleakala in Hawaii, the University of Hawaii 2.2 m, Gemini-North, and Keck I telescopes on Mauna Kea, the Danish 1.54 m telescope at La Silla, and the Isaac Newton Telescope on La Palma. An antisolar dust tail is observed from August 2010 through February 2011, while a dust trail aligned with the object's orbit plane is also observed from December 2010 through August 2011. Assuming typical phase darkening behavior, P/La Sagra is seen to increase in brightness by >1 mag between August 2010 and December 2010, suggesting that dust production is ongoing over this period. These results strongly suggest that the observed activity is cometary in nature (i.e., driven by the sublimation of volatile material), and that P/La Sagra is therefore the most recent main-belt comet to be discovered. We find an approximate absolute magnitude for the nucleus of H_R=17.9+/-0.2 mag, corresponding to a nucleus radius of ~0.7 km, assuming an albedo of p=0.05. Using optical spectroscopy, we find no evidence of sublimation products (i.e., gas emission), finding an upper limit CN production rate of Q_CN<6x10^23 mol/s, from which we infer an H2O production rate of Q_H2O<10^26 mol/s. Numerical simulations indicate that P/La Sagra is dynamically stable for >100 Myr, suggesting that it is likely native to its current location and that its composition is likely representative of other objects in the same region of the main belt, though the relatively close proximity of the 13:6 mean-motion resonance with Jupiter and the (3,-2,-1) three-body mean-motion resonance with Jupiter and Saturn mean that dynamical instability on larger timescales cannot be ruled out.

preprint2009arXiv

Detection of Earth-impacting asteroids with the next generation all-sky surveys

We have performed a simulation of a next generation sky survey's (Pan-STARRS 1) efficiency for detecting Earth-impacting asteroids. The steady-state sky-plane distribution of the impactors long before impact is concentrated towards small solar elongations (Chesley and Spahr, 2004) but we find that there is interesting and potentially exploitable behavior in the sky-plane distribution in the months leading up to impact. The next generation surveys will find most of the dangerous impactors (>140m diameter) during their decade-long survey missions though there is the potential to miss difficult objects with long synodic periods appearing in the direction of the Sun, as well as objects with long orbital periods that spend much of their time far from the Sun and Earth. A space-based platform that can observe close to the Sun may be needed to identify many of the potential impactors that spend much of their time interior to the Earth's orbit. The next generation surveys have a good chance of imaging a bolide like 2008TC3 before it enters the atmosphere but the difficulty will lie in obtaining enough images in advance of impact to allow an accurate pre-impact orbit to be computed.

preprint2009arXiv

The Thousand Asteroid Light Curve Survey

We present the results of our Thousand Asteroid Light Curve Survey (TALCS) conducted with the Canada-France-Hawaii Telescope in September 2006. Our untargeted survey detected 828 Main Belt asteroids to a limiting magnitude of g'~22.5 corresponding to a diameter range of 0.4 km <= D <= 10 km. Of these, 278 objects had photometry of sufficient quality to perform rotation period fits. We debiased the observations and light curve fitting process to determine the true distribution of rotation periods and light curve amplitudes of Main Belt asteroids. We confirm a previously reported excess in the fraction of fast rotators but find a much larger excess of slow rotating asteroids (~15% of our sample). A few percent of objects in the TALCS size range have large light curve amplitudes of ~1 mag. Fits to the debiased distribution of light curve amplitudes indicate that the distribution of triaxial ellipsoid asteroid shapes is proportional to the square of the axis-ratio, (b/a)^2, and may be bi-modal. Finally, we find six objects with rotation periods that may be less than 2 hours with diameters between 400 m and 1.5 km, well above the break-up limit for a gravitationally-bound aggregate. Our debiased data indicate that this population represents <4% of the Main Belt in the 1-10 km size range.