Search

The proper elements of asteroids are obtained from the instantaneous orbital elements by removing periodic oscillations produced by gravitational interactions with planets. They are unchanging in time, at least if chaotic dynamics and non-gravitational forces could be ignored, and can therefore be used to identify fragments of major collisions (asteroid families) that happened eons ago. Here we present a new catalog of proper elements for 1.25 million main belt asteroids. We explain the methodology, evaluate uncertainties, and discuss how the new catalog can be used to identify asteroid families. A systematic search for families yielded 153 cases not reported in Nesvorn\'y at al. (2015) -- 17 of these cases were identified in various other publications, 136 cases are new discoveries. There are now 274 families in the asteroid belt in total (plus a handful of families in the resonant Hilda population). We analyzed several compact families in detail. The new family around the middle belt asteroid (9332) 1990SB1 (9 members) is the youngest family found so far (estimated formation only 16-17 kyr ago). New families (1217) Maximiliana, (6084) Bascom, (10164) Akusekijima and (70208) 1999RX33 all formed 0.5-2.5 Myr ago. The (2110) Moore-Sitterly family is a close pair of relatively large bodies, 2110 and 44612, and 15 small members all located sunwards from 2110 and 44612, presumably a consequence of the Yarkovsky drift over the estimated family age (1.2-1.5 Myr). A systematic characterization of the new asteroid families is left for future work., Comment: ApJS, in press

Studies of micrometeorites in mid-Ordovician limestones and Earth's impact craters indicate that our planet witnessed a massive infall of ordinary L chondrite material 466 million years (My) ago (Heck et al. 2017, Schmieder & Kring 2020, Kenkmann 2021) that may have been at the origin of the first major mass extinction event (Schmitz et al. 2019). The breakup of a large asteroid in the main belt is the likely cause of this massive infall. In modern times, material originating from this breakup still dominates meteorite falls (>20% of all falls) (Swindle et al. 2014). Here, we provide spectroscopic observations and dynamical evidence that the Massalia collisional family is the only plausible source of this catastrophic event and of the most abundant class of meteorites falling on Earth today. It is suitably located in the inner belt, at low-inclination orbits, which corresponds to the observed distribution of L-chondrite-like near-Earth objects (NEOs) and of interplanetary dust concentrated at 1.4 degrees (Sykes 1990, Reach et al. 1997)., Comment: 35 pages, 11 pages, under revision

Jupiter Trojans (JTs) librate about the Lagrangian stationary centers L4 and L5 associated with this planet on a typically small-eccentricity and moderate-inclination heliocentric orbits. The physical and orbital properties of JTs provide important clues about the dynamical evolution of the giant planets in the early Solar System, as well as populations of planetesimals in their source regions. Here we use decade long observations from the Catalina Sky Survey (station G96) to determine the bias-corrected orbital and magnitude distributions of JTs. We distinguish the background JT population, filling smoothly the long-term stable orbital zone about L4 and L5 points, and collisional families. We find that the cumulative magnitude distribution of JTs (the background population in our case) has a steep slope for $H\leq 9$, followed with a moderately shallow slope till $H\simeq 14.5$, beyond which the distribution becomes even shallower. At $H=15$ we find a local power-law exponent $0.38\pm 0.01$. We confirm the asymmetry between the magnitude limited background populations in L4 and L5 clouds characterized by a ratio $1.45\pm 0.05$ for $H<15$. Our analysis suggests an asymmetry in the inclination distribution of JTs, with the L4 population being tighter and the L5 population being broader. We also provide a new catalog of the synthetic proper elements for JTs with an updated identification of statistically robust families (9 at L4, and 4 at L5). The previously known Ennomos family is found to consist of two, overlapping Deiphobus and Ennomos families., Comment: 46 pages, 36 figures, 4 tables, accepted for publication in The Astronomical Journal

We use observations from the Catalina Sky Survey (CSS) to determine the bias-corrected population of small members in four very young families down to sizes equivalent to several hundred meters. Using the most recent catalog of known asteroids, we identified members from four young families for which the population has grown appreciably over recent times. A large fraction of these bodies have also been detected by CSS. We used synthetic populations of asteroids, with their magnitude distribution controlled by a small number of parameters, as a template for the bias-corrected model of these families. Applying the known detection probability of the CSS observations, we could adjust these model parameters to match the observed (biased) populations in the young families. In the case of three families, Datura, Adelaide, and Rampo, we find evidence that the magnitude distribution transitions from steep to shallow slopes near $300$ to $400$ meters. Conversely, the Hobson family population may be represented by a single power-law model. The Lucascavin family has a limited population; no new members have been discovered over the past two decades. We consider a model of parent body rotational fission with the escaping secondary tidally split into two components (thereby providing three members within this family). In support of this idea, we find that no other asteroid with absolute magnitude $H\leq 18.3$ accompanies the known three members in the Lucascavin family. A similar result is found for the archetypal asteroid pair Rheinland--Kurpfalz., Comment: 32 pages, 27 figures, accepted for publication in Astronomy and Astrophysics

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 $\nu_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 $\theta \simeq 0^\circ$ for $a<2.5$~au and $\theta \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., Comment: AJ

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

Asteroid pairs, two objects that are not gravitationally bound to one another, but share a common origin, have been discovered in the Main belt and Hungaria populations. Such pairs are of major interest, as the study of their evolution under a variety of dynamical influences can indicate the time since the pair was created. To date, no asteroid pairs have been found in the Jovian Trojans, despite the presence of several binaries and collisional families in the population. The search for pairs in the Jovian Trojan population is of particular interest, given the importance of the Trojans as tracers of planetary migration during the Solar system's youth. Here we report a discovery of the first pair, (258656) 2002~ES$_{76}$ and 2013~CC$_{41}$, in the Jovian Trojans. The two objects are approximately the same size and are located very close to the L4 Lagrange point. Using numerical integrations, we find that the pair is at least $360$~Myr old, though its age could be as high as several Gyrs. The existence of the (258656) 2002~ES$_{76}$--2013~CC$_{41}$ pair implies there could be many such pairs scattered through the Trojan population. Our preferred formation mechanism for the newly discovered pair is through the dissociation of an ancient binary system, triggered by a sub-catastrophic impact, but we can not rule out rotation fission of a single object driven by YORP torques. A by-product of our work is an up-to-date catalog of Jovian Trojan proper elements, which we have made available for further studies., Comment: 15 figs. Accepted for publication in MNRAS

We present radiation hydrodynamic simulations in which binary planets form by close encounters in a system of several super-Earth embryos. The embryos are embedded in a protoplanetary disk consisting of gas and pebbles and evolve in a region where the disk structure supports convergent migration due to Type I torques. As the embryos accrete pebbles, they become heated and thus affected by the thermal torque (Ben\'{i}tez-Llambay et al. 2015) and the hot-trail effect (Chrenko et al. 2017) which excites orbital eccentricities. Motivated by findings of Eklund & Masset (2017), we assume the hot-trail effect operates also vertically and reduces the efficiency of inclination damping. Non-zero inclinations allow the embryos to become closely packed and also vertically stirred within the convergence zone. Subsequently, close encounters of two embryos assisted by the disk gravity can form transient binary planets which quickly dissolve. Binary planets with a longer lifetime $\sim$$10^{4}$ yr form in 3-body interactions of a transient pair with one of the remaining embryos. The separation of binary components generally decreases in subsequent encounters and due to pebble accretion until the binary merges, forming a giant planet core. We provide an order-of-magnitude estimate of the expected occurrence rate of binary planets, yielding one binary planet per $\simeq$$2$--$5\times10^{4}$ planetary systems. Therefore, although rare, the binary planets may exist in exoplanetary systems and they should be systematically searched for., Comment: Accepted for publication in ApJ, 18 pages, 10 figures

In order to fully understand the shapes of asteroids families in the 3-dimensional space of the proper elements $(a_{\rm p}, e_{\rm p}, \sin I_{\rm p})$ it is necessary to compare observed asteroids with N-body simulations. To this point, we describe a rigorous yet simple method which allows for a selection of the observed asteroids, assures the same size-frequency distribution of synthetic asteroids, accounts for a background population, and computes a $\chi^2$ metric. We study the Eos family as an example, and we are able to fully explain its non-isotropic features, including the distribution of pole latitudes $\beta$. We confirm its age $t = (1.3\pm0.3)\,{\rm Gyr}$; while this value still scales with the bulk density, it is verified by a Monte-Carlo collisional model. The method can be applied to other populous families (Flora, Eunomia, Hygiea , Koronis, Themis, Vesta, etc.)., Comment: 9 pages, 12 figures

Protoplanets of Super-Earth sizes may get trapped in convergence zones for planetary migration and form gas giants there. These growing planets undergo accretion heating, which triggers a hot-trail effect that can reverse migration directions, increase eccentricities and prevent resonant captures (Chrenko et al. 2017). We study populations of embryos accreting pebbles using Fargo-Thorin 2D hydrocode. We find that embryos in a disk with high surface density ($\Sigma_0 = 990\,{\rm g}\,{\rm cm}^{-2}$) undergo `unsuccessful' two-body encounters which do not lead to a merger. Only when a 3rd protoplanet arrives to the convergence zone, three-body encounters lead to mergers. For a low-viscosity disk ($\nu = 5\times10^{13}\,{\rm cm}^2\,{\rm s}^{-1}$) a massive coorbital is a possible outcome, for which a pebble isolation develops and the coorbital is stabilised. For more massive protoplanets ($5\,M_\oplus$), the convergence radius is located further out, in the ice-giant zone. After a series of encounters, there is an evolution driven by a dynamical torque of a tadpole region, which is systematically repeated several times, until the coorbital configuration is disrupted and planets merge. This may be a pathway how to solve the problem that coorbitals often form in simulations but they are not observed in nature. In contrast, the joint evolution of 120 low-mass protoplanets ($0.1\,M_\oplus$) reveals completely different dynamics. The evolution is no longer smooth, but rather a random walk. This is because the spiral arms, developed in the gas disk due to Lindblad resonances, overlap with each other and affect not only a single protoplanet but several in the surroundings. Our hydrodynamical simulations may have important implications for N-body simulations of planetary migration that use simplified torque prescriptions and are thus unable to capture protoplanet dynamics in its full glory., Comment: 15 pages, 25 figures, Astronomy & Astrophysics, accepted

We present the first release of the Czech Variable star catalogue that currently contains 1228 stars whose variability was discovered by 60 Czech observers. The catalogue contains confirmed variable stars of various types, but also candidates. We give precise coordinates, cross identification with other catalogues, information about constellation, variability type, brightness, light elements, name of the discoverer and year of discovery. In eighty-eight percent of stars the variability type is estimated, for more than 60 % of the stars the light ephemerides are given., Comment: 42 pages, 4 tables, 2 figures, Published in OEJV

We present physical models, i.e. convex shapes, directions of the rotation axis, and sidereal rotation periods, of 18 asteroids out of which 10 are new models and 8 are refined models based on much larger data sets than in previous work. The models were reconstructed by the lightcurve inversion method from archived publicly available lightcurves and our new observations with BlueEye600 robotic observatory. One of the new results is the shape model of asteroid (1663)~van~den~Bos with the rotation period of 749\,hr, which makes it the slowest rotator with known shape. We describe our strategy for target selection that aims at fast production of new models using the enormous potential of already available photometry stored in public databases. We also briefly describe the control software and scheduler of the robotic observatory and we discuss the importance of building a database of asteroid models for studying asteroid physical properties in collisional families.

Context: Planetary embryos can continue to grow by pebble accretion until they become giant planet cores. Simultaneously, these embryos mutually interact and also migrate due to torques arising from the protoplanetary disk. Aims: Our aim is to investigate how pebble accretion alters the orbital evolution of embryos undergoing the Type-I migration. In particular, we study whether they establish resonant chains, whether these chains are prone to instabilities and if giant planet cores form through embryo merging, thus occurring more rapidly than by pebble accretion alone. Methods: For the first time, we perform self-consistent global-scale radiative hydrodynamic simulations of a two-fluid protoplanetary disk consisting of gas and pebbles, the latter being accreted by embedded embryos. Accretion heating, along with other radiative processes, is accounted for to correctly model the Type-I migration. Results: We track the evolution of four super-Earth-like embryos, initially located in a region where the disk structure allows for a convergent migration. Generally, embryo merging is facilitated by rapidly increasing embryo masses and breaks the otherwise oligarchic growth. Moreover, we find that the orbital eccentricity of each embryo is considerably excited ($\simeq0.03$) due to the presence of an asymmetric underdense lobe of gas, a so-called `hot trail', produced by accretion heating of the embryo's vicinity. Eccentric orbits lead the embryos to frequent close encounters and make resonant locking more difficult. Conclusions: Embryo merging typically produces one massive core ($\gtrsim 10\,\mathrm{M_{E}}$) in our simulations, orbiting near $10\,\mathrm{AU}$. Pebble accretion is naturally accompanied by occurrence of eccentric orbits which should be considered in future efforts to explain the structure of exoplanetary systems., Comment: Accepted for publication in A&A, 25 pages, 17 figures

In this work, we analyze Jovian Trojans in the space of suitable resonant elements and we identify clusters of possible collisional origin by two independent methods: the hierarchical clustering and a so-called "randombox". Compared to our previous work (Bro\v{z} and Rozehnal 2011), we study a twice larger sample. Apart from Eurybates, Ennomos and $1996\,\rm RJ$ families, we have found three more clusters --- namely families around asteroids (20961)~Arkesilaos, (624)~Hektor in the $L_4$ libration zone and (247341)~$2001\,\rm UV_{209}$ in $L_5$. The families fulfill our stringent criteria, i.e. a high statistical significance, an albedo homogeneity and a steeper size-frequency distribution than that of background. In order to understand their nature, we simulate their long term collisional evolution with the Boulder code (Morbidelli et al. 2009) and dynamical evolution using a modified SWIFT integrator (Levison and Duncan, 1994). Within the framework of our evolutionary model, we were able to constrain the the age of the Hektor family to be either 1 to 4 Gyr or, less likely, 0.1 to 2.5 Gyr, depending on initial impact geometry. Since (624) Hektor itself seems to be a bilobed--shape body with a satellite (Marchis et al. 2014), i.e. an exceptional object, we address its association with the D--type family and we demonstrate that the moon and family could be created during a single impact event. We simulated the cratering event using a Smoothed Particle Hydrodynamics (SPH, Benz and Asphaug, 1994). This is also the first case of a family associated with a D--type parent body., Comment: accepted to MNRAS

In this work, we construct a new model for the collisional evolution of the main asteroid belt. Our goals are to test the scaling law of Benz and Asphaug (1999) and ascertain if it can be used for the whole belt. We want to find initial size-frequency distributions (SFDs) for the considered six parts of the belt (inner, middle, 'pristine', outer, Cybele zone, high-inclination region) and to verify if the number of synthetic asteroid families created during the simulation matches the number of observed families as well. We used new observational data from the WISE satellite (Masiero et al., 2011) to construct the observed SFDs. We simulate mutual collisions of asteroids with a modified version of the Boulder code (Morbidelli et al., 2009), where the results of hydrodynamic (SPH) simulations of Durda et al. (2007) and Benavidez et al. (2012) are included. Because material characteristics can significantly affect breakups, we created two models - for monolithic asteroids and for rubble-piles. To explain the observed SFDs in the size range D = 1 to 10 km we have to also account for dynamical depletion due to the Yarkovsky effect. The assumption of (purely) rubble-pile asteroids leads to a significantly worse fit to the observed data, so that we can conclude that majority of main-belt asteroids are rather monolithic. Our work may also serve as a motivation for further SPH simulations of disruptions of smaller targets (with a parent body size of the order of 1 km)., Comment: Accepted for publication in Icarus

We study K-type asteroids in the broad surroundings of the Eos family because they seem to be intimately related, according to their colours measured by the Sloan Digital Sky Survey. Such `halos' of asteroid families have been rarely used as constraints for dynamical studies to date. We explain its origin as bodies escaping from the family `core' due to the Yarkovsky semimajor-axis drift and interactions with gravitational resonances, mostly with the 9/4 mean-motion resonance with Jupiter at 3.03 AU. Our N-body dynamical model allows us to independently estimate the age of the family 1.5 to 1.9 Gyr. This is approximately in agreement with the previous age estimate by Vokrouhlick\'y et al. (2006) based on a simplified model (which accounts only for changes of semimajor axis). We can also constrain the geometry of the disruption event which had to occur at the true anomaly f =~ 120 to 180 deg., Comment: submitted to Icarus

We study orbital and physical properties of Trojan asteroids of Jupiter. We try to discern all families previously discussed in literature, but we conclude there is only one significant family among Trojans, namely the cluster around asteroid (3548) Eurybates. It is the only cluster, which has all of the following characteristics: (i) it is clearly concentrated in the proper-element space; (ii) size-frequency distribution is different from background asteroids; (iii) we have a reasonable collisional/dynamical model of the family. Henceforth, we can consider it as a real collisional family. We also report a discovery of a possible family around the asteroid (4709) Ennomos, composed mostly of small asteroids. The asteroid (4709) Ennomos is known to have a very high albedo $p_V \simeq 0.15$, which may be related to a hypothetical cratering event which exposed ice (Fern\'andez et al. 2003). The relation between the collisional family and the exposed surface of the parent body is a unique opportunity to study the physics of cratering events. However, more data are needed to confirm the existence of this family and its relationship with Ennomos.

Asteroids residing in the first-order mean motion resonances with Jupiter hold important information about the processes that set the final architecture of giant planets. Here we revise current populations of objects in the J2/1 (Hecuba-gap group), J3/2 (Hilda group) and J4/3 (Thule group) resonances. The number of multi-opposition asteroids found is 274 for J2/1, 1197 for J3/2 and 3 for J4/3. Using both hierarchical clustering technique and colour identification we characterise a collisionally-born asteroid family around the object (1911) Schubart in the J3/2 resonance. There is also a looser cluster around the largest asteroid (153) Hilda. Using N-body numerical simulations we prove that the Yarkovsky effect (infrared thermal emission from the surface of asteroids) causes a systematic drift in eccentricity for resonant asteroids, while their semimajor axis is almost fixed due to the strong coupling with Jupiter. This is a different mechanism from main belt families, where the Yarkovsky drift affects basically the semimajor axis. We use the eccentricity evolution to determine the following ages: (1.7+-0.7) Gyr for the Schubart family and >~4 Gyr for the Hilda family. We also find that collisionally-born clusters in the J2/1 resonance would efficiently dynamically disperse. The steep size distribution of the stable population inside this resonance could thus make sense if most of these bodies are fragments from an event older than ~1 Gyr. Finally, we test stability of resonant populations during Jupiter's and Saturn's crossing of their mutual mean motion resonances. In particular we find primordial objects in the J3/2 resonance were efficiently removed from their orbits when Jupiter and Saturn crossed their 1:2 mean motion resonance.

Tens of thousands of sparse-in-time lightcurves from astrometric projects are publicly available. We investigate these data and use them in the lightcurve inversion method to derive new asteroid models. By having a greater number of models with known physical properties, we can gain a better insight into the nature of individual objects and into the whole asteroid population. We use sparse photometry from selected observatories from the AstDyS database, either alone or in combination with dense lightcurves, to determine new asteroid models by the lightcurve inversion method. We present 80 new asteroid models derived from combined data sets where sparse photometry is taken from the AstDyS database and dense lightcurves are from the Uppsala Asteroid Photometric Catalogue (UAPC) and from several individual observers. For 18 asteroids, we present updated shape solutions based on new photometric data. The addition of new models increases the total number of models derived by the lightcurve inversion method to ~200. We also present a simple statistical analysis of physical properties of asteroids where we look for possible correlations between various physical parameters with an emphasis on the spin vector. We present the observed and de-biased distributions of ecliptic latitudes with respect to different size ranges of asteroids as well as a simple theoretical model of the latitude distribution and then compare its predictions with the observed distributions. From this analysis we find that the latitude distribution of small asteroids (D < 30 km) is clustered towards ecliptic poles and can be explained by the YORP thermal effect while the latitude distribution of larger asteroids (D > 60 km) exhibits an evident excess of prograde rotators, probably of primordial origin., Comment: 16 pages, 12 figures, 6 tables, accepted for publication in Astronomy & Astrophysics

The Yarkovsky effect is a thermal radiation force which causes objects to undergo semimajor axis drift and spin up/down as a function of their spin, orbit, and material properties. This mechanism can be used to (i) deliver asteroids (and meteoroids) with diameter D < 20 km from their parent bodies in the main belt to chaotic resonance zones capable of transporting this material to Earth-crossing orbits, (ii) disperse asteroid families, with drifting bodies jumping or becoming trapped in mean-motion and secular resonances within the main belt, and (iii) modify the rotation rates of asteroids a few km in diameter or smaller enough to explain the excessive number of very fast and very slow rotators among the small asteroids. Accordingly, we suggest that nongravitational forces, which produce small but meaningful effects on asteroid orbits and rotation rates over long timescales, should now be considered as important as collisions and gravitational perturbations to our overall understanding of asteroid evolution.