Search

Born as ice-rich planetesimals, cometary nuclei were gravitationally scattered onto their current orbits in the Kuiper Belt and the Oort Cloud during the giant planets' dynamical instability in the early stages of our Solar System's history. Here, we model the thermal evolution of planetesimals during and after the giant planet instability. We couple an adapted thermal evolution model to orbital trajectories provided by \textit{N}-body simulations to account for the planetesimals' orbital evolution, a parameter so far neglected by previous thermal evolution studies. Our simulations demonstrate intense thermal processing in all planetesimal populations, concerning mainly the hyper-volatile ice content. Unlike previous predictions, we show that hyper-volatile survival was possible in a significant number of planetesimals of the Kuiper Belt and the Oort Cloud. Planetesimals ejected into the interstellar space proved to be the most processed, while planetesimals ending in the Oort Cloud were the least processed population. We show that processing differences between populations are a direct consequence of their orbital evolution patterns, and that they provide a natural explanation for the observed variability in the abundance ratios of CO on cometary populations and on the recent observations of long-distance CO-driven activity on inbound Long-period Comets., Comment: 32 pages, 8 figures

Near-sun sky twilight observations allow for the detection of asteroid interior to the orbit of Venus (Aylos), the Earth (Atiras), and comets. We present the results of observations with the Palomar 48-inch telescope (P48)/Zwicky Transient Facility (ZTF) camera in 30 s r-band exposures taken during evening astronomical twilight from 2019 Sep 20 to 2022 March 7 and during morning astronomical twilight sky from 2019 Sep 21 to 2022 Sep 29. More than 46,000 exposures were taken in evening and morning astronomical twilight within 31 to 66 degrees from the Sun with an r-band limiting magnitude between 18.1 and 20.9. The twilight pointings show a slight seasonal dependence in limiting magnitude and ability to point closer towards the Sun, with limiting magnitude slightly improving during summer. In total, the one Aylo, (594913) 'Ayl\'o'chaxnim, and 4 Atiras, 2020 OV1, 2021 BS1, 2021 PB2, and 2021 VR3, were discovered in evening and morning twilight observations. Additional twilight survey discoveries also include 6 long-period comets: C/2020 T2, C/2020 V2, C/2021 D2, C/2021 E3, C/2022 E3, and C/2022 P3, and two short-period comets: P/2021 N1 and P/2022 P2 using deep learning comet detection pipelines. The P48/ZTF twilight survey also recovered 11 known Atiras, one Aylo, three short-period comes, two long-period comets, and one interstellar object. Lastly, the Vera Rubin Observatory will conduct a twilight survey starting in its first year of operations and will cover the sky within 45 degrees of the Sun. Twilight surveys such as those by ZTF and future surveys will provide opportunities for discovering asteroids inside the orbits of Earth and Venus., Comment: 26 pages, 13 figures, 4 tables, accepted for publication in Icarus

The process leading to the formation of the terrestrial planet remains elusive. In a previous publication, we have shown that, if the first generation of planetesimals forms in a ring at about 1 AU and the gas disk's density peaks at the ring location, planetary embryos of a few martian masses can grow and remain in the ring. In this work, we extend our simulations beyond the gas-disk stage, covering 200 Myr and accounting for the phase of giant planet instability, assumed to happen at different times. About half of the simulations form a pair of Venus and Earth analogues and, independently, about 10% form a Mars analogue. We find that the timing of the giant planet instability affects statistically the terrestrial system's excitation state and the timing of the last giant impacts. Hence a late instability (about 60 to 100 Myr after the Solar system's birth) is more consistent with a late Moon-formation time, as suggested by radioactive chronometers. However, the late veneer mass (LVM: mass accreted after the last giant impact) of Earth-sized planets suffering a giant impact after 80 My is usually an order of magnitude lower than the value inferred from geochemistry. In addition, the final angular momentum deficit (AMD) of the terrestrial planets tends to be too high. We tested the effect on the final AMD of the generation of debris during collisions and found that it is too small to change these conclusions. We argue that the best-case scenario is that the Moon-forming event occurred between 50 and 80 My, possibly just following the giant planet instability., Comment: Accepted for publication in Icarus

Understanding the origin of bright shooting stars and their meteorite samples is among the most ancient astronomy-related questions that at larger scales has human consequences [1-3]. As of today, only ${\sim}\,6\%$ of meteorite falls have been firmly linked to their sources (Moon, Mars, and asteroid (4) Vesta [4-6]). Here, we show that ${\sim}\,70\%$ of meteorites originate from three recent breakups of $D > 30\,{\rm km}$ asteroids that occurred 5.8, 7.5 and less than ${\sim}\,40$ million years ago. These breakups, including the well-known Karin family [7], took place in the prominent yet old Koronis and Massalia families and are at the origin of the dominance of H and L ordinary chondrites among meteorite falls. These young families distinguish themselves amidst all main belt asteroids by having a uniquely high abundance of small fragments. Their size-frequency distribution remains steep for a few tens of millions of years, exceeding temporarily the production of metre-sized fragments by the largest old asteroid families (e.g., Flora, Vesta). Supporting evidence includes the existence of associated dust bands [8-10], the cosmic-ray exposure ages of H-chondrite meteorites [11,12], or the distribution of pre-atmospheric orbits of meteorites [13-15]., Comment: 69 pages, 24 figures

We present a new crater chronology for Jupiter's Trojan asteroids. This tool can be used to interpret the collisional history of the bodies observed by NASA's Lucy mission. The Lucy mission will visit a total of six Trojan asteroids: Eurybates, Polymele, Orus, Leucus, and the near equal mass binary Patroclus-Menoetius. In addition, Eurybates and Polymele each have a small satellite. Here we present a prediction of Trojan cratering based on current models of how the Solar System and the objects themselves evolved. We give particular emphasis to the time lapsed since their implantation into stable regions near Jupiter's Lagrangian L4 and L4 points. We find that cratering on Trojans is generally dominated by mutual collisions, with the exception of a short period of time (~10 Myr) after implantation, in which cometary impacts may have been significant. For adopted crater scaling laws, we find that the overall spatial density of craters on Trojans is significantly lower than that of Main Belt asteroids on surfaces with similar formation ages. We also discuss specific predictions for similar-sized Eurybates and Orus, and the binary system Patroclus-Menoetius., Comment: Published on AJ

We report on a study of the mounds that dominate the appearance of Kuiper Belt Object (KBO) (486958) Arrokoth's larger lobe, named Wenu. We compare the geological context of these mounds, measure and intercompare their shapes, sizes/orientations, reflectance, and colors. We find the mounds are broadly self-similar in many respects and interpret them as the original building blocks of Arrokoth. It remains unclear why these building blocks are so similar in size, and this represents a new constrain and challenge for solar system formation models. We then discuss the interpretation of this interpretation., Comment: 24 pages, 8 figures

We studied impact processes by means of smoothed-particle hydrodynamics (SPH) simulations. The method was applied to modeling formation of main-belt families during the cometary bombardment (either early or late, ${\sim}\,3.85\,{\rm Gy}$ ago). If asteroids were bombarded by comets, as predicted by the Nice model, hundreds of asteroid families (catastrophic disruptions of diameter $D \ge 100\,{\rm km}$ bodies) should have been created, but the observed number is only 20. Therefore we computed a standard set of 125 simulations of collisions between representative $D = 100\,{\rm km}$ asteroids and high-speed icy projectiles (comets). According to our results, the largest remnant mass $M_{\rm lr}$ is similar as in low-speed collisions, due to appropriate scaling with the effective strength $Q_{\rm eff}$, but the largest fragment mass $M_{\rm lf}$ exhibits systematic differences - it is typically smaller for craterings and bigger for super-catastrophic events. This trend does not, however, explain the non-existence of old families. The respective parametric relations can be used in other statistical (Monte-Carlo) models to better understand collisions between asteroidal and cometary populations.

Massive cores of the giant planets are thought to have formed in a gas disk by accretion of pebble-size particles whose accretional cross-section is enhanced by aerodynamic gas drag [1][2]. A commonly held view is that the terrestrial planet system formed later (30-200 Myr after the dispersal of the gas disk) by giant collisions of tens of roughly Mars-size protoplanets [3]. Here we propose, instead, that the terrestrial planets formed earlier by gas-driven convergent migration of protoplanets toward $\sim\!1\,{\rm au}$ (related ref. [4] invoked a different process to concentrate planetesimals). To investigate situations in which convergent migration occurs, we developed a radiation-hydrodynamic model with realistic opacities [5][6] to determine the thermal structure of the gas and pebble disks in the terrestrial planet zone. We find that protoplanets rapidly grow by mutual collisions and pebble accretion, and gain orbital eccentricities by gravitational scattering and the hot-trail effect [7][8]. The orbital structure of the terrestrial planet system is well reproduced in our simulations, including its tight mass concentration at 0.7-1 au and the small sizes of Mercury and Mars. The early-stage protosolar disk temperature exceeds 1500 K inside 0.4 au implying that Mercury grew in a highly reducing environment, next to the evaporation lines of iron and silicates, influencing Mercury's bulk composition [9]. In a late-stage cold gas disk, accretion of icy/hydrated pebbles would contribute to Earth's water budget., Comment: original submission; final vers. is available at NA website

Very young asteroid families may record processes that accompanied their formation in the most pristine way. This makes analysis of this special class particularly interesting. We studied the very young Adelaide family in the inner part of the main belt. This cluster is extremely close to the previously known Datura family in the space of proper orbital elements and their ages overlap. As a result, we investigated the possibility of a causal relationship between the two families. We identified Adelaide family members in the up-to-date catalogue of asteroids. By computing their proper orbital elements we inferred the family structure. Backward orbital integration of selected members allowed us to determine the age of the family. The largest fragment (525) Adelaide, an S-type asteroid about $10$ km in size, is accompanied by 50 sub-kilometre fragments. This family is a typical example of a cratering event. The very tiny extent in the semi-major axis minimises chances that some significant mean motion resonances influence the dynamics of its members, though we recognise that part of the Adelaide family is affected by weak, three-body resonances. Weak chaos is also produced by distant encounters with Mars. Simultaneous convergence of longitude of node for the orbits of six selected members to that of (525) Adelaide constrains the Adelaide family age to $536\pm 12$ kyr (formal solution). While suspiciously overlapping with the age of the Datura family, we find it unlikely that the formation events of the two families are causally linked. In all likelihood, the similarity of their ages is just a coincidence., Comment: 15 pages, 9 figures, 2 tables, accepted for publication in Astronomy and Astrophysics

We combine several constraints provided by the crater records on Arrokoth and the worlds of the Pluto system to compute the size-frequency distribution (SFD) of the crater production function for craters with diameter D<10km. For this purpose, we use a Kuiper belt objects (KBO) population model calibrated on telescopic surveys, that describes also the evolution of the KBO population during the early Solar System. We further calibrate this model using the crater record on Pluto, Charon and Nix. Using this model, we compute the impact probability on Arrokoth, integrated over the age of the Solar System. This probability is then used together with other observational constraints to determine the slope of the crater-production function on Arrokoth. In addition, we use our Kuiper belt model also to compare the impact rates and velocities of KBOs on Arrokoth with those on Charon, integrated over the crater retention ages of their respective surfaces. This allows us to establish a relationship between the spatial density of sub-km craters on Arrokoth and of D~20km craters on Charon. Together, all these considerations suggest the crater production function on these worlds has a cumulative power law slope of -1.5

A recent examination of K2 lightcurves indicates that ~15% of Jupiter Trojans have very slow rotation (spin periods Ps>100 h). Here we consider the possibility that these bodies formed as equal-size binaries in the massive outer disk at ~20-30 au. Prior to their implantation as Jupiter Trojans, tight binaries tidally evolved toward a synchronous state with Ps~Pb, where Pb is the binary orbit period. They may have been subsequently dissociated by impacts and planetary encounters with at least one binary component retaining its slow rotation. Surviving binaries on Trojan orbits would continue to evolve by tides and spin-changing impacts over 4.5 Gyr. To explain the observed fraction of slow rotators, we find that at least ~15-20% of outer disk bodies with diameters 15

The New Horizons spacecraft's encounter with the cold classical Kuiper belt object (486958) Arrokoth (formerly 2014 MU69) revealed a contact-binary planetesimal. We investigate how it formed, finding it is the product of a gentle, low-speed merger in the early Solar System. Its two lenticular lobes suggest low-velocity accumulation of numerous smaller planetesimals within a gravitationally collapsing, solid particle cloud. The geometric alignment of the lobes indicates the lobes were a co-orbiting binary that experienced angular momentum loss and subsequent merger, possibly due to dynamical friction and collisions within the cloud or later gas drag. Arrokoth's contact-binary shape was preserved by the benign dynamical and collisional environment of the cold classical Kuiper belt, and so informs the accretion processes that operated in the early Solar System., Comment: Published in Science 28 Feb 2020 (First release 13 Feb 2020)

We characterize the radii and masses of the star and planets in the Kepler-59 system, as well as their orbital parameters. The star parameters are determined through a standard spectroscopic analysis, resulting in a mass of $1.359\pm 0.155\,M_\odot$ and a radius of $1.367\pm 0.078\,R_\odot$. The planetary radii obtained are $1.5\pm 0.1\,R_\oplus$ for the inner and $2.2\pm 0.1\,R_\oplus$ for the outer planet. The orbital parameters and the planetary masses are determined by the inversion of Transit Timing Variations (TTV) signals. For this, we consider two different data sets, one provided by Holczer et al. 2016, with TTVs only for the planet Kepler-59c, and the other provided by Rowe et al. 2015, with TTVs signals for both planets. The inversion method is carried out by applying an algorithm of Bayesian inference (MultiNest) combined with an efficient N-body integrator (Swift). For each of the data sets, two possible solutions are found, both having the same probability according to their corresponding Bayesian evidences. All four solutions appear to be indistinguishable within their 2-$\sigma$ uncertainties. Nevertheless, statistical analyses show that the solutions from Rowe et al. 2015 data better characterize the data. The first and second solutions identify masses of $5_{-2}^{+4}~M_{\mathrm{\oplus}}$ and $4.6_{-2.0}^{+3.6}~M_{\mathrm{\oplus}}$, and $3.0^{+0.8}_{-0.8}~M_{\mathrm{\oplus}}$ and $2.6^{+1.9}_{-0.8}~M_{\mathrm{\oplus}}$ for the inner and outer planet, respectively. This points to a system with an inner super-Earth and an outer mini-Neptune. Dynamical studies show the planets have almost co-planar orbits with small eccentricities ($e<0.1$), close but not into the 3:2 mean motion resonance. Stability analysis indicates that this configuration is stable over million years of evolution., Comment: 10 pages, 7 figures

We perform dynamical fits to the Transit Timing Variations (TTVs) of Kepler-419b. The TTVs from 17 Kepler quarters are obtained from Holczer et al (2016, ApJS 225,9). The dynamical fits are performed using the MultiNest Bayesian inference tool, coupled to an efficient symplectic N-body integrator. We find that the existing TTV data alone are able to uniquely constrain the planetary masses of Kepler-419b and c. Our estimates are in a good agreement with previous mass determinations that combined different techniques and observations, such as TTVs, radial velocity measurements, and the photoeccentric effect. As expected, however, our mass estimates have larger uncertainty. We study the global stability of the system within the parameters uncertainties to discard possible unstable solutions. We conclude that our method applied to the Holczer et al. data can provide reliable determinations of the planetary parameters. In the fore-coming work, we will use Holczer et al. (2016b) data to determine planetary parameters for a large set of Kepler systems.

We report on our study of asteroidal breakups, i.e. fragmentations of targets, subsequent gravitational reaccumulation and formation of small asteroid families. We focused on parent bodies with diameters $D_{\rm pb} = 10$ km. Simulations were performed with a smoothed-particle hydrodynamics (SPH) code combined with an efficient N-body integrator. We assumed various projectile sizes, impact velocities and impact angles (125 runs in total). Resulting size-frequency distributions are significantly different from scaled-down simulations with $D_{\rm pb} = 100$ km targets (Durda et al., 2007). We derive new parametric relations describing fragment distributions, suitable for Monte-Carlo collisional models. We also characterize velocity fields and angular distributions of fragments, which can be used as initial conditions for N-body simulations of small asteroid families. Finally, we discuss a number of uncertainties related to SPH simulations., Comment: Accepted to Icarus

The timeline of the lunar bombardment in the first Gy of the Solar System remains unclear. Some basin-forming impacts occurred 3.9-3.7Gy ago. Many other basins formed before, but their exact ages are not precisely known. There are two possible interpretations of the data: in the cataclysm scenario there was a surge in the impact rate approximately 3.9Gy ago, while in the accretion tail scenario the lunar bombardment declined since the era of planet formation and the latest basins formed in its tail-end. Here, we revisit the work of Morbidelli et al.(2012) that examined which scenario could be compatible with both the lunar crater record in the 3-4Gy period and the abundance of highly siderophile elements (HSE) in the lunar mantle. We use updated numerical simulations of the fluxes of impactors. Under the traditional assumption that the HSEs track the total amount of material accreted by the Moon since its formation, we conclude that only the cataclysm scenario can explain the data. The cataclysm should have started ~3.95Gy ago. However we show that HSEs could have been sequestered from the lunar mantle due to iron sulfide exsolution during magma ocean crystallization, followed by mantle overturn. Based on the hypothesis that the lunar magma ocean crystallized about 100-150My after Moon formation, and therefore that HSEs accumulated in the lunar mantle only after this time, we show that the bombardment in the 3-4Gy period can be explained in the accretion tail scenario. This hypothesis would also explain why the Moon appears so depleted in HSEs relative to the Earth. We also extend our analysis of the cataclysm and accretion tail scenarios to the case of Mars. The accretion tail scenario requires a global resurfacing event on Mars ~4.4Gy ago, possibly associated with the formation of the Borealis basin, and it is consistent with the HSE budget of the planet., Comment: in press in Icarus

Dynamical models of the asteroid delivery from the main belt suggest that the current impact flux of diameter D>10 km asteroids on the Earth is 0.5-1 per Gyr. Studies of the Near-Earth Asteroid (NEA) population find a much higher flux, with ~7 D>10-km asteroid impacts per Gyr. Here we show that this problem is rooted in the application of impact probability of small NEAs (1.5 per Gyr per object), whose population is well characterized, to large NEAs. In reality, large NEAs evolve from the main belt by different escape routes, have a different orbital distribution, and lower impact probabilities (0.8+/-0.3 per Gyr per object) than small NEAs. In addition, we find that the current population of two D>10 km NEAs (Ganymed and Eros) is a slight fluctuation over the long term average of 1.1+/-0.5 D>10 km NEAs in a steady state. These results have important implications for our understanding of the occurrence of the K/T-scale impacts on the terrestrial worlds., Comment: AJ, in press

Comets are icy objects that orbitally evolve from the trans-Neptunian region (the Kuiper belt and beyond) into the inner Solar System, where they are heated by solar radiation and become active due to sublimation of water ice. Here we perform end-to-end simulations in which cometary reservoirs are formed in the early Solar System and evolved over 4.5 Gyr. The gravitational effects of Planet 9 (P9), hypothesized to circle the Sun on a wide orbit, are included in some of our simulations. Different models are considered for comets to be active, including a simple assumption that comets remain active for Np(q) perihelion passages with perihelion distance q<2.5 au. The orbital distribution and number of active comets produced in our model is compared to observations. The orbital distribution of ecliptic comets (ECs) is well reproduced in models with Np(2.5)=500 and without P9. With P9, the inclination distribution of model ECs is wider than the observed one. We find that the known Halley-type comets (HTCs) have a nearly isotropic inclination distribution (with only a slight preference for prograde orbits). In our model, the HTCs appear to be an extension of the population of returning Oort-cloud comets (OCCs) to shorter orbital periods. The inclination distribution of model HTCs becomes broader with increasing Np, but the existing observational data are not good enough to constrain Np from orbital fits. Np(2.5)>1000 is required to obtain a steady-state population of large active HTCs that is consistent with observations. To fit the ratio of the returning-to-new OCCs, by contrast, our model implies that Np(2.5)<10, possibly because the detected long-period comets are smaller and much easier to disrupt than observed HTCs.

The Transit Timing Variations (TTVs) are deviations of the measured mid-transit times from the exact periodicity. One of the most interesting causes of TTVs is the gravitational interaction between planets. Here we consider a case of two planets in a mean motion resonance (orbital periods in a ratio of small integers). This case is important because the resonant interaction can amplify the TTV effect and allow planets to be detected more easily. We develop an analytic model of the resonant dynamics valid for small orbital eccentricities and use it to derive the principal TTV terms. We find that a resonant system should show TTV terms with two basic periods (and their harmonics). The resonant TTV period is proportional (m/M_*)^(-2/3), where m and M_* are the planetary and stellar masses. For m=10^(-4) M_*, for example, the TTV period exceeds the orbital period by ~2 orders of magnitude. The amplitude of the resonant TTV terms scales linearly with the libration amplitude. The ratio of the TTV amplitudes of two resonant planets is inversely proportional to the ratio of their masses. These and other relationships discussed in the main text can be used to aid the interpretation of TTV observations., Comment: The Astrophysical Journal

Asteroid families are groups of minor bodies produced by high-velocity collisions. After the initial dispersions of the parent bodies fragments, their orbits evolve because of several gravitational and non-gravitational effects,such as diffusion in mean-motion resonances, Yarkovsky and YORP effects, close encounters of collisions, etc. The subsequent dynamical evolution of asteroid family members may cause some of the original fragments to travel beyond the conventional limits of the asteroid family. Eventually, the whole family will dynamically disperse and no longer be recognizable. A natural question that may arise concerns the timescales for dispersion of large families. In particular, what is the oldest still recognizable family in the main belt? Are there any families that may date from the late stages of the Late Heavy Bombardment and that could provide clues on our understanding of the primitive Solar System? In this work, we investigate the dynamical stability of seven of the allegedly oldest families in the asteroid main belt. Our results show that none of the seven studied families has a nominally mean estimated age older than 2.7 Gyr, assuming standard values for the parameters describing the strength of the Yarkovsky force. Most "paleo-families" that formed between 2.7 and 3.8 Gyr would be characterized by a very shallow size-frequency distribution, and could be recognizable only if located in a dynamically less active region (such as that of the Koronis family). V-type asteroids in the central main belt could be compatible with a formation from a paleo-Eunomia family., Comment: 9 pages, 5 figures, 5 tables. Accepted for publication in MNRAS

Ceres is the largest and most massive body in the asteroid main belt. Observational data from the Dawn spacecraft reveal the presence of at least two impact craters about 280~km in diameter on the Ceres surface, that could have expelled a significant number of fragments. Yet, standard techniques for identifying dynamical asteroid families have not detected any Ceres family. In this work, we argue that linear secular resonances with Ceres deplete the population of objects near Ceres. Also, because of the high escape velocity from Ceres, family members are expected to be very dispersed, with a considerable fraction of km-sized fragments that should be able to reach the pristine region of the main belt, the area between the 5J:-2A and 7J:-3A mean-motion resonances, where the observed number of asteroids is low. Rather than looking for possible Ceres family members near Ceres, here we propose to search in the pristine region. We identified 156 asteroids whose taxonomy, colors, albedo could be compatible with being fragments from Ceres. Remarkably, most of these objects have inclinations near that of Ceres itself., Comment: 12 pages, 6 figures, 1 table. Accepted for publication in MNRAS

Themis is an old and statistically robust asteroid family populating the outer main belt, and resulting from a catastrophic collision that took place 2.5$\pm$1.0 Gyr ago. Within the old Themis family a young sub-family, Beagle, formed less than 10 Myr ago, has been identified. We present the results of a spectroscopic survey in the visible and near infrared range of 22 Themis and 8 Beagle families members. The Themis members investigated exhibit a wide range of spectral behaviors, while the younger Beagle family members look spectrally bluer with a smaller spectral slope variability. The best meteorite spectral analogues found for both Themis and Beagle families members are carbonaceous chondrites having experienced different degrees of aqueous alteration, prevalently CM2 but also CV3 and CI, and some of them are chondrite samples being unusual or heated. We extended the spectral analysis including the data available in the literature on Themis and Beagle families members, and we looked for correlations between spectral behavior and physical parameters using the albedo and size values derived from the WISE data. The analysis of this larger sample confirm the spectral diversity within the Themis family and that Beagle members tend to be bluer and to have an higher albedo. The differences between the two family may be partially explained by space weathering processes, which act on these primitive surfaces in a similar way than on S-type asteroids, i.e. producing reddening and darkening. However we see several Themis members having albedos and spectral slopes similar to the young Beagle members. Alternative scenarios are proposed including heterogeneity in the parent body having a compositional gradient with depth, and/or the survival of projectile fragments having a different composition than the parent body., Comment: Manuscript pages: 40; Figures: 15 ; Tables: 4 Icarus (2016),in press

The solar system has changed dramatically since its birth, and so did our understanding of it. A considerable research effort has been invested in the past decade in an attempt to reconstruct the solar system history, including the earliest stages some 4.5 billion years ago. The results indicate how several processes, such as planetary migration and dynamical instabilities, acted to relax the orbital spacing of the outer planets, and provided the needed perturbation to explain the present planetary orbits that are not precisely circular and coplanar. Here we highlight this work and illustrate the key results in a computer simulation that unifies several recently developed theories. The emerging view represents another step away from the initial perception of the solar system as part of unchanging heavens.

The 2:1 mean-motion resonance with Jupiter harbours two distinct groups of asteroids. The short-lived population is known to be a transient group sustained in steady state by the Yarkovsky semimajor axis drift. The long-lived asteroids, however, can exhibit dynamical lifetimes comparable to $4\,\mathrm{Gyr}$. They reside near two isolated islands of the phase space denoted $\mathrm{A}$ and $\mathrm{B}$, with an uneven population ratio $\mathrm{B}/\mathrm{A} \simeq 10$. The orbits of $\mathrm{A}$-island asteroids are predominantly highly inclined, compared to island $\mathrm{B}$. The size-frequency distribution is steep but the orbital distribution lacks any evidence of a collisional cluster. These observational constraints are somewhat puzzling and therefore the origin of the long-lived asteroids has not been explained so far. With the aim to provide a viable explanation, we first update the resonant population and revisit its physical properties. Using an $N$-body model with seven planets and the Yarkovsky effect included, we demonstrate that the dynamical depletion of island $\mathrm{A}$ is faster, in comparison with island $\mathrm{B}$. Then we investigate (i) the survivability of primordial resonant asteroids and (ii) capture of the population during planetary migration, following a recently described scenario with an escaping fifth giant planet and a jumping-Jupiter instability. We also model the collisional evolution of the resonant population over past $4\,\mathrm{Gyr}$. Our conclusion is that the long-lived group was created by resonant capture from a narrow part of hypothetical outer main-belt family during planetary migration. Primordial asteroids surviving the migration were probably not numerous enough to substantially contribute to the observed population., Comment: Accepted for publication in MNRAS. 20 pages, 17 figures

Asteroids formed in a dynamically quiescent disk but their orbits became gravitationally stirred enough by Jupiter to lead to high-speed collisions. As a result, many dozen large asteroids have been disrupted by impacts over the age of the Solar System, producing groups of fragments known as asteroid families. Here we explain how the asteroid families are identified, review their current inventory, and discuss how they can be used to get insights into long-term dynamics of main belt asteroids. Electronic tables of the membership for 122 notable families are reported on the Planetary Data System node., Comment: Asteroids IV chapter

The Nice model of the dynamical instability and migration of the giant planets can explain many properties of the present Solar System, and can be used to constrain its early architecture. In the jumping-Jupiter version of the Nice model, required from the terrestrial planet constraint and dynamical structure of the asteroid belt, Jupiter has encounters with an ice giant. Here we study the survival of the Galilean satellites in the jumping-Jupiter model. This is an important concern because the ice-giant encounters, if deep enough, could dynamically perturb the orbits of the Galilean satellites, and lead to implausible results. We performed numerical integrations where we tracked the effect of planetary encounters on the Galilean moons. We considered three instability cases from Nesvorny & Morbidelli (2012) that differed in the number and distribution of encounters. We found that in one case, where the number of close encounters was relatively small, the Galilean satellite orbits were not significantly affected. In the other two, the orbital eccentricities of all moons were excited by encounters, Callisto's semimajor axis changed, and, in a large fraction of trials, the Laplace resonance of the inner three moons was disrupted. The subsequent evolution by tides damps eccentricities and can recapture the moons in the Laplace resonance. A more important constraint is represented by the orbital inclinations of the moons, which can be excited during the encounters and not appreciably damped by tides. We find that one instability case taken from Nesvorny & Morbidelli (2012) clearly fails this constraint. This shows how the regular satellites of Jupiter can be used to set limits on the properties of encounters in the jumping-Jupiter model, and help us to better understand how the early Solar System evolved., Comment: The Astronomical Journal, in press

The irregular satellites of outer planets are thought to have been captured from heliocentric orbits. The exact nature of the capture process, however, remains uncertain. We examine the possibility that irregular satellites were captured from the planetesimal disk during the early Solar System instability when encounters between the outer planets occurred (Nesvorny, Vokrouhlicky & Morbidelli 2007, AJ 133; hereafter NVM07). NVM07 already showed that the irregular satellites of Saturn, Uranus and Neptune were plausibly captured during planetary encounters. Here we find that the current instability models present favorable conditions for capture of irregular satellites at Jupiter as well, mainly because Jupiter undergoes a phase of close encounters with an ice giant. We show that the orbital distribution of bodies captured during planetary encounters provides a good match to the observed distribution of irregular satellites at Jupiter. The capture efficiency for each particle in the original transplanetary disk is found to be (1.3-3.6)x10^-8. This is roughly enough to explain the observed population of jovian irregular moons. We also confirm NVM07's results for the irregular satellites of Saturn, Uranus and Neptune., Comment: Submitted to ApJ, revised version

Dawson and Murray-Clay (2012) pointed out that the inner part of the cold population in the Kuiper belt (that with semi major axis a<43.5 AU) has orbital eccentricities significantly smaller than the limit imposed by stability constraints. Here, we confirm their result by looking at the orbital distribution and stability properties in proper element space. We show that the observed distribution could have been produced by the slow sweeping of the 4/7 mean motion resonance with Neptune that accompanied the end of Neptune's migration process. The orbital distribution of the hot Kuiper belt is not significantly affected in this process, for the reasons discussed in the main text. Therefore, the peculiar eccentricity distribution of the inner cold population can not be unequivocally interpreted as evidence that the cold population formed in-situ and was only moderately excited in eccentricity; it can simply be the signature of Neptune's radial motion, starting from a moderately eccentric orbit. We discuss how this agrees with a scenario of giant planet evolution following a dynamical instability and, possibly, with the radial transport of the cold population., Comment: in press in Icarus

Evidence in the Solar System suggests that the giant planets underwent an epoch of radial migration that was very rapid, with an e-folding timescale shorter than 1~Myr. It is probable that the cause of this migration was that the giant planets experienced an orbital instability that caused them to encounter each other, resulting in radial migration. Several works suggest that this dynamical instability occurred `late', long after all the planets had formed and the solar nebula had dissipated. Assuming that the terrestrial planets had already formed, then their orbits would have been affected by the migration of the giant planets. As a result, how did the orbits of the terrestrial planets change? And can we use this migration to obtain information on the primordial orbits of the terrestrial planets? We directly model a large number of terrestrial planet systems and their response to giant planet migration. We study the change in the Angular Momentum Deficit (AMD) of the terrstrials. We conclude that the primordial AMD should have been lower than ~70\% of the current value, but higher than 10\%. We find that a scenario with five giant planets better satisfies the orbital constraints of the terrestrial planets. We predict that Mars was initially on an eccentric and inclined orbit while the orbits of Mercury, Venus and Earth were more circular and coplanar. The lower primordial dynamical excitement and the peculiar partitioning between planets impose new constraints for terrestrial planet formation simulations., Comment: Accepted for publication in MNRAS

Previous works have identified families halos by an analysis in proper elements domains, or by using Sloan Digital Sky Survey-Moving Object Catalog data, fourth release (SDSS-MOC4) multi-band photometry to infer the asteroid taxonomy, or by a combination of the two methods. The limited number of asteroids for which geometric albedo was known until recently discouraged in the past the extensive use of this additional parameter, which is however of great importance in identifying an asteroid taxonomy. The new availability of geometric albedo data from the Wide-field Infrared Survey Explorer (WISE) mission for about 100,000 asteroids significantly increased the sample of objects for which such information, with some errors, is now known. In this work we proposed a new method to identify families halos in a multi-domain space composed by proper elements, SDSS-MOC4 (a*,i-z) colors, and WISE geometric albedo for the whole main belt (and the Hungaria and Cybele orbital regions). Assuming that most families were created by the breakup of an undifferentiated parent body, they are expected to be homogeneous in colors and albedo. The new method is quite effective in determining objects belonging to a family halo, with low percentages of likely interlopers, and results that are quite consistent in term of taxonomy and geometric albedo of the halo members., Comment: 23 pages, 18 figures, 6 tables. Accepted for publication in MNRAS

In the Nice model, the late heavy bombardment (LHB) is related to an orbital instability of giant planets which causes a fast dynamical dispersion of a transneptunian cometary disk. We study effects produced by these hypothetical cometary projectiles on main-belt asteroids. In particular, we want to check whether the observed collisional families provide a lower or an upper limit for the cometary flux during the LHB. We present an updated list of observed asteroid families as identified in the space of synthetic proper elements by the hierarchical clustering method, colour data, albedo data and dynamical considerations and we estimate their physical parameters. We selected 12 families which may be related to the LHB according to their dynamical ages. We then used collisional models and N-body orbital simulations to gain insight into the long-term dynamical evolution of synthetic LHB families over 4 Gyr. We account for the mutual collisions, the physical disruptions of comets, the Yarkovsky/YORP drift, chaotic diffusion, or possible perturbations by the giant-planet migration. Assuming a "standard" size-frequency distribution of primordial comets, we predict the number of families with parent-body sizes D_PB >= 200 km which seems consistent with observations. However, more than 100 asteroid families with D_PB >= 100 km should be created at the same time which are not observed. This discrepancy can be nevertheless explained by the following processes: i) asteroid families are efficiently destroyed by comminution (via collisional cascade), ii) disruptions of comets below some critical perihelion distance (q <~ 1.5 AU) are common. Given the freedom in the cometary-disruption law, we cannot provide stringent limits on the cometary flux, but we can conclude that the observed distribution of asteroid families does not contradict with a cometary LHB., Comment: accepted in Astronomy and Astrophysics

We analyze the distribution of extrasolar planets (both confirmed and Kepler candidates) according to their orbital periods P and planetary radii R. Among confirmed planets, we find compelling evidence for a paucity of bodies with 3 < R < 10 R_\oplus, where R_\oplus in the Earth's radius, and P < 2-3 days. We have christened this region a "sub-Jovian Pampas". The same trend is detected in multiplanet Kepler candidates. Although approximately 16 Kepler single-planet candidates inhabit this Pampas, at least 7 are probable false positives (FP). This last number could be significantly higher if the ratio of FP is higher than 10%, as suggested by recent studies. In a second part of the paper we analyze the distribution of planets in the (P,R) plane according to stellar metallicities. We find two interesting trends: (i) a lack of small planets (R < 4 R_\oplus) with orbital periods P < 5 days in metal-poor stars, and (ii) a paucity of sub-Jovian planets (4 R_\oplus < R < 8 R_\oplus) with P < 100 days, also around metal-poor stars. Although all these trends are preliminary, they appear statistically significant and deserve further scrutiny. If confirmed, they could represent important constraints on theories of planetary formation and dynamical evolution., Comment: Accepted in ApJ

This paper explores the possibility that the progenitors of the small satellites of Pluto got captured in the Pluto-Charon system from the massive heliocentric planetesimal disk in which Pluto was originally embedded into. We find that, if the dynamical excitation of the disk is small, temporary capture in the Pluto-Charon system can occur with non-negligible probability, due to the dynamical perturbations exerted by the binary nature of the Pluto-Charon pair. However, the captured objects remain on very elliptic orbits and the typical capture time is only 100 years. In order to explain the origin of the small satellites of Pluto, we conjecture that some of these objects got disrupted during their Pluto-bound phase by a collision with a planetesimal of the disk. This could have generated a debris disk, which damped under internal collisional evolution, until turning itself into an accretional disk that could form small satellites on circular orbits, co-planar with Charon. Unfortunately, we find that objects large enough to carry a sufficient amount of mass to generate the small satellites of Pluto have collisional lifetimes orders of magnitude longer than the capture time. Thus, this scenario cannot explain the origin of the small satellites of Pluto, which remains elusive.

Aims: We present 11 high-precision photometric transit observations of the transiting super-Earth planet GJ1214b. Combining these data with observations from other authors, we investigate the ephemeris for possible signs of transit timing variations (TTVs) using a Bayesian approach. Methods: The observations were obtained using telescope-defocusing techniques, and achieve a high precision with random errors in the photometry as low as 1mmag per point. To investigate the possibility of TTVs in the light curve, we calculate the overall probability of a TTV signal using Bayesian methods. Results: The observations are used to determine the photometric parameters and the physical properties of the GJ1214 system. Our results are in good agreement with published values. Individual times of mid-transit are measured with uncertainties as low as 10s, allowing us to reduce the uncertainty in the orbital period by a factor of two. Conclusions: A Bayesian analysis reveals that it is highly improbable that the observed transit times is explained by TTV, when compared with the simpler alternative of a linear ephemeris., Comment: Submitted to A&A

Exoplanets show a pile-up of Jupiter-size planets in orbits with a 3-day period. A fraction of these hot Jupiters have retrograde orbits with respect to the parent star's rotation. To explain these observations we performed a series of numerical integrations of planet scattering followed by the tidal circularization. We considered planetary systems having 3 and 4 planets initially. We found that the standard Kozai migration is an inefficient mechanism for the formation of hot Jupiters. Our results show the formation of two distinct populations of hot Jupiters. The inner population of hot Jupiters with semimajor axis a < 0.03 AU formed in the systems where no planetary ejections occurred. This group contained a significant fraction of highly inclined and retrograde orbits, with distributions largely independent of the initial setup. However, our follow-up integrations showed that this populations was transient with most planets falling inside the Roche radius of the star in <1 Gyr. The outer population of hot Jupiters formed in systems where at least one planet was ejected. This population survived the effects of tides over >1 Gyr. The semimajor axis distribution of Population II fits nicely the observed 3-day pile-up. The inclination distribution of the outer hot planets depends on the number of planets in the initial systems and the 4-planet case showed a larger proportion (up to 10%), and a wider spread in inclination values. As the later results roughly agrees with observations, this may suggest that the planetary systems with observed hot Jupiters were originally rich in the number of planets, some of which were ejected. In a broad perspective, our work therefore hints on an unexpected link between the hot Jupiters and recently discovered free floating planets., Comment: submitted to ApJ

We model the long-term evolution of the Hilda collisional family located in the 3/2 mean-motion resonance with Jupiter. Its eccentricity distribution evolves mostly due to the Yarkovsky/YORP effect and assuming that: (i) impact disruption was isotropic, and (ii) albedo distribution of small asteroids is the same as for large ones, we can estimate the age of the Hilda family to be $4_{-1}^{+0}\,{\rm Gyr}$. We also calculate collisional activity in the J3/2 region. Our results indicate that current collisional rates are very low for a 200\,km parent body such that the number of expected events over Gyrs is much smaller than one. The large age and the low probability of the collisional disruption lead us to the conclusion that the Hilda family might have been created during the Late Heavy Bombardment when the collisions were much more frequent. The Hilda family may thus serve as a test of orbital behavior of planets during the LHB. We tested the influence of the giant-planet migration on the distribution of the family members. The scenarios that are consistent with the observed Hilda family are those with fast migration time scales $\simeq 0.3\,{\rm Myr}$ to $3\,{\rm Myr}$, because longer time scales produce a family that is depleted and too much spread in eccentricity. Moreover, there is an indication that Jupiter and Saturn were no longer in a compact configuration (with period ratio $P_{\rm S}/P_{\rm J} > 2.09$) at the time when the Hilda family was created.

This document covers five broad topics from theory of small bodies: planetesimal formation, cosmochemistry, thermal evolution, collisions and dynamics. Each of these topics is described in a separate section, where we prioritize selected main issues over completeness. The text points out the principal unresolved problems in each area, and suggests ways how progress can be made. This includes support for new code development, observations and experimental work that can be used to constrain theory, new research directions, and studies of cross-over regimes where the system's behavior is determined by several competing processes. The suggested development areas are placed in the context of NASA space exploration., Comment: http://www.lpi.usra.edu/sbag/

We report a transit timing study of the transiting exoplanetary system HD 189733. In total we observed ten transits in 2006 and 2008 with the 2.6-m Nordic Optical Telescope, and two transits in 2007 with the 4.2-m William Herschel Telescope. We used Markov-Chain Monte Carlo simulations to derive the system parameters and their uncertainties, and our results are in a good agreement with previously published values. We performed two independent analyses of transit timing residuals to place upper mass limits on putative perturbing planets. The results show no evidence for the presence of planets down to 1 Earth mass near the 1:2 and 2:1 resonance orbits, and planets down to 2.2 Earth masses near the 3:5 and 5:3 resonance orbits with HD 189733b. These are the strongest limits to date on the presence of other planets in this system., Comment: 10 pages, 4 figures, accepted by MNRAS

Using 8m-10m class telescopes and their Adaptive Optics (AO) systems, we conducted a long-term adaptive optics campaign initiated in 2003 focusing on four binary asteroid systems: (130) Elektra, (283) Emma, (379) Huenna, and (3749) Balam. The analysis of these data confirms the presence of their asteroidal satellite. We did not detect any additional satellite around these systems even though we have the capability of detecting a loosely-bound fragment (located at 1/4 x RHill) ~40 times smaller in diameter than the primary. The orbits derived for their satellites display significant eccentricity, ranging from 0.1 to 0.9, suggesting a different origin. Based on AO size estimate, we show that (130) Elektra and (283) Emma, G-type and P-type asteroids respectively, have a significant porosity (30-60% considering CI-CO meteorites as analogs) and their satellite's eccentricities (e~0.1) are possibly due to excitation by tidal effects. (379) Huenna and (3749) Balam, two loosely bound binary systems, are most likely formed by mutual capture. (3749) Balam's possible high bulk density is similar to (433) Eros, another S-type asteroid, and should be poorly fractured as well. (379) Huenna seems to display both characteristics: the moonlet orbits far away from the primary in term of stability (20% x RHill), but the primary's porosity is significant (30-60%)., Comment: 60 pages, 7 tables, 4 figures, in press to Icarus

The recent discovery of the first V-type asteroid in the middle belt, (21238) 1995WV7, located at ~2.54 AU, raises the question of whether it came from (4) Vesta or not. In this paper, we present spectroscopic observations indicating the existence of another V-type asteroid at ~2.53 AU, (40521) 1999RL95, and we investigate the possibility that these two asteroids evolved from the Vesta family to their present orbits by drifting in semi-major axis due to the Yarkovsky effect. The main problem with this scenario is that the asteroids need to cross the 3/1 mean motion resonance with Jupiter, which is highly unstable. Combining numerical simulations of the orbital evolution, that include the Yarkovsky effect, with Monte Carlo models, we compute the probability of an asteroid of given diameter D to evolve from the Vesta family and to cross over the 3/1 resonance, reaching a stable orbit in the middle belt. Our results indicate that an asteroid like (21238) 1995WV7 has a low probability of having evolved through this mechanism due to its large size (~5 km). However, the mechanism might explain the orbit of smaller bodies like (40521) 1999RL95 (~3 km), provided that we assume that the Vesta family formed > 3.5 Gy ago. We estimate that about 10% or more of the V-type bodies with D>1 km may come from the Vesta family by crossing over the 3/1 resonance. The remaining 90% must have a different origin., Comment: Submitted to Icarus - 4 tables, 8 figures

Among the many new irregular satellites that have been discovered in the last five years, at least six are in the so-called Kozai resonance. Due to solar perturbations, the argument of pericenter of a satellite usually precesses from 0 to 360 degrees. However, at inclinations higher than 39.3 degrees and lower than 140.7 degrees a new kind of behavior occurs for which the argument of pericenter oscillates around +/-90 degrees. In this work we will concentrate on the orbital history of the saturnian satellite S/2000 S5 Kiviuq, one of the satellites currently known to be in such resonance Kiviuq's orbit is very close to the separatrix of the Kozai resonance. Due to perturbations from the other jovian planets, it is expected that orbits near the Kozai separatrix may show significant chaotic behavior. This is important because chaotic diffusion may transfer orbits from libration to circulation, and vice versa. To identify chaotic orbits we used two well-known methods: the Frequency Analysis Method (Laskar 1990) and Maximum Lyapunov Exponents (Benettin et al. 1980). Our results show that the Kozai resonance is crossed by a web of secondary resonances, whose arguments involve combinations of the argument of pericenter, the argument of the Great Inequality, longitude of the node, and other terms related to the secular frequencies g5, g6, and s6. Many test orbits whose precession period is close to the period of the Great Inequality (883 yrs), or some of its harmonics, are trapped by these secondary resonances, and show significant chaotic behavior. Planetary migration, by moving the locations of these secondary resonances, may have depleted an original population of Kozai resonators., Comment: 47 pages, 18 figures, submitted to the Astronomical Journal

We present a path forward on a long-standing issue concerning the flux of small and slow meteoroids, which are believed to be the dominant portion of the incoming meteoric mass flux into the Earth's atmosphere. Such a flux, which is predicted by dynamical dust models of the Zodiacal Cloud, is not evident in ground-based radar observations. For decades this was attributed to the fact that the radars used for meteor observations lack the sensitivity to detect this population, due to the small amount of ionization produced by slow-velocity meteors. Such a hypothesis has been challenged by the introduction of meteor head echo (HE) observations with High Power and Large Aperture radars, in particular the Arecibo 430 MHz radar. Janches et al. developed a probabilistic approach to estimate the detectability of meteors by these radars and initially showed that, with the current knowledge of ablation and ionization, such particles should dominate the detected rates by one to two orders of magnitude compared to the actual observations. In this paper, we include results in our model from recently published laboratory measurements, which showed that (1) the ablation of Na is less intense covering a wider altitude range; and (2) the ionization probability, Beta ip, for Na atoms in the air is up to two orders of magnitude smaller for low speeds than originally believed. By applying these results and using a somewhat smaller size of the HE radar target we offer a solution that reconciles these observations with model predictions.