The gravity field of a giant planet is typically our best window into its interior structure and composition. Through comparison of a model planet's calculated gravitational potential with the observed potential, inferences can be made about interior quantities, including possible composition and the existence of a core. Necessarily, a host of assumptions go into such calculations, making every inference about a giant planet's structure strongly model dependent. In this work we present a more general picture by setting Saturn's gravity field, as measured during the \emph{Cassini} Grand Finale, as a likelihood function driving a Markov-Chain Monte Carlo exploration of the possible interior density profiles. The result is a posterior distribution of the interior structure that is not tied to assumed composition, thermal state, or material equations of state. Constraints on interior structure derived in this Bayesian framework are necessarily less informative, but are also less biased and more general. These empirical and probabilistic constraints on the density structure are our main data product which we archive for continued analysis. We find that the outer half of Saturn's radius is relatively well constrained, and we interpret our findings as suggesting a significant metal enrichment, in line with atmospheric abundances from remote sensing. As expected, the inner half of Saturn's radius is less well-constrained by gravity, but we generally find solutions that include a significant density enhancement, which can be interpreted as a core, although this core is often lower in density and larger in radial extent than typically found by standard models. This is consistent with a dilute core and/or composition gradients.

Seismology of the gas giants holds the potential to resolve long-standing questions about their internal structure and rotation state. We construct a family of Saturn interior models constrained by the gravity field and compute their adiabatic mode eigenfrequencies and corresponding Lindblad and vertical resonances in Saturn's C ring, where more than twenty waves with pattern speeds faster than the ring mean motion have been detected and characterized using high-resolution Cassini Visual and Infrared Mapping Spectrometer (VIMS) stellar occultation data. We present identifications of the fundamental modes of Saturn that appear to be the origin of these observed ring waves, and use their observed pattern speeds and azimuthal wavenumbers to estimate the bulk rotation period of Saturn's interior to be $10{\rm h}\, 33{\rm m}\, 38{\rm s}^{+1{\rm m}\, 52{\rm s}}_{-1{\rm m}\, 19{\rm s}}$ (median and 5%/95% quantiles), significantly faster than Voyager and Cassini measurements of periods in Saturn's kilometric radiation, the traditional proxy for Saturn's bulk rotation period. The global fit does not exhibit any clear systematics indicating strong differential rotation in Saturn's outer envelope.

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.

We advance the modeling of rubble-pile solid bodies by re-examining the tidal breakup of comet Shoemaker-Levy 9, an event that occurred during a 1.33 Jupiter radii encounter with that planet in July 1992. Tidal disruption of the comet nucleus led to a chain of sub-nuclei about 100-1000 m in diameter; these went on to collide with the planet two years later (Chodas & Yeomans 1996). They were intensively studied prior to and during the collisions, making SL9 the best natural benchmark for physical models of small body disruption. For the first time in the study of this event, we use numerical codes treating rubble-piles as collections of polyhedra (Korycansky & Asphaug 2009). This introduces forces of dilatation and friction, and inelastic response. As in our previous studies (Asphaug & Benz 1994,1996) we conclude that the progenitor must have been a rubble-pile, and we obtain approximately the same pre-breakup diameter (about 1.5 km) in our best fits to the data. We find that the inclusion of realistic fragment shapes leads to grain locking and dilatancy, so that even in the absence of friction or other dissipation we find that disruption is overall more difficult than in our spheres-based simulations. We constrain the comet's bulk density at about 300-400 kg/m^3, half that of our spheres-based predictions and consistent with recent estimates derived from spacecraft observations.

Numerical simulations, based on the core-nucleated accretion model, are presented for the formation of Jupiter at 5.2 AU in 3 primordial disks with three different assumed values of the surface density of solid particles. The grain opacities in the envelope of the protoplanet are computed using a detailed model that includes settling and coagulation of grains and that incorporates a recalculation of the grain size distribution at each point in time and space. We generally find lower opacities than the 2% of interstellar values used in previous calculations [Hubickyj, O., Bodenheimer, P., Lissauer, J. J., 2005. Icarus 179, 415--431; Lissauer, J. J., Hubickyj, O., D'Angelo, G., Bodenheimer, P., 2009. Icarus 199, 338-350]. These lower opacities result in more rapid heat loss from and more rapid contraction of the protoplanetary envelope. For a given surface density of solids, the new calculations result in a substantial speedup in formation time as compared with those previous calculations. Formation times are calculated to be 1.0, 1.9, and 4.0 Myr, and solid core masses are found to be 16.8, 8.9, and 4.7 Earth masses, for solid surface densities, sigma, of 10, 6, and 4 grams per squared centimeter, respectively. For sigma=10 and sigma=6 g/cm^2, respectively, these formation times are reduced by more than 50% and more than 80% compared with those in a previously published calculation with the old approximation to the opacity.

We have computed the size distribution of silicate grains in the outer radiative region of the envelope of a protoplanet evolving according to the scenario of Pollack et al. (1996). Our computation includes grain growth due to Brownian motion and overtake of smaller grains by larger ones. We also include the input of new grains due to the breakup of planetesimals in the atmosphere. We follow the procedure of Podolak (2003), but have speeded it up significantly. This allows us to test the sensitivity of the code to various parameters. We have also made a more careful estimate of the resulting grain opacity. We find that the grain opacity is of the order of $10^{-2}\ \mathrm{cm^2 g^{-1}}$ throughout most of the outer radiative zone as Hubickyj et al. (2005) assumed for their low opacity case, but near the outer edge of the envelope, the opacity can increase to $\sim{1} \mathrm{cm^2 g^{-1}}$. We discuss the effect of this on the evolution of the models.