skip to main content

Title: Theory of Figures to the Seventh Order and the Interiors of Jupiter and Saturn
Abstract Interior modeling of Jupiter and Saturn has advanced to a state where thousands of models are generated that cover the uncertainty space of many parameters. This approach demands a fast method of computing their gravity field and shape. Moreover, the Cassini mission at Saturn and the ongoing Juno mission delivered gravitational harmonics up to J 12 . Here we report the expansion of the theory of figures, which is a fast method for gravity field and shape computation, to the seventh order (ToF7), which allows for computation of up to J 14 . We apply three different codes to compare the accuracy using polytropic models. We apply ToF7 to Jupiter and Saturn interior models in conjunction with CMS-19 H/He equation of state. For Jupiter, we find that J 6 is best matched by a transition from an He-depleted to He-enriched envelope at 2–2.5 Mbar. However, the atmospheric metallicity reaches 1 × solar only if the adiabat is perturbed toward lower densities, or if the surface temperature is enhanced by ∼14 K from the Galileo value. Our Saturn models imply a largely homogeneous-in-Z envelope at 1.5–4 × solar atop a small core. Perturbing the adiabat yields metallicity profiles with extended, more » heavy-element-enriched deep interior (diffuse core) out to 0.4 R Sat , as for Jupiter. Classical models with compact, dilute, or no core are possible as long as the deep interior is enriched in heavy elements. Including a thermal wind fitted to the observed wind speeds, representative Jupiter and Saturn models are consistent with all observed J n values. « less
; ; ; ; ; ; ; ;
Award ID(s):
Publication Date:
Journal Name:
The Planetary Science Journal
Page Range or eLocation-ID:
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract The observed correlation between outer giant planets and inner super-Earths is emerging as an important constraint on planet formation theories. In this study, we focus on Kepler-167, which is currently the only system known to contain both inner transiting super-Earths and a confirmed outer transiting gas giant companion beyond 1 au. Using long-term radial velocity monitoring, we measure the mass of the gas giant Kepler-167e ( P = 1071 days) to be 1.01 − 0.15 + 0.16 M J , thus confirming it as a Jupiter analog. We refit the Kepler photometry to obtain updated radii for all four planets. Using a planetary structure model, we estimate that Kepler-167e contains 66 ± 19 M ⊕ of solids and is significantly enriched in metals relative to its solar-metallicity host star. We use these new constraints to explore the broader question of how systems like Kepler-167 form in the pebble accretion framework for giant planet core formation. We utilize simple disk evolution models to demonstrate that more massive and metal-rich disks, which are the most favorable sites for giant planet formation, can also deliver enough solids to the inner disk to form systems of super-Earths. We use these same models tomore »constrain the nature of Kepler-167's protoplanetary disk and find that it likely contained ≳300 M ⊕ of dust and was ≳40 au in size. These values overlap with the upper end of the observed dust mass and size distributions of Class 0 and I disks and are also consistent with the observed occurrence rate of Jupiter analogs around Sun-like stars.« less
  2. Abstract Isotope variations of nucleosynthetic origin among solar system solid samples are well documented, yet the origin of these variations is still uncertain. The observed variability of 54 Cr among materials formed in different regions of the protoplanetary disk has been attributed to variable amounts of presolar, chromium-rich oxide (chromite) grains, which exist within the meteoritic stardust inventory and most likely originated from some type of supernova explosion. To investigate if core-collapse supernovae (CCSNe) could be the site of origin of these grains, we analyze yields of CCSN models of stars with initial masses 15, 20, and 25 M ⊙ , and solar metallicity. We present an extensive abundance data set of the Cr, Mg, and Al isotopes as a function of enclosed mass. We find cases in which the explosive C ashes produce a composition in good agreement with the observed 54 Cr/ 52 Cr and 53 Cr/ 52 Cr ratios as well as the 50 Cr/ 52 Cr ratios. Taking into account that the signal at atomic mass 50 could also originate from 50 Ti, the ashes of explosive He burning also match the observed ratios. Addition of material from the He ashes (enriched in Al and Crmore »relative to Mg to simulate the make-up of chromite grains) to the solar system’s composition may reproduce the observed correlation between Mg and Cr anomalies, while material from the C ashes does not present significant Mg anomalies together with Cr isotopic variations. In all cases, nonradiogenic, stable Mg isotope variations dominate over the variations expected from 26 Al.« less
  3. Context. Dynamical models of Solar System evolution have suggested that the so-called P- and D-type volatile-rich asteroids formed in the outer Solar System beyond Neptune’s orbit and may be genetically related to the Jupiter Trojans, comets, and small Kuiper belt objects (KBOs). Indeed, the spectral properties of P- and D-type asteroids resemble that of anhydrous cometary dust. Aims. We aim to gain insights into the above classes of bodies by characterizing the internal structure of a large P- and D-type asteroid. Methods. We report high-angular-resolution imaging observations of the P-type asteroid (87) Sylvia with the Very Large Telescope Spectro-Polarimetric High-contrast Exoplanet REsearch (SPHERE) instrument. These images were used to reconstruct the 3D shape of Sylvia. Our images together with those obtained in the past with large ground-based telescopes were used to study the dynamics of its two satellites. We also modeled Sylvia’s thermal evolution. Results. The shape of Sylvia appears flattened and elongated (a/b ~1.45; a/c ~1.84). We derive a volume-equivalent diameter of 271 ± 5 km and a low density of 1378 ± 45 kg m −3 . The two satellites orbit Sylvia on circular, equatorial orbits. The oblateness of Sylvia should imply a detectable nodal precession which contrastsmore »with the fully-Keplerian dynamics of its two satellites. This reveals an inhomogeneous internal structure, suggesting that Sylvia is differentiated. Conclusions. Sylvia’s low density and differentiated interior can be explained by partial melting and mass redistribution through water percolation. The outer shell should be composed of material similar to interplanetary dust particles (IDPs) and the core should be similar to aqueously altered IDPs or carbonaceous chondrite meteorites such as the Tagish Lake meteorite. Numerical simulations of the thermal evolution of Sylvia show that for a body of such a size, partial melting was unavoidable due to the decay of long-lived radionuclides. In addition, we show that bodies as small as 130–150 km in diameter should have followed a similar thermal evolution, while smaller objects, such as comets and the KBO Arrokoth, must have remained pristine, which is in agreement with in situ observations of these bodies. NASA Lucy mission target (617) Patroclus (diameter ≈140 km) may, however, be differentiated.« less
  4. The slow neutron-capture process (s-process) efficiency in low-mass AGB stars (1.5 < M/M⊙ < 3) critically depends on how mixing processes in stellar interiors are handled, which is still affected by considerable uncertainties. In this work, we compute the evolution and nucleosynthesis of low-mass AGB stars at low metallicities using the MESA stellar evolution code. The combined data set includes models with initial masses Mini/M⊙=2 and 3 for initial metallicities Z=0.001 and 0.002. The nucleosynthesis was calculated for all relevant isotopes by post-processing with the NuGrid mppnp code. Using these models, we show the impact of the uncertainties affecting the main mixing processes on heavy element nucleosynthesis, such as convection and mixing at convective boundaries. We finally compare our theoretical predictions with observed surface abundances on low-metallicity stars. We find that mixing at the interface between the He-intershell and the CO-core has a critical impact on the s-process at low metallicities, and its importance is comparable to convective boundary mixing processes under the convective envelope, which determine the formation and size of the 13C-pocket. Additionally, our results indicate that models with very low to no mixing below the He-intershell during thermal pulses, and with a 13C-pocket size of at leastmore »∼3 × 10−4 M⊙, are strongly favored in reproducing observations. Online access to complete yield data tables is also provided.« less

    The hadron-quark phase transition in quantum chromodynamics has been suggested as an alternative explosion mechanism for core-collapse supernovae. We study the impact of three different hadron-quark equations of state (EoS) with first-order (DD2F_SF, STOS-B145) and second-order (CMF) phase transitions on supernova dynamics by performing 97 simulations for solar- and zero-metallicity progenitors in the range of $14\tt {-}100\, \text{M}_\odot$. We find explosions only for two low-compactness models (14 and $16\, \text{M}_\odot$) with the DD2F_SF EoS, both with low explosion energies of ${\sim }10^{50}\, \mathrm{erg}$. These weak explosions are characterized by a neutrino signal with several minibursts in the explosion phase due to complex reverse shock dynamics, in addition to the typical second neutrino burst for phase-transition-driven explosions. The nucleosynthesis shows significant overproduction of nuclei such as 90Zr for the $14\hbox{-} \text{M}_\odot$ zero-metallicity model and 94Zr for the $16\hbox{-}\text{M}_\odot$ solar-metallicity model, but the overproduction factors are not large enough to place constraints on the occurrence of such explosions. Several other low-compactness models using the DD2F_SF EoS and two high-compactness models using the STOS EoS end up as failed explosions and emit a second neutrino burst. For the CMF EoS, the phase transition never leads to a second bounce and explosion. For allmore »three EoS, inverted convection occurs deep in the core of the protocompact star due to anomalous behaviour of thermodynamic derivatives in the mixed phase, which heats the core to entropies up to 4kB/baryon and may have a distinctive gravitational-wave signature, also for a second-order phase transition.

    « less