Planet formation models: the interplay with the planetesimal disc

Context. According to the sequential accretion model (or core-nucleated accretion model), giant planet formation is based first on the formation of a solid core which, when massive enough, can gravitationally bind gas from the nebula to form the envelope. The most critical part of the model is the formation time of the core: to trigger the accretion of gas, the core has to grow up to several Earth masses before the gas component of the protoplanetary disc dissipates. Aims: We calculate planetary formation models including a detailed description of the dynamics of the planetesimal disc, taking into account both gas drag and excitation of forming planets. Methods: We computed the formation of planets, considering the oligarchic regime for the growth of the solid core. Embryos growing in the disc stir their neighbour planetesimals, exciting their relative velocities, which makes accretion more difficult. Here we introduce a more realistic treatment for the evolution of planetesimals' relative velocities, which directly impact on the formation timescale. For this, we computed the excitation state of planetesimals, as a result of stirring by forming planets, and gas-solid interactions. Results: We find that the formation of giant planets is favoured by the accretion of small planetesimals, as their random velocities are more easily damped by the gas drag of the nebula. Moreover, the capture radius of a protoplanet with a (tiny) envelope is also larger for small planetesimals. However, planets migrate as a result of disc-planet angular momentum exchange, with important consequences for their survival: due to the slow growth of a protoplanet in the oligarchic regime, rapid inward type I migration has important implications on intermediate-mass planets that have not yet started their runaway accretion phase of gas. Most of these planets are lost in the central star. Surviving planets have masses either below 10 M⊕ or above several Jupiter masses. Conclusions: To form giant planets before the dissipation of the disc, small planetesimals (~0.1 km) have to be the major contributors of the solid accretion process. However, the combination of oligarchic growth and fast inward migration leads to the absence of intermediate-mass planets. Other processes must therefore be at work to explain the population of extrasolar planets that are presently known.

Gas composition of the main volatile elements in protoplanetary discs and its implication for planet formation

Context. Direct observations of gaseous exoplanets reveal that their gas envelope has a higher C/O ratio than that of the host star (e.g., Wasp 12-b). This has been explained by considering that the gas phase of the disc could be inhomogeneous, exceeding the stellar C/O ratio in regions where these planets formed; but few studies have considered the drift of the gas and planet migration. Aims. We aim to derive the gas composition in planets through planet formation to evaluate if the formation of giant planets with an enriched C/O ratio is possible. The study focusses on the effects of different processes on the C/O ratio, such as the disc evolution, the drift of gas, and planet migration. Methods. We used our previous models for computing the chemical composition, together with a planet formation model, to which we added the composition and drift of the gas phase of the disc, which is composed of the main volatile species H2O, CO, CO2, NH3, N2, CH3OH, CH4, and H2S, H2 and He. The study focusses on the region where ice lines are present and influence the C/O ratio of the planets. Results. Modelling shows that the condensation of volatile species as a function of radial distance allows for C/O enrichment in specific parts of the protoplanetary disc of up to four times the solar value. This leads to the formation of planets that can be enriched in C/O in their envelope up to three times the solar value. Planet migration, gas phase evolution and disc irradiation enables the evolution of the initial C/O ratio that decreases in the outer part of the disc and increases in the inner part of the disc. The total C/O ratio of the planets is governed by the contribution of ices accreted, suggesting that high C/O ratios measured in planetary atmospheres are indicative of a lack of exchange of material between the core of a planet and its envelope or an observational bias. It also suggests that the observed C/O ratio is not representative of the total C/O ratio of the planet.

The Doppler shadow of WASP-3b. A tomographic analysis of Rossiter-McLaughlin observations

Context. Hot-Jupiter planets must form at large separations from their host stars where the temperatures are cool enough for their cores to condense. They then migrate inwards to their current observed orbital separations. Different theories of how this migration occurs lead to varying distributions of orbital eccentricity and the alignment between the rotation axis of the star and the orbital axis of the planet. Aims: The spin-orbit alignment of a transiting system is revealed via the Rossiter-McLaughlin effect, which is the anomaly present in the radial velocity measurements of the rotating star during transit due to the planet blocking some of the starlight. In this paper we aim to measure the spin-orbit alignment of the WASP-3 system via a new way of analysing the Rossiter-McLaughlin observations. Methods: We apply a new tomographic method for analysing the time variable asymmetry of stellar line profiles caused by the Rossiter-McLaughlin effect. This new method eliminates the systematic error inherent in previous methods used to analyse the effect. Results: We find a value for the projected stellar spin rate of v sin i = 13.9 ± 0.03 km s-1 which is in agreement with previous measurements but has a much higher precision. The system is found to be well aligned, with ? = 5-5+6° which favours an evolutionary history for WASP-3b involving migration through tidal interactions with a protoplanetary disc. From comparison with isochrones we put an upper limit on the age of the star of 2 Gyr.

Estudo da dinâmica de captura em discos proto-planetários

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Planet formation in a triple stellar system: implications of the third star's orbital inclination

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Theoretical models of planetary system formation: mass vs. semi-major axis

Context. Planet formation models have been developed during the past years to try to reproduce what has been observed of both the solar system and the extrasolar planets. Some of these models have partially succeeded, but they focus on massive planets and, for the sake of simplicity, exclude planets belonging to planetary systems. However, more and more planets are now found in planetary systems. This tendency, which is a result of radial velocity, transit, and direct imaging surveys, seems to be even more pronounced for low-mass planets. These new observations require improving planet formation models, including new physics, and considering the formation of systems. Aims: In a recent series of papers, we have presented some improvements in the physics of our models, focussing in particular on the internal structure of forming planets, and on the computation of the excitation state of planetesimals and their resulting accretion rate. In this paper, we focus on the concurrent effect of the formation of more than one planet in the same protoplanetary disc and show the effect, in terms of architecture and composition of this multiplicity. Methods: We used an N-body calculation including collision detection to compute the orbital evolution of a planetary system. Moreover, we describe the effect of competition for accretion of gas and solids, as well as the effect of gravitational interactions between planets. Results: We show that the masses and semi-major axes of planets are modified by both the effect of competition and gravitational interactions. We also present the effect of the assumed number of forming planets in the same system (a free parameter of the model), as well as the effect of the inclination and eccentricity damping. We find that the fraction of ejected planets increases from nearly 0 to 8% as we change the number of embryos we seed the system with from 2 to 20 planetary embryos. Moreover, our calculations show that, when considering planets more massive than ~5 M⊕, simulations with 10 or 20 planetary embryos statistically give the same results in terms of mass function and period distribution.

From stellar nebula to planetesimals

Context. Solar and extrasolar comets and extrasolar planets are the subject of numerous studies in order to determine their chemical composition and internal structure. In the case of planetesimals, their compositions are important as they govern in part the composition of future planets. Aims. The present works aims at determining the chemical composition of icy planetesimals, believed to be similar to present day comets, formed in stellar systems of solar chemical composition. The main objective of this work is to provide valuable theoretical data on chemical composition for models of planetesimals and comets, and models of planet formation and evolution. Methods. We have developed a model that calculates the composition of ices formed during the cooling of the stellar nebula. Coupled with a model of refractory element formation, it allows us to determine the chemical composition and mass ratio of ices to rocks in icy planetesimals throughout in the protoplanetary disc. Results. We provide relationships for ice line positions (for different volatile species) in the disc, and chemical compositions and mass ratios of ice relative to rock for icy planetesimals in stellar systems of solar chemical composition. From an initial homogeneous composition of the nebula, a wide variety of chemical compositions of planetesimals were produced as a function of the mass of the disc and distance to the star. Ices incorporated in planetesimals are mainly composed of H2O, CO, CO2, CH3OH, and NH3. The ice/rock mass ratio is equal to 1 ± 0.5 in icy planetesimals following assumptions. This last value is in good agreement with observations of solar system comets, but remains lower than usual assumptions made in planet formation models, taking this ratio to be of 2–3.

Global models of planet formation and evolution

Despite the strong increase in observational data on extrasolar planets, the processes that led to the formation of these planets are still not well understood. However, thanks to the high number of extrasolar planets that have been discovered, it is now possible to look at the planets as a population that puts statistical constraints on theoretical formation models. A method that uses these constraints is planetary population synthesis where synthetic planetary populations are generated and compared to the actual population. The key element of the population synthesis method is a global model of planet formation and evolution. These models directly predict observable planetary properties based on properties of the natal protoplanetary disc, linking two important classes of astrophysical objects. To do so, global models build on the simplified results of many specialized models that address one specific physical mechanism. We thoroughly review the physics of the sub-models included in global formation models. The sub-models can be classified as models describing the protoplanetary disc (of gas and solids), those that describe one (proto)planet (its solid core, gaseous envelope and atmosphere), and finally those that describe the interactions (orbital migration and N-body interaction). We compare the approaches taken in different global models, discuss the links between specialized and global models, and identify physical processes that require improved descriptions in future work. We then shortly address important results of planetary population synthesis like the planetary mass function or the mass-radius relationship. With these statistical results, the global effects of physical mechanisms occurring during planet formation and evolution become apparent, and specialized models describing them can be put to the observational test. Owing to their nature as meta models, global models depend on the results of specialized models, and therefore on the development of the field of planet formation theory as a whole. Because there are important uncertainties in this theory, it is likely that the global models will in future undergo significant modifications. Despite these limitations, global models can already now yield many testable predictions. With future global models addressing the geophysical characteristics of the synthetic planets, it should eventually become possible to make predictions about the habitability of planets based on their formation and evolution.

Discovery of a low-mass brown dwarf companion of the young nearby star G196-3.

A substellar-mass object in orbit at about 300 astronomical units from the young low-mass star G 196-3 was detected by direct imaging. Optical and infrared photometry and low- and intermediate-resolution spectroscopy of the faint companion, hereafter referred to as G 196-3B, confirm its cool atmosphere ?15 Jupiter masses. The separation and allow its mass to be estimated at 25?10 between the objects and their mass ratio suggest the fragmentation of a collapsing cloud as the most likely origin for G 196-3B, but alternatively it could have originated from a protoplanetary disc that has been dissipated. Whatever the formation process was, the young age of the primary star (about 100 million years) demonstrates that substellar companions can form on short time scales.

The mechanics of disc degeneration from an engineering and computational modelling viewpoint

The anatomy and microstructure of the spine and in particular the intervertebral disc are intimately linked to how they operate in vivo and how they distribute loads to the adjacent musculature and bony anatomy. The degeneration of the intervertebral discs may be characterised by a loss of hydration, loss of disc height, a granular texture and the presence of annular lesions. As such, degeneration of the intervertebral discs compromises the mechanical integrity of their components and results in adaption and modification in the mechanical means by which loads are distributed between adjacent spinal motion segments.

Expression of chondromodulin-1 in the temporomandibular joint condylar cartilage and disc

BACKGROUND: The temporomandibular joint (TMJ) cartilage consists of condylar cartilage and disc and undergoes continuous remodeling throughout post-natal life. To maintain the integrity of the TMJ cartilage, anti-angiogenic factors play an important role during the remodeling process. In this study, we investigated the expression of the anti-angiogenic factor, chondromodulin- 1 (ChM-1), in TMJ cartilage and evaluate its potential role in TMJ remodeling. METHODS: Eight TMJ specimens were collected from six 4-month-old Japanese white rabbits. Safranin-O staining was performed to determine proteoglycan content. ChM-1 expression in TMJ condylar cartilage and disc was determined by immunohistochemistry. Three human perforated disc tissue samples were collected for investigation of ChM-1 and vascular endothelial growth factor (VEGF) distribution in perforated TMJ disc. RESULTS: Safranin-O stained weakly in TMJ compared with tibial articular and epiphyseal cartilage. In TMJ, ChM-1 was expressed in the proliferative and hypertrophic zone of condylar cartilage and chondrocyte-like cells in the disc. No expression of ChM-1 was observed in osteoblasts and subchondral bone. ChM-1 and VEGF were both similarly expressed in perforated disc tissues. CONCLUSIONS: ChM-1 may play a role in the regulation of TMJ remodeling by preventing blood vessel invasion of the cartilage, thereby maintaining condylar cartilage and disc integrity.