Search

Small bodies are capable of delivering essential prerequisites for the development of life, such as volatiles and organics, to the terrestrial planets. For example, empirical evidence suggests that water was delivered to the Earth by hydrated planetesimals from distant regions of the Solar System. Recently, several morphologically inactive near-Earth objects (NEOs) were reported to experience significant nongravitational accelerations inconsistent with radiation-based effects, and possibly explained by volatile-driven outgassing. However, these "dark comets" display no evidence of comae in archival images, which are the defining feature of cometary activity. Here we report detections of nongravitational accelerations on seven additional objects previously classified as inactive (doubling the population) that could also be explainable by asymmetric mass loss. A detailed search of archival survey and targeted data rendered no detection of dust activity in any of these objects in individual or stacked images. We calculate dust production limits of $\sim10$, $0.1$, and $0.1$ kg s$^{-1}$ for 1998 FR$_{11}$, 2001 ME$_{1}$, and 2003 RM with these data, indicating little or no dust surrounding the objects during the observations. This set of dark comets reveals the delineation between two distinct populations: larger, "outer" dark comets on eccentric orbits that are end members of a continuum in activity level of comets, and smaller, "inner" dark comets on near-circular orbits that could signify a new population. These objects may trace various stages in the life cycle of a previously undetected, but potentially numerous, volatile-rich population that may have provided essential material to the Earth., Comment: 12 pages, 11 figures, Author's version of Proc. Natl. Acad. Sci. Vol. 121 No. 51 published at https://www.pnas.org/doi/10.1073/pnas.2406424121

In this paper, we investigate the nonprincipal axis (NPA) rotational state of 1I/`Oumuamua -- the first interstellar object discovered traversing the inner Solar System -- from its photometric light curve. Building upon Mashchenko (2019), we develop a model which incorporates NPA rotation and {Sun-induced, time-varying} outgassing torques to generate synthetic light curves of the object. The model neglects tidal forces, which are negligible compared to outgassing torques over the distances that `Oumuamua was observed. We implement an optimization scheme that incorporates the NPA rotation model to calculate the initial rotation state of the object. We find that an NPA rotation state with an average period of $\langle P \rangle\simeq7.34$ hr best reproduces the photometric data. The discrepancy between this period and previous estimates is due to continuous period modulation induced by outgassing torques in the rotational model, {as well as different periods being used}. The best fit to the October 2017 data does not reproduce the November 2017 data (although the later measurements are too sparse to fit). The light curve is consistent with no secular evolution of the angular momentum, somewhat in tension with the empirical correlations between nuclear spin-up and cometary outgassing. The complex rotation of `Oumuamua may be {the result of primordial rotation about the smallest principal axis} if (i) the object experienced hypervolatile outgassing and (ii) our idealized outgassing model is accurate., Comment: 22 pages, 8 figures, 1 animation. Accepted to the Planetary Science Journal. The animation can be found on YouTube (https://youtu.be/f5YEAMTvIeo) and in the online publication by PSJ (when available)

We present the detection of false positive double-lined spectroscopic binaries candidates (SB2) using medium-resolution survey (MRS) spectra from the one time-domain field of LAMOST data release 10 (DR10). The secondary component in all these binaries has near zero radial velocity and solar-like spectral lines. Highly likely this is light from the semi-transparent clouds illuminated by the full Moon. However we also suspect that partially this contamination can be caused by a solar light reflected from the surface of low-orbital artificial satellites launched in the beginning of 2022. We found several possible contaminant candidates using archival orbital data. We propose measures to reduce risk of such contamination for the future observations and methods to find it in archived ones., Comment: accepted in MNRAS letters 2023 July 07

We report statistically significant detections of non-radial nongravitational accelerations based on astrometric data in the photometrically inactive objects 1998 KY$_{26}$, 2005 VL$_1$, 2016 NJ$_{33}$, 2010 VL$_{65}$, 2016 RH$_{120}$, and 2010 RF$_{12}$. The magnitudes of the nongravitational accelerations are greater than those typically induced by the Yarkovsky effect and there is no radiation-based, non-radial effect that can be so large. Therefore, we hypothesize that the accelerations are driven by outgassing, and calculate implied H$_2$O production rates for each object. We attempt to reconcile outgassing induced acceleration with the lack of visible comae or photometric activity via the absence of surface dust and low levels of gas production. Although these objects are small and some are rapidly rotating, surface cohesive forces are stronger than the rotational forces and rapid rotation alone cannot explain the lack of surface debris. It is possible that surface dust was removed previously, perhaps via outgassing activity that increased the rotation rates to their present day value. We calculate dust production rates of order $\sim10^{-4}$ g s$^{-1}$ in each object assuming that the nuclei are bare, within the upper limits of dust production from a sample stacked image of 1998 KY$_{26}$ of $\dot{M}_{\rm Dust}<0.2$ g s$^{-1}$. This production corresponds to brightness variations of order $\sim0.0025\%$, which are undetectable in extant photometric data. We assess the future observability of each of these targets, and find that the orbit of 1998 KY$_{26}$ -- which is also the target for the extended Hayabusa2 mission -- exhibits favorable viewing geometry before 2025., Comment: 21 pages, 9 figures, 3 tables, accepted for publication in the Planetary Science Journal

2009 MS9 is a trans-Neptunian object (TNO) whose perihelion brings it close to the distance where some long period comets are seen to become active. Knowing this, and the fact that this object appears to brighten in excess of it's predicted nucleus brightness suggests that 2009 MS9 has a delayed onset of activity brought on by the sublimation of a species more volatile than water. In this paper we characterize 2009 MS9's physical properties and investigate potential outgassing through composite images, sublimation models, and measurements of spectral reflectivity. We find that deep composite images of the object at various epochs along its orbit show no evidence of dust yet place sensitive limits to the dust production. We estimate the nucleus radius to be 11.5 km $\pm 3.5$ km using thermal IR modeling from NEOWISE data and use this and data pre-perihelion to estimate a geometric albedo of 0.25. We compare a CO-sublimation activity model to its post perihelion heliocentric light curve and find this data supports an active fractional area of $5 \times 10^{-6}$ assuming 2 $\mu$m sized grains and other typical comet parameters. The spectral reflectivity of the surface materials obtained with the Gemini Observatory and CFHT at different epochs shows a reddening spectral slope. We compare the physical properties of 2009 MS9 to both TNO and comet populations, and speculate that 2009 MS9's reddening may be due to the buildup of a dust mantle on the surface and could be an explanation of why TNOs exhibit a color bimodality., Comment: 25 pages, 8 figures, Accepted to PSJ

In this paper, we present an open source software ($\texttt{SAMUS}$) which simulates constant-density, constant-viscosity liquid bodies subject to tidal forces for a range of assumed viscosites and sizes. This software solves the Navier-Stokes equations on a finite-element mesh, incorporating the centrifugal, Coriolis, self-gravitational, and tidal forces. The primary functionality is to simulate the deformation of minor bodies under the influence of tidal forces. It may therefore be used to constrain the composition and physical structure of bodies experiencing significant tidal forces, such as 99942 Apophis and 1I/`Oumuamua. We demonstrate that $\texttt{SAMUS}$ will be useful to constrain the material properties of Apophis during its near-Earth flyby in 2029. Depending on the material properties, Apophis may experience an area change of up to 0.5\%, with similar effects on the photometric brightness. We also apply $\texttt{SAMUS}$ to constrain the material dynamic viscosity of 1I/`Oumuamua, the first interstellar object discovered traversing the inner Solar System. `Oumuamua experienced a close approach to the Sun at perihelion ($q\simeq 0.25$ au) during which there were significant tidal forces that may have caused deformation of the body. This deformation could have lead to observable changes in the photometric light curve based on the material properties. The application of $\texttt{SAMUS}$ to produce synthetic observations which incorporate tidal deformation effects demonstrate that no deformation -- an infinite dynamic viscosity -- best reproduces the photometric data. While these results indicate that `Oumuamua did not experience significant tidal deformation, a sophisticated model incorporating non-principal axis rotation is necessary to conclusively analyze both `Oumuamua and Apophis., Comment: 26 pages, 15 figures, 5 tables. Accepted by AAS Planetary Science Journal. Split into two papers based on reviewer and editor suggestions. Comments welcome. Publicly available software at https://github.com/astertaylor/Oumuamua

Observatory end-to-end science operations is the overall process starting with a scientific question, represented by a proposal requesting observing time, and ending with the analysis of observation data addressing that question, and including all the intermediate steps needed to plan, schedule, obtain, and process these observations. Increasingly complex observing facilities demand a highly efficient science operations approach and at the same time be user friendly to the astronomical user community and enable the highest possible scientific return. Therefore, this process is supported by a collection of tools. In this paper, we describe the overall end-to-end process and its implementation for the three upcoming extremely large telescopes (ELTs), ESO's ELT, the Thirty Meter Telescope (TMT), and the Giant Magellan Telescope (GMT)., Comment: 38 pages, 20 figures. Revised submission to the SPIE Journal of Astronomical Telescopes, Instruments, and Systems

As a consequence of the large (and growing) number of near-Earth objects discovered, some of them are lost before their orbit can be firmly established to ensure long-term recovery. A fraction of these present non-negligible chances of impact with the Earth. We present a method of targeted observations that allowed us to eliminate that risk by obtaining deep images of the area where the object would be, should it be on a collision orbit. 2006 QV89 was one of these objects, with a chance of impact with the Earth on 2019 September 9. Its position uncertainty (of the order of 1 degree) and faintness (below V$\sim$24) made it a difficult candidate for a traditional direct recovery. However, the position of the virtual impactors could be determined with excellent accuracy. In July 2019, the virtual impactors of 2006 QV89 were particularly well placed, with a very small uncertainty region, and an expected magnitude of V$<$26. The area was imaged using ESO's Very Large Telescope, in the context of the ESA/ESO collaboration on Near-Earth Objects, resulting in very constraining a non-detection. This resulted in the elimination of the virtual impactor, even without effectively recovering 2006 QV89, indicating that it did not represent a threat. This method of deep non-detection of virtual impactors demonstrated a large potential to eliminate the threat of other-wise difficult to recover near-Earth objects, Comment: Accepted for publication in Astronomy & Astrophysics

The interstellar traveler, 2I/Borisov, is the first clearly active extrasolar comet, ever detected in our Solar system. We obtained high-resolution interferometric observations of 2I/Borisov with the Atacama Large Millimeter/submillimeter Array (ALMA), and multi-color optical observations with the Very Large Telescope (VLT) to gain a comprehensive understanding of the dust properties of this comet. We found that the dust coma of 2I/Borisov consists of compact "pebbles" of radii exceeding ~1 mm, suggesting that the dust particles have experienced compaction through mutual impacts during the bouncing collision phase in the protoplanetary disk. We derived a dust mass loss rate of >= 200 kg/s and a dust-to-gas ratio >=3. Our long term monitoring of 2I/Borisov with VLT indicates a steady dust mass loss with no significant dust fragmentation and/or sublimation occurring in the coma. We also detected emissions from carbon monoxide gas (CO) with ALMA and derived the gas production rate of Q(CO) (3.3+/-0.8)x10^{26} mole/s. We found that the CO/H$_2$O mixing ratio of 2I/Borisov changed drastically before and after perihelion, indicating the heterogeneity of the cometary nucleus, with components formed at different locations beyond the volatile snow-line with different chemical abundances. Our observations suggest that 2I/Borisov's home system, much like our own system, experienced efficient radial mixing from the innermost parts of its protoplanetary disk to beyond the frost line of CO., Comment: 12 pages, 4 figures

The effect of satellite constellations on observations in the visible and IR domains is estimated, considering 18 constellations in development by SpaceX, Amazon, OneWeb, and others, with over 26,000 satellites, constituting a representative distribution. This study uses a series of simplifications and assumptions to obtain conservative, order-of-magnitude estimates of the effects. The number of illuminated satellites from the constellations above the horizon ranges from ~1600 right after sunset, decreasing to 1100 at the end of astronomical twilight, most of them (~85%) close to the horizon (< 30deg). The large majority of these satellites will be too faint to be seen with the naked eye: at astronomical twilight, 110 brighter than mag 5. Most of them (~95%) will be close to the horizon. The number of naked-eye satellites plummets as the Sun reaches 30-40 deg below the horizon, depending on the latitude and season. The light trails caused by satellites would ruin a small fraction (below the 1% level) of exposures using narrow to normal field imaging or spectroscopic techniques in the visible and near IR during the first and last hours of the night. Similarly, the thermal emission of the satellite would affect only a negligible fraction of thermal IR observations. However, wide-field exposures, as well as long medium-field exposures,would be affected at the 3% level during the first and last hours of the night. Furthermore, ultra-wide imaging exposures on a very large telescope (eg NSF's Rubin Observatory, LSST), would be significantly affected, with 30 to 40% of such exposures being compromised during the first and last hours of the night. Coordination between the astronomical community, satellites companies, and government agencies is therefore critical to minimise and mitigate the effect on astronomical observations, in particular on survey telescopes., Comment: 12 pages, 11 figures, 2 tables; accepted for publication in Astronomy and Astrophysics

Comet 49P/Arend-Rigaux is a well known low-activity Jupiter Family comet. Despite the low activity, we have witnessed outgassing activity in 1992, 2004, and 2012. In 2012 a broad tail-like feature (PA$\sim270^\circ, \sim2.3\times10^5$ km) and a narrow jet-like feature (PA$\sim180^\circ, \sim9.3\times10^4$ km) were seen simultaneously. Using Finson-Probstein (FP) dust dynamical models we determine: grain sizes released in each event; duration of activity; when activity peaked; and velocity of the dust particles, allowing us to make comparisons between the events. We find that the tail feature in 2012 is similar to the tail in 1992 with large grains (40-4000 $\mu$m) peaking in activity near perihelion with a long outgassing duration greater than 150 days. The jet feature from 2012, however, is more similar to the 2004 event which we model with small grains (1-8 $\mu$m) with a short duration of activity ($\sim$1 month). The main difference between these two features is that the 2004 event occurs prior to perihelion, while the 2012 event is post-perihelion. We use the grain sizes from the FP models to constrain ice sublimation models. Between 1985 and 2018 we cover 6 apparitions with 26 nights of our own observations plus data from the literature and the Minor Planet Center, which together, allow us to model the heliocentric light curve. We find that the models are consistent with H$_2$O ice sublimation as the volatile responsible for driving activity over most of the active phases and a combination of H$_2$O and CO$_2$ ices are responsible for driving activity near perihelion. We measure the fractional active area over time for H$_2$O and discover that the activity decreases from an average active area of $\sim3\%$ to $\sim0.2\%$. This secular decrease in activity implies that the comet is becoming depleted of volatiles and is in the process of transitioning to a dormant or dead state., Comment: 18 pages, 8 figures, Accepted for publication by Icarus

We report on the catastrophic disintegration of P/2016 G1 (PANSTARRS), an active asteroid, in April 2016. Deep images over three months show that object is constituted by a central concentration of fragments surrounded by an elongated coma, and presents previously unreported sharp arc-like and a narrow linear features. The morphology and evolution of these characteristics independently point toward a brief event on 2016 March 6. The arc and the linear feature can be reproduced by large particles on a ring, moving at ~2.5 m/s. The expansion of the ring defines a cone with a ~40deg half-opening. We propose that G1 was hit by a small object which caused its (partial or total) disruption, and that the ring corresponds to large fragments ejected during the final stages of the crater formation.}

On 2019 January 5 a streamer associated with the 4--10 km main-belt asteroid (6478)~Gault was detected by the ATLAS sky survey, a rare discovery of activity around a main-belt asteroid. Archival data from ATLAS and Pan-STARRS1 show the trail in early December 2018, but not between 2010 and January 2018. The feature has significantly changed over one month, perfectly matching predictions of pure dust dynamical evolution and changes in observing geometry for a short release of dust around 2018 October 28. Follow-up observations with HST show a second narrow trail corresponding to a brief release of dust on 2018 December 30. Both releases occurred with negligible velocity. We find the dust grains to be fairly large, with power-law size distributions in the $10^{-5} - 10^{-3}$~m range and power-law indices of $\sim -1.5$. Three runs of ground-based data find a signature of $\sim 2\,\rm h$ rotation, close to the rotational limit, suggesting that the activity is the result of landslides or reconfigurations after YORP spin-up., Comment: 9 pages, 4 figures, 1 table

We show that `Oumuamua's excited spin could be in a high energy LAM state, which implies that its shape could be far from the highly elongated shape found in previous studies. CLEAN and ANOVA algorithms are used to analyze `Oumuamua's lightcurve using 818 observations over 29.3~days. Two fundamental periodicities are found at frequencies (2.77$\pm$0.11) and (6.42$\pm$0.18)~cycles/day, corresponding to (8.67$\pm$0.34)~h and (3.74$\pm$0.11)~h, respectively. The phased data show that the lightcurve does not repeat in a simple manner, but approximately shows a double minimum at 2.77~cycles/day and a single minimum at 6.42~cycles/day. This is characteristic of an excited spin state. `Oumuamua could be spinning in either the long (LAM) or short (SAM) axis mode. For both, the long axis precesses around the total angular momentum vector with an average period of (8.67$\pm$0.34)~h. For the three LAMs we have found, the possible rotation periods around the long axis are 6.58, 13.15, or 54.48~h, with 54.48~h being the most likely. `Oumuamua may also be nutating with respective periods of half of these values. We have also found two possible SAM states where `Oumuamua oscillates around the long axis with possible periods at 13.15 and 54.48~h, the latter as the most likely. In this case any nutation will occur with the same periods. Determination of the spin state, the amplitude of the nutation, the direction of the TAMV, and the average total spin period may be possible with a direct model fit to the lightcurve. We find that `Oumuamua is "cigar-shaped"', if close to its lowest rotational energy, and an extremely oblate spheroid if close to its highest energy state for its total angular momentum., Comment: 8 pages, 3 figures

Isotopic ratios in comets are critical to understanding the origin of cometary material and the physical and chemical conditions in the early solar nebula. Comet C/2015 ER61 (PANSTARRS) underwent an outburst with a total brightness increase of 2 magnitudes on the night of 2017 April 4. The sharp increase in brightness offered a rare opportunity to measure the isotopic ratios of the light elements in the coma of this comet. We obtained two high-resolution spectra of C/2015 ER61 with UVES/VLT on the nights of 2017 April 13 and 17. At the time of our observations, the comet was fading gradually following the outburst. We measured the nitrogen and carbon isotopic ratios from the CN violet (0,0) band and found that $^{12}$C/$^{13}$C=100 $\pm$ 15, $^{14}$N/$^{15}$N=130 $\pm$ 15. In addition, we determined the $^{14}$N/$^{15}$N ratio from four pairs of NH$_2$ isotopolog lines and measured $^{14}$N/$^{15}$N=140 $\pm$ 28. The measured isotopic ratios of C/2015 ER61 do not deviate significantly from those of other comets., Comment: 4 pages, 4 figures, accepted to be published by A&A

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

Comet 103P/Hartley~2 was observed on Nov. 1-6, 2010, coinciding with the fly-by of the space probe EPOXI. The goal was to connect the large scale phenomena observed from the ground, with those at small scale observed from the spacecraft. The comet showed strong activity correlated with the rotation of its nucleus, also observed by the spacecraft. We report here the characterization of the solid component produced by this activity, via observations of the emission in two spectral regions where only grain scattering of the solar radiation is present. We show that the grains produced by this activity had a lifetime of the order of 5 hours, compatible with the spacecraft observations of the large icy chunks. Moreover, the grains produced by one of the active regions have a very red color. This suggests an organic component mixed with the ice in the grains., Comment: 11 pages, 7 figures, Icarus in press

Asteroid (175706) 1996 FG3 is a binary asteroid and the baseline target for the proposed MarcoPolo-R sample-return mission. We present thermal IR photometry obtained with the ESO VLT+VISIR together with optical photometry obtained with the ESO NTT+EFOSC2 . An absolute visual magnitude H_V = 17.833 \pm 0.024 and phase parameter G = -0.041 \pm 0.005 is derived. The Near-Earth Asteroid Thermal Model (NEATM) has been fitted to the measured fluxes to derive a geometric visual albedo p_v = 0.046 \pm 0.014, effective diameter at the observed aspect D_eff = 1.68 \pm 0.25 km, and beaming parameter {\eta} = 1.15 for phase angle {\alpha} = 11.7{\deg}. The Advanced Thermophysical Model (ATPM) has been fitted to the measured fluxes to derive a more accurate effective diameter D_eff = 1.71 \pm 0.07 km and albedo p_v = 0.044 \pm 0.004. Based on the ATPM results, assuming the same albedo for primary and secondary, we derive a primary mean spherical diameter D_p = (1.69 +0.18/-0.12) km, secondary diameter D_s = 0.51 \pm 0.03 km, and a secondary orbital semi-major axis a = (2.8 +1.7/-0.7) km. A low surface thermal inertia {\Gamma} = 120 \pm 50 J m^-2 s^-1/2 K^-1 was also derived, suggesting a dusty surface and raising questions as to the binary formation mechanism of this asteroid. These physical properties are used to predict a Yarkovsky drift in semi-major axis of (-60 +31/-45) m yr-1., Comment: 24 pages, 7 Figures. Accepted by MNRAS

As a prerequisite for cosmological studies using adaptive optics techniques, we have begun to identify and characterize faint sources in the vicinity of bright stars at high Galactic latitudes. The initial phase of this work has been a program of K_s imaging conducted with SOFI at the ESO NTT. From observations of 42 southern fields evenly divided between the spring and autumn skies, we have identified 391 additional stars and 1589 galaxies lying at separations 60" from candidate guide stars in the magnitude range 9.0 R 12.4. When analyzed as a "discrete deep field" with 131 arcmin^2 area, our dataset gives galaxy number counts that agree with those derived previously over the range 16 K_s 20.5. This consistency indicates that in the aggregate, our fields should be suitable for future statistical studies. We provide our source catalogue as a resource for users of large telescopes in the southern hemisphere., Comment: 10 pages, 7 figures, accepted by A&A; Table 3 is available at http://www.rzg.mpg.de/~ajb/data.html pending upload to CDS