System Giant Planets Research Articles

The Solar System could have formed in a low-viscosity disc: A dynamical study from giant planet migration to the Nice model

Context. In the context of low-viscosity protoplanetary discs (PPDs), the formation scenarios of the Solar System should be revisited. In particular, the Jupiter-Saturn pair has been shown to lock in the 2:1 mean motion resonance while migrating generally inwards, making the Grand Tack scenario impossible. Aims. We explore what resonant chains of multiple giant planets can form in a low-viscosity disc, and whether these configurations can evolve into forming the Solar System in the post gas disc phase. Methods. We used hydrodynamical simulations with the code FARGOCA to study the migration of the giant planets in a disc with viscosity parameter of α = 10−4. After a transition phase to a gas-less configuration, we studied the stability of the obtained resonant chains through their interactions with a disc of leftover planetesimals by performing N-body simulations using rebound. Results. The gaps opened by giant planets are wider and deeper for lower viscosity, reducing the damping effect of the disc. Thus, when planets enter a resonance, the resonant angle remains closer to circulation, making the chain weaker. Exploring numerous configurations, we found five stable resonant chains of four or five planets. In a thin (cold) PPD, the four giant planets revert their migration and migrate outwards. After disc dispersal, under the influence of a belt of planetesimals, some resonant chains undergo an instability phase while others migrate smoothly over a billion years. For three of our resonant chains, about ~1% of the final configurations pass the four criteria to fit the Solar System. The most successful runs are obtained for systems formed in a cold PPD with a massive planetesimal disc. Conclusions. This work provides a fully consistent study of the dynamical history of the Solar System’s giant planets, from the protoplanetary disc phase up to the giant planet instability. Although building resonant configurations is difficult in low-viscosity discs, we find it possible to reproduce the Solar System from a cold, low-viscosity protoplanetary disc.

Uranus’s Influence on Neptune’s Exterior Mean-motion Resonances

Neptune’s external mean-motion resonances play an important role in sculpting the observed population of trans-Neptunian objects (TNOs). The population of scattering TNOs is known to “stick” to Neptune's resonances while evolving in semimajor axis (a), though simulations show that resonance sticking is less prevalent at a ≳ 200–250 au. Here we present an extensive numerical exploration of the strengths of Neptune's resonances for scattering TNOs with perihelion distances q = 33 au. We show that the drop-off in resonance sticking for the large a scattering TNOs is not a generic feature of scattering dynamics but can instead be attributed to the specific configuration of Neptune and Uranus in our solar system. In simulations with just Uranus removed from the giant planet system, Neptune's resonances are strong in the scattering population out to at least ∼300 au. Uranus and Neptune are near a 2:1 period ratio, and the variations in Neptune's orbit resulting from this near-resonance are responsible for destabilizing Neptune's resonances for high-e TNO orbits beyond the ∼20:1 resonance at a ≈ 220 au. Direct interactions between Uranus and the scattering population are responsible for slightly weakening Neptune's closer-in resonances. In simulations where Neptune and Uranus are placed in their mutual 2:1 resonance, we see almost no stable libration of scattering particles in Neptune's external resonances. Our results have important implications for how the strengths of Neptune's distant resonances varied during the epoch of planet migration when the Neptune–Uranus period ratio was evolving. These strength variations likely affected the distant scattering, resonant, and detached TNO populations.

Relative Occurrence Rate between Hot and Cold Jupiters as an Indicator to Probe Planet Migration

We propose a second-order statistic parameter ε, the relative occurrence rate between hot Jupiters (HJs) and cold Jupiters (CJs) (ε = η HJ/η CJ), to probe the migration of gas giants. Since the planet occurrence rate is the combined outcome of the formation and migration processes, a joint analysis of HJ and CJ frequency may shed light on the dynamical evolution of giant planet systems. We first investigate the behavior of ε as the stellar mass changes observationally. Based on the occurrence rate measurements of HJs (η HJ) from the Transiting Exoplanet Survey Satellite survey and CJs (η CJ) from the California Legacy Survey, we find a tentative trend (97% confidence) that ε drops when the stellar mass rises from 0.8 to 1.4 M ⊙, which can be explained by different giant planet growth and disk migration timescales around different stars. We carry out planetesimal and pebble accretion simulations, both of which can reproduce the results of η HJ, η CJ, and ε. Our findings indicate that the classical core accretion + disk migration model can explain the observed decreasing trend of ε. We propose two ways to increase the significance of the trend and verify the anticorrelation. Future works are required to better constrain ε, especially for M dwarfs and for more massive stars.

Long-term double synchronization in close-in gas giant planets

ABSTRACT Hot Jupiters, orbiting their host stars at extremely close distances, undergo tidal evolution, with some being engulfed by their stars due to angular momentum exchanges induced by tidal forces. However, achieving double synchronization can prolong their survival. Using the mesa stellar evolution code, combined with the magnetic braking model of Matt et al. (2015), we calculate 25 000 models with different metallicity and study how to attain the conditions that trigger the long-term double synchronization. Our results indicate that massive planets orbiting stars with lower convective turnover time are easier to achieve long-term double synchronization. The rotation angular velocity at the equilibrium point (Ωsta) is almost equal to orbital angular velocity of planet (n) for the majority of the main sequence lifetime if a system has undergone a long-term double synchronization, regardless of their state at this moment. We further compared our results with known parameters of giant planetary systems and found that those systems with larger planetary masses and lower convective turnover time seem to be less sensitive to changes in the tidal quality factor $Q^{\prime }_{_*}$. We suggest that for systems that fall on the state of Ωsta ≈ n, regardless of their current state, the synchronization will persist for a long time if orbital synchronization occurs at any stage of their evolution. Our results can be applied to estimate whether a system has experienced long-term double synchronization in the past or may experience it in the future.

Sweeping secular resonances and giant planet inclinations in transition discs

ABSTRACT The orbits of some warm Jupiters are highly inclined (20°–50°) to those of their exterior companions. Comparable misalignments are inferred between the outer and inner portions of some transition discs. These large inclinations may originate from planet–planet and planet–disc secular resonances that sweep across interplanetary space as parent discs disperse. The maximum factor by which a seed mutual inclination can be amplified is of the order of the square root of the angular momentum ratio of the resonant pair. We identify those giant planet systems (e.g. Kepler-448 and Kepler-693) that may have crossed a secular resonance, and estimate the required planet masses and semimajor axes in transition discs needed to warp their innermost portions (e.g. in CQ Tau). Passage through an inclination secular resonance could also explain the hypothesized large mutual inclinations in apsidally-orthogonal warm Jupiter systems (e.g. HD 147018).

Age Distribution of Exoplanet Host Stars: Chemical and Kinematic Age Proxies from GAIA DR3

The GAIA space mission is impacting astronomy in many significant ways by providing a uniform, homogeneous, and precise data set for over 1 billion stars and other celestial objects in the Milky Way and beyond. Exoplanet science has greatly benefited from the unprecedented accuracy of the stellar parameters obtained from GAIA. In this study, we combine photometric, astrometric, and spectroscopic data from the most recent Gaia DR3 to examine the kinematic and chemical age proxies for a large sample of 2611 exoplanets hosting stars whose parameters have been determined uniformly. Using spectroscopic data from the Radial Velocity Spectrometer on board GAIA, we show that stars hosting massive planets are metal-rich and α-poor in comparison to stars hosting small planets. The kinematic analysis of the sample reveals that stellar systems with small planets and those with giant planets differ in key aspects of galactic space velocity and orbital parameters, which are indicative of age. We find that the galactic orbital parameters have a statistically significant difference of 0.06 kpc for Z max and 0.03 for eccentricity, respectively. Furthermore, we estimated the stellar ages of the sample using the MIST-MESA isochrone models. The ages and their proxies for the planet-hosting stars indicate that the hosts of giant planetary systems are younger when compared to the population of stars harboring small planets. These age trends are also consistent with the chemical evolution of the galaxy and the formation of giant planets from the core-accretion process.

Large Interferometer for Exoplanets: VIII. Where Is the Phosphine? Observing Exoplanetary PH3 with a Space-Based Mid-Infrared Nulling Interferometer.

Phosphine could be a key molecule in the understanding of exotic chemistry that occurs in (exo)planetary atmospheres. While phosphine has been detected in the Solar System's giant planets, it has not been observed in exoplanets to date. In the exoplanetary context, however, it has been theorized to be a potential biosignature molecule. The goal of our study was to identify which illustrative science cases for PH3 chemistry are observable with a space-based mid-infrared nulling interferometric observatory like the Large Interferometer for Exoplanets (LIFE) concept. We identified a representative set of scenarios for PH3 detections in exoplanetary atmospheres that vary over the whole dynamic range of the LIFE mission. We used chemical kinetics and radiative transfer calculations to produce forward models of these informative, prototypical observational cases for LIFEsim, our observation simulator software for LIFE. In a detailed, yet first order approximation, it takes a mission like LIFE: (i) about 1 h to find phosphine in a warm giant around a G star at 10 pc, (ii) about 10 h in H2 or CO2 dominated temperate super-Earths around M star hosts at 5 pc, (iii) and even in 100 h it seems very unlikely that phosphine would be detectable in a Venus-Twin with extreme PH3 concentrations at 5 pc. Phosphine in concentrations previously discussed in the literature is detectable in 2 out of the 3 cases, and it is detected about an order of magnitude faster than in comparable cases with James Webb Space Telescope. We show that there is a significant number of objects accessible for these classes of observations. These results will be used to prioritize the parameter range for the next steps with more detailed retrieval simulations. They will also inform timely questions in the early design phase of a mission like LIFE and guide the community by providing easy-to-scale first estimates for a large part of detection space of such a mission.

Breakdown of planetary systems in embedded clusters

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.

Precise radial velocities of giant stars

Context. Radial velocity surveys of evolved stars allow us to probe a higher stellar mass range, on average, compared to main-sequence samples. Hence, differences between the planet populations around the two target classes can be caused by either the differing stellar mass or stellar evolution. To properly disentangle the effects of both variables, it is important to characterize the planet population around giant stars as accurately as possible. Aims. Our goal is to investigate the giant planet occurrence rate around evolved stars and determine its dependence on stellar mass, metallicity, and orbital period. Methods. We combine data from three different radial velocity surveys targeting giant stars: the Lick giant star survey, the radial velocity program EXoPlanets aRound Evolved StarS (EXPRESS), and the Pan-Pacific Planet Search (PPPS), yielding a sample of 482 stars and 37 planets. We homogeneously rederived the stellar parameters of all targets and accounted for varying observational coverage, precision and stellar noise properties by computing a detection probability map for each star via injection and retrieval of synthetic planetary signals. We then computed giant planet occurrence rates as a function of period, stellar mass, and metallicity, corrected for incompleteness. Results. Our findings agree with previous studies that found a positive planet-metallicity correlation for evolved stars and identified a peak in the giant planet occurrence rate as a function of stellar mass, but our results place it at a slightly smaller mass of (1.68 ± 0.59) M⊙. The period dependence of the giant planet occurrence rate seems to follow a broken power-law or log-normal distribution peaking at (718 ± 226) days or (797 ± 455) days, respectively, which roughly corresponds to 1.6 AU for a 1 M⊙ star and 2.0 AU for a 2 M⊙ star. This peak could be a remnant from halted migration around intermediate-mass stars, caused by stellar evolution, or an artifact from contamination by false positives. The completeness-corrected global occurrence rate of giant planetary systems around evolved stars is 10.7%−1.6%+2.2% for the entire sample, while the evolutionary subsets of RGB and HB stars exhibit 14.2%−2.7%+4.1% and 6.6%−1.3%+2.1%, respectively. However, both subsets have different stellar mass distributions and we demonstrate that the stellar mass dependence of the occurrence rate suffices to explain the apparent change of occurrence with the evolutionary stage.

A Geoscientific Review on CO and CO2 Ices in the Outer Solar System

Ground-based telescopes and space exploration have provided outstanding observations of the complexity of icy planetary surfaces. This work presents our review of the varying nature of carbon dioxide (CO2) and carbon monoxide (CO) ices from the cold traps on the Moon to Pluto in the Kuiper Belt. This review is organized into five parts. First, we review the mineral physics (e.g., rheology) relevant to these environments. Next, we review the radiation-induced chemical processes and the current interpretation of spectral signatures. The third section discusses the nature and distribution of CO2 in the giant planetary systems of Jupiter and Saturn, which are much better understood than the satellites of Uranus and Neptune, discussed in the subsequent section. The final sections focus on Pluto in comparison to Triton, having mainly CO, and a brief overview of cometary materials. We find that CO2 ices exist on many of these icy bodies by way of magnetospheric influence, while intermixing into solid ices with CH4 (methane) and N2 (nitrogen) out to Triton and Pluto. Such radiative mechanisms or intermixing can provide a wide diversity of icy surfaces, though we conclude where further experimental research of these ices is still needed.

Transformative Planetary Science with the US ELT Program

The proposed US Extremely Large Telescope (ELT) Program would secure national open access to at least 25% of the observing time on the Thirty Meter Telescope in the north and the Giant Magellan Telescope in the south. ELTs would advance solar system science via exceptional angular resolution, sensitivity, and advanced instrumentation. ELT contributions would include the study of interstellar objects, giant planet systems and ocean worlds, the formation of the solar system traced through small objects in the asteroid and Kuiper belts, and the active support of planetary missions. We recommend that (1) the US ELT Program be listed as critical infrastructure for solar system science, that (2) some support from NASA be provided to ensure mission support capabilities, and that (3) the US ELT Program expand solar-system community participation in development, planning, and operations.

Icy Exomoons Evidenced by Spallogenic Nuclides in Polluted White Dwarfs

Abstract We present evidence that excesses in Be in polluted white dwarfs (WDs) are the result of accretion of icy exomoons that formed in the radiation belts of giant exoplanets. Here we use excess Be in the white dwarf GALEX J2339–0424 as an example. We constrain the parent body abundances of rock-forming elements in GALEX J2339–0424 and show that the overabundance of beryllium in this WD cannot be accounted for by differences in diffusive fluxes through the WD outer envelope nor by chemical fractionations during typical rock-forming processes. We argue instead that the Be was produced by energetic proton irradiation of ice mixed with rock. We demonstrate that the MeV proton fluence required to form the high Be/O ratio in the accreted parent body is consistent with irradiation of ice in the rings of a giant planet within its radiation belt, followed by accretion of the ices to form a moon that is later accreted by the WD. The icy moons of Saturn serve as useful analogs. Our results provide an estimate of spallogenic nuclide excesses in icy moons formed by rings around giant planets in general, including those in the solar system. While excesses in Be have been detected in two polluted WDs to date, including the WD described here, we predict that excesses in the other spallogenic elements Li and B, although more difficult to detect, should also be observed, and that such detections would also indicate pollution by icy exomoons formed in the ring systems of giant planets.

Auroral emissions from Uranus and Neptune.

Uranus and Neptune possess highly tilted/offset magnetic fields whose interaction with the solar wind shapes unique twin asymmetric, highly dynamical, magnetospheres. These radiate complex auroral emissions, both reminiscent of those observed at the other planets and unique to the ice giants, which have been detected at radio and ultraviolet (UV) wavelengths to date. Our current knowledge of these radiations, which probe fundamental planetary properties (magnetic field, rotation period, magnetospheric processes, etc.), still mostly relies on Voyager 2 radio, UV and in situ measurements, when the spacecraft flew by each planet in the 1980s. These pioneering observations were, however, limited in time and sampled specific solar wind/magnetosphere configurations, which significantly vary at various timescales down to a fraction of a planetary rotation. Since then, despite repeated Earth-based observations at similar and other wavelengths, only the Uranian UV aurorae have been re-observed at scarce occasions by the Hubble Space Telescope. These observations revealed auroral features radically different from those seen by Voyager 2, diagnosing yet another solar wind/magnetosphere configuration. Perspectives for the in-depth study of the Uranian and Neptunian auroral processes, with implications for exoplanets, include follow-up remote Earth-based observations and future orbital exploration of one or both ice giant planetary systems. This article is part of a discussion meeting issue 'Future exploration of ice giant systems'.

On the origin of the eccentricity dichotomy displayed by compact super-Earths: dynamical heating by cold giants

ABSTRACT Approximately half of the planets discovered by NASA’s Kepler mission are in systems where just a single planet transits its host star, and the remaining planets are observed to be in multiplanet systems. Recent analyses have reported a dichotomy in the eccentricity distribution displayed by systems where a single planet transits compared with that displayed by the multiplanet systems. Using N-body simulations, we examine the hypothesis that this dichotomy has arisen because inner systems of super-Earths are frequently accompanied by outer systems of giant planets that can become dynamically unstable and perturb the inner systems. Our initial conditions are constructed using a subset of the known Kepler five-planet systems as templates for the inner systems, and systems of outer giant planets with masses between those of Neptune and Saturn that are centred on orbital radii 2 ≤ ap ≤ 10 au. The parameters of the outer systems are chosen so that they are always below an assumed radial velocity detection threshold of 3 m s−1. The results show an inverse relation between the mean eccentricities and the multiplicities of the systems. Performing synthetic transit observation of the final systems reveals dichotomies in both the eccentricity and multiplicity distributions that are close to being in agreement with the Kepler data. Hence, understanding the observed orbital and physical properties of the compact systems of super-Earths discovered by Kepler may require holistic modelling that couples the dynamics of both inner and outer systems of planets during and after the epoch of formation.

Resonance locking in giant planets indicated by the rapid orbital expansion of Titan

Tidal effects in planetary systems are the main driver in the orbital migration of natural satellites. They result from physical processes occurring deep inside celestial bodies, whose effects are rarely observable from surface imaging. For giant planet systems, the tidal migration rate is determined by poorly understood dissipative processes in the planet, and standard theories suggest an orbital expansion rate inversely proportional to the power 11/2 in distance, implying little migration for outer moons such as Saturn's largest moon, Titan. Here, we use two independent measurements obtained with the Cassini spacecraft to measure Titan's orbital expansion rate. We find Titan migrates away from Saturn at 11.3 $\pm$ 2.0 cm/year, corresponding to a tidal quality factor of Saturn of Q $\simeq$ 100, and a migration timescale of roughly 10 Gyr. This rapid orbital expansion suggests Titan formed significantly closer to Saturn and has migrated outward to its current position. Our results for Titan and five other moons agree with the predictions of a resonance locking tidal theory, sustained by excitation of inertial waves inside the planet. The associated tidal expansion is only weakly sensitive to orbital distance, motivating a revision of the evolutionary history of Saturn's moon system. The resonance locking mechanism could operate in other systems such as stellar binaries and exoplanet systems, and it may allow for tidal dissipation to occur at larger orbital separations than previously believed.

Do Metal-rich Stars Make Metal-rich Planets? New Insights on Giant Planet Formation from Host Star Abundances* †

Abstract The relationship between the compositions of giant planets and their host stars is of fundamental interest in understanding planet formation. The solar system giant planets are enhanced above solar composition in metals, both in their visible atmospheres and bulk compositions. A key question is whether the metal enrichment of giant exoplanets is correlated with that of their host stars. Thorngren et al. showed that in cool (T eq < 1000 K) giant exoplanets, the total heavy-element mass increases with total M p and the heavy-element enrichment relative to the parent star decreases with total M p . In their work, the host star metallicity was derived from literature [Fe/H] measurements. Here we conduct a more detailed and uniform study to determine whether different host star metals (C, O, Mg, Si, Fe, and Ni) correlate with the bulk metallicity of their planets, using correlation tests and Bayesian linear fits. We present new host star abundances of 19 cool giant planet systems, and combine these with existing host star data for a total of 22 cool giant planet systems (24 planets). Surprisingly, we find no clear correlation between stellar metallicity and planetary residual metallicity (the relative amount of metal versus that expected from the planet mass alone), which is in conflict with common predictions from formation models. We also find a potential correlation between residual planet metals and stellar volatile-to-refractory element ratios. These results provide intriguing new relationships between giant planet and host star compositions for future modeling studies of planet formation.

Prevalence of Fibonacci numbers in orbital period ratios in solar planetary and satellite systems and in exoplanetary systems

It is shown that orbital period ratios of successive planets in the Solar System, of satellites in giant planet systems and of exoplanets in exoplanetary systems are preferentially closer to irreducible fractions formed with Fibonacci numbers between 1 and 8 than to other fractions, in a ratio of approximately 60% vs 40%. Furthermore, if sets of minor planets are chosen with gradually smaller inclinations and eccentricities, the proximity to Fibonacci fractions of their period ratios with Jupiter or Mars’ period tends to increase. Finally, a simple model explains why the resonance of the form $\frac{P_{1}}{P_{2}} = \frac{p}{p+q}$ , with $P_{1}$ and $P_{2}$ orbital periods of planets or satellites and $p$ and $q$ small integers, are stronger and more commonly observed for $p$ and $( p+q )$ being both small Fibonacci numbers than for other cases.

Planet–planet scattering as the source of the highest eccentricity exoplanets

Most giant exoplanets discovered by radial velocity surveys have much higher eccentricities than those in the solar system. The planet–planet scattering mechanism has been shown to match the broad eccentricity distribution, but the highest-eccentricity planets are often attributed to Kozai-Lidov oscillations induced by a stellar companion. Here we investigate whether the highly eccentric exoplanet population can be produced entirely by scattering. We ran 500 N-body simulations of closely packed giant-planet systems that became unstable under their own mutual perturbations. We find that the surviving bound planets can have eccentricities up to e > 0.99, with a maximum of 0.999017 in our simulations. This suggests that there is no maximum eccentricity that can be produced by planet–planet scattering. Importantly, we find that extreme eccentricities are not extremely rare; the eccentricity distribution for all giant exoplanets with e > 0.3 is consistent with all planets concerned being generated by scattering. Our results show that the discovery of planets with extremely high eccentricities does not necessarily signal the action of the Kozai-Lidov mechanism.

Layered semi-convection and tides in giant planet interiors

Context. Recent Juno observations have suggested that the heavy elements in Jupiter could be diluted throughout a large fraction of its gaseous envelope, providing a stabilising compositional gradient over an extended region of the planet. This could trigger layered semi-convection, which, in the context of giant planets more generally, may explain Saturn’s luminosity excess and play a role in causing the abnormally large radii of some hot Jupiters. In giant planet interiors, it could take the form of density staircases, which are convective layers separated by thin stably stratified interfaces. In addition, the efficiency of tidal dissipation is known to depend strongly on the planetary internal structure. Aims. We aim to study the resulting tidal dissipation when internal waves are excited in a region of layered semi-convection by tidal gravitational forcing due to other bodies (such as moons in giant planet systems, or stars in hot Jupiter systems). Methods. We adopt a local Cartesian model with a background layered density profile subjected to an imposed tidal forcing, and we compute the viscous and thermal dissipation rates numerically. We consider two sets of boundary conditions in the vertical direction: periodic boundaries and impenetrable, stress-free boundaries, with periodic conditions in the horizontal directions in each case. These models are appropriate for studying the forcing of short-wavelength tidal waves in part of a region of layered semi-convection, and in an extended envelope containing layered semi-convection, respectively. Results. We find that the rates of tidal dissipation can be enhanced in a region of layered semi-convection compared to a uniformly convective medium, where the latter corresponds with the usual assumption adopted in giant planet interior models. In particular, a region of layered semi-convection possesses a richer set of resonances, allowing enhanced dissipation for a wider range of tidal frequencies. The details of these results significantly depend on the structural properties of the layered semi-convective regions. Conclusions. Layered semi-convection could contribute towards explaining the high tidal dissipation rates observed in Jupiter and Saturn, which have not yet been fully explained by theory. Further work is required to explore the efficiency of this mechanism in global models.

Formation of planetary systems by pebble accretion and migration: growth of gas giants

Giant planets migrate though the protoplanetary disc as they grow their solid core and attract their gaseous envelope. Previously, we have studied the growth and migration of an isolated planet in an evolving disc. Here, we generalise such models to include the mutual gravitational interaction between a high number of growing planetary bodies. We have investigated how the formation of planetary systems depends on the radial flux of pebbles through the protoplanetary disc and on the planet migration rate. Our N-body simulations confirm previous findings that Jupiter-like planets in orbits outside the water ice line originate from embryos starting out at 20–40 AU when using nominal type-I and type-II migration rates and a pebble flux of approximately 100–200 Earth masses per million years, enough to grow Jupiter within the lifetime of the solar nebula. The planetary embryos placed up to 30 AU migrate into the inner system (rP < 1AU). There they form super-Earths or hot and warm gas giants, producing systems that are inconsistent with the configuration of the solar system, but consistent with some exoplanetary systems. We also explored slower migration rates which allow the formation of gas giants from embryos originating from the 5–10 AU region, which are stranded exterior to 1 AU at the end of the gas-disc phase. These giant planets can also form in discs with lower pebbles fluxes (50–100 Earth masses per Myr). We identify a pebble flux threshold below which migration dominates and moves the planetary core to the inner disc, where the pebble isolation mass is too low for the planet to accrete gas efficiently. In our model, giant planet growth requires a sufficiently high pebble flux to enable growth to out-compete migration. An even higher pebble flux produces systems with multiple gas giants. We show that planetary embryos starting interior to 5 AU do not grow into gas giants, even if migration is slow and the pebble flux is large. These embryos instead grow to just a few Earth masses, the mass regime of super-Earths. This stunted growth is caused by the low pebble isolation mass in the inner disc and is therefore independent of the pebble flux. Additionally, we show that the long-term evolution of our formed planetary systems can naturally produce systems with inner super-Earths and outer gas giants as well as systems of giant planets on very eccentric orbits.