We extend the results of planetary formation synthesis by computing the long-term evolution of synthetic systems from the clearing of the gas disk into the dynamical evolution phase. We use the symplectic integrator SyMBA to numerically integrate the orbits of planets for 100 Ma, using populations from previous studies as initial conditions.We show that within the populations studied, mass and semi-major axis distributions experience only minor changes from post-formation evolution. We also show that, depending upon their initial distribution, planetary eccentricities can statistically increase or decrease as a result of gravitational interactions. We find that planetary masses and orbital spacings provided by planet formation models do not result in eccentricity distributions comparable to observed exoplanet eccentricities, requiring other phenomena such as e.g. stellar fly-bys to account for observed eccentricities.

In the standard model of core accretion, the cores of the giant planets form by the accretion of planetesimals. In this scenario, the evolution of the planetesimal population plays an important role in the formation of massive cores. Recently, we studied the role of planetesimal fragmentation in the in situ formation of a giant planet. However, the exchange of angular momentum between the planet and the gaseous disk causes the migration of the planet in the disk. In this new work, we incorporate the migration of the planet and globally study the role of planet migration in the formation of a massive core when the population of planetesimals evolves by planet accretion, migration due to the nebular drag, and fragmentation due to planetesimal collisions.

Planet formation models have been developed during the last years in order to try to reproduce the observations of both the solar system, and the extrasolar planets. Some of these models have partially succeeded, focussing however on massive planets, and for the sake of simplicity excluding planets belonging to planetary systems. However, more and more planets are now found in planetary systems. This tendency, which is a result of both radial velocity, transit and direct imaging surveys, seems to be even more pronounced for low mass planets. These new observations require the improvement of planet formation models, including new physics, and considering the formation of systems. 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 global architecture and composition of this multiplicity. We use a 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. We show that the masses and semi-major axis 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.