Search

Context. (469219) Kamo`oalewa is a small near-Earth asteroid, which is currently a quasi-satellite of the Earth. Lightcurve measurements also reveal a rotation period of only about 30 minutes. This asteroid has been selected as the target of the Tianwen-2 sample-return mission of the China National Space Administration. Aims. The first goal of this paper is to observe and improve the orbit determination of (469219) Kamo`oalewa, and better determine the Yarkovsky effect acting on it. The second goal is to estimate the thermal inertia of the asteroid, taking advantage of an improved Yarkovsky effect determination. Methods. Our observational campaign imaged the asteroid from the Loiano Astronomical Station and from the Calar Alto Observatory, in March 2024. We also accurately re-measured a precovery detection from the Sloan Digital Sky Survey from 2004. New astrometry was later used in a 7-dimensional orbit determination, aimed at estimating both the orbital elements and the Yarkovsky effect. Thermal inertia is later studied by using the ASTERIA, a new method that is suitable to estimate thermal inertia of small asteroids. Results. We detected a semi-major axis drift of $(-67.35 \pm 4.70) \times 10^{-4}$ au My$^{-1}$ due to the Yarkovsky effect, with a high signal-to-noise ratio of 14. The new orbit solution also significantly reduced the position uncertainty for the arrival of the Tianwen-2 spacecraft. By using different models for the physical parameters of Kamo`oalewa, the ASTERIA model estimated the thermal inertia at $\Gamma = 150^{+90}_{-45}$ J m$^{-2}$ K$^{-1}$ s$^{-1/2}$ or $\Gamma = 181^{+95}_{-60}$ J m$^{-2}$ K$^{-1}$ s$^{-1/2}$., Comment: Accepted for publication in Astronomy & Astrophysics

Asteroid Didymos, recently targeted by the NASA DART mission, is also planned to be visited by the ESA Hera mission. The main goal of the DART mission was to impact Dimorphos, the small satellite of Didymos, which was accomplished in September 2022. This collision altered the Didymos-Dimorphos system, generating a notable quantity of ejecta that turned Dimorphos into an active asteroid, with some ejecta potentially settling on the surfaces of both components. This prompts the investigation into the extent of post-impact surface alterations on these bodies compared to their original states. The purpose of this study is to evaluate the pre-impact thermal inertia of Didymos independently. We employed ASTERIA, an alternative to conventional thermophysical modeling, to estimate the surface thermal inertia of Didymos. The approach is based on a model-to-measurement comparison of the Yarkovsky effect-induced drift on the orbital semi-major axis. These results, alongside existing literature, enable an evaluation of the impact-induced alterations in Didymos's thermal inertia. Our nominal estimate with a constant thermal inertia model stands at $\Gamma = 211_{-55}^{+81}$ J m$^{-2}$ K$^{-1}$ s$^{-1/2}$, while assuming it varies with the heliocentric distance with an exponent of $-0.75$ thermal inertia of Didymos is found to be $258_{-63}^{+94}$ J m$^{-2}$ K$^{-1}$ s$^{-1/2}$. Subsequent verification confirmed that this result is robust against variations in unknown physical parameters. The thermal inertia estimates for Didymos align statistically with values reported in the literature, derived from both pre- and post-impact data. The forthcoming Hera mission will provide an opportunity to corroborate these findings further. Additionally, our results support the hypothesis that the thermal inertia of near-Earth asteroids is generally lower than previously expected., Comment: 11 pages, 2 figures, 2 tables

Thermal inertia estimates are available for a limited number of a few hundred objects, and the results are practically solely based on thermophysical modeling (TPM). We present a novel thermal inertia estimation method, Asteroid Thermal Inertia Analyzer (ASTERIA). The core of the ASTERIA model is the Monte Carlo approach, based on the Yarkovsky drift detection. We validate our model on asteroid Bennu plus ten well-characterized near-Earth asteroids (NEAs) for which a good estimation of the thermal inertia from the TPM exists. The tests show that the ASTERIA provides reliable results consistent with the literature values. The new method is independent from the TPM, allowing an independent verification of the results. As the Yarkovsky effect is more pronounced in small asteroids, the noteworthy advantage of the ASTERIA compared to the TPM is the ability to work with smaller asteroids for which TPM typically lacks the input data. We used the ASTERIA to estimate the thermal inertia of 38 NEAs, with 31 of them being sub-km asteroids. Twenty-nine objects in our sample are characterized as Potentially Hazardous Asteroids. On the limitation side, the ASTERIA is somewhat less accurate than the TPM. The applicability of our model is limited to NEAs, as the Yarkovsky effect is yet to be detected in main-belt asteroids. However, we can expect a significant increase in high-quality measurements of the input parameters relevant to the ASTERIA with upcoming surveys. This will surely increase the reliability of the results generated by the ASTERIA and widen the model's applicability., Comment: 28 pages, 7 figures, 11 tables; Accepted for publication in The Planetary Science Journal

Asteroids smaller than about 100 meters are observed to rotate very fast, with periods often much shorter than the critical limit of 2.2 h. Some of these super-fast rotators can also achieve a very large semi-major axis drift induced by the Yarkovsky effect, that in turn, is determined by internal and surface physical properties. We consider the small super-fast rotating near-Earth asteroid 2016 GE1. This object rotates in just 34 seconds, and a large Yarkovsky effect has been determined from astrometry. Here we aim to constrain the thermal inertia of the surface of this extreme object. We used a recently developed statistical method to determine the thermal properties of near-Earth asteroids. The method is based on the comparison between the observed and the modelled Yarkovsky effect, and the thermal conductivity (inertia) is determined by a Monte Carlo approach. Parameters of the Yarkovsky effect model are either fixed if their uncertainty is negligible, modelled with a Gaussian distribution of the errors if they are measured, or deduced from general properties of the population of near-Earth asteroids when they are unknown. Using a well-established orbit determination procedure, we determined the Yarkovsky effect on 2016 GE1, and verified a significant semi-major axis drift rate. Using a statistical method, we showed that this semi-major axis drift rate could be explained only by low thermal inertia values below 100 J m$^{-2}$ K$^{-1}$ s$^{-1/2}$: namely, 90\% of the probability density function of the model outcomes is contained at values smaller than 100 J m$^{-2}$ K$^{-1}$ s$^{-1/2}$. We propose two possible interpretations for the extremely low values: a high porosity or a cracked surface, or a thin layer of fine regolith on the surface. Though this seems unexpected in either case, it opens up the possibility of a subclass of low thermal inertia, super-fast rotating asteroids., Comment: Accepted for publication on A&A

We present the physical and dynamical properties of the recently discovered active asteroid (248370) 2005 QN173 (aka 433P). From our observations, we derived two possible rotation period solutions of 2.7 and 4.1 hours. The corresponding light curve amplitudes computed after correcting for the effect of coma are 0.28 and 0.58 mag, respectively. Both period solutions are shorter than the critical rotation limit computed for a strengthless triaxial ellipsoid, suggesting that rotation mass shedding should at least partly be responsible for the observed activity. We confirm that the activity level is fading further, but at a very modest rate of only 0.006 mag/day, still also compatible with sublimation-driven activity. We found that 248370 likely belongs to the Themis asteroid family, making it a fourth main-belt comet associated with this group. Orbital characteristics of 248370 are also consistent with its origin in the young 288P cluster of asteroids. The 288P cluster is associated with its namesake main-belt comet, providing an exciting possibility for a comparative analysis of intriguing main-belt comets 248370 and 288P., Comment: Accepted for publication in MNRAS

Collisions are one of the key processes shaping planetary systems. Asteroid families are outcomes of such collisions still identifiable across our solar system. The families provide a unique view of catastrophic disruption phenomena and have been in the focus of planetary scientists for more than a century. Most of them are located in the main belt, a ring of asteroids between Mars and Jupiter. Here we review the basic properties of the families, discuss some recent advances, and anticipate future challenges. This review pays more attention to dynamic aspects such as family identification, age determination, and long-term evolution. The text, however, goes beyond that. Especially, we cover the details of young families that see the major advances in the last years, and we anticipate it will develop even faster in the future. We also discuss the relevance of asteroid families for water-ice content in the asteroid belt and our current knowledge on links between families and main-belt comets., Comment: Review paper to appear in CeMDA's topical collection on "Main Belt Dynamics"

For studies of the long-term evolution of small Solar System objects, it is fundamental to add the Yarkovsky and Yarkovsky-O'Keefe-Radzievskii-Paddack (YORP) effects in the dynamical model. Still, implementations of these effects in publicly available N-body codes are often lacking, or the effects are implemented using significantly simplified models. In this paper, we present an implementation of the coupled Yarkovsky/YORP effects in the mercury and OrbFit $N$-body codes. Along with these two effects, we also included the effects of non-destructive collisions and rotationally induced breakups to model the asteroid spin state properly. Given the stochastic nature of many incorporated effects, the software is suitable for statistical dynamical studies. Here we explained the scientific aspect of the implementation, while technical details will be made freely available along with the source codes., Comment: 13 pages, 5 figures

Near-Earth asteroid (469219) Kamo'oalewa (aka 2016 HO3) is an Earth co-orbital and a potential space mission target. Its short-term dynamics is characterized by a periodic switching between quasi-satellite and horseshoe configurations. Due to its small diameter of only about 36 meters, the Yarkovsky effect may play a significant role in the long-term dynamics. In this work, we addressed this issue by studying the changes in the long-term motion of Kamo'oalewa caused by the Yarkovsky effect. We used an estimation of the magnitude of the Yarkovsky effect assuming different surface compositions and introduced the semi-major axis drift by propagating orbits of test particles representing the clones of the nominal orbit. Our simulations showed that the Yarkovsky effect may cause Kamo'oalewa to exit from the Earth co-orbital region a bit faster when compared to a purely gravitational model. Nevertheless, it still could remain an Earth companion for at least 0.5 My in the future. Our results imply that Kamo'oalewa is the most stable Earth's co-orbital object known so far, not only from a short-term perspective but also on long time scales., Comment: 14 pages, 9 figures, accepted for publication in The Astronomical Journal

Asteroid families with ages younger than $1$ Myr offer an interesting possibility of studying the outcomes of asteroid disruptions that are little modified by subsequent evolutionary processes. We analyze a very young asteroid family associated with (18777) Hobson in the central part of the main belt. We aim at (i) understanding its peculiar size distribution, and (ii) setting an upper limit on the characteristic dispersal velocity at subkilometer sizes corresponding to the smallest visible Hobson members. We identified the Hobson family using an up-to-date asteroid catalog. A significant increase in the number of its known members allowed us to study their size distribution and compare it with computer simulations of catastrophic disruptions. Backward orbital integrations of the heliocentric orbits allowed us to confirm the previously suggested age of Hobson and helped to estimate limits of the ejection speed. The Hobson family has an unusual size distribution: two nearly equal-size bodies, followed by a population of smaller asteroids, whose distribution takes a characteristic power law. There are two possibilities to explain these data. Either a canonical impact onto a single parent body, requiring fine-tuned impact conditions that have not been studied so far, or an unconventional model for the parent object of the Hobson family, namely a binary with $\simeq 7-9$ km primary and a $\simeq 2.5$ km secondary. In the latter case, the primary was disrupted, leaving behind the largest remnant (18777) Hobson and a suite of subkilometer asteroids. The second largest asteroid, (57738) 2001 UZ160, is the nearly intact satellite of the parent binary. The excellent convergence of nominal orbits of Hobson members sets an upper limit of $\simeq (10-20)$ m s$^{-1}$ for the initial dispersal velocity of the known members, which is consistent with both formation models., Comment: 17 pages, 12 figures, 1 table, accepted for publication in Astronomy and Astrophysics

Context: Asteroids with a diameter of up to a few dozen meters may spin very fast and complete an entire rotation within a few minutes. These small and fast-rotating bodies are thought to be monolithic objects because the gravitational force due to their small size is not strong enough to counteract the strong centripetal force caused by the fast rotation. Additionally, it is not clear whether the fast spin prevents dust and small particles (regolith) from being kept on their surface. Aims: We develop a model for constraining the thermal conductivity of the surface of the small, fast-rotating near-Earth asteroids. This model may suggest whether regolith is likely present on these objects. Methods: Our approach is based on the comparison of the measured Yarkovsky drift and a predicted value using a theoretical model that depends on the orbital, physical and thermal parameters of the object. The necessary parameters are either deduced from statistical distribution derived for near-Earth asteroids population or determined from observations with associated uncertainty. With this information, we performed Monte Carlo simulations and produced a probability density distribution for the thermal conductivity. Results: Applying our model to the superfast rotator asteroid (499998) 2011 PT, we find that the measured Yarkovsky drift can only be achieved when the thermal conductivity $K$ of the surface is low. The resulting probability density function for the conductivity is bimodal, with two most likely values being around 0.0001 and 0.005 W m$^{-1}$ K$^{-1}$. Based on this, we find that the probability that $K$ is lower than 0.1 W m$^{-1}$ K$^{-1}$ is at least 95\%. This low thermal conductivity might indicate that the surface of 2011 PT is covered with a thermal insulating layer, composed of a regolith-like material similar to lunar dust., Comment: Accepted for publication in A&A. 13 pages, 7 figures, 2 tables

The Karma asteroid family is a group of primitive asteroids in the middle part of the main belt, just at the outer edge of the 3J:1A mean-motion resonance. We obtained the list of the family members with 317 asteroids and estimated that it was formed by the catastrophic disruption of a parent body that was between 34 and 41 km in diameter. Based on the V-shape method, age of the Karma family is estimated to be about 137 Myr. A detailed dynamical map of the region combined with numerical simulations allowed us to reconstruct the long-term dynamical evolution of the family, and to identify the mechanisms responsible for this evolution. The numerical simulations successfully reproduced the main features in the orbital distribution of the family members but also showed that some regions of the Karma family could be missing. A more detailed analysis revealed that these regions likely consist of very dark objects, fainter than absolute magnitude H = 17, that have not yet been detected. Based on the obtained results, we concluded that magnitude-frequency distribution of family members up to H = 16 mag is neither affected by dynamical erosion nor observational incompleteness and therefore represents the result of collisional grinding of the original family population. Finally, we found that the Karma family have been supplying some asteroids to the near-Earth region via the 3J:1A resonance. Currently, there should about 10 family members larger than 1 km in diameter, orbiting in the near-Earth space., Comment: 12 pages, 13 Postscript figures, accepted for publication in MNRAS

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

Observational data suggest that a belt of planetesimals is expected close to the snow line in protoplanetary disks. Assuming there is such a belt in TRAPPIST-1 system, we examine possibilities of water delivery to the planets via planetesimals from the belt. The study is accomplished by numerical simulations of dynamical evolution of a hypothetical planetesimal belt. Our results show that the inner part of the belt is dynamically unstable and planetesimals located in this region are quickly scattered away, with many of them entering the region around the planets. The main dynamical mechanism responsible for the instability are close encounters with the outermost planet Trappist-1h. A low-order mean-motion resonance 2:3 with Trappist-1h, located in the same region, also contributes to the objects transport. In our nominal model, the planets have received non-negligible amount of water, with the smallest amount of $15$\% of the current Earth's water amount (EWA) being delivered to the planet 1b, while the planets Trappist-1e and Trappist-1g have received more than $60$\% of the EWA. We have found that while the estimated efficiency of water transport to the planets is robust, the amount of water delivered to each planet may vary significantly depending on the initial masses and orbits of the planets. The estimated dynamical "half-lives" have shown that the impactors' source region should be emptied in less then 1~Myr. Therefore, the obtained results suggest that transport of planetesimals through the system preferably occurs during an early phase of the planetary system evolution., Comment: 12 pages, 8 figures; MNRAS, accepted for publication

In this work we investigate the orbital distribution of interstellar objects (ISOs), observable by the future wide-field National Science Foundation Vera C. Rubin Observatory (VRO). We generate synthetic population of ISOs and simulate their ephemerides over a period of 10 years, in order to select those which may be observed by the VRO, based on the nominal characteristics of this survey. We find that the population of the observable ISOs should be significantly biased in favor of retrograde objects. The intensity of this bias is correlated with the slope of the size-frequency distribution (SFD) of the population, as well as with the perihelion distances. Steeper SFD slopes lead to an increased fraction of the retrograde orbits, and also of the median orbital inclination. On the other hand, larger perihelion distances result in more symmetric distribution of orbital inclinations. We believe that this is a result of the Holetschek effects, which is already suggested to cause observational bias in orbital distribution of long period comets. The most important implication of our findings is that an excess of retrograde orbits depends on the sizes and the perihelion distances. Therefore, the prograde/retrograde orbits ratio and the median inclination of the discovered population could, in turn, be used to estimate the SFD of the underlying true population of ISOs.

Recent dynamical analyses suggest that some Jupiter family comets (JFCs) may originate in the main asteroid belt instead of the outer solar system. This possibility is particularly interesting given evidence that icy main-belt objects are known to be present in the Themis asteroid family. We report results from dynamical analyses specifically investigating the possibility that icy Themis family members could contribute to the observed population of JFCs. Numerical integrations show that such dynamical evolution is indeed possible via a combination of eccentricity excitation apparently driven by the nearby 2:1 mean-motion resonance with Jupiter, gravitational interactions with planets other than Jupiter, and the Yarkovsky effect. We estimate that, at any given time, there may be tens of objects from the Themis family on JFC-like orbits with the potential to mimic active JFCs from the outer solar system, although not all, or even any, may necessarily be observably active. We find that dynamically evolved Themis family objects on JFC-like orbits have semimajor axes between 3.15 au and 3.40 au for the vast majority of their time on such orbits, consistent with the strong role that the 2:1 mean-motion resonance with Jupiter likely plays in their dynamical evolution. We conclude that a contribution from the Themis family to the active JFC population is plausible, although further work is needed to better characterize this contribution., Comment: 12 pages, 6 figures; accepted for publication in AJ

A specially interesting subgroup of asteroid families are young ones, i.e. those formed less than 10 Myrs ago. They have attracted a lot of attention since the first such family has been discovered in 2002. About 20 young families have been discovered so far, but new discoveries are still very important. Here we report discovery of four new young asteroid families. Three of them are located in the inner asteroid belt, namely the (525) Adelaide, (6142) Tantawi and (18429) 1994AO1, while the family of (2258) Viipuri is located in the middle part of the asteroid belt.

In this note we have shown that a newly discovered comet P/2017 S5 (ATLAS), that moves around the Sun in an asteroid-like orbit, is a member of the Theobalda asteroid family.

We report on the results of a systematic search for associated asteroid families for all active asteroids known to date. We find that 10 out of 12 main-belt comets (MBCs) and 5 out of 7 disrupted asteroids are linked with known or candidate families, rates that have ~0.1% and ~6% probabilities, respectively, of occurring by chance, given an overall family association rate of 37% for asteroids in the inner solar system. We find previously unidentified family associations between 238P/Read and the candidate Gorchakov family, 311P/PANSTARRS and the candidate Behrens family, 324P/La Sagra and the Alauda family, 354P/LINEAR and the Baptistina family, P/2013 R3-B (Catalina-PANSTARRS) and the Mandragora family, P/2015 X6 (PANSTARRS) and the Aeolia family, P/2016 G1 (PANSTARRS) and the Adeona family, and P/2016 J1-A/B (PANSTARRS) and the Theobalda family. All MBCs with family associations belong to families that contain asteroids with primitive taxonomic classifications and low average reported albedos (pV_avg < 0.10), while disrupted asteroids with family associations belong to families that contain asteroids that span wider ranges of taxonomic types and average reported albedos (0.06 < pV_avg < 0.25). These findings are consistent with MBC activity being closely correlated to composition (i.e., whether an object is likely to contain ice), while disrupted asteroid activity is not as sensitive to composition. Given our results, we describe a sequence of processes by which the formation of young asteroid families could lead to the production of present-day MBCs., Comment: 28 pages, 30 figures, accepted by AJ

We report the discovery of a new asteroid family among the dark asteroids residing in the Phocaea region the Tamara family. We make use of available physical data to separate asteroids in the region according to their surface reflectance properties, and establish the membership of the family. We determine the slope of the cumulative magnitude distribution of the family, and find it to be significantly steeper than the corresponding slope of all the asteroids in the Phocaea region. This implies that sub-kilometer dark Phocaeas are comparable in number to bright S-type objects, shedding light on an entirely new aspect of the composition of small Phocaea asteroids. We then use the Yarkovsky V-shape based method and estimate the age of the family to be 264$\pm$43 Myr. Finally, we carry out numerical simulations of the dynamical evolution of the Tamara family. The results suggest that up to 50 Tamara members with absolute magnitude H<19.4 may currently be found in the near-Earth region. Despite their relatively small number in the near-Earth space, the rate of Earth impacts by small, dark Phocaeas is non-negligible., Comment: Accepted for publication in AJ

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., Comment: Accepted by ApJ Letters, Feb. 17, 2017

In this work we explore dynamical perturbations induced by the massive asteroids Ceres and Vesta on main-belt asteroids through secular resonances. First we determine the location of the linear secular resonances with Ceres and Vesta in the main belt, using a purely numerical technique. Then we use a set of numerical simulations of fictitious asteroids to investigate the importance of these secular resonances in the orbital evolution of main-belt asteroids. We found, evaluating the magnitude of the perturbations in the proper elements of the test particles, that in some cases the strength of these secular resonances is comparable to that of known non-linear secular resonances with the giant planets. Finally we explore the asteroid families that are crossed by the secular resonances we studied, and identified several cases where the latter seem to play an important role in their post-impact evolution.

We consider the role of the dwarf planet Ceres on the secular dynamics of the asteroid main belt. Specifically, we examine the post impact evolution of asteroid families due to the interaction of their members with the linear nodal secular resonance with Ceres. First, we find the location of this resonance and identify which asteroid families are crossed by its path. Next, we summarize our results for three asteroid families, namely (1726) Hoffmeister, (1128) Astrid and (1521) Seinajoki which have irregular distributions of their members in the proper elements space, indicative of the effect of the resonance. We confirm this by performing a set of numerical simulations, showcasing that the perturbing action of Ceres through its linear nodal secular resonance is essential to reproduce the actual shape of the families., Comment: 9 pages, 4 figures; to appear in Proceedings of IAU Symposium 318

Here we report on the significant role of a so far overlooked dynamical aspect, namely a secular resonance between the dwarf planet Ceres and other asteroids. We demonstrate that this type of secular resonance can be the dominant dynamical factor in certain regions of the main asteroid belt. Specifically, we performed a dynamical analysis of the asteroids belonging to the (1726) Hoffmeister family. To identify which dynamical mechanisms are actually at work in this part of the main asteroid belt, i.e. to isolate the main perturber(s), we study the evolution of this family in time. The study is accomplished using numerical integrations of test particles performed within different dynamical models. The obtained results reveal that the post-impact evolution of the Hoffmeister asteroid family is a direct consequence of the nodal secular resonance with Ceres. This leads us to the conclusion that similar effects must exist in other parts of the asteroid belt. In this respect, the obtained results shed light on an important and entirely new aspect of the long-term dynamics of small bodies. Ceres' fingerprint in asteroid dynamics, expressed through the discovered secular resonance effect, completely changes our understanding of the way in which perturbations by Ceres-like objects affect the orbits of nearby bodies., Comment: Accepted for publication in the ApJL

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., Comment: 8 pages, 4 figures; accepted for publication in ApJ Letters

We present the results of our search for a dynamical family around the active asteroid P/2012F5 (Gibbs). By applying the hierarchical clustering method, we discover an extremely compact 9-body cluster associated with P/2012F5. The statistical significance of this newly discovered Gibbs cluster is estimated to be >99.9%, strongly suggesting that its members share a common origin. The cluster is located in a dynamically cold region of the outer main-belt at a proper semi-major axis of about 3.005 AU, and all members are found to be dynamically stable over very long time-scales. Backward numerical orbital integrations show that the age of the cluster is only 1.5 $\pm$ 0.1 Myr. Taxonomic classifications are unavailable for most of the cluster members, but SDSS spectrophotometry available for two cluster members indicate that both appear to be $Q$-type objects. We also estimate a lower limit of the size of the parent body to be about 10 km, and find that the impact event which produced the Gibbs cluster is intermediate between a cratering and a catastrophic collision. In addition, we search for new main-belt comets in the region of the Gibbs cluster by observing seven asteroids either belonging to the cluster, or being very close in the space of orbital proper elements. However, we do not detect any convincing evidence of the presence of a tail or coma in any our targets. Finally, we obtain optical images of P/2012F5, and find absolute R-band and V-band magnitudes of $H_R$ = 17.0 $\pm$ 0.1 mag and $H_V$ = 17.4 $\pm$ 0.1 mag, respectively, corresponding to an upper limit on the diameter of the P/2012F5 nucleus of about 2 km.

The number of asteroids with accurately determined orbits increases fast. The catalogs of asteroid physical observations have also increased, although the number of objects is still smaller than in the orbital catalogs. We developed a new approach to the asteroid family classification by combining the Hierarchical Clustering Method (HCM) with a method to add new members to existing families. This procedure makes use of the much larger amount of information contained in the proper elements catalogs, with respect to classifications using also physical observations for a smaller number of asteroids. Our work is based on the large catalog of the high accuracy synthetic proper elements (available from AstDyS). We first identify a number of core families; to these we attribute the next layer of smaller objects. Then, we remove all the family members from the catalog, and reapply the HCM to the rest. This gives both halo families which extend the core families and new independent families, consisting mainly of small asteroids. These two cases are discriminated by another step of attribution of new members and by merging intersecting families. By using information from absolute magnitudes, we take advantage of the larger size range in some families to analyze their shape in the proper semimajor axis vs. inverse diameter plane. This leads to a new method to estimate the family age (or ages). The results from the previous steps are then analyzed, using also auxiliary information on physical properties including WISE albedos and SDSS color indexes. This allows to solve some difficult cases of families overlapping in the proper elements space but generated by different collisional events. We analyze some examples of cratering families (Massalia, Vesta, Eunomia) which show internal structures, interpreted as multiple collisions. We also discuss why Ceres has no family.

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

Here we show an example of a young asteroid cluster located in a dynamically stable region, which was produced by partial disruption of a primitive body about 30 km in size. We estimate its age to be only 1.9 +/- 0.3 Myr, thus its post-impact evolution should have been very limited. The large difference in size between the largest object and the other cluster members means that this was a cratering event. The parent body had a large orbital inclination, and was subject to collisions with typical impact speeds higher by a factor of 2 than in the most common situations encountered in the main belt. For the first time we have at disposal the observable outcome of a very recent event to study high-speed collisions involving primitive asteroids, providing very useful constraints to numerical simulations of these events and to laboratory experiments., Comment: Accepted for publication by MNRAS

In this paper we apply different methods to examine the possibility that a small group of 24 asteroids dynamically linked to main-belt comet P/2006 VW139, recently discovered by the Pan-STARRS1 survey telescope, shares a common physical origin. By applying the Hierarchical Clustering and Backward Integration methods, we find strong evidence that 11 of these asteroids form a sub-group which likely originated in a recent collision event, and that this group includes P/2006 VW139. The objects not found to be part of the 11-member sub-group, which we designate as the P/2006 VW139 family, were either found to be dynamically unstable, or these are likely interlopers which should be expected due to the close proximity of the Themis family. As we demonstrated, statistical significance of P/2006 VW139 family is >99 per cent. We determine the age of the family to be 7.5 +/- 0.3 Myr, and estimate the diameter of the parent body to be about 11 km. Results show that the family is produced by an impact which can be best characterized as a transition from catastrophic to cratering regime. The dynamical environment of this family is studied as well, including the identification of the most influential mean motion and secular resonances in the region. Our findings make P/2006 VW139 now the second main-belt comet to be dynamically associated with a young asteroid family, a fact with important implications for the origin and activation mechanism of such objects., Comment: Accepted for publication in MNRAS

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

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

We present a new classification of families identified among the population of high-inclination asteroids. We computed synthetic proper elements for a sample of 18,560 numbered and multi-opposition objects having sine of proper inclination greater than 0.295. We considered three zones at different heliocentric distances (inner, intermediate and outer region) and used the standard approach based on the Hierarchical Clustering Method (HCM) to identify families in each zone. In doing so, we used slightly different approach with respect to previously published methodologies, to achieve a more reliable and robust classification. We also used available SDSS color data to improve membership and identify likely family interlopers. We found a total of 38 families, as well as a significant number of clumps and clusters deserving further investigation., Comment: Accepted by Icarus

The (778) Theobalda asteroid family attracted little attention so far, but our study shows that it is important in several aspects. In this paper we investigate the origin and evolution of Theobalda family. Firstly, we identify the family as a statistically relevant group in the space of synthetic proper elements. Using the hierarchical clustering method and adopted cut-off velocity of d_{cutoff}=85 m/s we found that Theobalda family currently consists of 128 members. This family is located in the outer belt, near proper semi-major axis a_{p}~3.175 au. This region is crossed by several three-body mean motion resonances which give rise to significant chaotic zones. Consequently, the majority of family members reside on chaotic orbits. Using two independent methods, chaotic chronology and backward integration, we found Theobalda family to be only 6.9 +/- 2.3 Myr old. We have also estimated, that the family was likely produced by the cratering impact on a parent body of diameter D_{PB}~78 +/- 9 km., Comment: Accepted for publication in MNRAS, 10 pages, 12 figures; http://onlinelibrary.wiley.com/doi/10.1111/j.1365-2966.2010.17051.x/full

We present a transport model that describes the orbital diffusion of asteroids in chaotic regions of the 3-D space of proper elements. Our goal is to use a simple random-walk model to study the evolution and derive accurate age estimates for dynamically complex asteroid families. To this purpose, we first compute local diffusion coefficients, which characterize chaotic diffusion in proper eccentricity (e_p) and inclination (I_p), in a selected phase-space region. Then, a Monte-Carlo-type code is constructed and used to track the evolution of random walkers (i.e. asteroids), by coupling diffusion in (e_p,I_p) with a drift in proper semi-major axis (a_p) induced by the Yarkovsky/YORP thermal effects. We validate our model by applying it to the family of (490) Veritas, for which we recover previous estimates of its age (~8.3 Myr). Moreover, we show that the spreading of chaotic family members in proper elements space is well reproduced in our random-walk simulations. Finally, we apply our model to the family of (3556) Lixiaohua, which is much older than Veritas and thus much more affected by thermal forces. We find the age of the Lixiaohua family to be 155+/-36 Myr., Comment: 11 pages, 11 figures, Accepted for publication in MNRAS

The celestial phenomenon have always been a source of wonder and interest to people, even as long ago as the ancient Egyptians. While the ancient Egyptians did not know all the things about astronomy that we do now, they had a good understanding of the some celestial phenomenon. The achievements in astronomy of ancient Egyptians are relatively well known, but we know very little about the people who made these achievements. The goal of this paper is to bring some light on the life of Senenmut, the chief architect and astronomer during the reign of Queen Hatshepsut., Comment: to appear in Publication of Astronomical Observatory of Belgrade

We presented here the orbital parameters for five visual binary stars calculated by using the new method which we named Sector Grid Search. Orbital parameters were obtained for the following stars: WDS 00152+2722 = ADS 195, WDS 02202+2949 = ADS 1780, WDS 11550$-$5606 = HIP 58106, WDS 16256$-$2327 = ADS 10049 and WDS 16256$-$2327 = ADS 10045. In addition, their masses, dynamical parallaxes and ephemerides were calculated as well., Comment: 9 pages, 5 figures, 4 tables

We present here the orbital elements for the pair BC of visual binary star WDS 18253+4846 (ADS 11344 = HD 170109) and ephemerides for the period 2008-2014. As the orbital elements of pair AB already existed we analyzed mass ratio between these two pairs and concluded that components C is the biggest one in this triple system., Comment: Accepted for publication in "The Observatory" v3 - references added