Showing 1–35 of 35 results for author: Ali-Dib, M

Empirical stability criteria for 3D hierarchical triple systems I: Circumbinary planets

Authors: Nikolaos Georgakarakos, Siegfried Eggl, Mohamad Ali-Dib, Ian Dobbs-Dixon

Abstract: In this work we revisit the problem of the dynamical stability of hierarchical triple systems with applications to circumbinary planetary orbits. We carry out more than 3 10^8 numerical simulations of planets between the size of Mercury and the lower fusion boundary (13 Jupiter masses) which revolve around the center of mass of a stellar binary over long timescales. For the first time, three dimen… ▽ More In this work we revisit the problem of the dynamical stability of hierarchical triple systems with applications to circumbinary planetary orbits. We carry out more than 3 10^8 numerical simulations of planets between the size of Mercury and the lower fusion boundary (13 Jupiter masses) which revolve around the center of mass of a stellar binary over long timescales. For the first time, three dimensional and eccentric planetary orbits are considered. We explore systems with a variety of binary and planetary mass ratios, binary and planetary eccentricities from 0 to 0.9 and orbital mutual inclinations ranging from 0 to 180 degrees. The simulation time is set to 10^6 planetary orbital periods. We classify the results of those long term numerical integrations into three categories: stable, unstable and mixed. We provide empirical expressions in the form of multidimensional, parameterized fits for the two borders that separate the three dynamical domains . In addition, we train a machine learning model on our data set in order to have an alternative tool of predicting the stability of circumbinary planets. Both the empirical fits and the machine learning model are tested against randomly generated circumbinary systems with very good results regarding the predictions of orbital stability. The empirical formulae are also applied to the Kepler and TESS circumbinary systems, confirming the stability of the planets in these systems. Finally, we present a REST API with a web based application for convenient access of our simulation data set. △ Less

Submitted 5 September, 2024; v1 submitted 21 April, 2024; originally announced April 2024.

Comments: Accepted for publication in AJ

Protoplanet collisions: new scaling laws from SPH simulations

Authors: Samuele Crespi, Mohamad Ali-Dib, Ian Dobbs-Dixon

Abstract: One common approach for solving collisions between protoplanets in simulations of planet formation is to employ analytical scaling laws. The most widely used one was developed by Leinhardt & Stewart (2012) from a catalog of ~ 180 N-body simulations of rubble-pile collisions. In this work, we use a new catalogue of more than 20,000 SPH simulations to test the validity and the prediction capability… ▽ More One common approach for solving collisions between protoplanets in simulations of planet formation is to employ analytical scaling laws. The most widely used one was developed by Leinhardt & Stewart (2012) from a catalog of ~ 180 N-body simulations of rubble-pile collisions. In this work, we use a new catalogue of more than 20,000 SPH simulations to test the validity and the prediction capability of Leinhardt & Stewart (2012) scaling laws. We find that these laws overestimate the fragmentation efficiency in the merging regime and they are not able to properly reproduce the collision outcomes in the super-catastrophic regime. In the merging regime, we also notice a significant dependence between the collision outcome, in terms of the largest remnant mass, and the relative mass of the colliding protoplanets. Here, we present a new set of scaling laws that are able to better predict the collision outcome in all regimes and it is also able to reproduce the observed dependence on the mass ratio. We compare our new scaling laws against a machine learning approach and obtain similar prediction efficiency. △ Less

Submitted 12 February, 2024; originally announced February 2024.

Comments: 8 pages, 4 figures, accepted for publication in A&A

Journal ref: A&A 685, A86 (2024)

Physics simulation capabilities of LLMs

Authors: Mohamad Ali-Dib, Kristen Menou

Abstract: [Abridged abstract] Large Language Models (LLMs) can solve some undergraduate-level to graduate-level physics textbook problems and are proficient at coding. Combining these two capabilities could one day enable AI systems to simulate and predict the physical world. We present an evaluation of state-of-the-art (SOTA) LLMs on PhD-level to research-level computational physics problems. We conditio… ▽ More [Abridged abstract] Large Language Models (LLMs) can solve some undergraduate-level to graduate-level physics textbook problems and are proficient at coding. Combining these two capabilities could one day enable AI systems to simulate and predict the physical world. We present an evaluation of state-of-the-art (SOTA) LLMs on PhD-level to research-level computational physics problems. We condition LLM generation on the use of well-documented and widely-used packages to elicit coding capabilities in the physics and astrophysics domains. We contribute $\sim 50$ original and challenging problems in celestial mechanics (with REBOUND), stellar physics (with MESA), 1D fluid dynamics (with Dedalus) and non-linear dynamics (with SciPy). Since our problems do not admit unique solutions, we evaluate LLM performance on several soft metrics: counts of lines that contain different types of errors (coding, physics, necessity and sufficiency) as well as a more "educational" Pass-Fail metric focused on capturing the salient physical ingredients of the problem at hand. As expected, today's SOTA LLM (GPT4) zero-shot fails most of our problems, although about 40\% of the solutions could plausibly get a passing grade. About $70-90 \%$ of the code lines produced are necessary, sufficient and correct (coding \& physics). Physics and coding errors are the most common, with some unnecessary or insufficient lines. We observe significant variations across problem class and difficulty. We identify several failure modes of GPT4 in the computational physics domain. Our reconnaissance work provides a snapshot of current computational capabilities in classical physics and points to obvious improvement targets if AI systems are ever to reach a basic level of autonomy in physics simulation capabilities. △ Less

Submitted 2 September, 2024; v1 submitted 4 December, 2023; originally announced December 2023.

Comments: Accepted for publication in Physica Scripta. Abridged abstract. 15 pages + appendix, 1 figure. Comments are welcome

The impermanent fate of massive stars in AGN disks

Authors: Mohamad Ali-Dib, Douglas N. C. Lin

Abstract: Stars are likely to form or to be captured in AGN disks. Their mass reaches an equilibrium when their rate of accretion is balanced by that of wind. If the exchanged gas is well mixed with the stellar core, this metabolic process would indefinitely sustain an "immortal" state on the main sequence (MS) and pollute the disk with He byproducts. This theoretical extrapolation is inconsistent with the… ▽ More Stars are likely to form or to be captured in AGN disks. Their mass reaches an equilibrium when their rate of accretion is balanced by that of wind. If the exchanged gas is well mixed with the stellar core, this metabolic process would indefinitely sustain an "immortal" state on the main sequence (MS) and pollute the disk with He byproducts. This theoretical extrapolation is inconsistent with the super-solar α element and Fe abundances inferred from the broad emission lines in active AGNs with modest He concentration. We show this paradox can be resolved with a highly-efficient retention of the He ashes or the suppression of chemical blending. The latter mechanism is robust in the geometrically-thin, dense, sub-pc regions of the disk where the embedded-stars' mass is limited by the gap-formation condition. These stars contain a radiative zone between their mass-exchange stellar surface and the nuclear-burning core. Insulation of the core lead to the gradual decrease of its H fuel and the stars' equilibrium masses. These stars transition to their post-main-sequence (PostMS) tracks on a chemical evolution time scale of a few Myr. Subsequently, the triple-α and α-chain reactions generate α and Fe byproducts which are released into their natal disks. These PostMS stars also undergo core collapse, set off type II supernova, and leave behind a few solar-mass residual black holes or neutron stars △ Less

Submitted 18 September, 2023; v1 submitted 8 September, 2023; originally announced September 2023.

Comments: 17 pages, 7 figures, Accepted for publication in MNRAS

Vanadium oxide and a sharp onset of cold-trapping on a giant exoplanet

Authors: Stefan Pelletier, Björn Benneke, Mohamad Ali-Dib, Bibiana Prinoth, David Kasper, Andreas Seifahrt, Jacob L. Bean, Florian Debras, Baptiste Klein, Luc Bazinet, H. Jens Hoeijmakers, Aurora Y. Kesseli, Olivia Lim, Andres Carmona, Lorenzo Pino, Núria Casasayas-Barris, Thea Hood, Julian Stürmer

Abstract: The abundance of refractory elements in giant planets can provide key insights into their formation histories. Due to the Solar System giants' low temperatures, refractory elements condense below the cloud deck limiting sensing capabilities to only highly volatile elements. Recently, ultra-hot giant exoplanets have allowed for some refractory elements to be measured showing abundances broadly cons… ▽ More The abundance of refractory elements in giant planets can provide key insights into their formation histories. Due to the Solar System giants' low temperatures, refractory elements condense below the cloud deck limiting sensing capabilities to only highly volatile elements. Recently, ultra-hot giant exoplanets have allowed for some refractory elements to be measured showing abundances broadly consistent with the solar nebula with titanium likely condensed out of the photosphere. Here we report precise abundance constraints of 14 major refractory elements on the ultra-hot giant planet WASP-76b that show distinct deviations from proto-solar, and a sharp onset in condensation temperature. In particular, we find nickel to be enriched, a possible sign of the accretion of a differentiated object's core during the planet's evolution. Elements with condensation temperatures below 1,550 K otherwise closely match those of the Sun before sharply transitioning to being strongly depleted above 1,550 K, well explained by nightside cold-trapping. We further unambiguously detect vanadium oxide on WASP-76b, a molecule long hypothesized to drive atmospheric thermal inversions, and also observe a global east-west asymmetry in its absorption signals. Overall, our findings indicate that giant planets have a mostly stellar-like refractory elemental content and suggest that temperature sequences of hot Jupiter spectra can show abrupt transitions wherein a mineral species is either present, or completely absent if a cold-trap exists below its condensation temperature. △ Less

Submitted 14 June, 2023; originally announced June 2023.

Comments: Published online in Nature on June 14, 2023

JWST/NIRCam discovery of the first Y+Y brown dwarf binary: WISE J033605.05$-$014350.4

Authors: Per Calissendorff, Matthew De Furio, Michael Meyer, Loïc Albert, Christian Aganze, Mohamad Ali-Dib, Daniella C. Bardalez Gagliuffi, Frederique Baron, Charles A. Beichman, Adam J. Burgasser, Michael C. Cushing, Jacqueline Kelly Faherty, Clémence Fontanive, Christopher R. Gelino, John E. Gizis, Alexandra Z. Greenbaum, J. Davy Kirkpatrick, Sandy K. Leggett, Frantz Martinache, David Mary, Mamadou N'Diaye, Benjamin J. S. Pope, Thomas L Roellig, Johannes Sahlmann, Anand Sivaramakrishnan , et al. (3 additional authors not shown)

Abstract: We report the discovery of the first brown dwarf binary system with a Y dwarf primary, WISE J033605.05$-$014350.4, observed with NIRCam on JWST with the F150W and F480M filters. We employed an empirical point spread function binary model to identify the companion, located at a projected separation of 84 milliarcseconds, position angle of 295 degrees, and with contrast of 2.8 and 1.8 magnitudes in… ▽ More We report the discovery of the first brown dwarf binary system with a Y dwarf primary, WISE J033605.05$-$014350.4, observed with NIRCam on JWST with the F150W and F480M filters. We employed an empirical point spread function binary model to identify the companion, located at a projected separation of 84 milliarcseconds, position angle of 295 degrees, and with contrast of 2.8 and 1.8 magnitudes in F150W and F480M, respectively. At a distance of 10$\,$pc based on its Spitzer parallax, and assuming a random inclination distribution, the physical separation is approximately 1$\,$au. Evolutionary models predict for that an age of 1-5 Gyr, the companion mass is about 4-12.5 Jupiter masses around the 7.5-20 Jupiter mass primary, corresponding to a companion-to-host mass fraction of $q=0.61\pm0.05$. Under the assumption of a Keplerian orbit the period for this extreme binary is in the range of 5-9 years. The system joins a small but growing sample of ultracool dwarf binaries with effective temperatures of a few hundreds of Kelvin. Brown dwarf binaries lie at the nexus of importance for understanding the formation mechanisms of these elusive objects, as they allow us to investigate whether the companions formed as stars or as planets in a disk around the primary. △ Less

Submitted 29 March, 2023; originally announced March 2023.

Comments: 8 pages, 3 figures, 1 table. Accepted for publication in Astrophysical Journal Letters

Dust processing in protoplanetary envelopes as the origin of hot minerals in comets

Authors: Mohamad Ali-Dib

Abstract: Crystalline silicates are found in a large number of comets. These pose a long-standing conundrum for solar system formation models as they can only be created in the inner hot disk at temperatures higher than 800 K, and there is no obvious mechanism to transport them out into the comets formation region. Here we propose that these particles could have formed inside the hydrostatic envelopes surro… ▽ More Crystalline silicates are found in a large number of comets. These pose a long-standing conundrum for solar system formation models as they can only be created in the inner hot disk at temperatures higher than 800 K, and there is no obvious mechanism to transport them out into the comets formation region. Here we propose that these particles could have formed inside the hydrostatic envelopes surrounding young protoplanets still embedded in the protoplanetary disk. Using a simplified 1D model we investigate the thermal structure of these envelopes, and find that for core masses ranging from 0.08 to 1.5 M_Earth, located anywhere between 1 and 30 AU, the temperature and pressure at the base of the envelopes are high enough to quickly vaporize silicate particles of various sizes. Moreover, if the grain abundance is at least solar, these envelopes become fully convective, allowing for dust ejection across the Bondi radius back into the disk. Amorphous silicates are hence thermally processed into crystalline particles in these envelopes, and then transported back to disk through convective diffusion to be finally incorporated into the cometary building blocks. △ Less

Submitted 4 January, 2023; originally announced January 2023.

Comments: Accepted for publication in MNRAS. 4 pages, 2 figures

A machine-generated catalogue of Charon's craters and implications for the Kuiper belt

Authors: Mohamad Ali-Dib

Abstract: In this paper we investigate Charon's craters size distribution using a deep learning model. This is motivated by the recent results of Singer et al. (2019) who, using manual cataloging, found a change in the size distribution slope of craters smaller than 12 km in diameter, translating into a paucity of small Kuiper Belt objects. These results were corroborated by Robbins and Singer (2021), but o… ▽ More In this paper we investigate Charon's craters size distribution using a deep learning model. This is motivated by the recent results of Singer et al. (2019) who, using manual cataloging, found a change in the size distribution slope of craters smaller than 12 km in diameter, translating into a paucity of small Kuiper Belt objects. These results were corroborated by Robbins and Singer (2021), but opposed by Morbidelli et al. (2021), necessitating an independent review. Our MaskRCNN-based ensemble of models was trained on Lunar, Mercurian, and Martian crater catalogues and both optical and digital elevation images. We use a robust image augmentation scheme to force the model to generalize and transfer-learn into icy objects. With no prior bias or exposure to Charon, our model find best fit slopes of q =-1.47+-0.33 for craters smaller than 10 km, and q =-2.91+-0.51 for craters larger than 15 km. These values indicate a clear change in slope around 15 km as suggested by Singer et al. (2019) and thus independently confirm their conclusions. Our slopes however are both slightly flatter than those found more recently by Robbins and Singer (2021). Our trained models and relevant codes are available online on github.com/malidib/ACID . △ Less

Submitted 16 June, 2022; originally announced June 2022.

Comments: 16 pages, 2 figures, accepted for publication in Icarus

Using Artificial Intelligence and real galaxy images to constrain parameters in galaxy formation simulations

Authors: Andrea V. Macciò, Mohamad Ali-Dib, Pavle Vulanović, Hind Al Noori, Fabian Walter, Nico Krieger, Tobias Buck

Abstract: Cosmological galaxy formation simulations are still limited by their spatial/mass resolution and cannot model from first principles some of the processes, like star formation, that are key in driving galaxy evolution. As a consequence they still rely on a set of 'effective parameters' that try to capture the scales and the physical processes that cannot be directly resolved in the simulation. In t… ▽ More Cosmological galaxy formation simulations are still limited by their spatial/mass resolution and cannot model from first principles some of the processes, like star formation, that are key in driving galaxy evolution. As a consequence they still rely on a set of 'effective parameters' that try to capture the scales and the physical processes that cannot be directly resolved in the simulation. In this study we show that it is possible to use Machine Learning techniques applied to real and simulated images of galaxies to discriminate between different values of these parameters by making use of the full information content of an astronomical image instead of collapsing it into a limited set of values like size, or stellar/ gas masses. In this work we apply our method to the NIHAO simulations and the THINGS and VLA-ANGST observations of HI maps in nearby galaxies to test the ability of different values of the star formation density threshold $n$ to reproduce observed HI maps. We show that observations indicate the need for a high value of $n \gtrsim 80$ ,cm$^{-3}$ (although the exact numerical value is model-dependent), which has important consequences for the dark matter distribution in galaxies. Our study shows that with innovative methods it is possible to take full advantage of the information content of galaxy images and compare simulations and observations in an interpretable, non-parametric and quantitative manner. △ Less

Submitted 18 February, 2022; originally announced February 2022.

Comments: 8 pages, 5 figures, accepted for publication on MNRAS

The effect of late giant collisions on the atmospheres of protoplanets and the formation of cold sub-Saturns

Authors: Mohamad Ali-Dib, Andrew Cumming, Douglas N. C. Lin

Abstract: We investigate the origins of cold sub-Saturns (CSS), an exoplanetary population inferred from microlensing surveys. If confirmed, these planets would rebut a theorised gap in planets' mass distribution between those of Neptune and Jupiter caused by the rapid runaway accretion of super-critical cores. In an attempt to resolve this theoretical-observational disparity, we examine the outcomes of gia… ▽ More We investigate the origins of cold sub-Saturns (CSS), an exoplanetary population inferred from microlensing surveys. If confirmed, these planets would rebut a theorised gap in planets' mass distribution between those of Neptune and Jupiter caused by the rapid runaway accretion of super-critical cores. In an attempt to resolve this theoretical-observational disparity, we examine the outcomes of giant collisions between sub-critical protoplanets. Due to the secular interaction among protoplanets, these events may occur in rapidly depleting discs. We show that impactors ~ 5% the mass of near-runaway envelopes around massive cores can efficiently remove these envelopes entirely via a thermally-driven super-Eddington wind emanating from the core itself, in contrast with the stellar Parker winds usually considered. After a brief cooling phase, the merged cores resume accretion. But, the evolution timescale of transitional discs is too brief for the cores to acquire sufficiently massive envelopes to undergo runaway accretion despite their large combined masses. Consequently, these events lead to the emergence of CSS without their transformation into gas giants. We show that these results are robust for a wide range of disc densities, grain opacities and silicate abundance in the envelope. Our fiducial case reproduces CSS with heavy (>= 30 M_Earth) cores and less massive (a few M_Earth) sub-critical envelopes. We also investigate the other limiting cases, where continuous mergers of comparable-mass cores yield CSS with wider ranges of core-to-envelope mass ratios and envelope opacities. Our results indicate that it is possible for CSS and Uranus and Neptune to emerge within the framework of well studied processes and they may be more common than previously postulated. △ Less

Submitted 15 October, 2021; originally announced October 2021.

Comments: Accepted for publication in MNRAS, 17 pages, 9 figures

The rarity of very red TNOs in the scattered disk

Authors: Mohamad Ali-Dib, Michael Marsset, Wing-Cheung Wong, Rola Dbouk

Abstract: We investigate the origins of the photometrically Very Red and Less Red Trans-Neptunian Objects. We first reanalyse the dataset of Marsset et al. 2019 and find that, in addition to the known color-inclination correlation in hot TNOs, a similar trend exists for color-eccentricity. We show that VR TNOs are sharply constrained to eccentricities < 0.42 and inclinations < 21 deg, leading to a paucity o… ▽ More We investigate the origins of the photometrically Very Red and Less Red Trans-Neptunian Objects. We first reanalyse the dataset of Marsset et al. 2019 and find that, in addition to the known color-inclination correlation in hot TNOs, a similar trend exists for color-eccentricity. We show that VR TNOs are sharply constrained to eccentricities < 0.42 and inclinations < 21 deg, leading to a paucity of VR scattered disk and distant MMR objects. We then interpret these findings using N-body simulations accounting for Neptune's outward migration into a massless particles disk, and find that these observations are best reproduced with a LR-to-VR color transition line between 38 and 42 AU in the primordial disk, separating the objects' formation locations. For an initial surface density profile $Σ\propto 1/r^2$, a color transition around 38 AU is needed to explain the high abundance of VR plutinos but creates too many VR scattered disk objects, while a transition line around 42 AU seems to better reproduces the scattered disk colors but creates virtually no VR plutinos. Our simulations furthermore show that the rarity of VR particles at high eccentricity is possibly due to the absence of sweeping higher order MMRs, and secular resonances, beyond 42 AU. Inspecting individual populations, we show that the majority of VR SDOs originate as objects trapped in Neptune's second and third order MMRs. These then evolve due to diffusion, scattering, Kozai-Lidov cycles, and secular resonances into their current orbits. Future unbiased color surveys are crucial to better constrain the TNOs dynamical origins. △ Less

Submitted 9 April, 2021; originally announced April 2021.

Comments: Accepted for publication in The Astronomical Journal. 25 pages, 15 figures

Constraining protoplanetary disks with exoplanetary dynamics: Kepler-419 as an example

Authors: Mohamad Ali-Dib, Cristobal Petrovich

Abstract: We investigate the origins of Kepler-419, a peculiar system hosting two nearly coplanar and highly eccentric gas giants with apsidal orientations librating around anti-alignment, and use this system to place constraints on the properties of their birth protoplanetary disk. We follow the proposal by Petrovich, Wu, & Ali-Dib (2019) that these planets have been placed on these orbits as a natural res… ▽ More We investigate the origins of Kepler-419, a peculiar system hosting two nearly coplanar and highly eccentric gas giants with apsidal orientations librating around anti-alignment, and use this system to place constraints on the properties of their birth protoplanetary disk. We follow the proposal by Petrovich, Wu, & Ali-Dib (2019) that these planets have been placed on these orbits as a natural result of the precessional effects of a dissipating massive disk and extend it by using direct N-body simulations and models for the evolution of the gas disks, including photoevaporation. Based on a parameter space exploration, we find that in order to reproduce the system the initial disk mass had to be at least 95 M_Jup and dissipate on a timescale of at least 10^4 yr. This mass is consistent with the upper end of the observed disk masses distribution, and the dissipation timescale is consistent with photoevaporation models. We study the properties of such disks using simplified 1D thin disk models and show that they are gravitationally stable, indicating that the two planets must have formed via core accretion and thus prone to disk migration. We hence finally investigate the sensitivity of this mechanism to the outer planet's semi major axis, and find that the nearby 7:1, 8:1, and 9:1 mean-motion resonances can completely quench this mechanism, while even higher order resonances can also significantly affect the system. Assuming the two planets avoid these high order resonances and/or close encounters, the dynamics seems to be rather insensitive to planet c semi major axis, and thus orbital migration driven by the disk. △ Less

Submitted 14 September, 2020; originally announced September 2020.

Comments: 8 pages, 5 figures, accepted for publication in MNRAS

The Imprint of the Protoplanetary Disk in the Accretion of Super-Earth Envelopes

Authors: Mohamad Ali-Dib, Andrew Cumming, Douglas N. C. Lin

Abstract: Super-Earths are by far the most dominant type of exoplanet, yet their formation is still not well understood. In particular, planet formation models predict that many of them should have accreted enough gas to become gas giants. Here we examine the role of the protoplanetary disk in the cooling and contraction of the protoplanetary envelope. In particular, we investigate the effects of 1) the the… ▽ More Super-Earths are by far the most dominant type of exoplanet, yet their formation is still not well understood. In particular, planet formation models predict that many of them should have accreted enough gas to become gas giants. Here we examine the role of the protoplanetary disk in the cooling and contraction of the protoplanetary envelope. In particular, we investigate the effects of 1) the thermal state of the disk as set by the relative size of heating by accretion or irradiation, and whether its energy is transported by radiation or convection, and 2) advection of entropy into the outer envelope by disk flows that penetrate the Hill sphere, as found in 3D global simulations. We find that, at 5 and 1 AU, this flow at the level reported in the non-isothermal simulations where it penetrates only to ~ 0.3 times the Hill radius has little effect on the cooling rate since most of the envelope mass is concentrated close to the core, and far from the flow. On the other hand, at 0.1 AU, the envelope quickly becomes fully-radiative, nearly isothermal, and thus cannot cool down, stalling gas accretion. This effect is significantly more pronounced in convective disks, leading to envelope mass orders of magnitude lower. Entropy advection at 0.1 AU in either radiative or convective disks could therefore explain why super-Earths failed to undergo runaway accretion. These results highlight the importance of the conditions and energy transport in the protoplanetary disk for the accretion of planetary envelopes. △ Less

Submitted 25 March, 2020; v1 submitted 3 March, 2020; originally announced March 2020.

Comments: 10 pages, 6 figures, accepted for publication in MNRAS

LRP2020: Signposts of planet formation in protoplanetary disks

Authors: Nienke van der Marel, Ruobing Dong, Ralph Pudritz, James Wadsley, Aaron Boley, Eve Lee, Mohamad Ali-Dib, Brenda Matthews, Christian Marois, Henry Ngo

Abstract: Successful exoplanet surveys in the last decade have revealed that planets are ubiquitous throughout the Milky Way, and show a large diversity in mass, location and composition. At the same time, new facilities such as the Atacama Large Millimeter/submillimeter Array (ALMA) and optical/infrared facilities including Gemini/GPI have provided us with sharper images than ever before of protoplanetary… ▽ More Successful exoplanet surveys in the last decade have revealed that planets are ubiquitous throughout the Milky Way, and show a large diversity in mass, location and composition. At the same time, new facilities such as the Atacama Large Millimeter/submillimeter Array (ALMA) and optical/infrared facilities including Gemini/GPI have provided us with sharper images than ever before of protoplanetary disks around young stars, the birth cradles of planets. The high spatial resolution has revealed astonishing structures in disks, such as rings, gaps, asymmetries and spiral arms, and the enormous jump in sensitivity has provided the tools for both large, statistically relevant surveys and deep, sensitive molecular line studies. These observations have revolutionized our view of planet formation, disk formation and disk evolution, bringing model simulations and observations closer to the same level of detail, with many contributions from Canadian researchers on theoretical, observational and technological sides. The new results have inevitably led to a range of new questions, which require next generation instruments such as the Next Generation Very Large Array (ngVLA) and large scale optical infrared facilities. In this white paper we will discuss the current transformation in our understanding of planet formation and the next steps and challenges in connecting theory with exoplanet demographics and protoplanetary disk observations for Canadian research. △ Less

Submitted 30 September, 2019; originally announced October 2019.

Comments: White paper E013 submitted to the Canadian Long Range Plan 2020

Automated crater shape retrieval using weakly-supervised deep learning

Authors: Mohamad Ali-Dib, Kristen Menou, Alan P. Jackson, Chenchong Zhu, Noah Hammond

Abstract: Crater ellipticity determination is a complex and time consuming task that so far has evaded successful automation. We train a state of the art computer vision algorithm to identify craters in Lunar digital elevation maps and retrieve their sizes and 2D shapes. The computational backbone of the model is MaskRCNN, an "instance segmentation" general framework that detects craters in an image while s… ▽ More Crater ellipticity determination is a complex and time consuming task that so far has evaded successful automation. We train a state of the art computer vision algorithm to identify craters in Lunar digital elevation maps and retrieve their sizes and 2D shapes. The computational backbone of the model is MaskRCNN, an "instance segmentation" general framework that detects craters in an image while simultaneously producing a mask for each crater that traces its outer rim. Our post-processing pipeline then finds the closest fitting ellipse to these masks, allowing us to retrieve the crater ellipticities. Our model is able to correctly identify 87\% of known craters in the longitude range we hid from the network during training and validation (test set), while predicting thousands of additional craters not present in our training data. Manual validation of a subset of these "new" craters indicates that a majority of them are real, which we take as an indicator of the strength of our model in learning to identify craters, despite incomplete training data. The crater size, ellipticity, and depth distributions predicted by our model are consistent with human-generated results. The model allows us to perform a large scale search for differences in crater diameter and shape distributions between the lunar highlands and maria, and we exclude any such differences with a high statistical significance. The predicted test set catalogue and trained model are available here: https://github.com/malidib/Craters_MaskRCNN/. △ Less

Submitted 11 March, 2020; v1 submitted 20 June, 2019; originally announced June 2019.

Comments: 59 pages, 13 figures, Accepted for publication in Icarus

Limits on Protoplanet Growth by Accretion of Small Solids

Authors: Mohamad Ali-Dib, Christopher Thompson

Abstract: This paper identifies constraints on the growth of a small planetary core (0.3 M$_{\oplus}$) that accretes millimeter-sized pebbles from a gaseous disk. We construct time-dependent spherical envelope models that capture physical processes not included in existing global hydrodynamic simulations, including particle size evolution, dust transport, and realistic gas equations of state. We assume a lo… ▽ More This paper identifies constraints on the growth of a small planetary core (0.3 M$_{\oplus}$) that accretes millimeter-sized pebbles from a gaseous disk. We construct time-dependent spherical envelope models that capture physical processes not included in existing global hydrodynamic simulations, including particle size evolution, dust transport, and realistic gas equations of state. We assume a low enough disk density that pebbles are marginally coupled to the gas and are trapped efficiently near the core Bondi radius. Pebbles then drift rapidly enough to experience erosion by sandblasting, mutual collisions, and sublimation of water ice. We find that pebble fragmentation is more efficient than dust re-sticking. Therefore the large pebble accretion rate $\dot{M}_p$ needed to build a core of mass > M$_{\oplus}$ leads to a high envelope metallicity and grain opacity. Above $\dot{M}_p$ $\sim 10^{-7}$ M$_{\oplus}$/yr, and without other luminosity sources, convective motions expand near the Bondi radius. The warm, dusty, and turbulent envelope buffers the inward drift of pebble debris: given a turbulent concentration factor $f_{turb} \geq$ 1 near the lower convective boundary, the core growth rate is limited to $1 \times 10^{-7} f_{turb}$ M$_{\oplus}$/yr and the e-folding time 3$f_{turb}^{-1}$ Myr. The remainder of the solid debris is expelled as highly processed silicates. Pebble ice never reaches the core, and the envelope contains comparable amounts of H2/He and metals. We interpret our results using simpler steady models and semi-analytical estimates. Future global simulations incorporating the processes modelled here are needed to understand the influence of rotation and vertical disk structure. △ Less

Submitted 9 July, 2020; v1 submitted 25 February, 2019; originally announced February 2019.

Comments: 41 pages, 17 figures. Accepted for publication in the Astrophysical Journal. Expanded analytic treatment of threshold for refluxing. Corrected treatment of ice sublimation in the envelope. Effect added of heating due to the accretion of a fraction of pebble debris onto the core

Journal ref: The Astrophysical Journal, Volume 900, Number 2 (2020)

Secular transport during disk dispersal: the case of Kepler-419

Authors: Cristobal Petrovich, Yanqin Wu, Mohamad Ali-Dib

Abstract: Due to fortuitous circumstances, the two giant planets around Kepler-419 have well characterized 3-D orbits. They are nearly coplanar to each other; the inner one has a large eccentricity (~0.82); and the apses of the two orbits librate around anti-alignment. Such a state defies available proposals for large eccentricities. We argue that it is instead uniquely produced by a decaying protoplanetary… ▽ More Due to fortuitous circumstances, the two giant planets around Kepler-419 have well characterized 3-D orbits. They are nearly coplanar to each other; the inner one has a large eccentricity (~0.82); and the apses of the two orbits librate around anti-alignment. Such a state defies available proposals for large eccentricities. We argue that it is instead uniquely produced by a decaying protoplanetary disk. When the disk was massive, its precessional effect on the planets forced the two apses to center around an anti-aligned state. And as the disk is gradually eroded, the pair of planets are adiabatically transported to a new state where most of the eccentricity (or rather, the angular momentum deficit) is transferred to the inner planet, and the two apses are largely anti-aligned. During this transport, any initial mutual inclination may be reduced or enhanced; either may be compatible with the current constraints. So a primordial disk can drive up planet eccentricities both in resonant planet pairs (as has been shown for GJ 876) and in secularly interacting, non-resonant pairs. The mechanism discussed here may be relevant for forming hot Jupiters and for explaining the observed eccentricities of warm and cold giant planets. △ Less

Submitted 20 November, 2018; v1 submitted 10 September, 2018; originally announced September 2018.

Comments: accepted for publication in AJ, 14 pages, 7 figures

Lunar Crater Identification via Deep Learning

Authors: Ari Silburt, Mohamad Ali-Dib, Chenchong Zhu, Alan Jackson, Diana Valencia, Yevgeni Kissin, Daniel Tamayo, Kristen Menou

Abstract: Crater counting on the Moon and other bodies is crucial to constrain the dynamical history of the Solar System. This has traditionally been done by visual inspection of images, thus limiting the scope, efficiency, and/or accuracy of retrieval. In this paper we demonstrate the viability of using convolutional neural networks (CNNs) to determine the positions and sizes of craters from Lunar digital… ▽ More Crater counting on the Moon and other bodies is crucial to constrain the dynamical history of the Solar System. This has traditionally been done by visual inspection of images, thus limiting the scope, efficiency, and/or accuracy of retrieval. In this paper we demonstrate the viability of using convolutional neural networks (CNNs) to determine the positions and sizes of craters from Lunar digital elevation maps (DEMs). We recover 92% of craters from the human-generated test set and almost double the total number of crater detections. Of these new craters, 15% are smaller in diameter than the minimum crater size in the ground-truth dataset. Our median fractional longitude, latitude and radius errors are 11% or less, representing good agreement with the human-generated datasets. From a manual inspection of 361 new craters we estimate the false positive rate of new craters to be 11%. Moreover, our Moon-trained CNN performs well when tested on DEM images of Mercury, detecting a large fraction of craters in each map. Our results suggest that deep learning will be a useful tool for rapidly and automatically extracting craters on various Solar System bodies. We make our code and data publicly available at https://github.com/silburt/DeepMoon.git and https://doi.org/10.5281/zenodo.1133969 . △ Less

Submitted 12 November, 2018; v1 submitted 6 March, 2018; originally announced March 2018.

Comments: 41 pages, 5 figures, Accepted for publication in Icarus. Comments and requests are welcomed

Ejection of rocky and icy material from binary star systems: Implications for the origin and composition of 1I/`Oumuamua

Authors: Alan P. Jackson, Daniel Tamayo, Noah Hammond, Mohamad Ali-Dib, Hanno Rein

Abstract: In single star systems like our own Solar system, comets dominate the mass budget of bodies that are ejected into interstellar space, since they form further away and are less tightly bound. However 1I/`Oumuamua, the first interstellar object detected, appears asteroidal in its spectra and in its lack of detectable activity. We argue that the galactic budget of interstellar objects like 1I/`Oumuam… ▽ More In single star systems like our own Solar system, comets dominate the mass budget of bodies that are ejected into interstellar space, since they form further away and are less tightly bound. However 1I/`Oumuamua, the first interstellar object detected, appears asteroidal in its spectra and in its lack of detectable activity. We argue that the galactic budget of interstellar objects like 1I/`Oumuamua should be dominated by planetesimal material ejected during planet formation in circumbinary systems, rather than in single star systems or widely separated binaries. We further show that in circumbinary systems, rocky bodies should be ejected in comparable numbers to icy ones. This suggests that a substantial fraction of additional interstellar objects discovered in the future should display an active coma. We find that the rocky population, of which 1I/`Oumuamua seems to be a member, should be predominantly sourced from A-type and late B-star binaries. △ Less

Submitted 28 March, 2018; v1 submitted 12 December, 2017; originally announced December 2017.

Comments: 5 pages, 2 figures. accepted in MNRAS Letters, updated with reviewer comments

What is Neptune's D/H ratio really telling us about its water abundance ?

Authors: Mohamad Ali-Dib, Gunjan Lakhlani

Abstract: We investigate the deep water abundance of Neptune using a simple 2-component (core + envelope) toy model. The free parameters of the model are the total mass of heavy elements in the planet (Z), the mass fraction of Z in the envelope (f_env), and the D/H ratio of the accreted building blocks (D/H_build ). We systematically search the allowed parameter space on a grid and constrain it using Neptun… ▽ More We investigate the deep water abundance of Neptune using a simple 2-component (core + envelope) toy model. The free parameters of the model are the total mass of heavy elements in the planet (Z), the mass fraction of Z in the envelope (f_env), and the D/H ratio of the accreted building blocks (D/H_build ). We systematically search the allowed parameter space on a grid and constrain it using Neptune's bulk carbon abundance, D/H ratio, and interior structure models. Assuming solar C/O ratio and cometary D/H for the accreted building blocks forming the planet, we can fit all of median ~ 7%), and the rest the constraints if less than ~ 15% of Z is in the envelope (f_env is locked in a solid core. This model predicts a maximum bulk oxygen abundance in Neptune of 65 \times solar value. If we assume a C/O of 0.17, corresponding to clathrate-hydrates building blocks, we predict a maximum oxygen abundance of 200 times solar value with a median value of ~ 140. Thus, both cases lead to an oxygen abundance significantly lower than the preferred value of Cavalié et al. (2017) (~ 540 times solar), inferred from model dependent deep CO observations. Such high water abundances are excluded by our simple but robust model. We attribute this discrepancy to our imperfect understanding of either the interior structure of Neptune or the chemistry of the primordial protosolar nebula. △ Less

Submitted 26 January, 2018; v1 submitted 13 September, 2017; originally announced September 2017.

Comments: 6 pages, 5 figures, Accepted for publication in MNRAS

Possible formation pathways for the low density Neptune-mass planet HAT-P-26b

Authors: Mohamad Ali-Dib, Gunjan Lakhlani

Abstract: We investigate possible pathways for the formation of the low density Neptune-mass planet HAT-P-26b. We use two formation different models based on pebbles and planetesimals accretion, and includes gas accretion, disk migration and simple photoevaporation. The models tracks the atmospheric oxygen abundance, in addition to the orbital period, and mass of the forming planets, that we compare to HAT-… ▽ More We investigate possible pathways for the formation of the low density Neptune-mass planet HAT-P-26b. We use two formation different models based on pebbles and planetesimals accretion, and includes gas accretion, disk migration and simple photoevaporation. The models tracks the atmospheric oxygen abundance, in addition to the orbital period, and mass of the forming planets, that we compare to HAT-P-26b. We find that pebbles accretion can explain this planet more naturally than planetesimals accretion that fails completely unless we artificially enhance the disk metallicity significantly. Pebble accretion models can reproduce HAT-P-26b with either a high initial core mass and low amount of envelope enrichment through core erosion or pebbles dissolution, or the opposite, with both scenarios being possible. Assuming a low envelope enrichment factor as expected from convection theory and comparable to the values we can infer from the D/H measurements in Uranus and Neptune, our most probable formation pathway for HAT-P-26b is through pebble accretion starting around 10 AU early in the disk's lifetime. △ Less

Submitted 11 September, 2017; v1 submitted 22 May, 2017; originally announced May 2017.

Comments: 8 pages, 6 figures, accepted for publication in MNRAS

The origin of the occurrence rate profile of gas giants inside 100 days

Authors: Mohamad Ali-Dib, Anders Johansen, Chelsea X. Huang

Abstract: We investigate the origin of the period distribution of giant planets. We fit the bias-corrected distribution of gas-giant planets inside 300 days found by Santerne et al. (2016) using a planet formation model based on pebble accretion. We investigate two possible initial conditions: a linear distribution of planetary seeds, and seeds injected exclusively on the water and CO icelines. Our simulati… ▽ More We investigate the origin of the period distribution of giant planets. We fit the bias-corrected distribution of gas-giant planets inside 300 days found by Santerne et al. (2016) using a planet formation model based on pebble accretion. We investigate two possible initial conditions: a linear distribution of planetary seeds, and seeds injected exclusively on the water and CO icelines. Our simulations exclude the linear initial distribution of seeds with a high degree of confidence. Our bimodal model based on snowlines give a more reasonable fit to the data, with the discrepancies reducing significantly if we assume the water snowline to be a factor 3-10 less efficient at producing planets. This model moreover performs better on both the warm/hot Jupiters ratio and a Gaussian mixture model as comparison criteria. Our results hint that the gas-giant exoplanets population inside 300 days is more compatible with planets forming preferentially at special locations. △ Less

Submitted 2 June, 2017; v1 submitted 20 April, 2017; originally announced April 2017.

Comments: 8 pages, 6 figures, Accepted for publication in MNRAS. Comments welcome

Disentangling Hot Jupiters formation location from their chemical composition

Authors: Mohamad Ali-Dib

Abstract: We use a population synthesis model that includes pebbles and gas accretion, planetary migration, and a simplified chemistry scheme to study the formation of hot-Jupiters. Models have been proposed that these planets can either originate beyond the snowline and then move inward via disk migration, or form "in-situ" inside the snowline. The goal of this work is to verify which of these two scenario… ▽ More We use a population synthesis model that includes pebbles and gas accretion, planetary migration, and a simplified chemistry scheme to study the formation of hot-Jupiters. Models have been proposed that these planets can either originate beyond the snowline and then move inward via disk migration, or form "in-situ" inside the snowline. The goal of this work is to verify which of these two scenarios is more compatible with pebble accretion, and whether we can distinguish observationally between them via the resulting planetary C/O ratios and core masses. Our results show that, for solar system composition, the C/O ratios will vary but moderately between the two populations, since a significant amount of carbon and oxygen are locked up in refractories. In this case, we find a strong correlation between the carbon and oxygen abundances and core mass. The C/O ratio variations are more pronounced in the case where we assume that all carbon and oxygen are in volatiles. On average, Hot-Jupiters forming "in-situ" inside the snowline will have higher C/O ratios because they accrete less water ice. However, only Hot-Jupiters forming in-situ around stars with C/O=0.8 can have a C/O ratio higher than unity. We finally find that, even with fast pebble accretion, it is significantly easier to form Hot-Jupiters outside of the snowline, even if forming these "in-situ" is not impossible in the limit of the simplifying assumptions made. △ Less

Submitted 26 January, 2017; v1 submitted 9 November, 2016; originally announced November 2016.

Comments: Accepted for publication in MNRAS, 12 pages, 10 figures. Feedback is welcomed

A Machine Learns to Predict the Stability of Tightly Packed Planetary Systems

Authors: Daniel Tamayo, Ari Silburt, Diana Valencia, Kristen Menou, Mohamad Ali-Dib, Cristobal Petrovich, Chelsea X. Huang, Hanno Rein, Christa van Laerhoven, Adiv Paradise, Alysa Obertas, Norman Murray

Abstract: The requirement that planetary systems be dynamically stable is often used to vet new discoveries or set limits on unconstrained masses or orbital elements. This is typically carried out via computationally expensive N-body simulations. We show that characterizing the complicated and multi-dimensional stability boundary of tightly packed systems is amenable to machine learning methods. We find tha… ▽ More The requirement that planetary systems be dynamically stable is often used to vet new discoveries or set limits on unconstrained masses or orbital elements. This is typically carried out via computationally expensive N-body simulations. We show that characterizing the complicated and multi-dimensional stability boundary of tightly packed systems is amenable to machine learning methods. We find that training an XGBoost machine learning algorithm on physically motivated features yields an accurate classifier of stability in packed systems. On the stability timescale investigated ($10^7$ orbits), it is 3 orders of magnitude faster than direct N-body simulations. Optimized machine learning classifiers for dynamical stability may thus prove useful across the discipline, e.g., to characterize the exoplanet sample discovered by the upcoming Transiting Exoplanet Survey Satellite (TESS). This proof of concept motivates investing computational resources to train algorithms capable of predicting stability over longer timescales and over broader regions of phase space. △ Less

Submitted 21 November, 2016; v1 submitted 17 October, 2016; originally announced October 2016.

Comments: Accepted in ApJ letters. 7 pages, 3 figures

A pebbles accretion model with chemistry and implications for the solar system

Authors: Mohamad Ali-Dib

Abstract: We investigate the chemical composition of the solar system's giant planets atmospheres using a physical formation model with chemistry. The model incorporate disk evolution, pebbles and gas accretion, type I and II migration, simplified disk photoevaporation and solar system chemical measurements. We track the chemical compositions of the formed giant planets and compare them to the observed valu… ▽ More We investigate the chemical composition of the solar system's giant planets atmospheres using a physical formation model with chemistry. The model incorporate disk evolution, pebbles and gas accretion, type I and II migration, simplified disk photoevaporation and solar system chemical measurements. We track the chemical compositions of the formed giant planets and compare them to the observed values. Two categories of models are studied: with and without disk chemical enrichment via photoevaporation. Predictions for the Oxygen and Nitrogen abundances, core masses, and total amount of heavy elements for the planets are made for each case. We find that in the case without disk PE, both Jupiter and Saturn will have a small residual core and comparable total amounts of heavy elements in the envelopes. We predict oxygen abundances enrichments in the same order as carbon, phosphorus and sulfur for both planets. Cometary Nitrogen abundances does not allow to easily reproduce Jupiter's nitrogen observations. In the case with disk PE, less core erosion is needed to reproduce the chemical composition of the atmospheres, so both planets will end up with possibly more massive residual cores, and higher total mass of heavy elements. It is also significantly easier to reproduce Jupiter's Nitrogen abundance. No single was disk was found to form both Jupiter and Saturn with all their constraints in the case without photoevaporation. No model was able to fit the constraints on Uranus & Neptune, hinting toward a more complicated formation mechanism for these planets. The predictions of these models should be compared to the upcoming Juno measurements to better understand the origins of the solar system giant planets. △ Less

Submitted 10 October, 2016; v1 submitted 11 September, 2016; originally announced September 2016.

Comments: Accepted for publication in MNRAS, 21 pages, 15 figures, feedback is welcome

A protosolar nebula origin for the ices agglomerated by Comet 67P/Churyumov-Gerasimenko

Authors: O. Mousis, J. I. Lunine, A. Luspay-Kuti, T. Guillot, B. Marty, M. Ali-Dib, P. Wurz, K. Altwegg, A. Bieler, M. Hässig, M. Rubin, P. Vernazza, J. H. Waite

Abstract: The nature of the icy material accreted by comets during their formation in the outer regions of the protosolar nebula is a major open question in planetary science. Some scenarios of comet formation predict that these bodies agglomerated from crystalline ices condensed in the protosolar nebula. Concurrently, alternative scenarios suggest that comets accreted amorphous ice originating from the int… ▽ More The nature of the icy material accreted by comets during their formation in the outer regions of the protosolar nebula is a major open question in planetary science. Some scenarios of comet formation predict that these bodies agglomerated from crystalline ices condensed in the protosolar nebula. Concurrently, alternative scenarios suggest that comets accreted amorphous ice originating from the interstellar cloud or from the very distant regions of the protosolar nebula. On the basis of existing laboratory and modeling data, we find that the N$_2$/CO and Ar/CO ratios measured in the coma of the Jupiter family comet 67P/Churyumov-Gerasimenko by the ROSINA instrument aboard the European Space Agency's Rosetta spacecraft match those predicted for gases trapped in clathrates. If these measurements are representative of the bulk N$_2$/CO and Ar/CO ratios in 67P/Churyumov-Gerasimenko, it implies that the ices accreted by the comet formed in the nebula and do not originate from the interstellar medium, supporting the idea that the building blocks of outer solar system bodies have been formed from clathrates and possibly from pure crystalline ices. Moreover, because 67P/Churyumov-Gerasimenko is impoverished in Ar and N$_2$, the volatile enrichments observed in Jupiter's atmosphere cannot be explained solely via the accretion of building blocks with similar compositions and require an additional delivery source. A potential source may be the accretion of gas from the nebula that has been progressively enriched in heavy elements due to photoevaporation. △ Less

Submitted 29 April, 2016; originally announced April 2016.

Comments: The Astrophysical Journal Letters, in press

Methane clathrates in the Solar System

Authors: Olivier Mousis, Eric Chassefière, Nils G. Holm, Alexis Bouquet, Jack Hunter Waite, Wolf Dietrich Geppert, Sylvain Picaud, Yuri Aikawa, Mohamad Ali-Dib, Jean-Luc Charlou, Philippe Rousselot

Abstract: We review the reservoirs of methane clathrates that may exist in the different bodies of the Solar System. Methane was formed in the interstellar medium prior to having been embedded in the protosolar nebula gas phase. This molecule was subsequently trapped in clathrates that formed from crystalline water ice during the cooling of the disk and incorporated in this form in the building blocks of co… ▽ More We review the reservoirs of methane clathrates that may exist in the different bodies of the Solar System. Methane was formed in the interstellar medium prior to having been embedded in the protosolar nebula gas phase. This molecule was subsequently trapped in clathrates that formed from crystalline water ice during the cooling of the disk and incorporated in this form in the building blocks of comets, icy bodies, and giant planets. Methane clathrates may play an important role in the evolution of planetary atmospheres. On Earth, the production of methane in clathrates is essentially biological, and these compounds are mostly found in permafrost regions or in the sediments of continental shelves. On Mars, methane would more likely derive from hydrothermal reactions with olivine-rich material. If they do exist, martian methane clathrates would be stable only at depth in the cryosphere and sporadically release some methane into the atmosphere via mechanisms that remain to be determined. △ Less

Submitted 26 October, 2015; originally announced October 2015.

Journal ref: Astrobiology, Volume 15, Number 4, 2015

Nebular dead zone effects on the D/H ratio in chondrites and comets

Authors: Mohamad Ali-Dib, Rebecca G. Martin, Jean-Marc Petit, Olivier Mousis, Pierre Vernazza, Jonathan I. Lunine

Abstract: Comets and chondrites show non-monotonic behaviour of their Deuterium to Hydrogen (D/H) ratio as a function of their formation location from the Sun. This is difficult to explain with a classical protoplanetary disk model that has a decreasing temperature structure with radius from the Sun. We want to understand if a protoplanetary disc with a dead zone, a region of zero or low turbulence, can e… ▽ More Comets and chondrites show non-monotonic behaviour of their Deuterium to Hydrogen (D/H) ratio as a function of their formation location from the Sun. This is difficult to explain with a classical protoplanetary disk model that has a decreasing temperature structure with radius from the Sun. We want to understand if a protoplanetary disc with a dead zone, a region of zero or low turbulence, can explain the measured D/H values in comets and chondrites. We use time snapshots of a vertically layered disk model with turbulent surface layers and a dead zone at the midplane. The disc has a non-monotonic temperature structure due to increased heating from self-gravity in the outer parts of the dead zone. We couple this to a D/H ratio evolution model in order to quantify the effect of such thermal profiles on D/H enrichment in the nebula. We find that the local temperature peak in the disk can explain the diversity in the D/H ratios of different chondritic families. This disk temperature profile leads to a non-monotonic D/H enrichment evolution, allowing these families to acquire their different D/H values while forming in close proximity. The formation order we infer for these families is compatible with that inferred from their water abundances. However, we find that even for very young disks, the thermal profile reversal is too close to the Sun to be relevant for comets. △ Less

Submitted 2 August, 2015; originally announced August 2015.

Comments: 9 pages. Accepted for publication in A&A

Journal ref: A&A 583, A58 (2015)

The composition of the protosolar disk and the formation conditions for comets

Authors: K. Willacy, C. Alexander, M. Ali-Dib, C. Ceccarelli, S. B. Charnley, M. Doronin, Y. Ellinger, P. Gast, E. Gibb, S. N. Milam, O. Mousis, E. Pauzat, C. Tornow, E. S. Wirstrom, E Zicler

Abstract: Conditions in the protosolar nebula have left their mark in the composition of cometary volatiles, thought to be some of the most pristine material in the solar system. Cometary compositions represent the end point of processing that began in the parent molecular cloud core and continued through the collapse of that core to form the protosun and the solar nebula, and finally during the evolution o… ▽ More Conditions in the protosolar nebula have left their mark in the composition of cometary volatiles, thought to be some of the most pristine material in the solar system. Cometary compositions represent the end point of processing that began in the parent molecular cloud core and continued through the collapse of that core to form the protosun and the solar nebula, and finally during the evolution of the solar nebula itself as the cometary bodies were accreting. Disentangling the effects of the various epochs on the final composition of a comet is complicated. But comets are not the only source of information about the solar nebula. Protostellar disks around young stars similar to the protosun provide a way of investigating the evolution of disks similar to the solar nebula while they are in the process of evolving to form their own solar systems. In this way we can learn about the physical and chemical conditions under which comets formed, and about the types of dynamical processing that shaped the solar system we see today. This paper summarizes some recent contributions to our understanding of both cometary volatiles and the composition, structure and evolution of protostellar disks. △ Less

Submitted 8 July, 2015; originally announced July 2015.

Comments: To appear in Space Science Reviews. The final publication is available at Springer via http://dx.doi.org/10.1007/s11214-015-0167-6

Detection and characterization of Io's atmosphere from high-resolution 4-μm spectroscopy

Authors: E. Lellouch, M. Ali-Dib, K. -L. Jessup, A. Smette, H. -U. Käufl, F. Marchis

Abstract: We report on high-resolution and spatially-resolved spectra of Io in the 4.0 μm region, recorded with the VLT/CRIRES instrument in 2008 and 2010, which provide the first detection of the ν1 + ν3 band of SO2 in Io's atmosphere. Data are analyzed to constrain the latitudinal, longitudinal, and diurnal distribution of Io's SO2 atmosphere as well as its characteristic temperature. equatorial SO2 colum… ▽ More We report on high-resolution and spatially-resolved spectra of Io in the 4.0 μm region, recorded with the VLT/CRIRES instrument in 2008 and 2010, which provide the first detection of the ν1 + ν3 band of SO2 in Io's atmosphere. Data are analyzed to constrain the latitudinal, longitudinal, and diurnal distribution of Io's SO2 atmosphere as well as its characteristic temperature. equatorial SO2 column densities clearly show longitudinal asymmetry, but with a maximum of around 1.5e17 cm-2 at central meridian longitude L = 200-220 and a minimum of around 3e16 cm-2 at L = 285-300, the longitudinal pattern somewhat differs from earlier inferences from Ly α and thermal IR measurements. Within the accuracy of the measurements, no evolution of the atmospheric density from mid-2008 to mid-2010 can be distinguished. The decrease of the SO2 column density towards high latitude is apparent, and the typical latitudinal extent of the atmosphere found to be (+-) 40° at half-maximum. The data show moderate diurnal variations of the equatorial atmosphere, which is evidence for a partially sublimation-supported atmospheric component. Compared to local noon, factor of 2 lower densities are observed around 40° before and 80° after noon. Best-fit gas temperatures range from 150 to 220 K, with a weighted mean value of 170 (+-) 20 K, which should represent the column-weighted mean kinetic temperature of Io's atmosphere. Finally, although the data include clear thermal emission due to Pillan (in outburst in July 2008) and Loki, no detectable enhancements in the SO2 atmosphere above these volcanic regions are found, with an upper limit of 4e16 cm-2 at Pillan and 1e17 cm-2 at Loki. △ Less

Submitted 19 February, 2015; originally announced February 2015.

Comments: Accepted for publication in Icarus

New insights on Saturn's formation from its nitrogen isotopic composition

Authors: Olivier Mousis, Jonathan I. Lunine, Leigh N. Fletcher, Kathleen E. Mandt, Mohamad Ali-Dib, Daniel Gautier, Sushil Atreya

Abstract: The recent derivation of a lower limit for the $^{14}$N/$^{15}$N ratio in Saturn's ammonia, which is found to be consistent with the Jovian value, prompted us to revise models of Saturn's formation using as constraints the supersolar abundances of heavy elements measured in its atmosphere. Here we find that it is possible to account for both Saturn's chemical and isotopic compositions if one assum… ▽ More The recent derivation of a lower limit for the $^{14}$N/$^{15}$N ratio in Saturn's ammonia, which is found to be consistent with the Jovian value, prompted us to revise models of Saturn's formation using as constraints the supersolar abundances of heavy elements measured in its atmosphere. Here we find that it is possible to account for both Saturn's chemical and isotopic compositions if one assumes the formation of its building blocks at $\sim$45 K in the protosolar nebula, provided that the O abundance was $\sim$2.6 times protosolar in its feeding zone. To do so, we used a statistical thermodynamic model to investigate the composition of the clathrate phase that formed during the cooling of the protosolar nebula and from which the building blocks of Saturn were agglomerated. We find that Saturn's O/H is at least $\sim$34.9 times protosolar and that the corresponding mass of heavy elements ($\sim$43.1 \Mearth) is within the range predicted by semi-convective interior models. △ Less

Submitted 20 October, 2014; originally announced October 2014.

Comments: Accepted for publication in Astrophysical Journal Letters

The measured compositions of Uranus and Neptune from their formation on the CO iceline

Authors: M. Ali-Dib, O. Mousis, J. -M. Petit, J. I. Lunine

Abstract: The formation mechanisms of the ice giants Uranus and Neptune, and the origin of their elemental and isotopic compositions, have long been debated. The density of solids in the outer protosolar nebula is too low to explain their formation, and spectroscopic observations show that both planets are highly enriched in carbon, very poor in nitrogen, and the ices from which they originally formed might… ▽ More The formation mechanisms of the ice giants Uranus and Neptune, and the origin of their elemental and isotopic compositions, have long been debated. The density of solids in the outer protosolar nebula is too low to explain their formation, and spectroscopic observations show that both planets are highly enriched in carbon, very poor in nitrogen, and the ices from which they originally formed might had deuterium-to-hydrogen ratios lower than the predicted cometary value, unexplained properties observed in no other planets. Here we show that all these properties can be explained naturally if Uranus and Neptune both formed at the carbon monoxide iceline. Due to the diffusive redistribution of vapors, this outer region of the protosolar nebula intrinsically has enough surface density to form both planets from carbon-rich solids but nitrogen-depleted gas, in abundances consistent with their observed values. Water rich interiors originating mostly from transformed CO ices reconcile the D/H value of Uranus and Neptune's building blocks with the cometary value. Finally, Our scenario generalizes a well known hypothesis that Jupiter formed on an iceline (water snowline) for the two ice giants, and might be a first step towards generalizing this mechanism for other giant planets. △ Less

Submitted 9 July, 2014; originally announced July 2014.

Comments: The Astrophysical Journal (in press), 8 pages, 5 figures

Scientific rationale of Saturn's in situ exploration

Authors: O. Mousis, L. N. Fletcher, J. -P. Lebreton, P. Wurz, T. Cavalié, A. Coustenis, R. Courtin, D. Gautier, R. Helled, P. G. J. Irwin, A. D. Morse, N. Nettelmann, B. Marty, P. Rousselot, O. Venot, D. H. Atkinson, J. H. Waite, K. R. Reh, A. Simon-Miller, S. Atreya, N. André, M. Blanc, I. A. Daglis, G. Fischer, W. D. Geppert , et al. (25 additional authors not shown)

Abstract: Remote sensing observations meet some limitations when used to study the bulk atmospheric composition of the giant planets of our solar system. A remarkable example of the superiority of in situ probe measurements is illustrated by the exploration of Jupiter, where key measurements such as the determination of the noble gases abundances and the precise measurement of the helium mixing ratio have o… ▽ More Remote sensing observations meet some limitations when used to study the bulk atmospheric composition of the giant planets of our solar system. A remarkable example of the superiority of in situ probe measurements is illustrated by the exploration of Jupiter, where key measurements such as the determination of the noble gases abundances and the precise measurement of the helium mixing ratio have only been made available through in situ measurements by the Galileo probe. This paper describes the main scientific goals to be addressed by the future in situ exploration of Saturn placing the Galileo probe exploration of Jupiter in a broader context and before the future probe exploration of the more remote ice giants. In situ exploration of Saturn's atmosphere addresses two broad themes that are discussed throughout this paper: first, the formation history of our solar system and second, the processes at play in planetary atmospheres. In this context, we detail the reasons why measurements of Saturn's bulk elemental and isotopic composition would place important constraints on the volatile reservoirs in the protosolar nebula. We also show that the in situ measurement of CO (or any other disequilibrium species that is depleted by reaction with water) in Saturn's upper troposphere would constrain its bulk O/H ratio. We highlight the key measurements required to distinguish competing theories to shed light on giant planet formation as a common process in planetary systems with potential applications to most extrasolar systems. In situ measurements of Saturn's stratospheric and tropospheric dynamics, chemistry and cloud-forming processes will provide access to phenomena unreachable to remote sensing studies. Different mission architectures are envisaged, which would benefit from strong international collaborations. △ Less

Submitted 18 April, 2014; originally announced April 2014.

Comments: Submitted to Planetary and Space Science

Carbon-rich planet formation in a solar composition disk

Authors: Mohamad Ali-Dib, Olivier Mousis, Jean-Marc Petit, Jonathan I. Lunine

Abstract: The C--to--O ratio is a crucial determinant of the chemical properties of planets. The recent observation of WASP 12b, a giant planet with a C/O value larger than that estimated for its host star, poses a conundrum for understanding the origin of this elemental ratio in any given planetary system. In this paper, we propose a mechanism for enhancing the value of C/O in the disk through the transpor… ▽ More The C--to--O ratio is a crucial determinant of the chemical properties of planets. The recent observation of WASP 12b, a giant planet with a C/O value larger than that estimated for its host star, poses a conundrum for understanding the origin of this elemental ratio in any given planetary system. In this paper, we propose a mechanism for enhancing the value of C/O in the disk through the transport and distribution of volatiles. We construct a model that computes the abundances of major C and O bearing volatiles under the influence of gas drag, sublimation, vapor diffusion, condensation and coagulation in a multi--iceline 1+1D protoplanetary disk. We find a gradual depletion in water and carbon monoxide vapors inside the water's iceline with carbon monoxide depleting slower than water. This effect increases the gaseous C/O and decreases the C/H ratio in this region to values similar to those found in WASP 12b's day side atmosphere. Giant planets whose envelopes were accreted inside the water's iceline should then display C/O values larger than those of their parent stars, making them members of the class of so-called ``carbon-rich planets''. △ Less

Submitted 20 February, 2014; originally announced February 2014.

Comments: 8 pages, 4 figures, accepted for publication ApJ

Influence of the C/O ratio on titanium and vanadium oxides in protoplanetary disks

Authors: M. Ali-Dib, O. Mousis, G. S. Pekmezci, J. I. Lunine, N. Madhusudhan, J. -M. Petit

Abstract: Context. The observation of carbon-rich disks have motivated several studies questioning the influence of the C/O ratio on their gas phase composition in order to establish the connection between the metallicity of hot-Jupiters and that of their parent stars. Aims. We to propose a method that allows the characterization of the adopted C/O ratio in protoplanetary disks independently from the dete… ▽ More Context. The observation of carbon-rich disks have motivated several studies questioning the influence of the C/O ratio on their gas phase composition in order to establish the connection between the metallicity of hot-Jupiters and that of their parent stars. Aims. We to propose a method that allows the characterization of the adopted C/O ratio in protoplanetary disks independently from the determination of the host star composition. Titanium and vanadium chemistries are investigated because they are strong optical absorbers and also because their oxides are known to be sensitive to the C/O ratio in some exoplanet atmospheres. Methods. We use a commercial package based on the Gibbs energy minimization technique to compute the titanium and vanadium equilibrium chemistries in protoplanetary disks for C/O ratios ranging from 0.05 to 10. Our calculations are performed for pressures ranging from 1e-6 to 1e-2 bar, and for temperatures ranging from 50 to 2000 K. Results. We find that the vanadium nitride/vanadium oxide and titanium hydride/titanium oxide gas phase ratios strongly depend on the C/O ratio in the hot parts of disks (T > 1000 K). Our calculations suggest that, in these regions, these ratios can be used as tracers of the C/O value in protoplanetary disks. △ Less

Submitted 12 November, 2013; originally announced November 2013.

Comments: Accepted for publication in A&A