In the widely accepted core accretion model of planet formation, small rocky and/or icy bodies (planetesimals) accrete to form protoplanetary cores. Gas giant planets are believed to have solid cores that must reach a critical mass, ∼10 Earth masses (ME), after which there is rapid inflow of gas from the gas disk. In order to accrete the gas giants’ massive atmospheres, this step must occur within the gas disk’s lifetime (1 − 10 million years). Numerical simulations of solid body accretion in the outer Solar System are performed using two integrators. The goal of these simulations is to investigate the effects of important dynamical processes instead of specifically recreating the formation of the Solar System’s giant planets. The first integrator uses the Symplectic Massive Body Algorithm (SyMBA) with a modification to allow for planetesimal fragmentation. Due to computational constraints, this code has some physical limitations, specifically that the planetesimals themselves cannot grow, so protoplanets must be seeded in the simulations. The second integrator, the Lagrangian Integrator for Planetary Accretion and Dynamics (LIPAD), is more computationally expensive. However, its treatment of planetesimals allows for growth of potential giant planetary cores from a disk consisting only of planetesimals. Thus, this thesis’ preliminary simulations use the first integrator to explore a wider range of parameters while the main simulations use LIPAD to further investigate some specific processes. These simulations are the first use of LIPAD to study giant planet formation and they identify a few important dynamical processes affecting core formation. Without any fragmentation, cores tend to grow to ∼2ME. When planetesimal fragmentation is included, the resulting fragments are easier to accrete and larger cores are formed (∼4ME). But, in half of the runs, the fragments force the entire system to migrate towards the Sun. In other half, outward migration via scattering off a large number of planetesimal helps the protoplanets grow and survive. However, in a preliminary set of simulations including protoplanetary fragmentation, very few collisions are found to result in accretion so it is difficult for any cores to form. / Thesis (Master, Physics, Engineering Physics and Astronomy) -- Queen's University, 2012-08-20 14:48:39.443

<p dir="ltr">We present a model describing an impact jetting origin for the formation of chondrules, the mm– scale, igneous components of chondritic meteorites which originated during the first few million years of solar system history. The ubiquity of chondrules in both non-carbonaceous and carbonaceous chondrites suggests their formation persisted throughout the protoplanetary disk, but their formation mechanism is debated and largely unexplored in the outer disk.<b> </b>Using the iSALE2D shock physics code, we generate models of the process of impact jetting during mixed material (dunite and water ice) impacts that mimic accretionary impacts that form giant planet cores. We show that the process of impact jetting provides the conditions necessary to satisfy critical first-order constraints on chondrule characteristics (size, shape, thermal history). We then explore the implications of chondrule formation by impact jetting during the formation of giant planet cores by combining the original results with simulations of giant planet core accretion generated using a Lagrangian Integrator for Planetary Accretion and Dynamics (LIPAD) code.</p><p dir="ltr">The second closest Galilean satellite to Jupiter is Europa, an ocean world with an outer ice shell and subsurface water ocean encapsulating its rocky core. The surface of Europa is covered in double ridges. These features are defined by two topographic highs about 100 meters tall, with a central trough between them, which extend for hundreds of kilometers over the surface of the moon. Accurate models for the formation of features as prominent as double ridges will help to further constrain the interior structure and dynamics of the interior of the body. We use analytical and numerical finite element models to show that the incremental growth of an ice wedge within the ice shell can cause deformation matching the observed size and shape of observed double ridges on Europa. These models indicate that the total height and width of the ridges correspond to the depth of the wedge, so that deeper wedges create shorter and broader ridges. We consider different sources for the wedge material and ultimately argue in favor of local sources of liquid water within the ice shell.</p>

This thesis describes efforts to improve the realism of numerical models of giant planet formation and migration in an attempt to better understand these processes. A new approach has been taken to the modelling of accretion, designed to mimic reality by allowing gas to accumulate upon a protoplanetary surface. Implementing this treatment in three-dimensional self-gravity radiation hydrodynamics calculations provides an excellent model for planet growth, allowing an exploration of the factors that affect accretion. Moreover, these calculations have also been extended to investigate the migration of protoplanets through their parent discs as they grow. When focusing on the growth of non-migrating protoplanets, the models are performed using small sections of disc, enabling excellent resolution right down to the core; gas structures and flow can be resolved on scales from ~ 10^4 to 10^11 metres. Using radiative transfer, these models reveal the importance of opacity in determining the accretion rates. For the low mass protoplanets, equivalent in mass to a giant planet core (~ 10 M⊕), the accretion rates were found to increase by up to an order of magnitude for a factor of 100 reduction in the grain opacity of the parent circumstellar disc. However, even these low opacities lead to growth rates that are an order of magnitude slower than those obtained in locally-isothermal conditions. For high mass protoplanets (>~ 100M⊕), the accretion rates show very little dependence upon opacity. Nevertheless, the rates obtained using radiative transfer are still lower than those obtained in locally-isothermal models by a factor of ~2, due to the release of accretion energy as heat. Only high mass protoplanets are found to be capable of developing circumplanetary discs, and this ability is dependent upon the opacity, as are the scaleheights of such discs. However, their radial extents were found to be independent of the opacity and the protoplanet mass, all reaching ≈ RH/3, inline with analytic predictions. Migration is investigated using global models, ensuring a self-consistently evolved disc. Using locally-isothermal calculations, it was found that the capture radius of an accreting sink particle, used to model a protoplanet without a surface, must be small (<< RH) to yield migration timescales consistent with linear theory of Type I migration. In the low mass regime of Type I migration, accreting sinks with such small radii yield timescales consistent with those models in which a protoplanetary surface is used. However, for high mass protoplanets, undergoing Type II migration, the surface treatment leads to faster rates of migration, indicating the importance of a realistic accretion model. Using radiative transfer, with high opacities, leads to a factor of ~ 3 increase in the migration timescale of the lowest mass protoplanets, improving their chances of survival. As suitable gas giant progenitors, their survival is key to understanding the growth of giant planets. An unexpected result of the radiative transfer was a reduction in the migration timescale of high mass planets. This appears to be a result of the less thoroughly evacuated gaps created by planets in non-locally-isothermal discs, which affects the corotation torque.

Orientador: Ernesto Vieira Neto / Resumo: O estudo da formação e evolução do Sistema Solar é uma fonte de informação para entender sob quais condições a vida poderia surgir e evoluir. Nós apresentamos, nesta Tese de doutorado, um estudo numérico da fase final de acresção dos planetas gigantes do Sistema Solar durante e após a fase do disco de gás protoplanetário. Em nossas simulações, utilizamos um modelo recente e confiável para a formação de Urano e Netuno para esculpir as propriedades do disco trans-Netuniano original (Izidoro et al. , 2015a). Nós fizemos este estudo de uma maneira autoconsistente considerando os efeitos do gás e da evolução dos embriões planetários que formam Urano e Netuno por colisões gigantescas. Consideramos diferentes histórias de migração de Júpiter, devido a incerteza de como Júpiter migrou, durante a fase de gás. As nossas simulações permitiram obter pela primeira vez as propriedades orbitais do disco trans-Netuniano original. Então, calculamos o tempo de instabilidade dos planetas gigantes a partir de sistemas planetários que formam similares Urano e Netuno. Nossos resultados indicam fortemente que a instabilidade dos planetas gigantes acontecem cedo em até 500 milhões de anos e mais provável ainda ter acontecido em 136 milhões de anos após a dissipação do gás. Nós também realizamos simulações para discutir alguns efeitos dinâmicos que acontecem na região do cinturão de Kuiper. Estes efeitos acontecem quando Netuno esteve em alta excentricidade durante a instabilidade planetária. Para es... (Resumo completo, clicar acesso eletrônico abaixo) / Abstract: A study of the formation and evolution of the Solar System is a source of information for an understanding of what conditions life could arise and evolve. We present a numerical study of the final stage of accretion of the giant planets of the Solar System during and after the protoplanetary gas disc phase. In our simulations, we use a recent and reliable model for the formation of Uranus and Neptune to sculpt the properties of the original trans-Neptunian disk (Izidoro et al. , 2015a). We have done this study in a self-consistent way considering the effects of gas and the evolution of planetary embryos which form Uranus and Neptune by mutual giant collisions. We considered different Jupiter migration stories due to the uncertainty of how Jupiter’s migration was during the gas phase. Our simulations provide for the first time to obtain the orbital properties of the original trans-Neptunian disk. We then calculate the instability time of the giant planets from planetary systems which form similar Uranus and Neptune. Our results strongly indicate that the instability of the giant planets occurs early within 500 million years and even more likely to happen at 136 million years after gas dissipation. We also perform simulations to discuss some dynamical effects that happen in the Kuiper belt region. These effects happen when Neptune was in high eccentricity during planetary instability. For this problem, we use the simulations performed by Gomes et al. (2018) who investigated the... (Complete abstract click electronic access below) / Doutor