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.

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.

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/.

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.

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.

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 .

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.

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\'e 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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

We report on high-resolution and spatially-resolved spectra of Io in the 4.0 {\mu}m region, recorded with the VLT/CRIRES instrument in 2008 and 2010, which provide the first detection of the {\nu}1 + {\nu}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 {\alpha} 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{\deg} 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{\deg} before and 80{\deg} 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.

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.

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.

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.

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''.

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.