We investigated the formation and evolution of satellite systems in a cold, extended circumplanetary disc around a 10 $M_{\rm{Jupiter}}$ gas giant which was generated by gravitational instability. We use a population synthesis approach, in which we seed the disc with satellite embryos and let them migrate, accrete mass, collide, with subsequent generations being created until the disc is dissipated. In each run we choose randomly the dust-to-gas ratio, dispersion- and refilling time-scales, the number of embryos and their starting locations in a plausible range. The disc structure is the result of a 3D global disc SPH simulation \citep{Szulagyi16b}. We also investigate the effect of the planet's semi-major axis on the resulting satellite systems, taking 50 AU as the nominal case. In the nominal case we find that most satellites are close in mass to the Galilean ones, with a maximum at around 3 $M_{\rm{Earth}}$ and form typically on time-scales comparable to the dispersion time-scale. We also find that about 10 $M_{\rm{Earth}}$ worth of satellites migrate into the planet, polluting it with metals. The influence of a different planet semi-major axis is mainly felt through the change in disc size. We find that for the discs closer to the star, the satellites are lighter, formation time-scales are longer, and more moons are lost into the planet, as a smaller disc means easier migration for the massive satellites. The probability of detecting satellites like the ones in this work is very low ($\leq$ 3%) even with upcoming powerful telescopes like E-ELT.

Satellites of giant planets thought to form in gaseous circumplanetary disks (CPDs) during the late planet-formation phase, but it was unknown so far whether smaller mass planets, such as the ice giants could form such disks, thus moons there. We combined radiative hydrodynamical simulations with satellite population synthesis to investigate the question in the case of Uranus and Neptune. For both ice giants we found that a gaseous CPD is created at the end of their formation. The population synthesis confirmed that Uranian-like, icy, prograde satellite-system could form in these CPDs within a couple of $10^5$ years. This means that Neptune could have a Uranian-like moon-system originally that was wiped away by the capture of Triton. Furthermore, the current moons of Uranus can be reproduced by our model without the need for planet-planet impact to create a debris disk for the moons to grow. These results highlight that even ice giants -- that among the most common mass-category of exoplanets -- can also form satellites, opening a way to a potentially much larger population of exomoons than previously thought.

The satellites of Jupiter are thought to form in a circumplanetary disc. Here we address their formation and orbital evolution with a population synthesis approach, by varying the dust-to-gas ratio, the disc dispersal timescale and the dust refilling timescale. The circumplanetary disc initial conditions (density and temperature) are directly drawn from the results of 3D radiative hydrodynamical simulations. The disc evolution is taken into account within the population synthesis. The satellitesimals were assumed to grow via streaming instability. We find that the moons form fast, often within $10^4$ years, due to the short orbital timescales in the circumplanetary disc. They form in sequence, and many are lost into the planet due to fast type I migration, polluting Jupiter's envelope with typically 15 Earth-masses of metals. The last generation of moons can form very late in the evolution of the giant planet, when the disc has already lost more than the 99% of its mass. The late circumplanetary disc is cold enough to sustain water ice, hence not surprisingly the 85% of the moon population has icy composition. The distribution of the satellite-masses is peaking slightly above Galilean masses, up until a few Earth-masses, in a regime which is observable with the current instrumentation around Jupiter-analog exoplanets orbiting sufficiently close to their host stars. We also find that systems with Galilean-like masses occur in 20% of the cases and they are more likely when discs have long dispersion timescales and high dust-to-gas ratios.