An observation of Jupiter's tidal response is anticipated for the on-going Juno spacecraft mission. We combine self-consistent, numerical models of Jupiter's equilibrium tidal response with observed Doppler shifts from the Juno gravity science experiment to test the sensitivity of the spacecraft to tides raised by the Galilean satellites and the Sun. The concentric Maclaurin spheroid (CMS) method finds the equilibrium shape and gravity field of a rotating, liquid planet with the tide raised by a satellite, expanded in Love numbers ($k_{nm}$). We present improvements to CMS theory that eliminate an unphysical center of mass offset and study in detail the convergence behavior of the CMS approach. We demonstrate that the dependence of $k_{nm}$ with orbital distance is important when considering the combined tidal response for Jupiter. Conversely, the details of the interior structure have a negligible influence on $k_{nm}$, for models that match the zonal harmonics $J_2$, $J_4$ and $J_6$, already measured to high precision by Juno. As the mission continues, improved coverage of Jupiter's gravity field at different phases of Io's orbit is expected to yield an observed value for the degree-2 Love number ($k_{22}$) and potentially select higher--degree $k_{nm}$. We present a test of the sensitivity of the Juno Doppler signal to the calculated $k_{nm}$, which suggests the detectability of $k_{33}$, $k_{42}$ and $k_{31}$, in addition to $k_{22}$. A mismatch of a robust Juno observation with the remarkably small range in calculated Io equilibrium $k_{22}=0.58976\pm0.0001$, would indicate a heretofore uncharacterized dynamic contribution to the tides.

The Cassini spacecraft's Grand Finale orbits provided a unique opportunity to probe Saturn's gravity field and interior structure. Doppler measurements yielded unexpectedly large values for the gravity harmonics J_6, J_8, and J_10 that cannot be matched with planetary interior models that assume uniform rotation. Instead we present a suite of models that assume the planet's interior rotates on cylinders, which allows us to match all the observed even gravity harmonics. For every interior model, the gravity field is calculated self-consistently with high precision using the Concentric Maclaurin Spheroid (CMS) method. We present an acceleration technique for this method, which drastically reduces the computational cost, allows us to efficiently optimize model parameters, map out allowed parameter regions with Monte Carlo sampling, and increases the precision of the calculated J_2n gravity harmonics to match the error bars of the observations, which would be difficult without acceleration. Based on our models, Saturn is predicted to have a dense central core of 15-18 Earth masses and an additional 1.5-5 Earth masses of heavy elements in the envelope. Finally, we vary the rotation period in the planet's deep interior and determine the resulting oblateness, which we compare with the value from radio occultation measurements by the Voyager spacecraft. We predict a rotation period of 10:33:34 h +- 55s, which is in agreement with recent estimates derived from ring seismology.

The Juno spacecraft has measured Jupiter's low-order, even gravitational moments, $J_2$--$J_8$, to an unprecedented precision, providing important constraints on the density profile and core mass of the planet. Here we report on a selection of interior models based on ab initio computer simulations of hydrogen-helium mixtures. We demonstrate that a dilute core, expanded to a significant fraction of the planet's radius, is helpful in reconciling the calculated $J_n$ with Juno's observations. Although model predictions are strongly affected by the chosen equation of state, the prediction of an enrichment of $Z$ in the deep, metallic envelope over that in the shallow, molecular envelope holds. We estimate Jupiter's core to contain an 7--25 Earth mass of heavy elements. We discuss the current difficulties in reconciling measured $J_n$ with the equations of state, and with theory for formation and evolution of the planet.

This article provides an overview of how models of giant planet interiors are constructed. We review measurements from past space missions that provide constraints for the interior structure of Jupiter. We discuss typical three-layer interior models that consist of a dense central core and an inner metallic and an outer molecular hydrogen-helium layer. These models rely heavily on experiments, analytical theory, and first-principle computer simulations of hydrogen and helium to understand their behavior up to the extreme pressures ~10 Mbar and temperatures ~10,000 K. We review the various equations of state used in Jupiter models and compare them with shock wave experiments. We discuss the possibility of helium rain, core erosion and double diffusive convection may have important consequences for the structure and evolution of giant planets. In July 2016 the Juno spacecraft entered orbit around Jupiter, promising high-precision measurements of the gravitational field that will allow us to test our understanding of gas giant interiors better than ever before.

The accretion of a terrestrial body and differentiation of its silicate/oxide mantle from iron core provide abundant energy for heating its interior to temperatures much higher than the present day Earth. The consequences of differentiation on the structure and composition of planets are typically addressed considering only the interaction of molten iron with an immiscible `rocky' phase. We demonstrate that mixing in a representative system of liquid or solid MgO and liquid iron to a single homogeneous liquid occurs at sufficiently low temperature to be present in the aftermath of a giant impact. Applying the thermodynamic integration technique to density functional theory molecular dynamics simulations, we determine the solvus closure temperature for the Fe-MgO system for pressures up to 400 GPa. Solvus closure occurs at $\sim$4000 K at low pressure, and has a weak positive pressure dependence, such that its gradient with respect to depth is less steep than an adiabatic temperature profile. This predicts a new mode of core-mantle differentiation following the most energetic giant impacts, with exsolution of iron from the mixture beginning in the outer layers of the planet. We demonstrate that high-temperature equilibration results in delivery of nominally insoluble Mg-rich material to the early core. Since MgO is the least soluble major mantle component in iron at low temperatures, these results may represent an upper bound on temperature for mixing in terrestrial planets.

In anticipation of improved observational data for Jupiter's gravitational field from the Juno spacecraft, we predict the static tidal response for a variety of Jupiter interior models based on ab initio computer simulations of hydrogen-helium mixtures. We calculate hydrostatic-equilibrium gravity terms using the non-perturbative concentric Maclaurin Spheroid (CMS) method that eliminates lengthy expansions used in the theory of figures. Our method captures terms arising from the coupled tidal and rotational perturbations, which we find to be important for a rapidly-rotating planet like Jupiter. Our predicted static tidal Love number $k_2 = 0.5900$ is $\sim$10\% larger than previous estimates. The value is, as expected, highly correlated with the zonal harmonic coefficient $J_2$, and is thus nearly constant when plausible changes are made to interior structure while holding $J_2$ fixed at the observed value. We note that the predicted static $k_2$ might change due to Jupiter's dynamical response to the Galilean moons, and find reasons to argue that the change may be detectable, although we do not present here a theory of dynamical tides for highly oblate Jovian planets. An accurate model of Jupiter's tidal response will be essential for interpreting Juno observations and identifying tidal signals from effects of other interior dynamics in Jupiter's gravitational field.

We extend to three dimensions the Concentric Maclaurin Spheroid method for obtaining the self-consistent shape and gravitational field of a rotating liquid planet, to include a tidal potential from a satellite. We exhibit, for the first time, the important effect of the planetary rotation rate on tidal response of gas giants. Simulations of planets with fast rotation rates like those of Jupiter and Saturn, exhibit significant changes in calculated tidal love numbers $k_{nm}$ when compared with non-rotating bodies. A test model of Saturn fitted to observed zonal gravitational multipole harmonics yields $k_2=0.413$, consistent with a recent observational determination from Cassini astrometry data (Lainey et al., 2016). The calculated love number is robust under reasonable assumptions of interior rotation rate, satellite parameters, and details of Saturn's interior structure. The method is benchmarked against several published test cases.

The formation of the giant planets in our solar system, and likely a majority of giant exoplanets, is commonly explained by the accretion of nebular hydrogen and helium onto a large core of terrestrial-like composition. The fate of this core has important consequences for the evolution of the interior structure of the planet. It has recently been shown that H2O, MgO and SiO2 dissolve in liquid metallic hydrogen at high temperature and pressure. In this study, we perform ab initio calculations to study the solubility of an innermost metallic core. We find dissolution of iron to be strongly favored above 2000 K over the entire pressure range (0.4-4 TPa) considered. We compare with and summarize the results for solubilities on other probable core constituents. The calculations imply that giant planet cores are in thermodynamic disequilibrium with surrounding layers, promoting erosion and redistribution of heavy elements. Differences in solubility behavior between iron and rock may influence evolution of interiors, particularly for Saturn-mass planets. Understanding the distribution of iron and other heavy elements in gas giants may be relevant in understanding mass-radius relationships, as well as deviations in transport properties from pure hydrogen-helium mixtures.