Showing 1–44 of 44 results for author: Bolin, B

AT2019pim: A Luminous Orphan Afterglow from a Moderately Relativistic Outflow

Authors: Daniel A. Perley, Anna Y. Q. Ho, Michael Fausnaugh, Gavin P. Lamb, Mansi M. Kasliwal, Tomas Ahumada, Shreya Anand, Igor Andreoni, Eric Bellm, Varun Bhalerao, Bryce Bolin, Thomas G. Brink, Eric Burns, S. Bradley Cenko, Alessandra Corsi, Alexei V. Filippenko, Dmitry Frederiks, Adam Goldstein, Rachel Hamburg, Rahul Jayaraman, Peter G. Jonker, Erik C. Kool, Shrinivas Kulkarni, Harsh Kumar, Russ Laher , et al. (12 additional authors not shown)

Abstract: Classical gamma-ray bursts (GRBs) have two distinct emission episodes: prompt emission from ultra-relativistic ejecta and afterglow from shocked circumstellar material. While both components are extremely luminous in known GRBs, a variety of scenarios predict the existence of luminous afterglow emission with little or no associated high-energy prompt emission. We present AT 2019pim, the first secu… ▽ More Classical gamma-ray bursts (GRBs) have two distinct emission episodes: prompt emission from ultra-relativistic ejecta and afterglow from shocked circumstellar material. While both components are extremely luminous in known GRBs, a variety of scenarios predict the existence of luminous afterglow emission with little or no associated high-energy prompt emission. We present AT 2019pim, the first secure example of this phenomenon to be identified. Serendipitously discovered during follow-up observations of a gravitational-wave trigger and located in a contemporaneous TESS sector, it is hallmarked by a fast-rising (t ~ 2 hr), luminous (M_UV,peak ~ -24.4 mag) optical transient with accompanying luminous X-ray and radio emission. No gamma-ray emission consistent with the time and location of the transient was detected by Fermi-GBM or by Konus, placing strong limits on an accompanying GRB. We investigate several independent observational aspects of the afterglow in the context of constraints on relativistic motion and find all of them are consistent with an initial Lorentz factor of Gamma_0 ~ 30-50, significantly lower than in any well-observed GRB and consistent with the theoretically-predicted "dirty fireball" scenario in which the high-energy prompt emission is stifled by pair production. However, we cannot rule out a structured jet model in which only the line-of-sight material was ejected at low-Gamma, off-axis from a classical high-Gamma jet core. This event represents a milestone in orphan afterglow searches, demonstrating that luminous afterglows with weak or no detectable gamma-ray radiation exist in nature and can be discovered by high-cadence optical surveys. △ Less

Submitted 29 January, 2024; originally announced January 2024.

Comments: Submitted to MNRAS

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

Palomar discovery and initial characterization of naked-eye long period comet C/2022 E3 (ZTF)

Authors: B. T. Bolin, F. J. Masci, D. A. Duev, J. W. Milburn, J. N. Purdum, C. Avdellidou, Y. -C. Cheng, M. Delbo, C. Fremling, M. Ghosal, Z. -Y. Lin, C. M. Lisse, A. Mahabal, M. Saki

Abstract: Long-period comets are planetesimal remnants constraining the environment and volatiles of the protoplanetary disc. We report the discovery of hyperbolic long-period comet C/2022 E3 (ZTF), which has a perihelion $\sim$1.11 au, an eccentricity $\gtrsim$1 and an inclination $\sim$109$^{\circ}$, from images taken with the Palomar 48-inch telescope during morning twilight on 2022 Mar 2. Additionally,… ▽ More Long-period comets are planetesimal remnants constraining the environment and volatiles of the protoplanetary disc. We report the discovery of hyperbolic long-period comet C/2022 E3 (ZTF), which has a perihelion $\sim$1.11 au, an eccentricity $\gtrsim$1 and an inclination $\sim$109$^{\circ}$, from images taken with the Palomar 48-inch telescope during morning twilight on 2022 Mar 2. Additionally, we report the characterization of C/2022 E3 (ZTF) from observations taken with the Palomar 200-inch, the Palomar 60-inch, and the NASA Infrared Telescope Facility in early 2023 February to 2023 March when the comet passed within $\sim$0.28 au of the Earth and reached a visible magnitude of $\sim$5. We measure g-r = 0.70$\pm$0.01, r-i = 0.20$\pm$0.01, i-z = 0.06$\pm$0.01, z-J = 0.90$\pm$0.01, J-H = 0.38$\pm$0.01 and H-K = 0.15$\pm$0.01 colours for the comet from observations. We measure the A(0$^\circ$)f$ρ$ (0.8~$μ$m) in a 6500~km radius from the nucleus of 1483$\pm$40~cm, and CN, C$_3$, and C$_2$ production of 5.43$\pm0.11\times$10$^{25}$~mol/s, 2.01$\pm0.04\times$10$^{24}$, and 3.08$\pm0.5\times$10$^{25}$~mol/s, similar to other long period comets. We additionally observe the appearance of jet-like structures at a scale of $\sim$4,000 km in wide-field g-band images, which may be caused by the presence of CN gas in the near-nucleus coma. △ Less

Submitted 25 September, 2023; originally announced September 2023.

Comments: Accepted for publication in MNRAS:L, 9 pages, 6 figures, 2 tables

Astreaks: Astrometry of NEOs with trailed background stars

Authors: Kritti Sharma, Harsh Kumar, Harsh Choudhary, Varun Bhalerao, Vishwajeet Swain, Bryce Bolin, G. C. Anupama, Sudhanshu Barway, Simran Joharle, Vedant Shenoy

Abstract: The detection and accurate astrometry of fast-moving near-Earth objects (NEOs) has been a challenge for the follow-up community. Their fast apparent motion results in streaks in sidereal images, thus affecting the telescope's limiting magnitude and astrometric accuracy. A widely adopted technique to mitigate trailing losses is non-sidereal tracking, which transfers the streaking to background refe… ▽ More The detection and accurate astrometry of fast-moving near-Earth objects (NEOs) has been a challenge for the follow-up community. Their fast apparent motion results in streaks in sidereal images, thus affecting the telescope's limiting magnitude and astrometric accuracy. A widely adopted technique to mitigate trailing losses is non-sidereal tracking, which transfers the streaking to background reference stars. However, no existing publicly available astrometry software is configured to detect such elongated stars. We present Astreaks, a streaking source detection algorithm, to obtain accurate astrometry of NEOs in non-sidereal data. We validate the astrometric accuracy of Astreaks on 371 non-sidereally tracked images for 115 NEOs with two instrument set-ups of the GROWTH-India Telescope. The observed NEOs had V-band magnitude in the range [15, 22] with proper motion up to 140$^{\prime\prime}$/min, thus resulting in stellar streaks as high as 6.5$^\prime$ (582 pixels) in our data. Our method obtained astrometric solutions for all images with 100% success rate. The standard deviation in Observed-minus-Computed (O-C) residuals is 0.52$^{\prime\prime}$, with O-C residuals <2$^{\prime\prime}$(<1$^{\prime\prime}$) for 98.4% (84.4%) of our measurements. These are appreciable, given the pixel scale of $\sim$0.3$^{\prime\prime}$ and $\sim$0.7$^{\prime\prime}$ of our two instrument set-ups. This demonstrates that our modular and fully-automated algorithm helps improve the telescope system's limiting magnitude without compromising astrometric accuracy by enabling non-sidereal tracking on the target. This will help the NEO follow-up community cope with the accelerated discovery rates and improved sensitivity of the next-generation NEO surveys. Astreaks has been made available to the community under an open-source license. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 10 pages, 7 figures

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

Keck and Gemini spectral characterization of Lucy mission fly-by target (152830) Dinkinesh

Authors: B. T. Bolin, K. S. Noll, I. Caiazzo, C. Fremling, R. P. Binzel

Abstract: Recently, the inner main belt asteroid (152830) Dinkinesh was identified as an additional fly-by target for the Lucy mission. The heliocentric orbit and approximate absolute magnitude of Dinkinesh are known, but little additional information was available prior to its selection as a target. In particular, the lack of color spectrophotometry or spectra made it impossible to assign a spectral type t… ▽ More Recently, the inner main belt asteroid (152830) Dinkinesh was identified as an additional fly-by target for the Lucy mission. The heliocentric orbit and approximate absolute magnitude of Dinkinesh are known, but little additional information was available prior to its selection as a target. In particular, the lack of color spectrophotometry or spectra made it impossible to assign a spectral type to Dinkinesh from which its albedo could be estimated. We set out to remedy this knowledge gap by obtaining visible wavelength spectra with the Keck telescope on 2022 November 23 and with Gemini-South on 2022 December 27. The spectra measured with the Keck I/Low Resolution Imaging Spectrometer (LRIS) and the Gemini South/Gemini Multi-Object Spectrograph South (GMOS-S) are most similar to the average spectrum of S- and Sq-type asteroids. The most diagnostic feature is the $\approx$15$\pm$1$\%$ silicate absorption feature at $\approx$0.9-1.0~micron. Small S- and Sq-type asteroids have moderately high albedos ranging from 0.17-0.35. Using this albedo range for Dinkinesh in combination with measured absolute magnitude, it is possible to derive an effective diameter and surface brightness for this body. The albedo, size and surface brightness are important inputs required for planning a successful encounter by the Lucy spacecraft. △ Less

Submitted 7 April, 2023; v1 submitted 14 March, 2023; originally announced March 2023.

Comments: 7 pages, 1 figure. Accepted for publication in Icarus

Tuning the Legacy Survey of Space and Time (LSST) Observing Strategy for Solar System Science

Authors: Megan E. Schwamb, R. Lynne Jones, Peter Yoachim, Kathryn Volk, Rosemary C. Dorsey, Cyrielle Opitom, Sarah Greenstreet, Tim Lister, Colin Snodgrass, Bryce T. Bolin, Laura Inno, Michele T. Bannister, Siegfried Eggl, Michael Solontoi, Michael S. P. Kelley, Mario Jurić, Hsing Wen Lin, Darin Ragozzine, Pedro H. Bernardinelli, Steven R. Chesley, Tansu Daylan, Josef Ďurech, Wesley C. Fraser, Mikael Granvik, Matthew M. Knight , et al. (5 additional authors not shown)

Abstract: The Vera C. Rubin Observatory is expected to start the Legacy Survey of Space and Time (LSST) in early to mid-2025. This multi-band wide-field synoptic survey will transform our view of the solar system, with the discovery and monitoring of over 5 million small bodies.The final survey strategy chosen for LSST has direct implications on the discoverability and characterization of solar system minor… ▽ More The Vera C. Rubin Observatory is expected to start the Legacy Survey of Space and Time (LSST) in early to mid-2025. This multi-band wide-field synoptic survey will transform our view of the solar system, with the discovery and monitoring of over 5 million small bodies.The final survey strategy chosen for LSST has direct implications on the discoverability and characterization of solar system minor planets and passing interstellar objects. Creating an inventory of the solar system is one of the four main LSST science drivers. The LSST observing cadence is a complex optimization problem that must balance the priorities and needs of all the key LSST science areas. To design the best LSST survey strategy, a series of operation simulations using the Rubin Observatory scheduler have been generated to explore the various options for tuning observing parameters and prioritizations. We explore the impact of the various simulated LSST observing strategies on studying the solar system's small body reservoirs. We examine what are the best observing scenarios and review what are the important considerations for maximizing LSST solar system science. In general, most of the LSST cadence simulations produce +/-5% or less variations in our chosen key metrics, but a subset of the simulations significantly hinder science returns with much larger losses in the discovery and light curve metrics. △ Less

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

Comments: Accepted to ApJS, 103 pages (including references), 43 figures, 9 Tables. Videos will be available in the online journal formatted and published version of the paper [v2.0 submission corrects the author list metadata from the arxiv initial submission and updates the abstract]

Keck, Gemini, and Palomar 200-inch visible photometry of red and very-red Neptunian Trojans

Authors: B. T. Bolin, C. Fremling, A. Morbidelli, K. S. Noll, J. van Roestel, E. K. Deibert, M. Delbo, G. Gimeno, J. -E. Heo, C. M. Lisse, T. Seccull, H. Suh

Abstract: Neptunian Trojans (NTs), trans-Neptunian objects in 1:1 mean-motion resonance with Neptune, are generally thought to have been captured from the original trans-Neptunian protoplanetary disk into co-orbital resonance with the ice giant during its outward migration. It is possible, therefore, that the colour distribution of NTs is a constraint on the location of any colour transition zones that may… ▽ More Neptunian Trojans (NTs), trans-Neptunian objects in 1:1 mean-motion resonance with Neptune, are generally thought to have been captured from the original trans-Neptunian protoplanetary disk into co-orbital resonance with the ice giant during its outward migration. It is possible, therefore, that the colour distribution of NTs is a constraint on the location of any colour transition zones that may have been present in the disk. In support of this possible test, we obtained $g$, $r$, and $i$-band observations of 18 NTs, more than doubling the sample of NTs with known visible colours to 31 objects. Out of the combined sample, we found $\approx$4 objects with $g$-$i$ colours of $>$1.2 mags placing them in the very red (VR) category as typically defined. We find, without taking observational selection effects into account, that the NT $g$-$i$ colour distribution is statistically distinct from other trans-Neptunian dynamical classes. The optical colours of Jovian Trojans and NTs are shown to be less similar than previously claimed with additional VR NTs. The presence of VR objects among the NTs may suggest that the location of the red to VR colour transition zone in the protoplanetary disk was interior to 30-35 au. △ Less

Submitted 8 February, 2023; originally announced February 2023.

Comments: 8 pages, 3 figures, 3 tables, accepted for publication in MNRAS:Letters

The discovery and characterization of a kilometre sized asteroid inside the orbit of Venus

Authors: Bryce T. Bolin, T. Ahumada, P. van Dokkum, C. Fremling, M. Granvik, K. K. Hardegree-Ullman, Y. Harikane, J. N. Purdum, E. Serabyn, J. Southworth, C. Zhai

Abstract: Near-Earth asteroid population models predict the existence of bodies located inside the orbit of Venus. Despite searches up to the end of 2019, none had been found. We report discovery and follow-up observations of (594913) 'Ayló'chaxnim, an asteroid with an orbit entirely interior to Venus. (594913) 'Ayló'chaxnim has an aphelion distance of ~0.65 au, is ~2 km in diameter and is red in colour. Th… ▽ More Near-Earth asteroid population models predict the existence of bodies located inside the orbit of Venus. Despite searches up to the end of 2019, none had been found. We report discovery and follow-up observations of (594913) 'Ayló'chaxnim, an asteroid with an orbit entirely interior to Venus. (594913) 'Ayló'chaxnim has an aphelion distance of ~0.65 au, is ~2 km in diameter and is red in colour. The detection of such a large asteroid inside the orbit of Venus is surprising given their rarity according to near-Earth asteroid population models. As the first officially numbered and named asteroid located entirely within the orbit of Venus, we propose that the class of interior to Venus asteroids be referred to as 'Ayló'chaxnim asteroids. △ Less

Submitted 15 August, 2022; originally announced August 2022.

Comments: 7 pages, 8 figures, 1 table. Published in Monthly Notices of the Royal Astronomical Society: Letters (2022). arXiv admin note: text overlap with arXiv:2009.04125

Science Cases for the Keck Wide-Field Imager

Authors: J. Cooke, C. Angus, K. Auchettl, J. Bally, B. Bolin, S. Brough, J. N. Burchett, R. Foley, G. Foran, D. Forbes, J. Gannon, R. Hirai, G. G. Kacprzak, R. Margutti, C. Martinez-Lombilla, U. Mestric, A. Moller, A. Rest, J. Rhodes, R. M. Rich, F. Schussler, R. Wainscoat, J. Walawender, I. Wold, J. Zhang

Abstract: The Keck Wide-Field Imager (KWFI) is a proposed 1-degree diameter field of view UV-sensitive optical camera for Keck prime focus. KWFI will be the most powerful optical wide-field camera in the world and the only such 8m-class camera sensitive down to ~3000 A for the foreseeable future. Twenty science cases are described for KWFI compiled largely during 2019-2021, preceded by a brief discussion of… ▽ More The Keck Wide-Field Imager (KWFI) is a proposed 1-degree diameter field of view UV-sensitive optical camera for Keck prime focus. KWFI will be the most powerful optical wide-field camera in the world and the only such 8m-class camera sensitive down to ~3000 A for the foreseeable future. Twenty science cases are described for KWFI compiled largely during 2019-2021, preceded by a brief discussion of the instrument, components, and capabilities for context. △ Less

Submitted 24 July, 2022; originally announced July 2022.

Comments: 56 pages, 31 figures

Characterizing Sparse Asteroid Light Curves with Gaussian Processes

Authors: Christina Willecke Lindberg, Daniela Huppenkothen, R. Lynne Jones, Bryce T. Bolin, Mario Juric, V. Zach Golkhou, Eric C. Bellm, Andrew J. Drake, Matthew J. Graham, Russ R. Laher, Ashish A. Mahabal, Frank J. Masci, Reed Riddle, Kyung Min Shin

Abstract: In the era of wide-field surveys like the Zwicky Transient Facility and the Rubin Observatory's Legacy Survey of Space and Time, sparse photometric measurements constitute an increasing percentage of asteroid observations, particularly for asteroids newly discovered in these large surveys. Follow-up observations to supplement these sparse data may be prohibitively expensive in many cases, so to ov… ▽ More In the era of wide-field surveys like the Zwicky Transient Facility and the Rubin Observatory's Legacy Survey of Space and Time, sparse photometric measurements constitute an increasing percentage of asteroid observations, particularly for asteroids newly discovered in these large surveys. Follow-up observations to supplement these sparse data may be prohibitively expensive in many cases, so to overcome these sampling limitations, we introduce a flexible model based on Gaussian Processes to enable Bayesian parameter inference of asteroid time series data. This model is designed to be flexible and extensible, and can model multiple asteroid properties such as the rotation period, light curve amplitude, changing pulse profile, and magnitude changes due to the phase angle evolution at the same time. Here, we focus on the inference of rotation periods. Based on both simulated light curves and real observations from the Zwicky Transient Facility, we show that the new model reliably infers rotational periods from sparsely sampled light curves, and generally provides well-constrained posterior probability densities for the model parameters. We propose this framework as an intermediate method between fast, but very limited period detection algorithms and much more comprehensive, but computationally expensive shape modeling based on ray-tracing codes. △ Less

Submitted 24 November, 2021; originally announced November 2021.

Comments: 27 pages, 18 figures, accepted for publication in AJ, associated software available at https://github.com/dirac-institute/asterogap/tree/v0.1

Discovery and confirmation of the shortest gamma ray burst from a collapsar

Authors: Tomas Ahumada, Leo P. Singer, Shreya Anand, Michael W. Coughlin, Mansi M. Kasliwal, Geoffrey Ryan, Igor Andreoni, S. Bradley Cenko, Christoffer Fremling, Harsh Kumar, Peter T. H. Pang, Eric Burns, Virginia Cunningham, Simone Dichiara, Tim Dietrich, Dmitry S. Svinkin, Mouza Almualla, Alberto J. Castro-Tirado, Kishalay De, Rachel Dunwoody, Pradip Gatkine, Erica Hammerstein, Shabnam Iyyani, Joseph Mangan, Dan Perley , et al. (32 additional authors not shown)

Abstract: Gamma-ray bursts (GRBs) are among the brightest and most energetic events in the universe. The duration and hardness distribution of GRBs has two clusters, now understood to reflect (at least) two different progenitors. Short-hard GRBs (SGRBs; T90 <2 s) arise from compact binary mergers, while long-soft GRBs (LGRBs; T90 >2 s) have been attributed to the collapse of peculiar massive stars (collapsa… ▽ More Gamma-ray bursts (GRBs) are among the brightest and most energetic events in the universe. The duration and hardness distribution of GRBs has two clusters, now understood to reflect (at least) two different progenitors. Short-hard GRBs (SGRBs; T90 <2 s) arise from compact binary mergers, while long-soft GRBs (LGRBs; T90 >2 s) have been attributed to the collapse of peculiar massive stars (collapsars). The discovery of SN 1998bw/GRB 980425 marked the first association of a LGRB with a collapsar and AT 2017gfo/GRB 170817A/GW170817 marked the first association of a SGRB with a binary neutron star merger, producing also gravitational wave (GW). Here, we present the discovery of ZTF20abwysqy (AT2020scz), a fast-fading optical transient in the Fermi Satellite and the InterPlanetary Network (IPN) localization regions of GRB 200826A; X-ray and radio emission further confirm that this is the afterglow. Follow-up imaging (at rest-frame 16.5 days) reveals excess emission above the afterglow that cannot be explained as an underlying kilonova (KN), but is consistent with being the supernova (SN). Despite the GRB duration being short (rest-frame T90 of 0.65 s), our panchromatic follow-up data confirms a collapsar origin. GRB 200826A is the shortest LGRB found with an associated collapsar; it appears to sit on the brink between a successful and a failed collapsar. Our discovery is consistent with the hypothesis that most collapsars fail to produce ultra-relativistic jets. △ Less

Submitted 13 May, 2021; v1 submitted 11 May, 2021; originally announced May 2021.

Comments: submitted to Nature Astronomy

Tails: Chasing Comets with the Zwicky Transient Facility and Deep Learning

Authors: Dmitry A. Duev, Bryce T. Bolin, Matthew J. Graham, Michael S. P. Kelley, Ashish Mahabal, Eric C. Bellm, Michael W. Coughlin, Richard Dekany, George Helou, Shrinivas R. Kulkarni, Frank J. Masci, Thomas A. Prince, Reed Riddle, Maayane T. Soumagnac, Stéfan J. van der Walt

Abstract: We present Tails, an open-source deep-learning framework for the identification and localization of comets in the image data of the Zwicky Transient Facility (ZTF), a robotic optical time-domain survey currently in operation at the Palomar Observatory in California, USA. Tails employs a custom EfficientDet-based architecture and is capable of finding comets in single images in near real time, rath… ▽ More We present Tails, an open-source deep-learning framework for the identification and localization of comets in the image data of the Zwicky Transient Facility (ZTF), a robotic optical time-domain survey currently in operation at the Palomar Observatory in California, USA. Tails employs a custom EfficientDet-based architecture and is capable of finding comets in single images in near real time, rather than requiring multiple epochs as with traditional methods. The system achieves state-of-the-art performance with 99% recall, 0.01% false positive rate, and 1-2 pixel root mean square error in the predicted position. We report the initial results of the Tails efficiency evaluation in a production setting on the data of the ZTF Twilight survey, including the first AI-assisted discovery of a comet (C/2020 T2) and the recovery of a comet (P/2016 J3 = P/2021 A3). △ Less

Submitted 26 February, 2021; originally announced February 2021.

Time-series and Phasecurve Photometry of Episodically-Active Asteroid (6478) Gault in a Quiescent State Using APO, GROWTH, P200 and ZTF

Authors: Josiah N. Purdum, Zhong-Yi Lin, Bryce T. Bolin, Kritti Sharma, Philip I. Choi, Varun Bhalerao, Harsh Kumar, Robert Quimby, Joannes C. Van Roestel, Chengxing Zhai, Yanga R. Fernandez, Josef Hanuš, Carey M. Lisse, Dennis Bodewits, Christoffer Fremling, Nathan Ryan Golovich, Chen-Yen Hsu, Wing-Huen Ip, Chow-Choong Ngeow, Navtej S. Saini, Michael Shao, Yuhan Yao, Tomás Ahumada, Shreya Anand, Igor Andreoni , et al. (27 additional authors not shown)

Abstract: We observed Episodically Active Asteroid (6478) Gault in 2020 with multiple telescopes in Asia and North America and have found that it is no longer active after its recent outbursts at the end of 2018 and start of 2019. The inactivity during this apparation allowed us to measure the absolute magnitude of Gault of H_r = 14.63 +/- 0.02, G_r = 0.21 +/- 0.02 from our secular phasecurve observations.… ▽ More We observed Episodically Active Asteroid (6478) Gault in 2020 with multiple telescopes in Asia and North America and have found that it is no longer active after its recent outbursts at the end of 2018 and start of 2019. The inactivity during this apparation allowed us to measure the absolute magnitude of Gault of H_r = 14.63 +/- 0.02, G_r = 0.21 +/- 0.02 from our secular phasecurve observations. In addition, we were able to constrain Gault's rotation period using time-series photometric lightcurves taken over 17 hours on multiple days in 2020 August, September and October. The photometric lightcurves have a repeating $\lesssim$0.05 magnitude feature suggesting that (6478) Gault has a rotation period of ~2.5 hours and may have a semi-spherical or top-like shape, much like Near-Earth Asteroids Ryugu and Bennu. The rotation period of ~2.5 hours is near to the expected critical rotation period for an asteroid with the physical properties of (6478) Gault suggesting that its activity observed over multiple epochs is due to surface mass shedding from its fast rotation spun up by the Yarkovsky-O'Keefe-Radzievskii-Paddack effect. △ Less

Submitted 27 March, 2021; v1 submitted 25 February, 2021; originally announced February 2021.

Comments: 23 pages, 18 figures; Accepted by ApJ Letters

Initial Characterization of Active Transitioning Centaur, P/2019 LD2 (ATLAS), using Hubble, Spitzer, ZTF, Keck, APO and GROWTH Visible & Infrared Imaging and Spectroscopy

Authors: Bryce T. Bolin, Yanga R. Fernandez, Carey M. Lisse, Timothy R. Holt, Zhong-Yi Lin, Josiah N. Purdum, Kunal P. Deshmukh, James M. Bauer, Eric C. Bellm, Dennis Bodewits, Kevin B. Burdge, Sean J. Carey, Chris M. Copperwheat, George Helou, Anna Y. Q. Ho, Jonathan Horner, Jan van Roestel, Varun Bhalerao, Chan-Kao Chang, Christine Chen, Chen-Yen Hsu, Wing-Huen Ip, Mansi M. Kasliwal, Frank J. Masci, Chow-Choong Ngeow , et al. (21 additional authors not shown)

Abstract: We present visible and mid-infrared imagery and photometry of temporary Jovian co-orbital comet P/2019 LD$_2$ taken with HST/WFC3, Spitzer/IRAC, the GROWTH telescope network, visible spectroscopy from Keck/LRIS and archival ZTF observations taken between 2019 April and 2020 August. Our observations indicate that the nucleus of LD$_2$ has a radius between 0.2-1.8 km assuming a 0.08 albedo and a com… ▽ More We present visible and mid-infrared imagery and photometry of temporary Jovian co-orbital comet P/2019 LD$_2$ taken with HST/WFC3, Spitzer/IRAC, the GROWTH telescope network, visible spectroscopy from Keck/LRIS and archival ZTF observations taken between 2019 April and 2020 August. Our observations indicate that the nucleus of LD$_2$ has a radius between 0.2-1.8 km assuming a 0.08 albedo and a coma dominated by $\sim$100$μ$ m-scale dust ejected at $\sim$1 m/s speeds with a $\sim$1'' jet pointing in the SW direction. LD$_2$ experienced a total dust mass loss of $\sim$10$^8$ kg at a loss rate of $\sim$6 kg/s with Af$ρ$/cross-section varying between $\sim$85 cm/125 km$^2$ and $\sim$200 cm/310 km$^2$ from 2019 April 9 to 2019 Nov 8. If the increase in Af$ρ$/cross-section remained constant, it implies LD$_2$'s activity began $\sim$2018 November when within 4.8 au of the Sun, implying the onset of H$_2$O sublimation. We measure CO/CO$_2$ gas production of $\lesssim$10$^{27}$ mol/s /$\lesssim$10$^{26}$ mol/s from our 4.5 $μ$m Spitzer observations, $g$-$r$ = 0.59$\pm$0.03, $r$-$i$ = 0.18$\pm$0.05, $i$-$z$ = 0.01$\pm$0.07 from GROWTH observations, H$_2$O gas production of $\lesssim$80 kg/s scaling from our estimated $C_2$ production of $Q_{C_2}\lesssim$7.5$\times10^{24}$ mol/s from Keck/LRIS spectroscopy. We determine that the long-term orbit of LD$_2$ is similar to Jupiter family comets having close encounters with Jupiter within $\sim$0.5 Hill radius in the last $\sim$3 y, within 0.8 Hill radius in $\sim$9 y. Additionally, 78.8$\%$ of our orbital clones are ejected from the Solar System within $1 \times 10^{6}$ years having a dynamical half-life of 3.4 $\times 10^5$ years. △ Less

Submitted 5 January, 2021; v1 submitted 7 November, 2020; originally announced November 2020.

Comments: 28 pages, 8 figures, 2 Tables, accepted for publication in AJ

The Scientific Impact of the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) for Solar System Science

Authors: Vera C. Rubin Observatory LSST Solar System Science Collaboration, R. Lynne Jones, Michelle T. Bannister, Bryce T. Bolin, Colin Orion Chandler, Steven R. Chesley, Siegfried Eggl, Sarah Greenstreet, Timothy R. Holt, Henry H. Hsieh, Zeljko Ivezić, Mario Jurić, Michael S. P. Kelley, Matthew M. Knight, Renu Malhotra, William J. Oldroyd, Gal Sarid, Megan E. Schwamb, Colin Snodgrass, Michael Solontoi, David E. Trilling

Abstract: Vera C. Rubin Observatory will be a key facility for small body science in planetary astronomy over the next decade. It will carry out the Legacy Survey of Space and Time (LSST), observing the sky repeatedly in u, g, r, i, z, and y over the course of ten years using a 6.5 m effective diameter telescope with a 9.6 square degree field of view, reaching approximately r = 24.5 mag (5-σ depth) per visi… ▽ More Vera C. Rubin Observatory will be a key facility for small body science in planetary astronomy over the next decade. It will carry out the Legacy Survey of Space and Time (LSST), observing the sky repeatedly in u, g, r, i, z, and y over the course of ten years using a 6.5 m effective diameter telescope with a 9.6 square degree field of view, reaching approximately r = 24.5 mag (5-σ depth) per visit. The resulting dataset will provide extraordinary opportunities for both discovery and characterization of large numbers (10--100 times more than currently known) of small solar system bodies, furthering studies of planetary formation and evolution. This white paper summarizes some of the expected science from the ten years of LSST, and emphasizes that the planetary astronomy community should remain invested in the path of Rubin Observatory once the LSST is complete. △ Less

Submitted 14 September, 2020; originally announced September 2020.

Comments: White paper submitted to the 2020 Planetary Astronomy Decadal Survey (7 pages, 1 figure)

Optical follow-up of the neutron star-black hole mergers S200105ae and S200115j

Authors: Shreya Anand, Michael W. Coughlin, Mansi M. Kasliwal, Mattia Bulla, Tomás Ahumada, Ana Sagués Carracedo, Mouza Almualla, Igor Andreoni, Robert Stein, Francois Foucart, Leo P. Singer, Jesper Sollerman, Eric C. Bellm, Bryce Bolin, M. D. Caballero-García, Alberto J. Castro-Tirado, S. Bradley Cenko, Kishalay De, Richard G. Dekany, Dmitry A. Duev, Michael Feeney, Christoffer Fremling, Daniel A. Goldstein, V. Zach Golkhou, Matthew J. Graham , et al. (24 additional authors not shown)

Abstract: LIGO and Virgo's third observing run (O3) revealed the first neutron star-black hole (NSBH) merger candidates in gravitational waves. These events are predicted to synthesize r-process elements creating optical/near-IR "kilonova" (KN) emission. The joint gravitational-wave (GW) and electromagnetic detection of an NSBH merger could be used to constrain the equation of state of dense nuclear matter,… ▽ More LIGO and Virgo's third observing run (O3) revealed the first neutron star-black hole (NSBH) merger candidates in gravitational waves. These events are predicted to synthesize r-process elements creating optical/near-IR "kilonova" (KN) emission. The joint gravitational-wave (GW) and electromagnetic detection of an NSBH merger could be used to constrain the equation of state of dense nuclear matter, and independently measure the local expansion rate of the universe. Here, we present the optical follow-up and analysis of two of the only three high-significance NSBH merger candidates detected to date, S200105ae and S200115j, with the Zwicky Transient Facility (ZTF). ZTF observed $\sim$\,48\% of S200105ae and $\sim$\,22\% of S200115j's localization probabilities, with observations sensitive to KNe brighter than $-$17.5\,mag fading at 0.5\,mag/day in g- and r-bands; extensive searches and systematic follow-up of candidates did not yield a viable counterpart. We present state-of-the-art KN models tailored to NSBH systems that place constraints on the ejecta properties of these NSBH mergers. We show that with depths of $\rm m_{\rm AB}\approx 22$ mag, attainable in meter-class, wide field-of-view survey instruments, strong constraints on ejecta mass are possible, with the potential to rule out low mass ratios, high BH spins, and large neutron star radii. △ Less

Submitted 14 September, 2020; originally announced September 2020.

Preliminary estimates of the Zwicky Transient Facility 'Ayló'chaxnim asteroid population completeness

Authors: B. T. Bolin, T. Ahumada, P. van Dokkum, C. Fremling, K. K. Hardegree-Ullman, J. N. Purdum, E. Serabyn, J. Southworth

Abstract: Near-Earth asteroids (NEAs) are organized into five main classes: Amor, Apollo, Aten, Atira and 'Ayló'chaxnim. Asteroids belonging to the 'Ayló'chaxnim class are located entirely within the orbit of Venus making them difficult to detect by ground-based observatories. The first-known asteroid of this class, (594913) 'Ayló'chaxnim, was discovered by the Zwicky Transient Facility (ZTF) in 2020 Januar… ▽ More Near-Earth asteroids (NEAs) are organized into five main classes: Amor, Apollo, Aten, Atira and 'Ayló'chaxnim. Asteroids belonging to the 'Ayló'chaxnim class are located entirely within the orbit of Venus making them difficult to detect by ground-based observatories. The first-known asteroid of this class, (594913) 'Ayló'chaxnim, was discovered by the Zwicky Transient Facility (ZTF) in 2020 January during a twilight search for asteroids at small solar elongations that ran between September 2019 and January 2020. Due to its large diameter of $\sim$2 km, the discovery of (594913) 'Ayló'chaxnim is surprising because contemporary NEA population models predict a scarcity of asteroids of this size located inside the orbit of Venus. To compare the discovery of (594913) 'Ayló'chaxnim by ZTF with the predictions of NEA population models, we estimated the ZTF survey completeness at detecting 'Ayló'chaxnim asteroids and the number of 'Ayló'chaxnim asteroids expected to have been discovered by simulating observations of synthetic 'Ayló'chaxnim asteroids. We find that the 'Ayló'chaxnim population completeness of the survey is $\sim$18$\%$ and there is only a 5$\%$ probability that a single 'Ayló'chaxnim asteroid would have been discovered. Given the small chance for (594913) 'Ayló'chaxnim to have been discovered, its presence is either a statistical fluke or it implies that asteroid population models may need to be revised. △ Less

Submitted 19 January, 2023; v1 submitted 9 September, 2020; originally announced September 2020.

Comments: 9 pages, 2 figures. Accepted for publication in Icarus

Characterization of Temporarily-Captured Minimoon 2020 CD$_3$ by Keck Time-resolved Spectrophotometry

Authors: Bryce T. Bolin, Christoffer Fremling, Timothy R. Holt, Matthew J. Hankins, Tomás Ahumada, Shreya Anand, Varun Bhalerao, Kevin B. Burdge, Chris M. Copperwheat, Michael Coughlin, Kunal P. Deshmukh, Kishalay De, Mansi M. Kasliwal, Alessandro Morbidelli, Josiah N. Purdum, Robert Quimby, Dennis Bodewits, Chan-Kao Chang, Wing-Huen Ip, Chen-Yen Hsu, Russ R. Laher, Zhong-Yi Lin, Carey M. Lisse, Frank J. Masci, Chow-Choong Ngeow , et al. (20 additional authors not shown)

Abstract: We present time-resolved visible spectrophotometry of minimoon 2020 CD$_3$, the second asteroid known to become temporarily captured by the Earth-Moon system's gravitational field. The spectrophotometry was taken with Keck I/LRIS between wavelengths 434 nm and 912 nm in $B$, $g$, $V$, $R$, $I$ and RG850 filters as it was leaving the Earth-Moon system on 2020 March 23 UTC. The spectrophotometry of… ▽ More We present time-resolved visible spectrophotometry of minimoon 2020 CD$_3$, the second asteroid known to become temporarily captured by the Earth-Moon system's gravitational field. The spectrophotometry was taken with Keck I/LRIS between wavelengths 434 nm and 912 nm in $B$, $g$, $V$, $R$, $I$ and RG850 filters as it was leaving the Earth-Moon system on 2020 March 23 UTC. The spectrophotometry of 2020 CD$_3$ most closely resembles the spectra of V-type asteroids and some Lunar rock samples with a reddish slope of ~18$\%$/100 nm between 434 nm and 761 nm corresponding to colors of $g$-$r$ = 0.62$\pm$0.08, $r$-$i$ = 0.21 $\pm$ 0.06 and an absorption band at ~900 nm corresponding to $i$-$z$ = -0.54$\pm$0.10. Combining our measured 31.9$\pm$0.1 absolute magnitude with an albedo of 0.35 typical for V-type asteroids, we determine 2020 CD$_3$'s diameter to be ~0.9$\pm$0.1 m making it the first minimoon and one of the smallest asteroids to be spectrally studied. We use our time-series photometry to detect periodic lightcurve variations with a $<$10$^{-4}$ false alarm probability corresponding to a lightcurve period of ~573 s and a lightcurve amplitude of ~1 mag implying 2020 CD$_3$ possesses a $b/a$ axial ratio of ~2.5. In addition, we extend the observational arc of 2020 CD$_3$ to 37 days between 2020 February 15 UTC and 2020 March 23 UTC. From the improved orbital solution for 2020 CD$_3$, we estimate its likely duration of its capture to be ~2 y, and we measure the non-gravitation perturbation on its orbit due to radiation pressure with an area-to-mass ratio of 6.9$\pm$2.4$\times$10$^{-4}$ m$^2$/kg implying a density of 2.3$\pm$0.8 g/cm$^3$, broadly compatible with the densities of other meter-scale asteroids and Lunar rock. We searched for pre-discovery detections of 2020 CD$_3$ in the ZTF archive as far back as 2018 October, but were unable to locate any positive detections. △ Less

Submitted 12 August, 2020; originally announced August 2020.

Comments: 26 pages, 8 figures, 3 Tables, accepted for publication in ApJL

Kilonova Luminosity Function Constraints based on Zwicky Transient Facility Searches for 13 Neutron Star Mergers

Authors: Mansi M. Kasliwal, Shreya Anand, Tomas Ahumada, Robert Stein, Ana Sagues Carracedo, Igor Andreoni, Michael W. Coughlin, Leo P. Singer, Erik C. Kool, Kishalay De, Harsh Kumar, Mouza AlMualla, Yuhan Yao, Mattia Bulla, Dougal Dobie, Simeon Reusch, Daniel A. Perley, S. Bradley Cenko, Varun Bhalerao, David L. Kaplan, Jesper Sollerman, Ariel Goobar, Christopher M. Copperwheat, Eric C. Bellm, G. C. Anupama , et al. (78 additional authors not shown)

Abstract: We present a systematic search for optical counterparts to 13 gravitational wave (GW) triggers involving at least one neutron star during LIGO/Virgo's third observing run. We searched binary neutron star (BNS) and neutron star black hole (NSBH) merger localizations with the Zwicky Transient Facility (ZTF) and undertook follow-up with the Global Relay of Observatories Watching Transients Happen (GR… ▽ More We present a systematic search for optical counterparts to 13 gravitational wave (GW) triggers involving at least one neutron star during LIGO/Virgo's third observing run. We searched binary neutron star (BNS) and neutron star black hole (NSBH) merger localizations with the Zwicky Transient Facility (ZTF) and undertook follow-up with the Global Relay of Observatories Watching Transients Happen (GROWTH) collaboration. The GW triggers had a median localization of 4480 deg^2, median distance of 267 Mpc and false alarm rates ranging from 1.5 to 1e-25 per yr. The ZTF coverage had a median enclosed probability of 39%, median depth of 20.8mag, and median response time of 1.5 hr. The O3 follow-up by the GROWTH team comprised 340 UVOIR photometric points, 64 OIR spectra, and 3 radio. We find no promising kilonova (radioactivity-powered counterpart) and we convert the upper limits to constrain the underlying kilonova luminosity function. Assuming that all kilonovae are at least as luminous as GW170817 at discovery (-16.1mag), we calculate our joint probability of detecting zero kilonovae is only 4.2%. If we assume that all kilonovae are brighter than -16.6mag (extrapolated peak magnitude of GW170817) and fade at 1 mag/day (similar to GW170817), the joint probability of zero detections is 7%. If we separate the NSBH and BNS populations, the joint probability of zero detections, assuming all kilonovae are brighter than -16.6mag, is 9.7% for NSBH and 7.9% for BNS mergers. Moreover, <57% (<89%) of putative kilonovae could be brighter than -16.6mag assuming flat (fading) evolution, at 90% confidence. If we further account for the online terrestrial probability for each GW trigger, we find that <68% of putative kilonovae could be brighter than -16.6mag. Comparing to model grids, we find that some kilonovae must have Mej < 0.03 Msun or Xlan>1e-4 or phi>30deg to be consistent with our limits. (Abridged) △ Less

Submitted 19 June, 2020; originally announced June 2020.

Comments: Submitted to ApJ

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

Constraints on the Spin-Pole Orientation, Jet Morphology and Rotation of Interstellar Comet 2I/Borisov with Deep HST Imaging

Authors: Bryce T. Bolin, Carey M. Lisse

Abstract: We present high resolution, deep imaging of interstellar comet 2I/Borisov taken with the Hubble Space Telescope/Wide Field Camera 3 (HST/WFC3) on 2019 December 8 UTC and 2020 January 27 UTC (HST GO 16040, PI Bolin) before and after its perihelion passage in combination with HST/WFC3 images taken on 2019 October 12 UTC and 2019 November 16 UTC (HST GO/DD 16009, PI Jewitt) before its outburst and fr… ▽ More We present high resolution, deep imaging of interstellar comet 2I/Borisov taken with the Hubble Space Telescope/Wide Field Camera 3 (HST/WFC3) on 2019 December 8 UTC and 2020 January 27 UTC (HST GO 16040, PI Bolin) before and after its perihelion passage in combination with HST/WFC3 images taken on 2019 October 12 UTC and 2019 November 16 UTC (HST GO/DD 16009, PI Jewitt) before its outburst and fragmentation of March 2020, thus observing the comet in a relatively undisrupted state. We locate 1-2\arcsec~long (2,000 - 3,000 km projected length) jet-like structures near the optocenter of 2I that appear to change position angles from epoch to epoch. With the assumption that the jet is located near the rotational pole supported by its stationary appearance on $\sim$10-100 h time frames in HST images, we determine that 2I's pole points near $α$ = 322$\pm$10$^\circ$, $δ$ = 37$\pm$10$^\circ$ ($λ$ = 341$^\circ$, $β$ = 48$^\circ$) and may be in a simple rotation state. Additionally, we find evidence for possible periodicity in the HST time-series lightcurve on the time scale of $\sim$5.3 h with a small amplitude of $\sim$0.05 mag implying a lower limit on its $b/a$ ratio of $\sim$1.5 unlike the large $\sim$2 mag lightcurve observed for 1I/`Oumuamua. However, these small lightcurve variations may not be the result of the rotation of 2I's nucleus due to its dust-dominated light-scattering cross-section. Therefore, uniquely constraining the pre-Solar System encounter, pre-outburst rotation state of 2I may not be possible even with the resolution and sensitivity provided by HST observations. △ Less

Submitted 23 July, 2020; v1 submitted 13 December, 2019; originally announced December 2019.

Comments: 14 pages, 6 figures and 3 tables, accepted for publication in MNRAS on 23 July 2020

Pre-discovery Activity of New Interstellar Comet 2I/Borisov Beyond 5 AU

Authors: Quanzhi Ye, Michael S. P. Kelley, Bryce T. Bolin, Dennis Bodewits, Davide Farnocchia, Frank J. Masci, Karen J. Meech, Marco Micheli, Robert Weryk, Eric C. Bellm, Eric Christensen, Richard Dekany, Alexandre Delacroix, Matthew J. Graham, Shrinivas R. Kulkarni, Russ R. Laher, Ben Rusholme, Roger M. Smith

Abstract: Comet 2I/Borisov, the first unambiguous interstellar comet ever found, was discovered in August 2019 at $\sim3$ au from the Sun on its inbound leg. No pre-discovery detection beyond 3 au has yet been reported, mostly due to the comet's proximity to the Sun as seen from the Earth. Here we present a search for pre-discovery detections of comet Borisov using images taken by the Catalina Sky Survey (C… ▽ More Comet 2I/Borisov, the first unambiguous interstellar comet ever found, was discovered in August 2019 at $\sim3$ au from the Sun on its inbound leg. No pre-discovery detection beyond 3 au has yet been reported, mostly due to the comet's proximity to the Sun as seen from the Earth. Here we present a search for pre-discovery detections of comet Borisov using images taken by the Catalina Sky Survey (CSS), Pan-STARRS and Zwicky Transient Facility (ZTF), with a further comprehensive follow-up campaign being presented in \citet{Bolin2019}. We identified comet Borisov in ZTF images taken in May 2019 and use these data to update its orbit. This allowed us to identify the comet in images acquired as far back as December 2018, when it was 7.8 au from the Sun. The comet was not detected in November 2018 when it was 8.6 au from the Sun, possibly implying an onset of activity around this time. This suggests that the activity of the comet is either driven by a more volatile species other than H$_2$O, such as CO or CO$_2$, or by exothermic crystallization of amorphous ice. We derive the radius of the nucleus to be $<7$ km using the non-detection in November 2018, and estimate an area of $\sim0.5$---$10 \mathrm{km^2}$ has been active between December 2018 and September 2019, though this number is model-dependent and is highly uncertain. The behavior of comet Borisov during its inbound leg is observationally consistent with dynamically new comets observed in our solar system, suggesting some similarities between the two. △ Less

Submitted 23 December, 2019; v1 submitted 13 November, 2019; originally announced November 2019.

Comments: AJ in press

Characterization of the Nucleus, Morphology and Activity of Interstellar Comet 2I/Borisov by Optical and Near-Infrared GROWTH, Apache Point, IRTF, ZTF and Keck Observations

Authors: Bryce T. Bolin, Carey M. Lisse, Mansi M. Kasliwal, Robert Quimby, Hanjie Tan, Chris Copperwheat, Zhong-Yi Lin, Alessandro Morbidelli, Lyu Abe, Philippe Bendjoya, James Bauer, Kevin B. Burdge, Michael Coughlin, Christoffer Fremling, Ryosuke Itoh, Michael Koss, Frank J. Masci, Syota Maeno, Eric E. Mamajek, Federico Marocco, Katsuhiro Murata, Jean-Pierre Rivet, Michael L. Sitko, Daniel Stern, David Vernet , et al. (30 additional authors not shown)

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

Submitted 12 May, 2020; v1 submitted 30 October, 2019; originally announced October 2019.

Comments: 22 pages, 7 figures, 1 table. Accepted for publication in AJ on 12 May 2020

Maximizing LSST Solar System Science: Approaches, Software Tools, and Infrastructure Needs

Authors: Henry H. Hsieh, Michele T. Bannister, Bryce T. Bolin, Josef Durech, Siegfried Eggl, Wesley C. Fraser, Mikael Granvik, Michael S. P. Kelley, Matthew M. Knight, Rodrigo Leiva, Marco Micheli, Joachim Moeyens, Michael Mommert, Darin Ragozzine, Cristina A. Thomas

Abstract: The Large Synoptic Survey Telescope (LSST) is expected to increase known small solar system object populations by an order of magnitude or more over the next decade, enabling a broad array of transformative solar system science investigations to be performed. In this white paper, we discuss software tools and infrastructure that we anticipate will be needed to conduct these investigations and outl… ▽ More The Large Synoptic Survey Telescope (LSST) is expected to increase known small solar system object populations by an order of magnitude or more over the next decade, enabling a broad array of transformative solar system science investigations to be performed. In this white paper, we discuss software tools and infrastructure that we anticipate will be needed to conduct these investigations and outline possible approaches for implementing them. Feedback from the community or contributions to future updates of this work are welcome. Our aim is for this white paper to encourage further consideration of the software development needs of the LSST solar system science community, and also to be a call to action for working to meet those needs in advance of the expected start of the survey in late 2022. △ Less

Submitted 26 June, 2019; originally announced June 2019.

Comments: 82 pages

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

Towards Efficient Detection of Small Near-Earth Asteroids Using the Zwicky Transient Facility (ZTF)

Authors: Quanzhi Ye, Frank J. Masci, Hsing Wen Lin, Bryce Bolin, Chan-Kao Chang, Dmitry A. Duev, George Helou, Wing-Huen Ip, David L. Kaplan, Emily Kramer, Ashish Mahabal, Chow-Choong Ngeow, Avery J. Nielsen, Thomas A. Prince, Hanjie Tan, Ting-Shuo Yeh, Eric C. Bellm, Richard Dekany, Matteo Giomi, Matthew J. Graham, Shrinivas R. Kulkarni, Thomas Kupfer, Russ R. Laher, Ben Rusholme, David L. Shupe , et al. (1 additional authors not shown)

Abstract: We describe ZStreak, a semi-real-time pipeline specialized in detecting small, fast-moving near-Earth asteroids (NEAs) that is currently operating on the data from the newly-commissioned Zwicky Transient Facility (ZTF) survey. Based on a prototype originally developed by Waszczak et al. (2017) for the Palomar Transient Factory (PTF), the predecessor of ZTF, ZStreak features an improved machine-lea… ▽ More We describe ZStreak, a semi-real-time pipeline specialized in detecting small, fast-moving near-Earth asteroids (NEAs) that is currently operating on the data from the newly-commissioned Zwicky Transient Facility (ZTF) survey. Based on a prototype originally developed by Waszczak et al. (2017) for the Palomar Transient Factory (PTF), the predecessor of ZTF, ZStreak features an improved machine-learning model that can cope with the $10\times$ data rate increment between PTF and ZTF. Since its first discovery on 2018 February 5 (2018 CL), ZTF/ZStreak has discovered $45$ confirmed new NEAs over a total of 232 observable nights until 2018 December 31. Most of the discoveries are small NEAs, with diameters less than $\sim100$ m. By analyzing the discovery circumstances, we find that objects having the first to last detection time interval under 2 hr are at risk of being lost. We will further improve real-time follow-up capabilities, and work on suppressing false positives using deep learning. △ Less

Submitted 21 April, 2019; originally announced April 2019.

Comments: PASP in press

Multiple Outbursts of Asteroid (6478) Gault

Authors: Quanzhi Ye, Michael S. P. Kelley, Dennis Bodewits, Bryce Bolin, Zhong-Yi Lin, Eric C. Bellm, Richard Dekany, Dmitry A. Duev, Steven Groom, George Helou, Shrinivas R. Kulkarni, Thomas Kupfer, Frank J. Masci, Thomas A. Prince, Maayane T. Soumagnac

Abstract: Main-belt asteroid (6478) Gault unexpectedly sprouted two tails in late 2018 and early 2019, identifying it as a new active asteroid. Here we present observations obtained by the 1.2-m Zwicky Transient Facility survey telescope that provide detailed time-series coverage of the onset and evolution of Gault's activity. Gault exhibited two brightening events, with the first one starting on 2018 Oct.… ▽ More Main-belt asteroid (6478) Gault unexpectedly sprouted two tails in late 2018 and early 2019, identifying it as a new active asteroid. Here we present observations obtained by the 1.2-m Zwicky Transient Facility survey telescope that provide detailed time-series coverage of the onset and evolution of Gault's activity. Gault exhibited two brightening events, with the first one starting on 2018 Oct. 18$\pm5$ days and a second one starting on 2018 Dec. 24$\pm1$ days. The amounts of mass released are $2\times10^7$ kg and $1\times10^6$ kg, respectively. Based on photometric measurements, each event persisted for about a month. Gault's color has not changed appreciably over time, with a pre-outburst color of $g_\mathrm{PS1}-r_\mathrm{PS1}=0.50\pm0.04$ and $g_\mathrm{PS1}-r_\mathrm{PS1}=0.46\pm0.04$ during the two outbursts. Simulations of dust dynamics shows that the ejecta consists of dust grains of up to 10 $μ$m in size that are ejected at low velocities below $1~\mathrm{m~s^{-1}}$ regardless of particle sizes. This is consistent with non-sublimation-driven ejection events. The size distribution of the dust exhibits a broken power-law, with particles at 10--20 $μ$m following a power-law of $-2.5$ to $-3.0$, while larger particles follow a steeper slope of $-4.0$. The derived properties can be explained by either rotational excitation of the nucleus or a merger of a near-contact binary, with the latter scenario to be statistically more likely. △ Less

Submitted 13 March, 2019; originally announced March 2019.

Comments: ApJL in press. Animation of Figure 1 at https://youtu.be/_wuuabhe5B0. Data and codes that generate the figures of this paper are available at https://github.com/Yeqzids/activation_of_6478_gault

The Zwicky Transient Facility: Science Objectives

Authors: Matthew J. Graham, S. R. Kulkarni, Eric C. Bellm, Scott M. Adams, Cristina Barbarino, Nadejda Blagorodnova, Dennis Bodewits, Bryce Bolin, Patrick R. Brady, S. Bradley Cenko, Chan-Kao Chang, Michael W. Coughlin, Kishalay De, Gwendolyn Eadie, Tony L. Farnham, Ulrich Feindt, Anna Franckowiak, Christoffer Fremling, Avishay Gal-yam, Suvi Gezari, Shaon Ghosh, Daniel A. Goldstein, V. Zach Golkhou, Ariel Goobar, Anna Y. Q. Ho , et al. (92 additional authors not shown)

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

Submitted 5 February, 2019; originally announced February 2019.

Comments: 26 pages, 7 figures, Published in PASP Focus Issue on the Zwicky Transient Facility

The Zwicky Transient Facility: System Overview, Performance, and First Results

Authors: Eric C. Bellm, Shrinivas R. Kulkarni, Matthew J. Graham, Richard Dekany, Roger M. Smith, Reed Riddle, Frank J. Masci, George Helou, Thomas A. Prince, Scott M. Adams, C. Barbarino, Tom Barlow, James Bauer, Ron Beck, Justin Belicki, Rahul Biswas, Nadejda Blagorodnova, Dennis Bodewits, Bryce Bolin, Valery Brinnel, Tim Brooke, Brian Bue, Mattia Bulla, Rick Burruss, S. Bradley Cenko , et al. (91 additional authors not shown)

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

Submitted 5 February, 2019; originally announced February 2019.

Comments: Published in PASP Focus Issue on the Zwicky Transient Facility (https://dx.doi.org/10.1088/1538-3873/aaecbe). 21 Pages, 12 Figures

Journal ref: Publications of the Astronomical Society of the Pacific, Volume 131, Issue 995, pp. 018002 (2019)

Fast algorithms for slow moving asteroids: constraints on the distribution of Kuiper Belt Objects

Authors: Peter J. Whidden, J. Bryce Kalmbach, Andrew J. Connolly, R. Lynne Jones, Hayden Smotherman, Dino Bektesevic, Colin Slater, Andrew C. Becker, Željko Ivezić, Mario Jurić, Bryce Bolin, Joachim Moeyens, Francisco Förster, V. Zach Golkhou

Abstract: We introduce a new computational technique for searching for faint moving sources in astronomical images. Starting from a maximum likelihood estimate for the probability of the detection of a source within a series of images, we develop a massively parallel algorithm for searching through candidate asteroid trajectories that utilizes Graphics Processing Units (GPU). This technique can search over… ▽ More We introduce a new computational technique for searching for faint moving sources in astronomical images. Starting from a maximum likelihood estimate for the probability of the detection of a source within a series of images, we develop a massively parallel algorithm for searching through candidate asteroid trajectories that utilizes Graphics Processing Units (GPU). This technique can search over 10^10 possible asteroid trajectories in stacks of the order 10-15 4K x 4K images in under a minute using a single consumer grade GPU. We apply this algorithm to data from the 2015 campaign of the High Cadence Transient Survey (HiTS) obtained with the Dark Energy Camera (DECam). We find 39 previously unknown Kuiper Belt Objects in the 150 square degrees of the survey. Comparing these asteroids to an existing model for the inclination distribution of the Kuiper Belt we demonstrate that we recover a KBO population above our detection limit consistent with previous studies. Software used in this analysis is made available as an open source package. △ Less

Submitted 8 January, 2019; originally announced January 2019.

Comments: Accepted to 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

APO Time Resolved Color Photometry of Highly-Elongated Interstellar Object 1I/'Oumuamua

Authors: Bryce T. Bolin, Harold A. Weaver, Yanga R. Fernandez, Carey M. Lisse, Daniela Huppenkothen, R. Lynne Jones, Mario Juric, Joachim Moeyens, Charles A. Schambeau, Colin T. Slater, Zeljko Ivezic, Andrew J. Connolly

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

Submitted 29 January, 2018; v1 submitted 13 November, 2017; originally announced November 2017.

Comments: 13 pages, 6 figures, 1 table, accepted for publication in the Astrophysical Journal Letters, data are available at: https://github.com/dhuppenkothen/CometPeriodSearch and https://zenodo.org/record/1068392

Size-dependent modification of asteroid family Yarkovsky V-shapes

Authors: Bryce T. Bolin, Alessandro Morbidelli, Kevin J. Walsh

Abstract: The thermal properties of the surfaces of asteroids determine the magnitude of the drift rate cause by the Yarkovsky force. In the general case of Main Belt asteroids, the Yarkovsky force is indirectly proportional to the thermal inertia, $Γ$. Following the proposed relationship between $Γ$ and asteroid diameter $D$, we find that asteroids' Yarkovsky drift rates might have a more complex size depe… ▽ More The thermal properties of the surfaces of asteroids determine the magnitude of the drift rate cause by the Yarkovsky force. In the general case of Main Belt asteroids, the Yarkovsky force is indirectly proportional to the thermal inertia, $Γ$. Following the proposed relationship between $Γ$ and asteroid diameter $D$, we find that asteroids' Yarkovsky drift rates might have a more complex size dependence than previous thought, leading to a curved family V-shape boundary in semi-major axis, a, vs. 1/$D$ space. This implies that asteroids are drifting faster at larger sizes than previously considered decreasing on average the known ages of asteroid families. The V-Shape curvature is determined for $>$25 families located throughout the Main Belt to quantify the Yarkovsky size-dependent drift rate. We find that there is no correlation between family age and V-shape curvature. In addition, the V-shape curvature decreases for asteroid families with larger heliocentric distances suggesting that the relationship between $Γ$ and $D$ is weaker in the outer MB possibly due to homogenous surface roughness among family members. △ Less

Submitted 4 March, 2018; v1 submitted 11 October, 2017; originally announced October 2017.

Comments: 68 Pages, 69 Figures, 4 Tables, accepted for publication in Astronomy & Astrophysics. arXiv admin note: text overlap with arXiv:1708.03662

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

Initial velocity V-shapes of young asteroid families

Authors: Bryce T. Bolin, Kevin J. Walsh, Alessandro Morbidelli, Marco Delbo

Abstract: Ejection velocity fields of asteroid families are largely unconstrained due to the fact that members disperse relatively quickly on Myr time-scales by secular resonances and the Yarkovsky effect. The spreading of fragments in $a$ by the Yarkovsky effect is indistinguishable from the spreading caused by the initial ejection of fragments. By examining families $<$20 Myrs-old, we can use the V-shape… ▽ More Ejection velocity fields of asteroid families are largely unconstrained due to the fact that members disperse relatively quickly on Myr time-scales by secular resonances and the Yarkovsky effect. The spreading of fragments in $a$ by the Yarkovsky effect is indistinguishable from the spreading caused by the initial ejection of fragments. By examining families $<$20 Myrs-old, we can use the V-shape identification technique to separate family shapes that are due to the initial ejection velocity field and those that are due to the Yarkovsky effect. $<$20 Myr-old asteroid families provide an opportunity to study the velocity field of family fragments before they become too dispersed. Only the Karin family's initial velocity field has been determined and scales inversely with diameter, $D^{-1}$. We have applied the V-shape identification technique to constrain young families' initial ejection velocity fields by measuring the curvature of their fragments' V-shape correlation in semi-major axis, $a$, vs. $D^{-1}$ space. Curvature from a straight line implies a deviation from a scaling of $D^{-1}$. We measure the V-shape curvature of 11 young asteroid families including the \FYnospace, Aeolia, Brangane, Brasilia, Clarissa, Iannini, Karin, Konig, Koronis(2), Theobalda and Veritas asteroid families. We find that the majority of asteroid families have initial ejection velocity fields consistent with $\sim D^{-1}$ supporting laboratory impact experiments and computer simulations of disrupting asteroid parent bodies. △ Less

Submitted 26 September, 2017; v1 submitted 11 August, 2017; originally announced August 2017.

Comments: 68 Pages, 25 Figures, 2 Tables, accepted for publication in MNRAS

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

Yarkovsky V-shape identification of asteroid families

Authors: Bryce T. Bolin, Marco Delbo, Alessandro Morbidelli, Kevin J. Walsh

Abstract: There are only a few known main belt (MB) asteroid families with ages greater than 2 Gyr (Brož et al. 2013, Spoto et al. 2015). Estimates based on the family producing collision rate suggest that the lack of >2 Gyr-old families may be due to a selection bias in current techniques used to identify families. Family fragments disperse in their orbital elements, semi-major axis, $a$, eccentricity,… ▽ More There are only a few known main belt (MB) asteroid families with ages greater than 2 Gyr (Brož et al. 2013, Spoto et al. 2015). Estimates based on the family producing collision rate suggest that the lack of >2 Gyr-old families may be due to a selection bias in current techniques used to identify families. Family fragments disperse in their orbital elements, semi-major axis, $a$, eccentricity, $e$, and inclination, $i$, due to secular resonances, close encounters with massive asteroids and the non-gravitational Yarkovsky force. This causes the family fragments to be indistinguishable from the background of the main belt making them more difficult to identify with the hierarchical clustering method (HCM) with increasing family age. The discovery of the Eulalia and new Polana families in the inner belt relied on new techniques because Yarkovsky spreading made them too disperse to be identified using the classical HCM. The technique used to discover the new Polana and Eulalia families is modified to identify asteroid families by searching for correlations between $a$ and asteroid diameter, $D$, or absolute magnitude, $H$. A group of asteroids is identified as a collisional family if its boundary in the $a$ vs. 1/$D$ or $a$ vs. $H$ planes has a characteristic V-shape which is due to the size dependent Yarkovsky spreading. The V-shape boundary is identified with two separate techniques. The first identifies a border by measuring the drop between the number of objects inside and outside of the border. The second identifies the V-shape border by measuring the number density of objects in $a$ vs. 1/$D$,$H$ space. The V-shape techniques are demonstrated on known families Erigone, Vesta, Koronis, and families difficult to identify by HCM such as Flora, Baptistina, new Polana, Eulalia and Karin. Application of the technique in a search for >2 Gyr-old families throughout the MB is discussed. △ Less

Submitted 20 September, 2016; originally announced September 2016.

Comments: 74 pages, 17 figures, accepted for publication in Icarus

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)

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

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

Authors: Henry H. Hsieh, Olivier Hainaut, Bojan Novakovic, Bryce Bolin, Larry Denneau, Alan Fitzsimmons, Nader Haghighipour, Jan Kleyna, Rosita Kokotanekova, Pedro Lacerda, Karen J. Meech, Marco Micheli, Nick Moskovitz, Eva Schunova, Colin Snodgrass, Richard J. Wainscoat, Lawrence Wasserman, Adam Waszczak

Abstract: 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 t… ▽ More 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. △ Less

Submitted 15 January, 2015; originally announced January 2015.

Comments: 8 pages, 4 figures; accepted for publication in ApJ Letters

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

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

The Pan-STARRS Moving Object Processing System

Authors: Larry Denneau, Robert Jedicke, Tommy Grav, Mikael Granvik, Jeremy Kubica, Andrea Milani, Peter Veres, Richard Wainscoat, Daniel Chang, Francesco Pierfederici, N. Kaiser, K. C. Chambers, J. N. Heasley, Eugene. A. Magnier, P. A. Price, Jonathan Myers, Jan Kleyna, Henry Hsieh, Davide Farnocchia, Chris Waters, W. H. Sweeney, Denver Green, Bryce Bolin, W. S. Burgett, J. S. Morgan , et al. (19 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 28 February, 2013; originally announced February 2013.

Comments: 57 Pages, 26 Figures, 13 Tables