Search NASA⌕ Search

SEARCH · Search NASA

Results for “planetary system evolution”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Nonlinear Modeling of Macroscopic Plasma Dynamics Final Scientific Report May 1, 2006 – April 30, 2018

The project entitled, "Nonlinear Modeling of Macroscopic Plasma Dynamics," applied numerical computation to solve mathematical models of plasma behavior in laboratory systems. Plasma is the state of matter with sufficiently large internal energy density such that electrons dissociate from atoms to form gas of electrically charged particles. It occurs naturally in stars, planetary magnetospheres, and interstellar and intergalactic media. It is also used for industrial purposes and in lighting. Fusion energy science seeks to create and sustain plasma conditions where atomic nuclei fuse to create energy that can be captured for power production. Many laboratory experiments for this effort use magnetic field to contain plasma and to thermally insulate it from its surroundings. Possible plasma dynamics in magnetized systems include collective behavior over spatial scales that are comparable to the configuration size, and these dynamics often include evolution of the magnetic field that is used for confinement. Solving the mathematical models that describe these dynamics in specific experiments often requires computer simulations, as done in this project. This project analyzed plasma magnetic relaxation dynamics in the laboratory configurations known as the reversed-field pinch (RFP), the spheromak, and the spherical torus (ST). It also developed capabilities for modeling disruptive plasma dynamics in the tokamak, which has achieved the highest performance for magnetized plasma and is the design of the large international ITER experiment (https://www.iter.org) that is under construction and is expected to be the first to achieve self-heating conditions. Dynamic evolution to states of low energy, the aforementioned magnetic relaxation, is known to be important for generating the RFP and spheromak configurations. Nonlinear computational modeling with the Non-Ideal Magnetohydrodynamics with Rotation, Open Discussion code (NIMROD, https://nimrodteam.org) was applied to understand specific aspects of magnetic relaxation in RFPs, spheromaks, and STs. The RFP study considered how the nonlinear transfer of power among different fluctuation harmonics is altered due to the application of pulses of electric field, as done experimentally in the Madison Symmetric Torus RFP device at the University of Wisconsin-Madison. It also investigated the separation of electron and ion dynamics, known as two-fluid effects, including effects on the evolution of the plasma flow profile. The spheromak study considered the influence of these two-fluid effects in conditions relevant to the Sustained Spheromak Experiment (SSPX), formerly at Lawrence Livermore National Laboratory. Our study found that two-fluid effects enhance stability with respect to interchange in SSPX-relevant conditions, but having a complete model with two-fluid effects in separate electron and ion temperature equations is important. The study for STs focused on plasma startup in the Pegasus Toroidal Experiment at UW-Madison, which uses DC voltage applied through a small plasma gun. The simulations show how the resulting helical current stream relaxes through repetitive merging, forming current rings that accumulate into a tokamak-like configuration. Analytical theory was also developed to help explain behavior that was observed in the RFP and spheromak computations. Although the tokamak configuration has outperformed other magnetic confinement configurations, it has a tendency to undergo discharge-terminating dynamics that release energy onto plasma-facing surfaces and exert electromechanical forces on magnetic coils and structures. Our efforts have improved NIMROD's capabilities to model tokamak disruption by incorporating a resistive wall model that can be applied to axisymmetric, as well as asymmetric, distortions of the magnetic field configuration. Computations addressed whether mathematical conditions applied to plasma flow at the wall influence disruptive dynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

ORBITAL STABILITY OF MULTI-PLANET SYSTEMS: BEHAVIOR AT HIGH MASSES

In the coming years, high-contrast imaging surveys are expected to reveal the characteristics of the population of wide-orbit, massive, exoplanets. To date, a handful of wide planetary mass companions are known, but only one such multi-planet system has been discovered: HR 8799. For low mass planetary systems, multi-planet interactions play an important role in setting system architecture. In this paper, we explore the stability of these high mass, multi-planet systems. While empirical relationships exist that predict how system stability scales with planet spacing at low masses, we show that extrapolating to super-Jupiter masses can lead to up to an order of magnitude overestimate of stability for massive, tightly packed systems. We show that at both low and high planet masses, overlapping mean-motion resonances trigger chaotic orbital evolution, which leads to system instability. We attribute some of the difference in behavior as a function of mass to the increasing importance of second order resonances at high planet–star mass ratios. We use our tailored high mass planet results to estimate the maximum number of planets that might reside in double component debris disk systems, whose gaps may indicate the presence of massive bodies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

An Eccentric Massive Jupiter Orbiting a Subgiant on a 9.5-day Period Discovered in the Transiting Exoplanet Survey Satellite Full Frame Images

We report the discovery of TOI-172 b from the Transiting Exoplanet Survey Satellite (TESS) mission, a massive hot Jupiter transiting a slightly evolved G star with a 9.48-day orbital period. This is the first planet to be confirmed from analysis of only the TESS full frame images, because the host star was not chosen as a two-minute cadence target. From a global analysis of the TESS photometry and follow-up observations carried out by the TESS Follow-up Observing Program Working Group, TOI-172 (TIC 29857954) is a slightly evolved star with an effective temperature of T {sub eff} = 5645 ± 50 K, a mass of M {sub ⋆} = 1.128{sub −0.061}{sup +0.065} M {sub ⊙}, radius of R {sub ⋆} = 1.777{sub −0.044}{sup +0.047} R {sub ⊙}, a surface gravity of log g {sub ⋆} = 3.993{sub −0.028}{sup +0.027}, and an age of 7.4{sub −1.5}{sup +1.6} Gyr. Its planetary companion (TOI-172 b) has a radius of R {sub P} = 0.965{sub −0.029}{sup +0.032} R {sub J}, a mass of M {sub P} = 5.42{sub −0.20}{sup +0.22} M {sub J}, and is on an eccentric orbit (e=0.3806{sub −0.0090}{sup +0.0093}). TOI-172 b is one of the few known massive giant planets on a highly eccentric short-period orbit. Future study of the atmosphere of this planet and its system architecture offer opportunities to understand the formation and evolution of similar systems.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The homogeneous mixing of MgO and H 2 O at extreme conditions

Investigating water worlds presents a unique opportunity to understand the fundamental processes of planetary formation and evolution. One key aspect is characterizing the interactions between water and rock under the pressures and temperatures present within these worlds. Investigating the conditions for the homogeneous mixing of these materials is imperative to characterizing bulk properties and evolution of water-rich exoplanets. Here, in this study, we use density functional molecular dynamics simulations to study MgO-H 2 O mixtures at high pressure–temperature conditions where H 2 O occurs in solid, superionic, or liquid form. MgO, the representative rocky material, can be either solid or liquid. We start from 500 K at 120 GPa, increasing the temperature step by step up to 8000 K. By inspection, we determine the temperature at which MgO-HO homogeneously mix in our simulations. At 6000 K and 174 GPa is when we find the system to homogeneously mix. This heat-until-it-mixes approach provides us with an upper bound on the temperature for the mixing of MgO and H 2 O. We find that homogeneous mixing occurs at sufficiently low temperatures to be relevant for the collisional growth of a water-rich planet.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Theoretical studies of chemical reactions related to the formation and growth of polycyclic aromatic hydrocarbons (PAH) and molecular properties of their key intermediates (Final Progress Report)

The formation mechanisms of polycyclic aromatic hydrocarbons, (PAHs) – organic molecules carrying fused benzene rings – are of great interest to scientists and engineers due to their importance in combustion chemistry and astrochemistry. On Earth, PAHs are largely produced in incomplete combustion of fossil fuel and are considered as critical precursors to unwanted soot particles leading to combustion inefficiency and causing air pollution along with detrimental health effects. Simple PAH molecules initially formed in the gas phase, are further involved in a build-up process in combustion flames leading to larger PAH, bowl-shaped nanostructures, fullerenes, and solid-phase species including carbonaceous dust, graphene particles, and soot. In deep space, PAH and their derivatives are potential key intermediates and nucleation sites leading eventually to carbonaceous nanoparticles (“interstellar grains”). Therefore, the understanding of the key processes in the synthesis of PAHs along with their precursors and their degradation mechanisms in combustion systems and in interstellar, circumstellar, and planetary atmospheric environments will provide critical insights into how complex aromatic structures, carbonaceous nanoparticles, and fullerenes are formed and destroyed. Achieving this understanding is an important step in the development of the efficient combustion processes and of the ecofriendly devices with reduced environmental pollution as well as technological strategies for the production of hydrogen and solid carbon through thermal or plasma-assisted pyrolysis of natural gas and biomass. Also, the understanding of the key processes of PAH and soot growth will help in our comprehension of chemical evolution in the universe. Detailed information on the mechanisms and reliable rate constants of the key elementary chemical reactions involved in PAH formation and destruction processes and in inception of soot particles is often missing, with the main deficiencies being the absence of temperature- and pressure-dependent rate constants for the broad range of conditions occurring in various terrestrial and interstellar processes and the lack of data on the reaction products and their branching ratios. Complementary to experimental studies, these gaps in knowledge can be filled by using quantum chemical calculations of reaction potential energy surfaces providing us with accurate energies of reaction products, intermediates, and transition states, revealing the reaction mechanism, and giving the molecular properties required to compute rate constants for relevant reaction steps and product branching ratios using the RRKM-Master Equation (ME) method. Molecular dynamics (MD) simulations can be used in cases when a reaction rate cannot be properly described by statistical theories. During the terminal renewal project period we employed these ab initio/RRKM-ME and MD approaches to complete our studies on several key reactions relevant to the formation/growth of PAH and inception of soot particles including (1) the reaction mechanism and kinetics of the resonance stabilized fulvenallenyl radical with propargyl and C 3 H 4 isomers; (2) the reaction mechanism and kinetics for the C + indene and C 2 + styrene reactions producing naphthyl or azulenyl radicals in low-temperature environments; (3) the MD study of non-equilibrium dimerization of acepyrene and coronene and its radical. The information derived from our theoretical calculations contributed to a better fundamental understanding of the reaction mechanisms and provide missing critical kinetic data to improve combustion models of hydrocarbon fuels and astrochemical models of the growth of carbonaceous molecules and particles in cold molecular clouds, circumstellar envelopes, and planetary atmospheres.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Understanding the Impacts of Stellar Companions on Planet Formation and Evolution: A Survey of Stellar and Planetary Companions within 25 pc

We explore the impact of outer stellar companions on the occurrence rate of giant planets detected with radial velocities. We searched for stellar and planetary companions to a volume-limited sample of solar-type stars within 25 pc. Using adaptive optics imaging observations from the Lick 3 m and Palomar 200″ Telescopes, we characterized the multiplicity of our sample stars, down to the bottom of the main sequence. With these data, we confirm field star multiplicity statistics from previous surveys. We additionally combined three decades of radial velocity (RV) data from the California Planet Search with newly collected RV data from Keck/HIRES and the Automated Planet Finder/Levy Spectrometer to search for planetary companions in these same systems. Using an updated catalog of both stellar and planetary companions, as well as detailed injection/recovery tests to determine our sensitivity and completeness, we measured the occurrence rate of planets among the single- and multiple-star systems. We found that planets with masses in the range of 0.1–10 M {sub J} and with semimajor axes of 0.1–10 au have an occurrence rate of 0.18{sub −0.03}{sup +0.04} planets per star when they orbit single stars and an occurrence rate of 0.12 ± 0.04 planets per star when they orbit a star in a binary system. Breaking the sample down by the binary separation, we found that only one planet-hosting binary system had a binary separation <100 au, and none had a separation <50 au. These numbers yielded planet occurrence rates of 0.20{sub −0.06}{sup +0.07} planets per star for binaries with separation a {sub B} > 100 au and 0.04{sub −0.02}{sup +0.04} planets per star for binaries with separation a {sub B} < 100 au. The similarity in the planet occurrence rate around single stars and wide primaries implies that wide binary systems should actually host more planets than single-star systems, since they have more potential host stars. We estimated a system-wide planet occurrence rate of 0.3 planets per wide binary system for binaries with separations a {sub B} > 100 au. Finally, we found evidence that giant planets in binary systems have a different semimajor-axis distribution than their counterparts in single-star systems. The planets in the single-star sample had a significantly higher occurrence rate outside of 1 au than inside 1 au by nearly 4σ, in line with expectations that giant planets are most common near the snow line. However, the planets in the wide binary systems did not follow this distribution, but rather had equivalent occurrence rates interior and exterior to 1 au. This may point to binary-mediated planet migration acting on our sample, even in binaries wider than 100 au.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2019

Instabilities in the Early Solar System Due to a Self-gravitating Disk

Modern studies of the early solar system routinely invoke the possibility of an orbital instability among the giant planets triggered by gravitational interactions between the planets and a massive exterior disk of planetesimals. Previous works have suggested that this instability can be substantially delayed (about hundreds of megayears) after the formation of the giant planets. Bodies in the disk are typically treated in a semi-active manner, wherein their gravitational force on the planets is included, but interactions between the planetesimals are ignored. We perform N-body numerical simulations using GENGA, which makes use of GPUs to allow for the inclusion of all gravitational interactions between bodies. Although our simulated Kuiper Belt particles are more massive than the probable masses of real primordial Kuiper Belt objects, our simulations indicate that the self-stirring of the primordial Kuiper Belt is very important to the dynamics of the giant planet instability. We find that interactions between planetesimals dynamically heat the disk and typically prevent the outer solar system’s instability from being delayed by more than a few tens of megayears after giant planet formation. Longer delays occur in a small fraction of systems that have at least 3.5 au gaps between the planets and planetesimal disk. Our final planetary configurations match the solar system at a rate consistent with other previous works in most regards. Pre-instability heating of the disk typically yields final Jovian eccentricities comparable to the modern solar system’s value, which has been a difficult constraint to match in past works.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Northern Arabian Sea Circulation-Autonomous Research (NASCar): A Research Initiative Based on Autonomous Sensors

The Arabian Sea circulation is forced by strong monsoonal winds and is characterized by vigorous seasonally reversing currents, extreme differences in sea surface salinity, localized substantial upwelling, and widespread submesoscale thermohaline structures. Its complicated sea surface temperature patterns are important for the onset and evolution of the Asian monsoon. This article describes a program that aims to elucidate the role of upper-ocean processes and atmospheric feedbacks in setting the sea surface temperature properties of the region. The wide range of spatial and temporal scales and the difficulty of accessing much of the region with ships due to piracy motivated a novel approach based on state-of-the-art autonomous ocean sensors and platforms. The extensive data set that is being collected, combined with numerical models and remote sensing data, confirms the role of planetary waves in the reversal of the Somali Current system. These data also document the fast response of the upper equatorial ocean to monsoon winds through changes in temperature and salinity and the connectivity of the surface currents across the northern Indian Ocean. New observations of thermohaline interleaving structures and mixing in setting the surface temperature properties of the northern Arabian Sea are also discussed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Collaborative Project: Ocean-Atmosphere Interaction From Meso- to Planetary-Scale: Mechanisms, Parameterization, and Variability

The project focused on the dynamics and scale dependence of the boundary layer response, the coupled responses in the Pacific, and low frequency oceanic processes. Using theory and model simulations, we found coupling of sea surface temperatures associated with the ocean eddy field and the atmospheric boundary layer be dependent on scale, and captured in wave-number space by transfer functions. Investigations of coupled processes focused on the sea surface temperatures scales of eastern boundary upwelling system, and showed stability and eddy variance in the Kuroshio Extension to be predictable with lead times of several years through ocean Rossby waves and coupling with the atmosphere. In addition to these dynamical processes, we found that a significant part of the ocean evolution is associated with density compensating temperature and salinity anomalies, that, unlike Rossby waves, behave like passive tracers and follow the mean circulation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

SCExAO and GPI Y JH band photometry and integral field spectroscopy of the young brown dwarf companion to HD 1160

Here, we present high signal-to-noise ratio, precise Y JH photometry and Y band (0.957–1.120 μm) spectroscopy of HD 1160 B, a young substellar companion discovered from the Gemini NICI Planet Finding Campaign using the Subaru Coronagraphic Extreme Adaptive Optics instrument and the Gemini Planet Imager. HD 1160 B has typical mid-M dwarf-like infrared colors and a spectral type of M5.5${}_{-0.5}^{+1.0}$, where the blue edge of our Y band spectrum rules out earlier spectral types. Atmospheric modeling suggests HD 1160 B has an effective temperature of 3000–3100 K, a surface gravity of log g = 4–4.5, a radius of 1.55 ± 0.10 R J, and a luminosity of log L/L ⊙ = –2.76 ± 0.05. Neither the primary's Hertzspring–Russell diagram position nor atmospheric modeling of HD 1160 B show evidence for a subsolar metallicity. Interpretation of the HD 1160 B spectroscopy depends on which stellar system components are used to estimate the age. Considering HD 1160 A, B and C jointly, we derive an age of 80–125 Myr, implying that HD 1160 B straddles the hydrogen-burning limit (70–90 M J ). If we consider HD 1160 A alone, younger ages (20–125 Myr) and a brown dwarf-like mass (35–90 M J ) are possible. Interferometric measurements of the primary, a precise Gaia parallax, and moderate-resolution spectroscopy can better constrain the system's age and how HD 1160 B fits within the context of (sub)stellar evolution.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Hydrodynamics and Survivability during Post-main-sequence Planetary Engulfment

The engulfment of substellar bodies (SBs), such as brown dwarfs and planets, by giant stars is a possible explanation for rapidly rotating giants, lithium-rich giants, and the presence of SBs in close orbits around subdwarfs and white dwarfs. We perform three-dimensional hydrodynamical simulations of the flow in the vicinity of an engulfed SB. We model the SB as a rigid body with a reflective surface because it cannot accrete. This reflective boundary changes the flow morphology to resemble that of engulfed compact objects with outflows. We measure the drag coefficients for the ram-pressure and gravitational drag forces acting on the SB, and use them to integrate its trajectory inside the star. We find that engulfment can increase the luminosity of a 1 M ⊙ star by up to a few orders of magnitude. The time for the star to return to its original luminosity is up to a few thousand years when the star has evolved to ≈10 R ⊙ and up to a few decades at the tip of the red giant branch (RGB). No SBs can eject the envelope of a 1 M ⊙ star before it evolves to ≈10 R ⊙ if the orbit of the SB is the only energy source contributing to the ejection. In contrast, SBs as small as ≈10 MJup can eject the envelope at the tip of the RGB. The numerical framework we introduce here can be used to study planetary engulfment in a simplified setting that captures the physics of the flow at the scale of the SB.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Asteroids seen by JWST-MIRI: Radiometric size, distance, and orbit constraints

Infrared measurements of asteroids are crucial for the determination of physical and thermal properties of individual objects, and for understanding the small-body populations in the solar system as a whole. However, standard radiometric methods can only be applied if the orbit of an object is known, hence its position at the time of the observation. With JWST-MIRI observations the situation will change and many unknown, often very small, solar system objects will be detected. Later orbit determinations are difficult due to the faintness of the objects and the lack of dedicated follow-up concepts. We present MIRI observations of the outer-belt asteroid (10920) 1998 BC1 and an unknown object, detected in all nine MIRI bands in close apparent proximity to (10920). We developed a new method called STM-ORBIT to interpret the multi-band measurements without knowing the object’s true location. The power of the new technique is that it determines the most-likely heliocentric and observer-centric distance and phase angle ranges, allowing us to make a radiometric size estimate. The application to the MIRI fluxes of (10920) was used to validate the method. It leads to a confirmation of the known radiometric size-albedo solution, and puts constraints on the asteroid’s location and orbit in agreement with its true orbit. To back up the validation of the method, we obtained additional ground-based light curve observations of (10920), combined with Gaia data, which indicate a very elongated object (a/b ≥ 1.5), with a spin-pole at (λ, β) ecl = (178°, +81°), with an estimated error of about 20°, and a rotation period of 4.861191 ± 0.000015 h. A thermophysical study of all available JWST-MIRI and WISE measurements leads to a size of 14.5–16.5 km (diameter of an equal-volume sphere), a geometric albedo pV between 0.05 and 0.10, and a thermal inertia in the range 9–35 (best value 15) J m -2 s -0.5 K -1 . For the newly discovered MIRI object, the STM-ORBIT method revealed a size of 100–230 m. The new asteroid must be on a low-inclination orbit (0.7° < i < 2.0°) and it was located in the inner main-belt region during JWST observations. A beaming parameter η larger than 1.0 would push the size even below 100 meters, a main-belt regime that has escaped IR detections so far. This kind of MIRI observations can therefore contribute to formation and evolution studies via classical size-frequency studies, which are currently limited to objects larger than about one kilometer in size. We estimate that MIRI frames with pointings close to the ecliptic and short integration times of only a few seconds will always include a few asteroids; most of them will be unknown objects.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

The Strong Impact of Weak Horizontal Convergence on Continental Shallow Convection

Idealized large-eddy simulation (LES) is a basic tool for studying three-dimensional turbulence in the planetary boundary layer. LES is capable of providing benchmark solutions for parameterization development efforts. However, real small-scale atmospheric flows develop in heterogeneous and transient environments with locally varying vertical motions inherent to open multiscale interactive dynamical systems. These variations are often too subtle to detect them by state-of-the-art remote and in situ measurements, and are typically excluded from idealized simulations. The present study addresses the impact of weak [i.e., O(10-6) s-1] short-lived low-level large-scale convergence/divergence perturbations on continental shallow convection. The results show a strong response of shallow nonprecipitating convection to the applied weak large-scale dynamical forcing. Evolutions of CAPE, mean liquid water path, and cloud-top heights are significantly affected by the imposed convergence/divergence. In contrast, evolving cloud-base properties, such as the area coverage and mass flux, are only weakly affected. To contrast those impacts with microphysical sensitivity, the baseline simulations are perturbed assuming different observationally based cloud droplet number concentrations and thus different rainfall. For the tested range of microphysical perturbations, the imposed convergence/divergence provides significantly larger impact than changes in the cloud microphysics. Simulation results presented here provide a stringent test for convection parameterizations, especially important for large-scale models progressing toward resolving some nonhydrostatic effects.

Meteorology & Atmospheric Sciences↗
DOE OSTI2016

SPIN–ORBIT MISALIGNMENT AS A DRIVER OF THE KEPLER DICHOTOMY

During its five-year mission, the Kepler spacecraft has uncovered a diverse population of planetary systems with orbital configurations ranging from single-transiting planets to systems of multiple planets co-transiting the parent star. By comparing the relative occurrences of multiple to single-transiting systems, recent analyses have revealed a significant over-abundance of singles. Dubbed the “ Kepler Dichotomy,” this feature has been interpreted as evidence for two separate populations of planetary systems: one where all orbits are confined to a single plane, and a second where the constituent planetary orbits possess significant mutual inclinations, allowing only a single member to be observed in transit at a given epoch. In this work, we demonstrate that stellar obliquity, excited within the disk-hosting stage, can explain this dichotomy. Young stars rotate rapidly, generating a significant quadrupole moment, which torques the planetary orbits, with inner planets influenced more strongly. Given nominal parameters, this torque is sufficiently strong to excite significant mutual inclinations between planets, enhancing the number of single-transiting planets, sometimes through a dynamical instability. Furthermore, as hot stars appear to possess systematically higher obliquities, we predict that single-transiting systems should be relatively more prevalent around more massive stars. We analyze the Kepler data and confirm this signal to be present.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Direct Simulations of H–He Mixtures at Planetary Interior Conditions: Demixing, Insulator–Metal Transition and Miscibility Boundaries

Accurate knowledge of the electrical and thermal conductivities and structural properties of hydrogen–helium mixtures under thermodynamic conditions within and beyond the immiscibility range is very important to predict the thermal evolution and internal structure of gas giant planets like Jupiter and Saturn. Here, we propose a novel method to determine the immiscibility boundary accurately without the need for free energy calculations, while providing consistent insights into structural and transport properties of mixtures. We show with direct large-scale ab initio simulations that the insulator–metal transition (IMT) of the hydrogen subsystem is strongly affected by an admixture with a small fraction of helium and occurs at temperatures significantly higher than those of pure hydrogen. At pressures below 150 GPa, the IMT boundary is not related anymore to the H 2 subsystem dissociation, the system remains insulating even after the full dissociation of H 2 molecules and its transition to an H–He mixture. The offset of the IMT in the H–He mixture relative to the dissociation region in the hydrogen subsystem and the significant reduction of static electrical and thermal conductivity by a factor between two and a few thousand relative to pure hydrogen found in mixtures have consequences for Jupiter and Saturn’s thermal evolution, internal structure, and dynamo action, affecting a large fraction of the interior of both planets.

Helium↗
DOE OSTI2017

THE EFFECT OF ROTATION ON OSCILLATORY DOUBLE-DIFFUSIVE CONVECTION (SEMICONVECTION)

Oscillatory double-diffusive convection (ODDC, more traditionally called semiconvection) is a form of linear double-diffusive instability that occurs in fluids that are unstably stratified in temperature (Schwarzschild unstable), but stably stratified in chemical composition (Ledoux stable). This scenario is thought to be quite common in the interiors of stars and giant planets, and understanding the transport of heat and chemical species by ODDC is of great importance to stellar and planetary evolution models. Fluids unstable to ODDC have a tendency to form convective thermocompositional layers that significantly enhance the fluxes of temperature and chemical composition compared with microscopic diffusion. Although a number of recent studies have focused on studying properties of both layered and nonlayered ODDC, few have addressed how additional physical processes such as global rotation affect its dynamics. In this work, we study first how rotation affects the linear stability properties of rotating ODDC. Using direct numerical simulations, we then analyze the effect of rotation on properties of layered and nonlayered ODDC, and we study how the angle of the rotation axis with respect to the direction of gravity affects layering. We find that rotating systems can be broadly grouped into two categories based on the strength of rotation. The qualitative behavior in the more weakly rotating group is similar to nonrotating ODDC, but strongly rotating systems become dominated by vortices that are invariant in the direction of the rotation vector and strongly influence transport. We find that whenever layers form, rotation always acts to reduce thermal and compositional transport.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

The Scattering Outcomes of Kepler Circumbinary Planets: Planet Mass Ratio

Recent studies reveal that the free eccentricities of Kepler-34b and Kepler-413b are much larger than their forced eccentricities, implying that scattering events may take place in their formation. The observed orbital configuration of Kepler-34b cannot be well reproduced in disk-driven migration models, whereas a two-planet scattering scenario can play a significant role of shaping the planetary configuration. These studies indicate that circumbinary planets discovered by Kepler may have experienced scattering process. In this work, we extensively investigate the scattering outcomes of circumbinary planets focusing on the effects of planet mass ratio . We find that the planetary mass ratio and the the initial relative locations of planets act as two important parameters that affect the eccentricity distribution of the surviving planets. As an application of our model, we discuss the observed orbital configurations of Kepler-34b and Kepler-413b. We first adopt the results from the disk-driven models as the initial conditions, then simulate the scattering process that occurs in the late evolution stage of circumbinary planets. We show that the present orbital configurations of Kepler-34b and Kepler-413b can be well reproduced when considering a two unequal-mass planet ejection model. Our work further suggests that some of the currently discovered circumbinary single-planet systems may be survivors of original multiple-planet systems. The disk-driven migration and scattering events occurring in the late stage both play an irreplaceable role in sculpting the final systems.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

The possibility of hydrogen-water demixing in Uranus, Neptune, K2-18 b and TOI-270 d

Context. The internal structures of Uranus and Neptune remain unknown. In addition, sub-Neptunes are now thought to be the most common type of exoplanets. Improving our understanding of the physical processes that govern the interiors of such planets is therefore essential. Phase separation between planetary constituents may occur, in particular, hydrogen-water immiscibility in cold, water-rich intermediate-mass planets. Aims. We assess whether hydrogen-water demixing could occur in Uranus, Neptune, K2-18 b and TOI-270 d, and investigate its effect on the planetary evolution and inferred internal structure. Methods. We couple planetary evolution models with recent ab initio calculations of the hydrogen-water phase diagram, allowing for temperature shifts to account for uncertainties in miscibility gaps. Results. We find that demixing may occur and could lead to a complete depletion of water in the outermost regions of Uranus and Neptune. Temperature offsets of up to 1100 K lead to a depleted region comprising as much as 16% of the planet’s mass, and an increase in planetary radius by nearly 20%. For K2-18 b, our models suggest that hydrogen-water demixing is ongoing and may explain the absence of water features in its JWST spectrum. A temperature offset of 500 K is required to get a complete depletion of water in the atmosphere of K2-18b. TOI-270d may also have experienced hydrogen-water demixing. When applying a similar temperature offset on the phase diagram as for K2-18 b, we find a partial depletion of water in the atmosphere of TOI-270 d, consistent with JWST’s detection of water. Conclusions. Hydrogen-water immiscibility may play a key role in shaping the structure and evolution of both Solar System giant planets like Uranus and Neptune, and cold/temperate exoplanets such as K2-18 b and TOI-270 d. Accounting for such internal processes is crucial to accurately interpret atmospheric observations from current (e.g., JWST) and upcoming (e.g., ARIEL) missions.

planets and satellites: composition↗