Following its flyby and first imaging the Pluto-Charon binary, the New Horizons spacecraft visited the Kuiper-Belt-Object (KBO) (486958) 2014 MU69 (Arrokoth). Imaging showed MU69 to be a contact-binary, made of two individual lobes connected by a narrow neck, rotating at low spin period (15.92 h), and having high obliquity (~98 deg), similar to other KBO contact-binaries inferred through photometric observations. The origin of such peculiar configurations is puzzling, and all scenarios suggested for the origins of contact-binaries fail to reproduce such properties and their likely high frequency. Here we show that semi-secular perturbations operating only on ultra-wide (~0.1-0.4 Hill-radius) KBO-binaries can robustly lead to gentle, slow-speed binary mergers at arbitrarily high obliquities, but low rotational velocities, that can reproduce MU69's (and similar oblique contact binaries) characteristics. Using N-body simulations, we find that ~15% of all ultra-wide binaries with cosine-uniform inclination distribution are likely to merge through this process. Moreover, we find that such mergers are sufficiently gentle as to only slightly deform the KBO shape, and can produce the measured rotation speed of MU69. The semi-secular contact-binary formation channel not only explains the observed properties of MU69, but could also apply for other Kuiper/asteroid belt binaries, and for Solar/extra-solar moon systems.

Exomoons orbiting terrestrial or super-terrestrial exoplanets have not yet been discovered; their possible existence and properties are therefore still an unresolved question. Here we explore the collisional formation of exomoons through giant planetary impacts. We make use of smooth particle hydrodynamical (SPH) collision simulations and survey a large phase-space of terrestrial/super-terrestrial planetary collisions. We characterize the properties of such collisions, finding one rare case in which an exomoon forms through a graze&capture scenario, in addition to a few graze&merge or hit&run scenarios. Typically however, our collisions form massive circumplanetary discs, for which we use follow-up N-body simulations in order to derive lower-limit mass estimates for the ensuing exomoons. We investigate the mass, long-term tidal-stability, composition and origin of material in both the discs and the exomoons. Our giant-impact models often generate relatively iron-rich moons, that form beyond the synchronous radius of the planet, and would thus tidally evolve outward with stable orbits, rather than be destroyed. Our results suggest that it is extremely difficult to collisionally form currently-detectable exomoons orbiting super-terrestrial planets, through single giant impacts. It might be possible to form massive, detectable exomoons through several mergers of smaller exomoons, formed by multiple impacts, however more studies are required in order to reach a conclusion. Given the current observational initiatives, the search should focus primarily on more massive planet categories. However, about a quarter of the exomoons predicted by our models are approximately Mercury-mass or more, and are much more likely to be detectable given a factor 2 improvement in the detection capability of future instruments, providing further motivation for their development.

We make use of a new hybrid method to simulate the long-term, multiple-orbit disc formation through tidal disruptions of rocky bodies by white dwarfs, at high-resolution and realistic semi-major axis. We perform the largest-yet suite of simulations for dwarf and terrestrial planets, spanning four orders of magnitude in mass, various pericentre distances and semi-major axes between 3 AU and 150 AU. This large phase space of tidal disruption conditions has not been accessible through the use of previous codes. We analyse the statistical and structural properties of the emerging debris discs, as well as the ejected unbound debris contributing to the population of interstellar asteroids. Unlike previous tidal disruption studies of small asteroids which form ring-like structures on the original orbit, we find that the tidal disruption of larger bodies usually forms dispersed structures of interlaced elliptic eccentric annuli on tighter orbits. We characterize the (typically power-law) size-distribution of the ejected interstellar bodies as well as their composition, rotation velocities and ejection velocities. We find them to be sensitive to the depth (impact parameter) of the tidal disruption. Finally, we briefly discuss possible implications of our results in explaining the peculiar variability of Tabby's star, the origin of the transit events of ZTF J0139+5245 and the formation of a planetary core around SDSS J1228+1040.

We introduce a new hybrid method to perform high-resolution tidal disruption simulations, at arbitrary orbits. An SPH code is used to simulate tidal disruptions only in the immediate spatial domain of the star, namely, where the tidal forces dominate over gravity, and then during the fragmentation phase in which the emerging tidal stream may collapse under its own gravity to form fragments. Following each hydrodynamical simulation, an analytical treatment is then applied to instantaneously transfer each fragment back to the tidal sphere for its subsequent disruption, in an iterative process. We validate the hybrid model by comparing it to both an analytical impulse approximation model of single tidal disruptions, as well as full-scale SPH simulations spanning the entire disc formation. The hybrid simulations are essentially indistinguishable from the full-scale SPH simulations, while computationally outperforming their counterparts by orders of magnitude. Thereby our new hybrid approach uniquely enables us to follow the long-term formation and continuous tidal disruption of the planet/planetesimal debris, without the resolution and orbital configuration limitation of previous studies. In addition, we describe a variety of future directions and applications for our hybrid model, which is in principle applicable to any star, not merely white dwarfs

During the last stages of the terrestrial planet formation, planets grow mainly through giant-impacts with large planetary embryos. The Earth's Moon was suggested to form through one of these impacts. However, since the proto-Earth has experienced many giant-impacts, several moons are naturally expected to form through a sequence of multiple (including smaller scale) impacts. Each impact potentially forms a sub-Lunar mass moonlet that interacts gravitationally with the proto-Earth and possibly with previously-formed moonlets. Such interactions result in either moonlet-moonlet mergers, moonlet ejections or infall of moonlets on the Earth. The latter possibility, leading to low-velocity moonlet-Earth collisions is explored here for the first time. We make use of SPH simulations and consider a range of moonlet masses, collision impact-angles and initial proto-Earth rotation rates. We find that grazing/tidal-collisions are the most frequent and produce comparable fractions of accreted-material and debris. The latter typically clump in smaller moonlets that can potentially later interact with other moonlets. Other collision geometries are more rare. Head-on collisions do not produce much debris and are effectively perfect mergers. Intermediate impact angles result in debris mass-fractions in the range of 2-25% where most of the material is unbound. Retrograde collisions produce more debris than prograde collisions, whose fractions depend on the proto-Earth initial rotation rate. Moonfalls can slightly change the rotation-rate of the proto-Earth. Accreted moonfall material is highly localized, potentially explaining the isotopic heterogeneities in highly siderophile elements in terrestrial rocks, and possibly forming primordial super-continent topographic features. Our results can be used for simple scaling laws and applied to n-body studies of the formation of the Earth and Moon.

Studies suggest that the pollution of white dwarf (WD) atmospheres arises from the accretion of minor planets, but the exact properties of polluting material, and in particular the evidence for water in some cases are not yet understood. Previous works studied the water retention in minor planets around main-sequence and evolving host stars, in order to evaluate the possibility that water survives inside minor planets around WDs. However, all of these studies focused on small, comet-sized to moonlet-sized minor planets, when the inferred mass inside the convection zones of He-dominated WDs could actually also be compatible with much more massive minor planets. In this study we therefore explore for the first time, the water retention inside exo-planetary dwarf planets, or moderate-sized moons, with radii of the order of hundreds of kilometres. We now cover nearly the entire potential mass range of minor planets. The rest of the parameter space considered in this study is identical to that of our previous study, and also includes multiple WD progenitor star masses. We find that water retention in more massive minor planets is still affected by the mass of the WD progenitor, however not as much as when small minor planets were considered. We also find that water retention is now almost always greater than zero. On average, the detected water fraction in He-dominated WD atmospheres should be at least 5\%, irrespective of the assumed initial water composition, if it came from a single accretion event of an icy dwarf planet or moon. This finding also strengthens the possibility of WD habitability. To finalize our previous and current findings, we provide a code which may be freely used as a service to the community. The code calculates ice and water retention by interpolation, spanning the full mass range of both minor planets and their host stars.

Most studies suggest the pollution of white dwarf (WD) atmospheres arises from accretion of minor planets, but the exact properties of polluting material, and in particular the evidence for water in some cases are not yet understood. Here we study the water retention of small icy bodies in exo-solar planetary systems, as their respective host stars evolve through and off the main sequence and eventually become WDs. We explore for the first time a wide range of star masses and metallicities. We find that the mass of the WD progenitor star is of crucial importance for the retention of water, while its metallicity is relatively unimportant. We predict that minor planets around lower-mass WD progenitors would retain more water in general, and would do so at closer distances from the WD, compared with high-mass progenitors. The dependence of water retention on progenitor mass and other parameters has direct implications for the origin of observed WD pollution, and we discuss how our results and predictions might be tested in the future as more observations of WDs with long cooling-ages become available.

The New Horizons mission to the Kuiper Belt has recently revealed intriguing features on the surface of Charon, including a network of chasmata, cutting across or around a series of high topography features, conjoining to form a belt. It is proposed that this tectonic belt is a consequence of contraction/expansion episodes in the moon's evolution associated particularly with compaction, differentiation and geochemical reactions of the interior. The proposed scenario involves no need for solidification of a vast subsurface ocean and/or a warm initial state. This scenario is based on a new, detailed thermo-physical evolution model of Charon that includes multiple processes. According to the model, Charon experiences two contraction/expansion episodes in its history that may provide the proper environment for the formation of the tectonic belt. This outcome remains qualitatively the same, for several different initial conditions and parameter variations. The precise orientation of Charon's tectonic belt, and the cryovolcanic features observed south of the tectonic belt may have involved a planetary-scale impact, that occurred only after the belt had already formed.

Most observations of polluted white dwarf atmospheres are consistent with accretion of water depleted planetary material. Among tens of known cases, merely two cases involve accretion of objects that contain a considerable mass fraction of water. The purpose of this study is to investigate the relative scarcity of these detections. Based on a new and highly detailed model, we evaluate the retention of water inside icy minor planets during the high luminosity stellar evolution that follows the main sequence. Our model fully considers the thermal, physical, and chemical evolution of icy bodies, following their internal differentiation as well as water depletion, from the moment of their birth and through all stellar evolution phases preceding the formation of the white dwarf. We also account for different initial compositions and formation times. Our results show that previous studies have either underestimated or overestimated water retention. We also reaffirm that water can survive in a variety of circumstances and in great quantities, and therefore other possibilities are discussed in order to explain the infrequency of water detections. We find that the most likely explanation is that water does not survive the dissociation of tidally disrupted minor planets, and instead a direct impact onto the white dwarf is required. We show that this possibility is statistically plausible, and predict that if this hypothesis is correct, future observations of extremely water rich atmospheres will involve only helium dominated DB white dwarfs.

We develop a long-term 1-D evolution model for icy satellites that couples multiple processes: water migration and differentiation, geochemical reactions and silicate phase transitions, compaction by self-gravity, and ablation. The model further considers the following energy sources and sinks: tidal heating, radiogenic heating, geochemical energy released by serpentinization or absorbed by mineral dehydration, gravitational energy and insolation, and heat transport by conduction, convection, and advection. We apply the model to Enceladus, by guessing the initial conditions that would render a structure compatible with present-day observations, assuming the initial structure to have been homogeneous. Assuming the satellite has been losing water continually along its evolution, we postulate that it was formed as a more massive, more icy and more porous satellite, and gradually transformed into its present day state due to sustained long-term tidal heating. We consider several initial compositions and evolution scenarios and follow the evolution for the age of the Solar System, testing the present day model results against the available observational constraints. Our model shows the present configuration to be differentiated into a pure icy mantle, several tens of km thick, overlying a rocky core, composed of dehydrated rock at the center and hydrated rock in the outer part. For Enceladus, it predicts a higher rock/ice mass ratio than previously assumed and a thinner ice mantle, compatible with recent estimates based on gravity field measurements. Although, obviously, the model cannot be used to explain local phenomena, it sheds light on the internal structure invoked in explanations of localized features and activities.

Although many models have been proposed, the physical mechanisms responsible for the formation of low-mass brown dwarfs are poorly understood. The multiplicity properties and minimum mass of the brown-dwarf mass function provide critical empirical diagnostics of these mechanisms. We present the discovery via gravitational microlensing of two very low-mass, very tight binary systems. These binaries have directly and precisely measured total system masses of 0.025 Msun and 0.034 Msun, and projected separations of 0.31 AU and 0.19 AU, making them the lowest-mass and tightest field brown-dwarf binaries known. The discovery of a population of such binaries indicates that brown dwarf binaries can robustly form at least down to masses of ~0.02 Msun. Future microlensing surveys will measure a mass-selected sample of brown-dwarf binary systems, which can then be directly compared to similar samples of stellar binaries.

Brown dwarfs are important objects because they may provide a missing link between stars and planets, two populations that have dramatically different formation history. In this paper, we present the candidate binaries with brown dwarf companions that are found by analyzing binary microlensing events discovered during 2004 - 2011 observation seasons. Based on the low mass ratio criterion of q < 0.2, we found 7 candidate events, including OGLE-2004-BLG-035, OGLE-2004-BLG-039, OGLE-2007-BLG-006, OGLE-2007-BLG-399/MOA-2007-BLG-334, MOA-2011-BLG-104/OGLE-2011-BLG-0172, MOA-2011-BLG-149, and MOA-201-BLG-278/OGLE-2011-BLG-012N. Among them, we are able to confirm that the companions of the lenses of MOA-2011-BLG-104/OGLE-2011-BLG-0172 and MOA-2011-BLG-149 are brown dwarfs by determining the mass of the lens based on the simultaneous measurement of the Einstein radius and the lens parallax. The measured mass of the brown dwarf companions are (0.02 +/- 0.01) M_Sun and (0.019 +/- 0.002) M_Sun for MOA-2011-BLG-104/OGLE-2011-BLG-0172 and MOA-2011-BLG-149, respectively, and both companions are orbiting low mass M dwarf host stars. More microlensing brown dwarfs are expected to be detected as the number of lensing events with well covered light curves increases with new generation searches.