Showing 1–23 of 23 results for author: Malamud, U

A link between the size and composition of comets

Authors: James E. Robinson, Uri Malamud, Cyrielle Opitom, Hagai Perets, Jürgen Blum

Abstract: All cometary nuclei that formed in the early Solar System incorporated radionuclides and therefore were subject to internal radiogenic heating. Previous work predicts that if comets have a pebble-pile structure internal temperature build-up is enhanced due to very low thermal conductivity, leading to internal differentiation. An internal thermal gradient causes widespread sublimation and migration… ▽ More All cometary nuclei that formed in the early Solar System incorporated radionuclides and therefore were subject to internal radiogenic heating. Previous work predicts that if comets have a pebble-pile structure internal temperature build-up is enhanced due to very low thermal conductivity, leading to internal differentiation. An internal thermal gradient causes widespread sublimation and migration of either ice condensates, or gases released from amorphous ice hosts during their crystallisation. Overall, the models predict that the degree of differentiation and re-distribution of volatile species to a shallower near-surface layer depends primarily on nucleus size. Hence, we hypothesise that cometary activity should reveal a correlation between the abundance of volatile species and the size of the nucleus. To explore this hypothesis we have conducted a thorough literature search for measurements of the composition and size of cometary nuclei, compiling these into a unified database. We report a statistically significant correlation between the measured abundance of CO/H$_{2}$O and the size of cometary nuclei. We further recover the measured slope of abundance as a function of size, using a theoretical model based on our previous thermophysical models, invoking re-entrapment of outward migrating high volatility gases in the near-surface pristine amorphous ice layers. This model replicates the observed trend and supports the theory of internal differentiation of cometary nuclei by early radiogenic heating. We make our database available for future studies, and we advocate for collection of more measurements to allow more precise and statistically significant analyses to be conducted in the future. △ Less

Submitted 11 July, 2024; v1 submitted 22 March, 2024; originally announced March 2024.

Comments: 25 pages, 10 figures, published in MNRAS. Minor copyedit changes and updated links to data repositories

Uranus and Neptune as methane planets: producing icy giants from refractory planetesimals

Authors: Uri Malamud, Morris Podolak, Joshua Podolak, Peter Bodenheimer

Abstract: Uranus and Neptune are commonly considered ice giants, and it is often assumed that, in addition to a solar mix of hydrogen and helium, they contain roughly twice as much water as rock. This classical picture has led to successful models of their internal structure and has been understood to be compatible with the composition of the solar nebula during their formation (Reynolds and Summers 1965; P… ▽ More Uranus and Neptune are commonly considered ice giants, and it is often assumed that, in addition to a solar mix of hydrogen and helium, they contain roughly twice as much water as rock. This classical picture has led to successful models of their internal structure and has been understood to be compatible with the composition of the solar nebula during their formation (Reynolds and Summers 1965; Podolak and Cameron 1974; Podolak and Reynolds 1984; Podolak et al. 1995; Nettelmann et al. 2013). However, the dominance of water has been recently questioned (Teanby et al. 2020; Helled and Fortney 2020; Podolak et al. 2022). Planetesimals in the outer solar system are composed mainly of refractory materials, leading to an inconsistency between the icy composition of Uranus and Neptune and the ice-poor planetesimals they accreted during formation (Podolak et al. 2022). Here we elaborate on this problem, and propose a new potential solution. We show that chemical reactions between planetesimals dominated by organic-rich refractory materials and the hydrogen in gaseous atmospheres of protoplanets can form large amounts of methane 'ice'. Uranus and Neptune could thus be compatible with having accreted refractory-dominated planetesimals, while still remaining icy. Using random statistical computer models for a wide parameter space, we show that the resulting methane-rich internal composition could be a natural solution, giving a good match to the size, mass and moment of inertia of Uranus and Neptune, whereas rock-rich models appear to only work if a rocky interior is heavily mixed with hydrogen. Our model predicts a lower than solar hydrogen to helium ratio, which can be tested. We conclude that Uranus, Neptune and similar exoplanets could be methane-rich, and discuss why Jupiter and Saturn cannot. △ Less

Submitted 10 July, 2024; v1 submitted 19 March, 2024; originally announced March 2024.

Comments: Accepted for publication in Icarus

White dwarf systems: exoplanets and debris disks

Authors: Uri Malamud

Abstract: Although there is abundant and diverse observational evidence in support of white dwarf stars hosting planets or debris disks which form in the catastrophic destruction of various planetary bodies, the key processes that explain these observations are still being intensely investigated. The study of white dwarf planetary systems offers a unique perspective on exo-solar composition, that cannot be… ▽ More Although there is abundant and diverse observational evidence in support of white dwarf stars hosting planets or debris disks which form in the catastrophic destruction of various planetary bodies, the key processes that explain these observations are still being intensely investigated. The study of white dwarf planetary systems offers a unique perspective on exo-solar composition, that cannot be obtained by any other means. This chapter describes the various observational techniques that are used in order to find and characterize exo-planets and debris disks around white dwarfs. In turn, it discusses how to theoretically interpret these observations by surveying an array of various research tools and models currently employed in this field. △ Less

Submitted 21 April, 2024; v1 submitted 12 March, 2024; originally announced March 2024.

Comments: Preprint of a chapter for the 'Encyclopedia of Astrophysics' (Editor-in-Chief Ilya Mandel, Section Editor Dimitri Veras) to be published by Elsevier as a Reference Module

An Equation of State of CO for use in Planetary Modeling

Authors: Morris Podolak, Amit Levi, Allona Vazan, Uri Malamud

Abstract: Although carbon monoxide (CO) is an abundant molecule and may have great importance for planetary interiors, measurements of its properties are difficult due to its extreme volatility. We calculate the equation of state for CO over a range of temperature and density that is applicable to the conditions in planetary interiors. Previous experimental and theoretical studies cover only a limited tempe… ▽ More Although carbon monoxide (CO) is an abundant molecule and may have great importance for planetary interiors, measurements of its properties are difficult due to its extreme volatility. We calculate the equation of state for CO over a range of temperature and density that is applicable to the conditions in planetary interiors. Previous experimental and theoretical studies cover only a limited temperature-density range. Our calculations match these early results well, but now cover the full range of relevance. The method of calculation is based on the general-purpose quotidian equation of state described by More et al. (1988), which is here used in order to generate a freely downloadable look-up table to be used by the community. △ Less

Submitted 5 January, 2023; originally announced January 2023.

Comments: 12 pages, 5 figures. Accepted for publication in Icarus

Asynchronous accretion can mimic diverse white dwarf pollutants I: core and mantle fragments

Authors: Marc G. Brouwers, Amy Bonsor, Uri Malamud

Abstract: Polluted white dwarfs serve as astrophysical mass spectrometers - their photospheric abundances are used to infer the composition of planetary objects that accrete onto them. We show that due to asymmetries in the accretion process, the composition of the material falling onto a star may vary with time during the accretion of a single planetary body. Consequently, the instantaneous photospheric ab… ▽ More Polluted white dwarfs serve as astrophysical mass spectrometers - their photospheric abundances are used to infer the composition of planetary objects that accrete onto them. We show that due to asymmetries in the accretion process, the composition of the material falling onto a star may vary with time during the accretion of a single planetary body. Consequently, the instantaneous photospheric abundances of white dwarfs do not necessarily reflect the bulk composition of their pollutants, especially when their diffusion timescales are short. In particular, we predict that when an asteroid with an iron core tidally disrupts around a white dwarf, a larger share of its mantle is ejected, and that the core/mantle fraction of the accreting material varies with time during the event. Crucially, this implies that the core fraction of differentiated pollutants cannot be determined for white dwarfs with short diffusion timescales, which sample only brief episodes of longer accretion processes. The observed population of polluted white dwarfs backs up the proposed theory. More white dwarfs have accreted material with high Fe/Ca than low Fe/Ca relative to stellar abundance ratios, indicating the ejection of mantle material. Additionally, we find tentative evidence that the accretion rate of iron decreases more rapidly than that of magnesium or calcium, hinting at variability of the accreted composition. Further corroboration of the proposed theory will come from the up-coming analysis of large samples of young white dwarfs. △ Less

Submitted 9 November, 2022; originally announced November 2022.

Comments: Accepted for publication in MNRAS. Part one of a series of two papers. Comments and questions welcome

Asynchronous accretion can mimic diverse white dwarf pollutants II: water content

Authors: Marc G. Brouwers, Andrew M. Buchan, Amy Bonsor, Uri Malamud, Elliot Lynch Laura Rogers, Detlev Koester

Abstract: Volatiles, notably water, are key to the habitability of rocky planets. The presence of water in planetary material can be inferred from the atmospheric oxygen abundances of polluted white dwarfs, but this interpretation is often complex. We study the accretion process, and find that ices may sublimate and accrete before more refractory minerals reach the star. As a result, a white dwarf's relativ… ▽ More Volatiles, notably water, are key to the habitability of rocky planets. The presence of water in planetary material can be inferred from the atmospheric oxygen abundances of polluted white dwarfs, but this interpretation is often complex. We study the accretion process, and find that ices may sublimate and accrete before more refractory minerals reach the star. As a result, a white dwarf's relative photospheric abundances may vary with time during a single accretion event, and do not necessarily reflect the bulk composition of a pollutant. We offer two testable predictions for this hypothesis: 1. cooler stars will more often be inferred to have accreted wet pollutants, and 2. there will be rare occurrences of accretion events with inferred volatile levels far exceeding those of pristine comets. To observationally test these predictions, we statistically constrain the water content of white dwarf pollutants. We find that in the current sample, only three stars show statistically significant evidence of water at the 2$σ$ level, due to large typical uncertainties in atmospheric abundances and accretion states. In the future, an expanded sample of polluted white dwarfs with hydrogen-dominated atmospheres will allow for the corroboration of our theoretical predictions. Our work also shows the importance of interpreting pollutant compositions statistically, and emphasizes the requirement to reduce uncertainties on measured abundances to allow for statistically significant constraints on their water content. △ Less

Submitted 9 November, 2022; originally announced November 2022.

Comments: Accepted for publication in MNRAS. Part two of a series of two papers. Comments and questions welcome

Are there any pristine comets? Constraints from pebble structure

Authors: Uri Malamud, Wolf A. Landeck, Dorothea Bischoff, Christopher Kreuzig, Hagai B. Perets, Bastian Gundlach, Jurgen Blum

Abstract: We show that if comets (or any small icy planetesimals such as Kuiper belt objects) are composed of pebble piles, their internal radiogenic as well as geochemical heating results in considerably different evolutionary outcomes compared to similar past studies. We utilize a 1D thermo-physical evolution code, modified to include state-of-the-art empirical measurements of pebble thermal conductivity… ▽ More We show that if comets (or any small icy planetesimals such as Kuiper belt objects) are composed of pebble piles, their internal radiogenic as well as geochemical heating results in considerably different evolutionary outcomes compared to similar past studies. We utilize a 1D thermo-physical evolution code, modified to include state-of-the-art empirical measurements of pebble thermal conductivity and compression, the latter obtained through a new laboratory experiment presented here for the first time. Results indicate that due to the low pebble thermal conductivity, the peak temperatures attained during evolution are much higher than in any previous study given the same formation time. Assuming meteoritic radiogenic abundances, we find that only extremely small, sub-kilometre comets have the potential to retain the primordial, uniform and thermally unprocessed composition from which they formed. Comets with radii in excess of about 20 km are typically swept by rapid and energetically powerful aqueous hydration reactions. Across the full range of comet sizes and formation times, evolutions result in the processing and differentiation of various volatile species, and a radially heterogeneous nucleus stucture. Our computations however also indicate that the assumed fraction of radionuclides is a pivotal free parameter, because isotopic analyses of the only available cometary samples suggest that no 26Al was ever present in comet 81P/Wild 2. We show that if comets formed early in the protoplanetary disc (within 1-3 Myr), the radionuclide abundances indeed must be much smaller than those typically assumed based on meteoritic samples. We discuss the importance of our findings for the formation, present-day attributes and future research of comets. △ Less

Submitted 31 May, 2022; originally announced June 2022.

Random Models for Exploring Planet Compositions I: Uranus as an Example

Authors: Joshua Podolak, Uri Malamud, Morris Podolak

Abstract: Modeling the interior of a planet is difficult because the small number of measured parameters is insufficient to constrain the many variables involved in describing the interior structure and composition. One solution is to invoke additional constraints based on arguments about how the planet formed. However, a planet's actual structure and composition may hold clues to its formation which would… ▽ More Modeling the interior of a planet is difficult because the small number of measured parameters is insufficient to constrain the many variables involved in describing the interior structure and composition. One solution is to invoke additional constraints based on arguments about how the planet formed. However, a planet's actual structure and composition may hold clues to its formation which would be lost if this structure were not allowed by the initial assumptions. It is therefore interesting to explore the space of allowable compositions and structures in order to better understand which cosmogonic constraints are absolutely necessary. To this end, we describe a code for generating random, monotonic, density distributions, rho(r), that fit a given mass, radius, and moment of inertia. Integrating the equation of hydrostatic equilibrium gives the pressure, P(r), at each point in the body. We then provide three algorithms for generating a monotonic temperature distribution, T(r), and an associated composition that is consistent with the rho-P relation, and realistic equations of state. We apply this code to Uranus as a proof of concept, and show that the ratio of rock to water cannot be much larger than 2. △ Less

Submitted 2 March, 2022; originally announced March 2022.

Comments: 32 pages, 15 figures

A road-map to white dwarf pollution: Tidal disruption, eccentric grind-down, and dust accretion

Authors: Marc G. Brouwers, Amy Bonsor, Uri Malamud

Abstract: A significant fraction of white dwarfs show metal lines indicative of pollution with planetary material but the accretion process remains poorly understood. The main aim of this paper is to produce a road-map illustrating several potential routes for white dwarf pollution and to link these paths to observational outcomes. Our proposed main road begins with the tidal disruption of a scattered aster… ▽ More A significant fraction of white dwarfs show metal lines indicative of pollution with planetary material but the accretion process remains poorly understood. The main aim of this paper is to produce a road-map illustrating several potential routes for white dwarf pollution and to link these paths to observational outcomes. Our proposed main road begins with the tidal disruption of a scattered asteroid and the formation of a highly eccentric tidal disc with a wide range of fragment sizes. Accretion of these fragments by Poynting-Robertson (PR) drag alone is too slow to explain the observed rates. Instead, in the second stage, several processes including differential apsidal precession cause high-velocity collisions between the eccentric fragments. Large asteroids produce more fragments when they disrupt, causing rapid grind-down and generating short and intense bursts of dust production, whereas smaller asteroids grind down over longer periods of time. In the final stage, the collisionally produced dust circularises and accretes onto the white dwarf via drag forces. We show that optically thin dust accretion by PR drag produces large infrared (IR) excesses when the accretion rate exceeds 10^7 g/s. We hypothesise that around white dwarfs accreting at a high rate, but with no detected infrared excess, dust circularisation requires enhanced drag - for instance due to the presence of gas near the disc's pericentre. △ Less

Submitted 11 November, 2021; v1 submitted 14 October, 2021; originally announced October 2021.

Comments: Published in MNRAS

The entry geometry and velocity of planetary debris into the Roche sphere of a white dwarf

Authors: Dimitri Veras, Nikolaos Georgakarakos, Alexander J. Mustill, Uri Malamud, Tim Cunningham, Ian Dobbs-Dixon

Abstract: Our knowledge of white dwarf planetary systems predominately arises from the region within a few Solar radii of the white dwarfs, where minor planets break up, form rings and discs, and accrete onto the star. The entry location, angle and speed into this Roche sphere has rarely been explored but crucially determines the initial geometry of the debris, accretion rates onto the photosphere, and ulti… ▽ More Our knowledge of white dwarf planetary systems predominately arises from the region within a few Solar radii of the white dwarfs, where minor planets break up, form rings and discs, and accrete onto the star. The entry location, angle and speed into this Roche sphere has rarely been explored but crucially determines the initial geometry of the debris, accretion rates onto the photosphere, and ultimately the composition of the minor planet. Here we evolve a total of over 10^5 asteroids with single-planet N-body simulations across the giant branch and white dwarf stellar evolution phases to quantify the geometry of asteroid injection into the white dwarf Roche sphere as a function of planetary mass and eccentricity. We find that lower planetary masses increase the extent of anisotropic injection and decrease the probability of head-on (normal to the Roche sphere) encounters. Our results suggest that one can use dynamical activity within the Roche sphere to make inferences about the hidden architectures of these planetary systems. △ Less

Submitted 8 June, 2021; originally announced June 2021.

Comments: Accepted for publication in MNRAS

Circularization of tidal debris around white dwarfs: implications for gas production and dust variability

Authors: Uri Malamud, Evgeni Grishin, Marc Brouwers

Abstract: White dwarf (WD) pollution is thought to arise from the tidal disruption of planetary bodies. The initial fragment stream is extremely eccentric, while observational evidence suggest that discs are circular or nearly so. Here we propose a novel mechanism to bridge this gap and show that the fragments can rapidly circularise through dust or gas drag when they interact with a pre-existing compact di… ▽ More White dwarf (WD) pollution is thought to arise from the tidal disruption of planetary bodies. The initial fragment stream is extremely eccentric, while observational evidence suggest that discs are circular or nearly so. Here we propose a novel mechanism to bridge this gap and show that the fragments can rapidly circularise through dust or gas drag when they interact with a pre-existing compact disc. We assume that the tidal stream mainly consists of small cohesive fragments in the size range 10-1000 m, capable of resisting the WD tidal forces, whereas the compact discs span a wide mass range. We provide an analytical model, accompanied by N-body simulations, and find a large parameter space in fragment sizes and orbital separation that leads to full circularization. Partial circularization is possible for compact discs that are several orders of magnitudes less massive. We show that dust-induced circularization inherently produces gas as tidal fragments collisionally vaporize the pre-existing dust along their path. We show that ongoing gas production has a higher probability to occur during the early stages of tidal disruption events, resulting from the fact that smaller fragments are the first to circularize. Intermittent gas production however becomes more likely as the tidal stream matures. This could explain why only a small subset of systems with dusty compact discs also have an observed gaseous component. Additionally, the interaction yields fragment erosion by collisional shattering, sputtering, sublimation and possibly ram-pressure. Material scattered by the collisions might form a thin dusty halo that evolves through PR drag, in compatibility with observed infrared variability. △ Less

Submitted 18 January, 2021; v1 submitted 14 December, 2020; originally announced December 2020.

Origin of (2014) MU69-like Kuiper-belt contact binaries from wide binaries

Authors: Evgeni Grishin, Uri Malamud, Hagai B. Perets, Oliver Wandel, Christoph M. Schaefer

Abstract: 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 photometr… ▽ More 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. △ Less

Submitted 3 March, 2020; originally announced March 2020.

Comments: Submitted on 27 June 2019

Journal ref: Nature, 580, 463 (2020)

Tidal disruption of planetary bodies by white dwarfs II: Debris disc structure and ejected interstellar asteroids

Authors: Uri Malamud, Hagai Perets

Abstract: 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.… ▽ More 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. △ Less

Submitted 15 January, 2020; v1 submitted 27 November, 2019; originally announced November 2019.

Tidal disruption of planetary bodies by white dwarfs I: A hybrid SPH-analytical approach

Authors: Uri Malamud, Hagai Perets

Abstract: 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. Follow… ▽ More 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 △ Less

Submitted 15 January, 2020; v1 submitted 27 November, 2019; originally announced November 2019.

Collisional formation of massive exomoons of super-terrestrial exoplanets

Authors: Uri Malamud, Hagai B. Perets, Christoph Schaefer, Christoph Burger

Abstract: 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… ▽ More 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. △ Less

Submitted 28 January, 2020; v1 submitted 29 April, 2019; originally announced April 2019.

Moonfalls: Collisions between the Earth and its past moons

Authors: Uri Malamud, Hagai B. Perets, Christoph Schafer, Christoph Burger

Abstract: 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 potential… ▽ More 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. △ Less

Submitted 30 April, 2018; originally announced May 2018.

Post-main-sequence evolution of icy minor planets. III. water retention in dwarf planets and exo-moons and implications for white dwarf pollution

Authors: Uri Malamud, Hagai B. Perets

Abstract: 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… ▽ More 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. △ Less

Submitted 24 August, 2017; originally announced August 2017.

Post main sequence evolution of icy minor planets II: water retention and white dwarf pollution around massive progenitor stars

Authors: Uri Malamud, Hagai B. Perets

Abstract: 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 even… ▽ More 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. △ Less

Submitted 14 May, 2017; v1 submitted 4 April, 2017; originally announced April 2017.

Post main sequence evolution of icy minor planets: Implications for water retention and white dwarf pollution

Authors: Uri Malamud, Hagai B. Perets

Abstract: 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 retentio… ▽ More 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. △ Less

Submitted 14 May, 2017; v1 submitted 1 August, 2016; originally announced August 2016.

The Contraction/Expansion History of Charon with implication for its Planetary Scale Tectonic Belt

Authors: Uri Malamud, Hagai B. Perets, Gerald Schubert

Abstract: 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, differenti… ▽ More 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. △ Less

Submitted 14 May, 2017; v1 submitted 2 March, 2016; originally announced March 2016.

Comments: Comments are welcome

A 1-D evolutionary model for icy satellites, applied to Enceladus

Authors: Uri Malamud, Dina Prialnik

Abstract: 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 minera… ▽ More 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. △ Less

Submitted 25 October, 2015; originally announced October 2015.

Comments: submitted to Icarus

Microlensing Discovery of a Population of Very Tight, Very Low-mass Binary Brown Dwarfs

Authors: J. -Y. Choi, C. Han, A. Udalski, T. Sumi, B. S. Gaudi, A. Gould, D. P. Bennett, M. Dominik, J. -P. Beaulieu, Y. Tsapras, V. Bozza, F. Abe, I. A. Bond, C. S. Botzler, P. Chote, M. Freeman, A. Fukui, K. Furusawa, Y. Itow, C. H. Ling, K. Masuda, Y. Matsubara, N. Miyake, Y. Muraki, K. Ohnishi , et al. (99 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 20 March, 2013; v1 submitted 18 February, 2013; originally announced February 2013.

Comments: 7 pages, 5 figures, 3 tables, ApJ submitted

Microlensig Binaries with Candidate Brown Dwarf Companions

Authors: I. -G. Shin, C. Han, A. Gould, A. Udalski, T. Sumi, M. Dominik, J. -P. Beaulieu, Y. Tsapras, V. Bozza, M. K. Szymański, M. Kubiak, I. Soszyński, G. Pietrzyński, R. Poleski, K. Ulaczyk, P. Pietrukowicz, S. Kozłowski, J. Skowron, Ł. Wyrzykowski, F. Abe, D. P. Bennett, I. A. Bond, C. S. Botzler, M. Freeman, A. Fukui , et al. (130 additional authors not shown)

Abstract: 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.… ▽ More 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. △ Less

Submitted 2 October, 2012; v1 submitted 11 August, 2012; originally announced August 2012.

Comments: 10 pages, 9 figures, 4 tables