The intrinsic luminosity of Uranus is a factor of 10 less than that of Neptune, an observation that standard giant planetary evolution models, which assume negligible viscosity, fail to capture. Here we show that more than half of the interior of Uranus is likely to be in a solid state, and that thermal evolution models that account for this high viscosity region satisfy the observed faintness of Uranus by storing accretional heat deep in the interior. A frozen interior also explains the quality factor of Uranus required by the evolution of the orbits of its satellites.

We investigated the structure of the low density regions of the inner crust of neutron stars using the Hartree-Fock-Bogoliubov (HFB) model to predict the proton content $Z$ of the nuclear clusters and, together with the lattice spacing, the proton content of the crust as a function of the total baryonic density $\rho_b$. The exploration of the energy surface in the $(Z,\rho_b)$ configuration space and the search for the local minima require thousands of calculations. Each of them implies an HFB calculation in a box with a large number of particles, thus making the whole process very demanding. In this work, we apply a statistical model based on a Gaussian Process Emulator that makes the exploration of the energy surface ten times faster. We also present a novel treatment of the HFB equations that leads to an uncertainty on the total energy of $\approx 4$ keV per particle. Such a high precision is necessary to distinguish neighbour configurations around the energy minima.