In this review, the formation of the Oort Cloud is illuminated from several aspects. One is the history of the subject with an outline of the fundamental discoveries by Öpik [1], Oort [2] and Hills [3]. It is argued that the basic reason for judging Oort as the real discoverer is that he had access to observational data in the form of original orbits of long period comets. Further landmarks are identified, like the exploration of the role of the Galactic tide in the supply of observable comets by Heisler and Tremaine [4], the clarification of a synergy between tide and stars as the reason for a continued, efficient supply by Rickman et al. [5], the discovery by Kaib and Quinn [6] that inner core comets become observable due to planetary perturbations, disguised as new comets, and the demonstration of how Oort Cloud formation may work in the realm of the Nice Model by Brasser and Morbidelli [7]. It is finally argued that an important Galactic sculpting has occurred since the primordial Oort Cloud was formed by means of global shake-up events resulting from impulses imparted to the Sun by external perturbers like massive stars or Giant Molecular Clouds, and that this may be the real reason for the survival of an outer halo that reveals the existence of the Oort Cloud through the Oort spike.

We investigate the survivability of solar system-like planetary systems during close encounters in stellar associations using a suite of 1980 N-body simulations. Each system is based on one of the possible five-planet resonant configurations proposed to represent the initial solar system architecture and is systematically scaled in both planetary mass and orbital compactness to explore the parameter space of observed exoplanetary architectures. Simulations explore a range of stellar encounter scenarios drawn from four distinct cluster environments. Our results show that system survival depends critically on the interplay between planetary mass and orbital scale: compact configurations are more resistant to external perturbations, while increased planetary mass improves resilience only up to a threshold, beyond which internal instabilities dominate. No system whose planets are twice as massive as the ones in the solar system survives stellar encounters. Systems that are at least an order of magnitude more compact than the solar system remain stable under typical encounter conditions. These findings place strong constraints on the initial architectures of planetary systems that can endure stellar-dense birth environments.

This paper presents a study of the dynamical mechanisms behind the extraction of Jupiter–Saturn planetesimals into detached orbits belonging to the Oort Cloud or the sednoid population. In an embedded cluster, these are the global cluster tide and gravitational impulses acquired due to stellar encounters. Our study builds on a large set of numerical simulations of Oort Cloud and sednoid formation, whose statistical results are reported in an accompanying paper. The material is extensive, amounting to detailed data on the orbital evolutions of 288000 planetesimals. We have developed two graphical tools for our work: (1) diagrams showing the extraction history of sednoids and Oort Cloud comets for each of 72 solar template stars; and (2) time profiles of the evolution of semi-major axis and perihelion distance for individual planetesimals. These help us identify the main triggers behind individual extractions due to intrinsic differences between the tidal and stellar flyby triggering mechanisms. The time scale of an encounter-based extraction never exceeds the duration of the encounter, while the tide operates on two time scales — the orbital periods of the planetesimal around the star and the star around the cluster centre. These periods are very different for sednoids but often similar for Oort Cloud comets. We note the absence of evidence for synergy between the tide and the encounters or between subsequent encounters, likely due to the fact that each effect depends strongly on the geometrical circumstances, and those are rarely repeatable. The balance between tide and encounters as dominating trigger varies between different cluster types. For clusters which are strongly concentrated with the primordial gas remaining for a few Myr, we discover a previously unnoticed phenomenon, which we call tidal purging. In this case, nearly all Oort Cloud comets with semi-major axes exceeding 7000astronomical units are ejected from the solar system on hyperbolic orbits due to the non-conservative nature of the cluster tide. In a separate study of the end states reached by all planetesimals, we find that an overwhelming majority undergo hyperbolic ejections, remain in orbits similar to the initial ones, or get physically removed from the solar system due to collisions with a planet or the Sun.

In this paper we report on massive computer simulations aimed to clarify the formation of an Oort Cloud and the adjoining sednoid population during the earliest times of the solar system due to the perturbing effects of the Sun’s birth aggregate, for which we use different embedded cluster models. We define the sednoid population to have semi-major axes between 100 and 1600 astronomical units, and it can be seen as the innermost extension of the Oort Cloud. By also restricting the orbits to perihelion distances exceeding 60 astronomical units, the sednoid population can be regarded as dynamically inert with respect to both planetary and Galactic influences. It includes the objects (90377) Sedna, 2012 VP113 and (541132) Leleākūhonua. We start by creating 12 cluster models, using the NBODY6++GPU code, with three random realisations of four physically different models. In each of these we select six random solar template stars and equip those with a planetary system composed of Jupiter, Saturn and 4000 massless planetesimals representing a population of stray planetesimals in the Jupiter-Saturn zone with a total mass assumed to be about 20 Earth masses. From dynamical simulations of all these 72 variants of the early solar system including the cluster perturbations, we derive the total numbers and orbital characteristics of the sednoid and Oort Cloud populations, including the size distributions resulting from the limiting effects of solar nebula gas drag on extraction efficiency. We find that the predicted number of sednoids varies dramatically depending on the orbital evolution of the solar template star in the cluster. In general, numbers consistent with the existence of the observed sednoids is an expected outcome only for a birth aggregate with a strong central condensation and a relatively long lifetime of the residual gas, whereas other cluster models tend to fail to explain this. In the most successful, concentrated cluster models, we predict the existence of an as yet unseen population of at least 108 km-sized sednoids. However, we note that these cluster models produce an overheated inclination distribution at variance with the observed inclinations. Hence, we conclude that our study points to a possible, severe problem of the embedded cluster scenario in explaining the origin of the sednoids. The Oort Cloud predicted to originate in local, Jupiter-Saturn planetesimals by our simulations amounts to ∼10% of the estimated, current size of the entire Oort Cloud. This shows that such comets, likely characterised by Earth-like D/H ratios, may be a minority though a significant one. In an accompanying paper we discuss the extraction dynamics in some detail, based on the same simulations as here.

ABSTRACT We report the first simulations of planetary system dynamics as affected by an embedded cluster environment. Such environments are generally believed to be relevant for the large majority of newborn stars of solar type. Moreover, our cluster model is more realistic than in previous work. We focus on a giant planet system with five members, which represents a likely precursor of our solar system. Our main result is that the perturbing effects of close encounters with cluster stars trigger dynamical chaos leading to breakdown of the system with a significant probability, especially if the natal gas discs are short-lived and the clusters are highly concentrated. When breakdown occurs, all planets except Jupiter suffer a large risk of being ejected from the system or extracted into distant orbits with semimajor axes of hundreds or thousands of astronomical units. This is consistent with recent estimates of a large abundance of low-mass, free-floating planets. We demonstrate a possibility for Jupiter and Saturn to evolve into hot Jupiter orbits by tidal circularization during the chaotic evolution. Even so, the low occurrence rate of this outcome indicates that the real hot Jupiters in general have an origin unrelated to dynamical evolution in birth clusters.

Ariel, the Atmospheric Remote-sensing Infrared Exoplanet Large-survey, was adopted as the fourth medium-class mission in ESA's Cosmic Vision programme to be launched in 2029. During its 4-year mission, Ariel will study what exoplanets are made of, how they formed and how they evolve, by surveying a diverse sample of about 1000 extrasolar planets, simultaneously in visible and infrared wavelengths. It is the first mission dedicated to measuring the chemical composition and thermal structures of hundreds of transiting exoplanets, enabling planetary science far beyond the boundaries of the Solar System. The payload consists of an off-axis Cassegrain telescope (primary mirror 1100 mm x 730 mm ellipse) and two separate instruments (FGS and AIRS) covering simultaneously 0.5-7.8 micron spectral range. The satellite is best placed into an L2 orbit to maximise the thermal stability and the field of regard. The payload module is passively cooled via a series of V-Groove radiators; the detectors for the AIRS are the only items that require active cooling via an active Ne JT cooler. The Ariel payload is developed by a consortium of more than 50 institutes from 16 ESA countries, which include the UK, France, Italy, Belgium, Poland, Spain, Austria, Denmark, Ireland, Portugal, Czech Republic, Hungary, the Netherlands, Sweden, Norway, Estonia, and a NASA contribution.

Aims. We want to find the conditions under which planetary close encounters transform the orbits of small Solar System bodies from direct to retrograde, and vice versa. Methods. We derive analytical constraints on the orbital elements of the small body that allow direct-retrograde transitions at close encounters. We check the validity of the analytical constraints with numerical integrations of close encounters in the restricted, circular, three-dimensional three-body problem. Results. For bound orbits, inclination flips at close encounters are possible only for values of the Tisserand parameter, computed with respect to the planet actually encountered, which are within certain limits. We give an analytical expression for the probability per revolution of this transition, as function of the orbital parameters. We show how to identify, among the known asteroids and comets on direct orbits, those that can flip to retrograde motion due to an encounter with an outer planet. Conclusions. Inclination flips at planetary close encounters can be quantitatively characterized with the analytical theory of close encounters.

ABSTRACT Simulations of angular momentum transfer from a surrounding star cluster to distant planetesimals orbiting around a cluster member have been used since more than two decades to study the formation of small body populations in the outskirts of the Solar system (Oort Cloud and sednoids). We present a new model for these interactions, for the first time combining two features of earlier works: (1) a self-consistent treatment of cluster evolution based on N-body simulations, and (2) a treatment of circumstellar dynamics as resulting from a combination of a smooth tidal field representing the whole cluster and close encounters by individual cluster members. The model is expected to be both flexible, accurate, and efficient in terms of CPU time and hence a suitable tool when simulating large or long-lived clusters or when many independent runs are needed for statistical significance. We describe the model in detail and give examples of its outputs. These are of relevance not only for the perihelion extraction of small bodies but also for the stability properties of young or nascent planetary systems or protoplanetary discs.

Between 2014 December 31 and 2015 March 17, the OSIRIS cameras on Rosetta documented the growth of a 140 -m wide and 0.5 -m deep depression in the Hapi region on Comet 67P/Churyumov-Gerasimenko. This shallow pit is one of several that later formed elsewhere on the comet, all in smooth terrain that primarily is the result of airfall of coma particles. We have compiled observations of this region in Hapi by the microwave instrument MIRO on Rosetta, acquired during October and November 2014. We use thermophysical and radiative transfer models in order to reproduce the MIRO observations. This allows us to place constraints on the thermal inertia, diffusivity, chemical composition, stratification, extinction coefficients, and scattering properties of the surface material, and how they evolved during the months prior to pit formation. The results are placed in context through long-term comet nucleus evolution modelling. We propose that (1) MIRO observes signatures that are consistent with a solid-state greenhouse effect in airfall material; (2) CO2 ice is sufficiently close to the surface to have a measurable effect on MIRO antenna temperatures, and likely is responsible for the pit formation in Hapi observed by OSIRIS; (3) the pressure at the CO2 sublimation front is sufficiently strong to expel dust and water ice outwards, and to compress comet material inwards, thereby causing the near-surface compaction observed by CONSERT, SESAME, and groundbased radar, manifested as the 'consolidated terrain' texture observed by OSIRIS.

The goal of the Ariel space mission is to observe a large and diversified population of transiting planets around a range of host star types to collect information on their atmospheric composition. The planetary bulk and atmospheric compositions bear the marks of the way the planets formed: Ariel's observations will therefore provide an unprecedented wealth of data to advance our understanding of planet formation in our Galaxy. A number of environmental and evolutionary factors, however, can affect the final atmospheric composition. Here we provide a concise overview of which factors and effects of the star and planet formation processes can shape the atmospheric compositions that will be observed by Ariel, and highlight how Ariel's characteristics make this mission optimally suited to address this very complex problem.

We report on N-body simulations of embedded stellar clusters over a time span of 10 Myr from inception until several crossing times after gas expulsion. We focus our attention on the survivability of bound stellar clusters following gas expulsion. The final values of the bound fraction are significantly larger in case the initial structure is relatively open with a Plummer radius of approximate to 1 pc than for a much stronger degree of concentration. However, only one of 12 simulated clusters shows any promise of leaving behind a long-lived, bound system. We tentatively conclude that the predictions of our model are in rough agreement with the observed "infant mortality" of newborn stellar aggregates. We also find that relaxation effects due to close encounters, while present, are significantly curtailed by the rapid expansion of the clusters.

Ariel, the Atmospheric Remote-sensing Infrared Exoplanet Large-survey, was adopted as the fourth medium-class mission in ESA's Cosmic Vision programme to be launched in 2029. During its 4-year mission, Ariel will study what exoplanets are made of, how they formed and how they evolve, by surveying a diverse sample of about 1000 extrasolar planets, simultaneously in visible and infrared wavelengths. It is the first mission dedicated to measuring the chemical composition and thermal structures of hundreds of transiting exoplanets, enabling planetary science far beyond the boundaries of the Solar System. The payload consists of an off-axis Cassegrain telescope (primary mirror 1100 mm x 730 mm ellipse) and two separate instruments (FGS and AIRS) covering simultaneously 0.5-7.8 micron spectral range. The satellite is best placed into an L2 orbit to maximise the thermal stability and the field of regard. The payload module is passively cooled via a series of V-Groove radiators; the detectors for the AIRS are the only items that require active cooling via an active Ne JT cooler. The Ariel payload is developed by a consortium of more than 50 institutes from 16 ESA countries, which include the UK, France, Italy, Belgium, Poland, Spain, Austria, Denmark, Ireland, Portugal, Czech Republic, Hungary, the Netherlands, Sweden, Norway, Estonia, and a NASA contribution.

Bilobate comets—small icy bodies with two distinct lobes—are a common configuration among comets, but the factors shaping these bodies are largely unknown. Cometary nuclei, the solid centres of comets, erode by ice sublimation when they are sufficiently close to the Sun, but the importance of a comet's internal structure on its erosion is unclear. Here we present three-dimensional analyses of images from the Rosetta mission to illuminate the process that shaped the Jupiter-family bilobate comet 67P/Churyumov–Gerasimenko over billions of years. We show that the comet's surface and interior exhibit shear-fracture and fault networks, on spatial scales of tens to hundreds of metres. Fractures propagate up to 500 m below the surface through a mechanically homogeneous material. Through fracture network analysis and stress modelling, we show that shear deformation generates fracture networks that control mechanical surface erosion, particularly in the strongly marked neck trough of 67P/Churyumov–Gerasimenko, exposing its interior. We conclude that shear deformation shapes and structures the surface and interior of bilobate comets, particularly in the outer Solar System where water ice sublimation is negligible. The shape and internal structure of bilobate comet 67P is controlled by shear deformation inducing mechanically driven erosion along shear fracture networks, according to a 3D analysis of images from the Rosetta mission.

Aims. We investigate the surface distribution of the source regions of dust jets on comet 67P/Churyumov-Gerasimenko as a function of time. Methods. The dust jet source regions were traced by the comprehensive imaging data set provided by the OSIRIS scientific camera. Results. We show in detail how the projected footpoints of the dust jets and hence the outgassing zone would move in consonance with the sunlit belt. Furthermore, a number of source regions characterized by repeated jet activity might be the result of local topographical variations or compositional heterogeneities. Conclusions. The spatial and temporal variations in source regions of the dust jets are influenced significantly by the seasonal effect. The strong dependence on the solar zenith angle and local time could be related to the gas sublimation process driven by solar insolation on a surface layer of low thermal inertia.

We directly measured twenty overhanging cliffs on the surface of comet 67P/Churyumov-Gerasimenko extracted from the latest shape model and estimated the minimum tensile strengths needed to support them against collapse under the comet's gravity. We find extremely low strengths of around 1 Pa or less (1 to 5 Pa, when scaled to a metre length). The presence of eroded material at the base of most overhangs, as well as the observed collapse of two features and the implied previous collapse of another, suggests that they are prone to failure and that the true material strengths are close to these lower limits (although we only consider static stresses and not dynamic stress from, for example, cometary activity). Thus, a tensile strength of a few pascals is a good approximation for the tensile strength of the 67P nucleus material, which is in agreement with previous work. We find no particular trends in overhang properties either with size over the similar to 10-100 m range studied here or location on the nucleus. There are no obvious differences, in terms of strength, height or evidence of collapse, between the populations of overhangs on the two cometary lobes, suggesting that 67P is relatively homogenous in terms of tensile strength. Low material strengths are supportive of cometary formation as a primordial rubble pile or by collisional fragmentation of a small body (tens of km).