The late-stage formation of giant planetary systems is rich in interesting dynamical mechanisms. Previous simulations of three giant planets initially on quasi-circular and quasi-coplanar orbits in the gas disc have shown that highly mutually inclined configurations can be formed, despite the strong eccentricity and inclination damping exerted by the disc. Much attention has been directed to inclination-type resonance, asking for large eccentricities to be acquired during the migration of the planets. Here we show that inclination excitation is also present at small to moderate eccentricities in two-planet systems that have previously experienced an ejection or a merging and are close to resonant commensurabilities at the end of the gas phase. We perform a dynamical analysis of these planetary systems, guided by the computation of planar families of periodic orbits and the bifurcation of families of spatial periodic orbits. We show that inclination excitation at small to moderate eccentricities can be produced by (temporary) capture in inclination-type resonance and the possible proximity of the non-coplanar systems to spatial periodic orbits contributes to maintaining their mutual inclination over long periods of time.

Evidence of mutually inclined planetary orbits has been reported for giant planets these last years. Here we aim to study the impact of eccentric and inclined massive giant planets on the terrestrial planet formation process, and investigate whether it can possibly lead to the existence of inclined terrestrial planets. We have performed 126 simulations of the late-stage planetary accretion in eccentric and inclined giant planet systems. The physical and orbital parameters of the giant planet systems result from n-body simulations of three giant planets in the late stage of the gas disc, under the combined action of Type II migration and planet-planet scattering. Fourteen two- and three-planet configurations have been selected, with diversified masses, semi-major axes (resonant configurations or not), eccentricities and inclinations (including coplanar systems) at the dispersal of the gas disc. We have then followed the gravitational interactions of these systems with an inner disc of planetesimals and embryos (9 runs per system), studying in detail the final configurations of the formed terrestrial planets. While coplanar giant systems harbour several terrestrial planets, generally as massive as the Earth and mainly on low eccentric and low inclined orbits, terrestrial planets formed in systems with mutually inclined giant planets are usually fewer, less massive (<0.5 M_{\Earth}) and with larger eccentricities and inclinations. This work shows that terrestrial planets can form on stable inclined orbits through the classical accretion theory, even in coplanar giant planet systems emerging from the disc phase.

We aim to investigate the influence of the eccentricity and inclination damping due to planet-disc interactions on the final configurations of the systems, generalizing previous studies on the combined action of the gas disc and planet-planet scattering during the disc phase. Instead of the simplistic $K$-prescription, our n-body simulations adopt the damping formulae for eccentricity and inclination provided by the hydrodynamical simulations of our companion paper. We follow the evolution of $11000$ numerical experiments of three giant planets in the late stage of the gas disc, exploring different initial configurations, planetary mass ratios and disc masses. The dynamical evolutions of the planetary systems are studied along the simulations, with emphasis on the resonance captures and inclination-growth mechanisms. Most of the systems are found with small inclinations ($\le10^{\circ}$) at the dispersal of the disc. Even though many systems enter an inclination-type resonance during the migration, the disc usually damps the inclinations on a short timescale. Although the majority of the multiple systems in our results are quasi-coplanar, $\sim5\%$ of them end up with high mutual inclinations ($\ge10^{\circ}$). Half of these highly mutually inclined systems result from two- or three-body MMR captures, the other half being produced by orbital instability and/or planet-planet scattering. When considering the long-term evolution over $100$ Myr, destabilization of the resonant systems is common, and the percentage of highly mutually inclined systems still evolving in resonance drops to $30\%$. Finally, the parameters of the final system configurations are in very good agreement with the semi-major axis and eccentricity distributions in the observations, showing that planet-planet interactions during the disc phase could have played an important role in sculpting planetary systems.