Source author record

David J. Stevenson

David J. Stevenson appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

17works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

17 published item(s)

preprint2022arXiv

Revelations on Jupiter's Formation, Evolution and Interior: Challenges from Juno Results

The Juno mission has revolutionized and challenged our understanding of Jupiter. As Juno transitioned to its extended mission, we review the major findings of Jupiter's internal structure relevant to understanding Jupiter's formation and evolution. Results from Juno's investigation of Jupiter's interior structure imply that the planet has compositional gradients and is accordingly non-adiabatic, with a complex internal structure. These new results imply that current models of Jupiter's formation and evolution require a revision. In this paper, we discuss potential formation and evolution paths that can lead to an internal structure model consistent with Juno data, and the constraints they provide. We note that standard core accretion formation models, including the heavy-element enrichment during planetary growth is consistent with an interior that is inhomogeneous with composition gradients in its deep interior. However, such formation models typically predict that this region, which could be interpreted as a primordial dilute core, is confined to about 10% of Jupiter's total mass. In contrast, structure models that fit Juno data imply that this region contains 30% of the mass or more. One way to explain the origin of this extended region is by invoking a relatively long (about 2 Myrs) formation phase where the growing planet accretes gas and planetesimals delaying the runaway gas accretion. Alternatively, Jupiter's fuzzy core could be a result of a giant impact or convection post-formation. These novel scenarios require somewhat special and specific conditions. Clarity on the plausibility of such conditions could come from future high-resolution observations of planet-forming regions around other stars, from the observed and modeled architectures of extrasolar systems with giant planets, and future Juno data obtained during its extended mission.

preprint2022arXiv

The gravitational imprint of an interior-orbital resonance in Jupiter-Io

At mid-mission perijove 17, NASA's Juno mission has revealed a $7σ$ discrepancy between Jupiter's observed high-degree tidal response and the theoretical equilibrium tidal response, namely the Love number $k_{42}$. Here, we propose an interpretation for this puzzling disagreement based on an interior-orbital resonance between internal gravity waves trapped in Jupiter's dilute core and the orbital motion of Io. We use simple Jupiter models to calculate a fractional correction $Δk_{42}$ to the equilibrium tidal response that comes from the dynamical tidal response of a $g$-mode trapped in Jupiter's dilute core. Our results suggest that an extended dilute core ($r\gtrsim0.7R_J$) produces an interior-orbital resonance with Io that modifies Jupiter's tidal response in $Δk_{42}\sim-11\%$, allowing us to fit Juno's $k_{42}$. In our proposed self-consistent scenario, Jupiter's dilute core evolves in resonant locking with Io's orbital migration, which allows the interior-orbital resonance to persist over geological timescales. This scenario requires a dilute core that becomes smoother or shrinks over time, together with a $_4^2g_1$ mode ($\ell,m,n=4,2,1$) with resonant tidal dissipation reaching $Q_4\sim1000$. Jupiter's dilute core evolution path and the dissipation mechanism for the resonant $_4^2g_1$ mode are uncertain and motivate future analysis. No other alternative exists so far to explain the $7σ$ discrepancy in Juno $k_{42}$. Our proposed interior-orbital resonance can be tested by Juno observations of $k_{42}$ tides raised on Jupiter by Europa as obtained at the end of the extended mission (mid 2025), and by future seismological observations of Jupiter's $_4^2g_1$ mode oscillation frequency.

preprint2021arXiv

The lost meaning of Jupiter's high-degree Love numbers

NASA's Juno mission recently reported Jupiter's high-degree (degree $\ell$, azimuthal order $m$ $=4,2$) Love number $k_{42}=1.289\pm0.063$ ($1σ$), an order of magnitude above the hydrostatic $k_{42}$ obtained in a nonrotating Jupiter model. After numerically modeling rotation, the hydrostatic $k_{42}=1.743\pm0.002$ is still $7σ$ away from the observation, raising doubts about our understanding of Jupiter's tidal response. Here, we use first-order perturbation theory to explain the hydrostatic $k_{42}$ result analytically. We use a simple Jupiter equation of state ($n=1$ polytrope) to obtain the fractional change in $k_{42}$ when comparing a rotating model with a nonrotating model. Our analytical result shows that the hydrostatic $k_{42}$ is dominated by the tidal response at $\ell=m=2$ coupled into the spherical harmonic $\ell,m=4,2$ by the planet's oblate figure. The $\ell=4$ normalization in $k_{42}$ introduces an orbital factor $(a/s)^2$ into $k_{42}$, where $a$ is the satellite semimajor axis and $s$ is Jupiter's average radius. As a result, different Galilean satellites produce a different $k_{42}$. We conclude that high-degree tesseral Love numbers ($\ell> m$, $m\geq2$) are dominated by lower-degree Love numbers and thus provide little additional information about interior structure, at least when they are primarily hydrostatic. Our results entail important implications for a future interpretation of the currently observed Juno $k_{42}$. After including the coupling from the well-understood $\ell=2$ dynamical tides ($Δk_2 \approx -4\%$), Jupiter's hydrostatic $k_{42}$ requires an unknown dynamical effect to produce a fractional correction $Δk_{42}\approx-11\%$ in order to fit Juno's observation within $3σ$. Future work is required to explain the required $Δk_{42}$.

preprint2021arXiv

Thermodynamically Governed Interior Models of Uranus and Neptune

Interior models of Uranus and Neptune often assume discrete layers, but sharp interfaces are expected only if major constituents are immiscible. Diffuse interfaces could arise if accretion favored a central concentration of the least volatile constituents (also incidentally the most dense); compositional gradients arising in such a structure would likely inhibit convection. Currently, two lines of evidence suggest possible hydrogen-water immiscibility in ice giant interiors. The first arises from crude extrapolation of the experimental hydrogen-water critical curve to $\sim 3$ GPa (Bali et al. 2013). The data are obtained for an impure system containing silicates, though Uranus and Neptune could also be "dirty." Current ab initio models disagree (Soubiran & Militzer 2015), though hydrogen and water are difficult to model from first-principles quantum mechanics with the necessary precision. The second argument for hydrogen-water immiscibility in ice giants, outlined herein, invokes reasoning about the gravitational and magnetic fields. While consensus remains lacking, here we examine the immiscible case. Applying the resulting thermodynamic constraints, we find that Neptune models with envelopes containing a substantial water mole fraction, as much as $χ\gtrsim 0.1$ relative to hydrogen, can satisfy observations. In contrast, Uranus models appear to require $χ\lesssim 0.01$, potentially suggestive of fully demixed hydrogen and water. Enough gravitational potential energy would be available from gradual hydrogen-water demixing, to supply Neptune's present-day heatflow for roughly ten solar system lifetimes. Hydrogen-water demixing could slow Neptune's cooling rate by an order of magnitude; different hydrogen-water demixing states could account for the different heatflows of Uranus and Neptune.

preprint2019arXiv

Comparison of the deep atmospheric dynamics of Jupiter and Saturn in light of the Juno and Cassini gravity measurements

The nature and structure of the observed east-west flows on Jupiter and Saturn has been one of the longest-lasting mysteries in planetary science. This mystery has been recently unraveled due to the accurate gravity measurements provided by the Juno mission to Jupiter and the Grand Finale of the Cassini mission to Saturn. These two experiments, which coincidentally happened around the same time, allowed determination of the vertical and meridional profiles of the zonal flows on both planets. This paper reviews the topic of zonal jets on the gas giants in light of the new data from these two experiments. The gravity measurements not only allow the depth of the jets to be constrained, yielding the inference that the jets extend roughly 3000 and 9000 km below the observed clouds on Jupiter and Saturn, respectively, but also provide insights into the mechanisms controlling these zonal flows. Specifically, for both planets this depth corresponds to the depth where electrical conductivity is within an order of magnitude of 1 S/m, implying that the magnetic field likely plays a key role in damping the zonal flows.

preprint2016arXiv

Analysis of a Precambrian resonance-stabilized day length

During the Precambrian era, Earth's decelerating rotation would have passed a 21-hour period that would have been resonant with the semidiurnal atmospheric thermal tide. Near this point, the atmospheric torque would have been maximized, being comparable in magnitude but opposite in direction to the lunar torque, halting Earth's rotational deceleration, maintaining a constant day length, as detailed by Zahnle and Walker (1987). We develop a computational model to determine necessary conditions for formation and breakage of this resonant effect. Our simulations show the resonance to be resilient to atmospheric thermal noise but suggest a sudden atmospheric temperature increase like the deglaciation period following a possible "snowball Earth" near the end of the Precambrian would break this resonance; the Marinoan and Sturtian glaciations seem the most likely candidates for this event. Our model provides a simulated day length over time that resembles existing paleorotational data, though further data is needed to verify this hypothesis.

preprint2016arXiv

The Interior of Saturn

We review our current understanding of the interior structure and thermal evolution of Saturn, with a focus on recent results in the Cassini era. There has been important progress in understanding physical inputs, including equations of state of planetary materials and their mixtures, physical parameters like the gravity field and rotation rate, and constraints on Saturnian free oscillations. At the same time, new methods of calculation, including work on the gravity field of rotating fluid bodies, and the role of interior composition gradients, should help to better constrain the state of Saturn's interior, now and earlier in its history. However, a better appreciation of modeling uncertainties and degeneracies, along with a greater exploration of modeling phase space, still leave great uncertainties in our understanding of Saturn's interior. Further analysis of Cassini data sets, as well as precise gravity field measurements from the Cassini Grand Finale orbits, will further revolutionize our understanding of Saturn's interior over the next few years.

preprint2015arXiv

Large Impacts around a Solar Analog Star in the Era of Terrestrial Planet Formation

The final assembly of terrestrial planets occurs via massive collisions, which can launch copious clouds of dust that are warmed by the star and glow in the infrared. We report the real-time detection of a debris-producing impact in the terrestrial planet zone around a 35-million year-old solar analog star. We observed a substantial brightening of the debris disk at 3-5 μm, followed by a decay over a year, with quasi-periodic modulations of the disk flux. The behavior is consistent with the occurrence of a violent impact that produced vapor out of which a thick cloud of silicate spherules condensed that were ground into dust by collisions. These results demonstrate how the time domain can become a new dimension for the study of terrestrial planet formation.

preprint2015arXiv

Melting and Mixing States of the Earth's Mantle after the Moon-Forming Impact

The Earth's Moon is thought to have formed by an impact between the Earth and an impactor around 4.5 billion years ago. This impact could have been so energetic that it could have mixed and homogenized the Earth's mantle. However, this view appears to be inconsistent with geochemical studies that suggest that the Earth's mantle was not mixed by the impact. Another plausible outcome is that this energetic impact melted the whole mantle, but the extent of mantle melting is not well understood even though it must have had a significant effect on the subsequent evolution of the Earth's interior and atmosphere. To understand the initial state of the Earth's mantle, we perform giant impact simulations using smoothed particle hydrodynamics (SPH) for three different models: (a) standard: a Mars-sized impactor hits the proto-Earth, (b) fast-spinning Earth: a small impactor hits a rapidly rotating proto-Earth, and (c) sub-Earths: two half Earth-sized planets collide. We use two types of equations of state (MgSiO3 liquid and forsterite) to describe the Earth's mantle. We find that the mantle remains unmixed in (a), but it may be mixed in (b) and (c). The extent of mixing is most extensive in (c). Therefore, (a) is most consistent and (c) may be least consistent with the preservation of the mantle heterogeneity, while (b) may fall between. We determine that the Earth's mantle becomes mostly molten by the impact in all of the models. The choice of the equations of state does not affect these outcomes. Additionally, our results indicate that entropy gains of the mantle materials by a giant impact cannot be predicted well by the Rankine-Hugoniot equations. Moreover, we show that the mantle can remain unmixed on a Moon-forming timescale if it does not become mixed by the impact.

preprint2014arXiv

Investigation of the Initial State of the Moon-Forming Disk: Bridging SPH Simulations and Hydrostatic Models

According to the standard giant impact hypothesis, the Moon formed from a partially vaporized disk generated by a collision between the proto Earth and a Mars sized impactor. The initial structure of the disk significantly affects the Moon forming process, including the Moons mass, its accretion time scale, and its isotopic similarity to Earth. The dynamics of the impact event determines the initial structure of a nearly hydrostatic Moon forming disk. However, the hydrostatic and hydrodynamic models have been studied separately and their connection has not previously been well quantified. Here, we show the extent to which the properties of the disk can be inferred from Smoothed Particle Hydrodynamic (SPH) simulations. By using entropy, angular momentum and mass distributions of the SPH outputs as approximately conserved quantities, we compute the two dimensional disk structure. We investigate four different models: (a) standard, the canonical giant impact model, (b) fast spinning Earth, a collision between a fast spinning Earth and a small impactor, (c) sub Earths, a collision between two objects with half Earths mass, and (d) intermediate, a collision of two bodies whose mass ratio is 7:3. Our SPH calculations show that the initial disk has approximately uniform entropy. The disks of the fast spinning Earth and sub Earths cases are hotter and more vaporized (80-90% vapor) than the standard case (20%). The intermediate case falls between these values. In the highly vaporized cases, our procedure fails to establish a unique surface density profile of the disk because the disk is unstable according to the Rayleigh criterion. In these cases, we estimate non-unique disk models by conserving global quantities. We also develop a semi analytic model for the thermal structure of the disk, which requires only two inputs: the average entropy and the surface density of the disk.

preprint2013arXiv

Magnetically Controlled Circulation on Hot Extrasolar Planets

Through the process of thermal ionization, intense stellar irradiation renders Hot Jupiter atmospheres electrically conductive. Simultaneously, lateral variability in the irradiation drives the global circulation with peak wind speeds of order ~ km/s. In turn, the interactions between the atmospheric flows and the background magnetic field give rise to Lorentz forces that can act to perturb the flow away from its purely hydrodynamical counterpart. Using analytical theory and numerical simulations, here we show that significant deviations away from axisymmetric circulation are unstable in presence of a non-negligible axisymmetric magnetic field. Specifically, our results suggest that dayside-to-nightside flows, often obtained within the context of three-dimensional circulation models, only exist on objects with anomalously low magnetic fields, while the majority of highly irradiated exoplanetary atmospheres are entirely dominated by zonal jets.

preprint2013arXiv

Mass-Radius Relationships for Very Low Mass Gaseous Planets

Recently, the Kepler spacecraft has detected a sizable aggregate of objects, characterized by giant-planet-like radii and modest levels of stellar irradiation. With the exception of a handful of objects, the physical nature, and specifically the average densities, of these bodies remain unknown. Here, we propose that the detected giant planet radii may partially belong to planets somewhat less massive than Uranus and Neptune. Accordingly, in this work, we seek to identify a physically sound upper limit to planetary radii at low masses and moderate equilibrium temperatures. As a guiding example, we analyze the interior structure of the Neptune-mass planet Kepler-30d and show that it is acutely deficient in heavy elements, especially compared with its solar system counterparts. Subsequently, we perform numerical simulations of planetary thermal evolution and in agreement with previous studies, show that generally, 10 - 20 Earth mass, multi-billion year old planets, composed of high density cores and extended H/He envelopes can have radii that firmly reside in the giant planet range. We subject our results to stability criteria based on extreme ultraviolet radiation, as well as Roche-lobe overflow driven mass-loss and construct mass-radius relationships for the considered objects. We conclude by discussing observational avenues that may be used to confirm or repudiate the existence of putative low mass, gas-dominated planets.

preprint2013arXiv

Nonhydrostatic Effects and the Determination of Icy Satellites' Moment of Inertia

We compare the moment of inertia (MOI) of a simple hydrostatic, two layer body as determined by the Radau-Darwin Approximation (RDA) to its exact hydrostatic MOI calculated to first order in the parameter q = w^2R^3/GM, where w, R, and M are the spin angular velocity, radius, and mass of the body, and G is the gravitational constant. RDA is in error by less than 1% for many configurations of core sizes and layer densities congruent with those of solid bodies in the Solar System. We determine the error in the MOI of icy satellites calculated with the RDA due to nonhydrostatic effects by using a simple model in which the core and outer shell have slight degree 2 distortions away from their expected hydrostatic shapes. Since the hydrostatic shape has an associated stress of order pw^2R^2 (where p is density) it follows that the importance of nonhydrostatic effects scales with the dimensionless number s/pw^2R^2, where s is the nonhydrostatic stress. This highlights the likely importance of this error for slowly rotating bodies (e.g., Titan and Callisto) and small bodies (e.g., Saturn moons other than Titan). We apply this model to Titan, Callisto, and Enceladus and find that the RDA-derived MOI can be 10% greater than the actual MOI for nonhydrostatic stresses as small as ~0.1 bars at the surface or ~1 bar at the core-mantle boundary, while for Ganymede the stresses necessary to produce the same MOI errors are an order of magnitude greater due to its faster rotation. If satellites can reorient to the lowest energy state then RDA will always give an overestimate of the true MOI. Observations have shown that small nonhydrostatic gravity anomalies exist on Ganymede and Titan. We conclude that nonhydrostatic effects could be present to an extent that allows Callisto and Titan to be fully differentiated.

preprint2013arXiv

Stability of Ice/Rock Mixtures with Application to a Partially Differentiated Titan

Titan's moment of inertia, calculated assuming hydrostatic equilibrium from gravity field data obtained during the Cassini-Huygens mission, implies an internal mass distribution that may be incompatible with complete differentiation. This suggests that Titan may have a mixed ice/rock core, possibly consistent with slow accretion in a gas-starved disk, which may initially spare Titan from widespread ice melting and subsequent differentiation. A partially differentiated Titan, however, must still efficiently remove radiogenic heat over geologic time. We argue that compositional heterogeneity in the major Saturnian satellites indicates that Titan formed from planetesimals with disparate densities. The resulting compositional anomalies would quickly redistribute to form a vertical density gradient that would oppose thermal convection. We use elements of the theory of double-diffusive convection to create a parameterized model for the thermal evolution of ice/rock mixtures with a stabilizing compositional gradient. Simulations are performed for a wide range of initial conditions to account for large uncertainties in material properties and accretionary processes. Ultimately, for realistic density gradients, double-diffusive convection in the ice/rock interior can delay, but not prevent, ice melting and differentiation, even if a substantial fraction of potassium is leached from the rock component. Consequently, Titan is not partially differentiated.

preprint2012arXiv

Evolution of Ohmically Heated Hot Jupiters

We present calculations of thermal evolution of Hot Jupiters with various masses and effective temperatures under Ohmic dissipation. The resulting evolutionary sequences show a clear tendency towards inflated radii for effective temperatures that give rise to significant ionization of alkali metals in the atmosphere, compatible with the trend of the data. The degree of inflation shows that Ohmic dissipation, along with the likely variability in heavy element content can account for all of the currently detected radius anomalies. Furthermore, we find that in absence of a massive core, low-mass hot Jupiters can over-flow their Roche-lobes and evaporate on Gyr time-scales, possibly leaving behind small rocky cores.

preprint2011arXiv

Jupiter's Moment of Inertia: A Possible Determination by JUNO

The moment of inertia of a giant planet reveals important information about the planet's internal density structure and this information is not identical to that contained in the gravitational moments. The forthcoming Juno mission to Jupiter might determine Jupiter's normalized moment of inertia NMoI=C/MR^2 by measuring Jupiter's pole precession and the Lense-Thirring acceleration of the spacecraft (C is the axial moment of inertia, and M and R are Jupiter's mass and mean radius, respectively). We investigate the possible range of NMoI values for Jupiter based on its measured gravitational field using a simple core/envelope model of the planet assuming that J_2 and J_4 are perfectly known and are equal to their measured values. The model suggests that for fixed values of J_2 and J_4 a range of NMOI values between 0.2629 and 0.2645 can be found. The Radau-Darwin relation gives a NMoI value that is larger than the model values by less than 1%. A low NMoI of ~ 0.236, inferred from a dynamical model (Ward & Canup, 2006, ApJ, 640, L91) is inconsistent with this range, but the range is model dependent. Although we conclude that the NMoI is tightly constrained by the gravity coefficients, a measurement of Jupiter's NMoI to a few tenths of percent by Juno could provide an important constraint on Jupiter's internal structure. We carry out a simplified assessment of the error involved in Juno's possible determination of Jupiter's NMoI.

preprint2010arXiv

Inflating Hot Jupiters With Ohmic Dissipation

We present a new, magnetohydrodynamic mechanism for inflation of close-in giant extrasolar planets. The idea behind the mechanism is that current, which is induced through interaction of atmospheric winds and the planetary magnetic field, results in significant Ohmic dissipation of energy in the interior. We develop an analytical model for computation of interior Ohmic dissipation, with a simplified treatment of the atmosphere. We apply our model to HD209458b, Tres-4b and HD189733b. With conservative assumptions for wind speed and field strength, our model predicts a generated power that appears to be large enough to maintain the transit radii, opening an unexplored avenue towards solving a decade-old puzzle of extrasolar gas giant radius anomalies.