We examine the comparative thermal evolution of Jupiter and Saturn applying recent theoretical results for helium's immiscibility in fluid metallic hydrogen. The redistribution of helium in their interiors proceeds very differently for the two planets. We confirm that based on Jupiter's atmospheric helium depletion as observed in situ by the Galileo entry probe, Jupiter's interior helium has differentiated modestly, and we present models reconciling Jupiter's helium depletion, radius, and heat flow at the solar age. Jupiter's recently revised Bond albedo implies a lower intrinsic flux for the planet, accommodating more luminosity from helium differentiation such that mildly superadiabatic interiors can satisfy all constraints. The same phase diagram applied to the less massive Saturn predicts dramatic helium differentiation to the degree that Saturn inevitably forms a helium-rich shell or core, an outcome previously proposed by Stevenson & Salpeter and others. The luminosity from Saturn's helium differentiation is sufficient to extend its cooling time to the solar age, even for adiabatic interiors. This model predicts Saturn's atmospheric helium to be depleted to Y=0.07+/-0.01, corresponding to a He/H_2 mixing ratio 0.036+/-0.006. We also show that neon differentiation may have contributed to both planets' luminosity in the past. These results demonstrate that Jupiter and Saturn's thermal evolution can be explained self-consistently with a single physical model, and emphasize that nontrivial helium distributions should be considered in future models for Saturn's internal structure and dynamo.

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.

We present Direct Numerical Simulations of the transport of heat and heavy elements across a double-diffusive interface or a double-diffusive staircase, in conditions that are close to those one may expect to find near the boundary between the heavy-element rich core and the hydrogen-helium envelope of giant planets such as Jupiter. We find that the non-dimensional ratio of the buoyancy flux associated with heavy element transport to the buoyancy flux associated with heat transport lies roughly between 0.5 and 1, which is much larger than previous estimates derived by analogy with geophysical double-diffusive convection. Using these results in combination with a core-erosion model proposed by Guillot et al. (2004), we find that the entire core of Jupiter would be eroded within less than 1Myr assuming that the core-envelope boundary is composed of a single interface. We also propose an alternative model that is more appropriate in the presence of a well-established double-diffusive staircase, and find that in this limit a large fraction of the core could be preserved. These findings are interesting in the context of Juno's recent results, but call for further modeling efforts to better understand the process of core erosion from first principles.

Hydrogen and helium demix when sufficiently cool, and this bears on the evolution of all giant planets at large separations at or below roughly a Jupiter mass. We model the thermal evolution of Jupiter, including its evolving helium distribution following results of ab initio simulations for helium immiscibility in metallic hydrogen. After 4 Gyr of homogeneous evolution, differentiation establishes a thin helium gradient below 1 Mbar that dynamically stabilizes the fluid to convection. The region undergoes overstable double-diffusive convection (ODDC), whose weak heat transport maintains a superadiabatic temperature gradient. With a generic parameterization for the ODDC efficiency, the models can reconcile Jupiter's intrinsic flux, atmospheric helium content, and radius at the age of the solar system if the Lorenzen et al. H-He phase diagram is translated to lower temperatures. We cast the evolutionary models in an MCMC framework to explore tens of thousands of evolutionary sequences, retrieving probability distributions for the total heavy element mass, the superadiabaticity of the temperature gradient due to ODDC, and the phase diagram perturbation. The adopted SCvH-I equation of state favors inefficient ODDC such that a thermal boundary layer is formed, allowing the molecular envelope to cool rapidly while the deeper interior actually heats up over time. If the overall cooling time is modulated with an additional free parameter to imitate the effect of a colder or warmer EOS, the models favor those that are colder than SCvH-I. In this case the superadiabaticity is modest and a warming or cooling deep interior are equally likely.

Fingering convection (or thermohaline convection) is a weak yet important kind of mixing that occurs in stably-stratified stellar radiation zones in the presence of an inverse mean-molecular-weight gradient. Brown et al. (2013) recently proposed a new model for mixing by fingering convection, which contains no free parameter, and was found to fit the results of direct numerical simulations in almost all cases. Notably, however, they found that mixing was substantially enhanced above their predicted values in the few cases where large-scale gravity waves, followed by thermo-compositional layering, grew spontaneously from the fingering convection. This effect is well-known in the oceanographic context, and is attributed to the excitation of the so-called "collective instability". In this work, we build on the results of Brown et al. (2013) and of Traxler et al. (2011b) to determine the conditions under which the collective instability may be expected. We find that it is only relevant in stellar regions which have a relatively large Prandtl number (the ratio of the kinematic viscosity to the thermal diffusivity), $O(10^{-3})$ or larger. This implies that the collective instability cannot occur in main sequence stars, where the Prandtl number is always much smaller than this (except in the outer layers of surface convection zones where fingering is irrelevant anyway). It could in principle be excited in regions of high electron degeneracy, during He core flash, or in the interiors of white dwarfs. We discuss the implications of our findings for these objects, both from a theoretical and from an observational point of view.

Jupiter's atmosphere has been observed to be depleted in helium (Yatm~0.24), suggesting active helium sedimentation in the interior. This is accounted for in standard Jupiter structure and evolution models through the assumption of an outer, He-depleted envelope that is separated from the He-enriched deep interior by a sharp boundary. Here we aim to develop a model for Jupiter's inhomogeneous thermal evolution that relies on a more self-consistent description of the internal profiles of He abundance, temperature, and heat flux. We make use of recent numerical simulations on H/He demixing, and on layered (LDD) and oscillatory (ODD) double diffusive convection, and assume an idealized planet model composed of a H/He envelope and a massive core. A general framework for the construction of interior models with He rain is described. Despite, or perhaps because of, our simplifications made we find that self-consistent models are rare. For instance, no model for ODD convection is found. We modify the H/He phase diagram of Lorenzen et al. to reproduce Jupiter's atmospheric helium abundance and examine evolution models as a function of the LDD layer height, from those that prolong Jupiter's cooling time to those that actually shorten it. Resulting models that meet the luminosity constraint have layer heights of about 0.1-1 km, corresponding to ~10,-20,000 layers in the rain zone between ~1 and 3-4.5 Mbars. Present limitations and directions for future work are discussed, such as the formation and sinking of He droplets.

Asteroseismology of 1.0-2.0 Msun red giants by the Kepler satellite has enabled the first definitive measurements of interior rotation in both first ascent red giant branch (RGB) stars and those on the Helium burning clump. The inferred rotation rates are 10-30 days for the ~0.2Msun He degenerate cores on the RGB and 30-100 days for the He burning core in a clump star. Using the MESA code we calculate state-of-the-art stellar evolution models of low mass rotating stars from the zero-age main sequence to the cooling white dwarf (WD) stage. We include transport of angular momentum due to rotationally induced instabilities and circulations, as well as magnetic fields in radiative zones (generated by the Tayler-Spruit dynamo). We find that all models fail to predict core rotation as slow as observed on the RGB and during core He burning, implying that an unmodeled angular momentum transport process must be operating on the early RGB of low mass stars. Later evolution of the star from the He burning clump to the cooling WD phase appears to be at nearly constant core angular momentum. We also incorporate the adiabatic pulsation code, ADIPLS, to explicitly highlight this shortfall when applied to a specific Kepler asteroseismic target, KIC8366239. The MESA inlist adopted to calculate the models in this paper can be found at \url{https://authorea.com/1608/} (bottom of the document).

We substantially update the capabilities of the open source software package Modules for Experiments in Stellar Astrophysics (MESA), and its one-dimensional stellar evolution module, MESA Star. Improvements in MESA Star's ability to model the evolution of giant planets now extends its applicability down to masses as low as one-tenth that of Jupiter. The dramatic improvement in asteroseismology enabled by the space-based Kepler and CoRoT missions motivates our full coupling of the ADIPLS adiabatic pulsation code with MESA Star. This also motivates a numerical recasting of the Ledoux criterion that is more easily implemented when many nuclei are present at non-negligible abundances. This impacts the way in which MESA Star calculates semi-convective and thermohaline mixing. We exhibit the evolution of 3-8 Msun stars through the end of core He burning, the onset of He thermal pulses, and arrival on the white dwarf cooling sequence. We implement diffusion of angular momentum and chemical abundances that enable calculations of rotating-star models, which we compare thoroughly with earlier work. We introduce a new treatment of radiation-dominated envelopes that allows the uninterrupted evolution of massive stars to core collapse. This enables the generation of new sets of supernovae, long gamma-ray burst, and pair-instability progenitor models. We substantially modify the way in which MESA Star solves the fully coupled stellar structure and composition equations, and we show how this has improved MESA's performance scaling on multi-core processors. Updates to the modules for equation of state, opacity, nuclear reaction rates, and atmospheric boundary conditions are also provided. We describe the MESA Software Development Kit (SDK) that packages all the required components needed to form a unified and maintained build environment for MESA. [Abridged]