Search NASA⌕ Search

SEARCH · Search NASA

Results for “planetary dynamics”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Rare Event Sampling Improves Mercury Instability Statistics

Due to the chaotic nature of planetary dynamics, there is a non-zero probability that Mercury's orbit will become unstable in the future. Previous efforts have estimated the probability of this happening between 3 and 5 billion years in the future using a large number of direct numerical simulations with an N-body code, but were not able to obtain accurate estimates before 3 billion years in the future because Mercury instability events are too rare. In this paper we use a new rare-event sampling technique, Quantile Diffusion Monte Carlo (QDMC), to estimate that the probability of a Mercury instability event in the next 2 billion years is approximately 10 –4 in the REBOUND N-body code. We show that QDMC provides unbiased probability estimates at a computational cost of up to 100 times less than direct numerical simulation. QDMC is easy to implement and could be applied to many problems in planetary dynamics in which it is necessary to estimate the probability of a rare event.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Eject, crash, or survive: Using machine learning to predict orbital instability of exoplanetary systems

Astronomers throughout history, including titans like Kepler and Newton, have tackled planetary dynamics and orbital instability. Despite strides taken in research, understanding the evolution of planetary orbits remains an intricate, computationally expensive, and analytically unsolved problem. I apply machine learning classification methods to numerical simulations of planetary systems in order to predict the long-term fate of the planet - whether the planet remains in a stable orbit or not. My method uses the first 41.1 years (≤ 500 orbits) of data from a planet’s simulation to calculate 17 dynamically-motivated metrics; I trained my classifier on these features to predict a planet’s stability after 107 years. At 84.33%, my classifier was comparable in accuracy to pre-existing literature, despite using significantly less computational power than most other methods. In my research, I found that the standard deviation of eccentricity, mass ratios for neighboring planets, and semi-major axis ratio with the outer planet neighbor to be the most predictive features of instability. I propose reasons for the importance of these features, their role in planetary dynamics, as well as possible explanations for why some planets were misclassified. By understanding the important metrics of instability and reasons for misclassification, we can begin to understand more about system architectures, orbital motion and dynamics, and the formation and evolution of the exoplanetary systems. This is applicable in our own Solar System, and with exoplanet discovery missions such as TESS, this research becomes especially relevant in understanding the new exoplanetary systems we discover.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Eject, crash, or survive: Using machine learning to predict orbital instability of exoplanetary systems

Astronomers throughout history, including titans like Kepler and Newton, have tackled planetary dynamics and orbital instability. Despite strides taken in research, understanding the evolution of planetary orbits remains an intricate, computationally expensive, and analytically unsolved problem. I apply machine learning classification methods to numerical simulations of planetary systems in order to predict the long-term fate of the planet - whether the planet remains in a stable orbit or not. My method uses the first 41.1 years (≤ 500 orbits) of data from a planet’s simulation to calculate 17 dynamically-motivated metrics; I trained my classifier on these features to predict a planet’s stability after 107 years. At 84.33%, my classifier was comparable in accuracy to pre-existing literature, despite using significantly less computational power than most other methods. In my research, I found that the standard deviation of eccentricity, mass ratios for neighboring planets, and semi-major axis ratio with the outer planet neighbor to be the most predictive features of instability. I propose reasons for the importance of these features, their role in planetary dynamics, as well as possible explanations for why some planets were misclassified. By understanding the important metrics of instability and reasons for misclassification, we can begin to understand more about system architectures, orbital motion and dynamics, and the formation and evolution of the exoplanetary systems. This is applicable in our own Solar System, and with exoplanet discovery missions such as TESS, this research becomes especially relevant in understanding the new exoplanetary systems we discover.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Long-period Jovian Tilts the Orbits of Two sub-Neptunes Relative to Stellar Spin Axis in Kepler-129

We present the discovery of Kepler-129 d (P{sub d}=7.2{sub -0.3}{sup +0.4} yr, msini{sub d}=8.3{sub -0.7}{sup +1.1} M{sub Jup}, e{sub d}=0.15{sub -0.05}{sup +0.07}) based on six years of radial-velocity observations from Keck/HIRES. Kepler-129 also hosts two transiting sub-Neptunes: Kepler-129 b (P {sub b} = 15.79 days, r {sub b} = 2.40 ± 0.04 R {sub ⊕}) and Kepler-129 c (P {sub c} = 82.20 days, r {sub c} = 2.52 ± 0.07 R {sub ⊕}) for which we measure masses of m {sub b} < 20 M {sub ⊕} and m{sub c}=43{sub -12}{sup +13} M{sub ⊕}. Kepler-129 is a hierarchical system consisting of two tightly packed inner planets and a massive external companion. In such a system, two inner planets precess around the orbital normal of the outer companion, causing their inclinations to oscillate with time. Based on an asteroseismic analysis of Kepler data, we find tentative evidence that Kepler-129 b and c are misaligned with stellar spin axis by ≳38°, which could be torqued by Kepler-129 d if it is inclined by ≳19° relative to inner planets. Using N-body simulations, we provide additional constraints on the mutual inclination between Kepler-129 d and inner planets by estimating the fraction of time during which two inner planets both transit. The probability that two planets both transit decreases as their misalignment with Kepler-129 d increases. We also find a more massive Kepler-129 c enables the two inner planets to become strongly coupled and more resistant to perturbations from Kepler-129 d. The unusually high mass of Kepler-129 c provides a valuable benchmark for both planetary dynamics and interior structure, since the best-fit mass is consistent with this 2.5 R {sub ⊕} planet having a rocky surface.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Strength, plasticity, and spin transition of Fe-N compounds in planetary cores

Elastic and plastic properties of Fe-light element alloys and compounds are needed to determine the compositions and dynamics of planetary cores. Elastic strength and plastic deformation mechanisms and their relationship to electronic properties of ε-Fe 7 N 3 and γ’-Fe 4 N mixture were investigated by x-ray diffraction and x-ray emission spectroscopy in the diamond anvil cell from 1 bar up to 60 GPa. X-ray diffraction shows that ε-Fe 7 N 3 reaches a pressure of 15–20 GPa before undergoing bulk plasticity at a differential stress of 4.4–10.4 GPa. ε-Fe 7 N 3 is stronger than γ’-Fe 4 N and hcp-Fe which achieve a flow stress of 1.5–3.6 GPa at 10–15 GPa and 2–3 GPa at ~20 GPa, respectively. X-ray emission spectroscopy shows that a decrease in electronic spin moment begins before and completes after plastic flow onset for each nitride, suggesting that pressure-driven changes in electronic arrangement do not trigger a plastic response although they may modify the strength and plastic behavior of Fe-N compounds. Plastic deformation in ε-Fe 7 N 3 and hcp-Fe results in a preferred orientation of (0001) normal to maximum compression, while γ’-Fe 4 N develops a maximum in the (110). Furthermore, these observations may be combined with measurements of elasticity to model seismic properties of cores of small planetary bodies such as Mars, Mercury, and the Moon.

58 GEOSCIENCES↗
DOE OSTI2025

Experimental Pathways for Detecting Double Superionicity in Planetary Ices

The ice giant planets Uranus and Neptune are assumed to contain large amounts of planetary ices such as water, methane, and ammonia. The properties of mixtures of such ices at the extreme pressures and temperatures of planetary interiors are not yet well understood. Ab initio computer simulations have predicted that a number of ices exhibit a hydrogen superionic state and a doubly superionic state. Since the latter state has not yet been generated with experiments, we outline here two possible pathways for reaching and detecting such a state with dynamic compression experiments. Here, we suggest X-ray diffraction as the principal tool for detecting when the material becomes doubly superionic and the sublattice of one of the heavy nuclei melts. That would require a temperature of ~3500 K and pressures greater than ~200 GPa for H 3 NO 4 , which we use as an example material here. Such conditions can be reached with experiments that employ an initial shock that is followed by a ramp compression wave. Alternatively, one may use triple-shock compression because a single shock does not yield sufficiently high densities.

Neptune↗
DOE OSTI2024

Phase Transitions and Thermal Equation of State of Fe‐9wt.%Si Applied to the Moon and Mercury

Abstract Accurate knowledge of the phase transitions and thermoelastic properties of candidate iron alloys, such as Fe‐Si alloys, is essential for understanding the nature and dynamics of planetary cores. The phase diagrams of some Fe‐Si alloys between 1 atm and 16 GPa have been back‐extrapolated from higher pressures, but the resulting phase diagram of Fe 83.6 Si 16.4 (9 wt.% Si) is inconsistent with temperature‐induced changes in its electrical resistivity between 6 and 8 GPa. This study reports in situ synchrotron X‐ray diffraction (XRD) measurements on pre‐melted and powder Fe 83.6 Si 16.4 samples from ambient conditions to 60 GPa and 900 K using an externally heated diamond‐anvil cell. Upon compression at 300 K, the bcc phase persisted up to ∼38 GPa. The hcp phase appeared near 8 GPa in the pre‐melted sample, and near 17 GPa in the powder sample. The appearance of the hcp phase in the pre‐melted sample reconciles the reported changes in electrical resistivity of a similar sample, thus resolving the low‐pressure region of the phase diagram. The resulting high‐temperature Birch‐Murnaghan equation of state (EoS) and thermal EoS based on the Mie‐Gruneisen‐Debye model of the bcc and hcp structures are consistent with, and complement the literature data at higher pressures. The calculated densities based on the thermal EoS of Fe‐9wt.%Si indicate that both bcc and hcp phases agree with the reported core density estimates for the Moon and Mercury.

Berrada, Meryem↗
DOE OSTI2022

Preface for frontiers of magnetic reconnection research in heliophysical, astrophysical, and laboratory plasmas

Magnetic reconnection—the topological rearrangement of magnetic field—underlies many explosive phenomena across a wide range of natural and laboratory plasmas.3 It plays a pivotal role in electron and ion heating, particle acceleration to high energies, energy transport, and self-organization. Reconnection can have a complex relationship with turbulence at both large and small scales, leading to various effects that are only beginning to be understood. In heliophysics, magnetic reconnection plays a key role in solar flares, coronal mass ejections, coronal heating, solar wind dissipation, the interaction of interplanetary plasma with magnetospheres, dynamics of planetary magnetospheres such as magnetic substorms, and the heliospheric boundary with the interstellar medium. Here, the magnetic reconnection is integral to the solar and planetary dynamo processes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Hydrogen solubility in FeSi alloy phases at high pressures and temperatures

Abstract Light elements alloying with metallic Fe can change the properties and therefore play a key role in the structure and dynamics of planetary cores. Hydrogen and silicon are possible light elements in the rocky planets’ cores. However, hydrogen storage in Fe-Si alloy systems remains unclear at high pressures and high temperatures because of experimental difficulties. Taking advantage of pulsed laser heating combined with high-energy synchrotron X-ray diffraction, we studied reactions between FeSi and H in laser-heated diamond-anvil cells (LHDACs) up to 61.9 GPa and 3500 K. We found that under H-saturated conditions the amount of H alloying with FeSi (0.3 and <0.1 wt% for the B20 and B2 structures, respectively) is much smaller than that in pure Fe metal (>1.8 wt%). Our experiments also suggest that H remains in the crystal structure of FeSi alloy when recovered to 1 bar. Further density functional theory (DFT) calculations indicate that the low-H solubility likely results from the highly distorted interstitial sites in the B20 and B2 structures, which are not favorable for H incorporation. The recovery of H in the B20 FeSi crystal structure at ambient conditions could open up possibilities to understand geochemical behaviors of H during core formation in future experiments. The low-H content in FeSi alloys suggests that if a planetary core is Si-rich, Si can limit the ingassing of H into the Fe-rich core.

Geochemistry & Geophysics↗
DOE OSTI2024

BD+44°493: Chemo-dynamical Analysis and Constraints on Companion Planetary Masses from WIYN/NEID Spectroscopy

In this work, we present high-resolution (R ~ 100,000), high signal-to-noise ratio (S/N ~ 800) spectroscopic observations for the well-known, bright, extremely metal-poor, carbon-enhanced star BD+44°493. We determined chemical abundances and upper limits for 17 elements from WIYN/NEID data, complemented with 11 abundances redetermined from Subaru and Hubble data, using the new, more accurate, stellar atmospheric parameters calculated in this work. Our analysis suggests that BD+44°493 is a low-mass (0.83 M ⊙ ), old (12.1–13.2 Gyr) second-generation star likely formed from a gas cloud enriched by a single metal-free 20.5 M ⊙ Population III star in the early Universe. With a disk-like orbit, BD+44°493 does not appear to be associated with any major merger event in the early history of the Milky Way. From the precision radial-velocity NEID measurements (median absolute deviation = 16 m s –1 ), we were able to constrain companion planetary masses around BD+44°493 and rule out the presence of planets as small as $m\sin i=2$ M J out to periods of 100 days. This study opens a new avenue of exploration for the intersection between stellar archaeology and exoplanet science using NEID.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

A possible misaligned orbit for the young planet AU Mic c

ABSTRACT The AU Microscopii planetary system is only 24 Myr old, and its geometry may provide clues about the early dynamical history of planetary systems. Here, we present the first measurement of the Rossiter–McLaughlin effect for the warm sub-Neptune AU Mic c, using two transits observed simultaneously with the European Southern Observatory's (ESO's) Very Large Telescope (VLT)/Echelle SPectrograph for Rocky Exoplanets and Stable Spectroscopic Observations (ESPRESSO), CHaracterising ExOPlanet Satellite (CHEOPS), and Next-Generation Transit Survey (NGTS). After correcting for flares and for the magnetic activity of the host star, and accounting for transit-timing variations, we find the sky-projected spin–orbit angle of planet c to be in the range $\lambda _{\mathrm{c}}=67.8_{-49.0}^{+31.7}$ degrees (1$\sigma$). We examine the possibility that planet c is misaligned with respect to the orbit of the inner planet b ($\lambda _{\mathrm{b}}=-2.96_{-10.30}^{+10.44}$), and the equatorial plane of the host star, and discuss scenarios that could explain both this and the planet’s high density, including secular interactions with other bodies in the system or a giant impact. We note that a significantly misaligned orbit for planet c is in some degree of tension with the dynamical stability of the system, and with the fact that we see both planets in transit, though these arguments alone do not preclude such an orbit. Further observations would be highly desirable to constrain the spin–orbit angle of planet c more precisely.

Yu, H. (ORCID:0000000209716078)↗
DOE OSTI2026

Geostrophic to magnetically-damped turbulence in liquid metal rotating magnetoconvection

Understanding planetary core convection dynamics requires the study of convective flows in which the Coriolis and Lorentz forces attain a leading-order, so-called magnetostrophic balance. Experimental investigations of rotating magnetoconvection (RMC) in the magnetostrophic regime are therefore essential to broadly characterize the properties of local-scale planetary core flow. Toward this end, we present here the first thermovelocimetric measurements of magnetostrophic, liquid metal convection, which are made using liquid gallium as the working fluid, at moderate rotation rates (Ekman numbers 10 -4 ≤ Ek ≤ 10 -5 ) and in the presence of dynamically strong magnetic fields (Elsasser number Λ = 1). Complementary rotating convection (RC) experiments are performed at the same rotation rates to serve as reference cases. Our RMC velocity measurements adequately follow a geostrophic turbulent scaling for cases in which local-scale convective inertial forces exceed the Lorentz forces in the fluid bulk. In cases where Lorentz forces exceed local-scale inertia (local or convective interaction parameter N ℓ or N C ≳ 3), the root-mean-square RMC velocities are magnetically damped, yielding values below the geostrophic turbulent RC scaling prediction. An enhancement in heat transfer is observed, which we attribute to the increased coherence of vertically aligned magnetostrophic convective flow. Finally, our results and comparisons to dynamo simulations suggest that the system is in the strong-field regime, where convection-scale flows are magnetically damped with N C ≳ 3).

Convective turbulence↗
DOE OSTI2020

Dynamical Mass Estimates of the β Pictoris Planetary System through Gaussian Process Stellar Activity Modeling

Nearly 15 yr of radial velocity (RV) monitoring and direct imaging enabled the detection of two giant planets orbiting the young, nearby star β Pictoris. The δ Scuti pulsations of the star, which overwhelm planetary signals, need to be carefully suppressed. In this work, we independently revisit the analysis of the RV data following a different approach than available in the literature to model the activity of the star. We show that a Gaussian process (GP) with a stochastically driven damped harmonic oscillator kernel can model the δ Scuti pulsations. It provides similar results to parametric models but with a simpler framework, using only three hyperparameters. It also enables us to model poorly sampled RV data that were excluded from previous analyses, hence extending the RV baseline by nearly five years. Altogether, the orbit and mass of both planets can be constrained from RV only, which was not possible with the parametric modeling. To characterize the system more accurately, we also perform a joint fit of all available relative astrometry and RV data. Our orbital solutions for β Pic b favor a low eccentricity of 0.029$_{−0.024}^{+0.061}$ and a relatively short period of 21.1$_{−0.8}^{+2.0}$ yr. The orbit of β Pic c is eccentric with 0.206$_{−0.063}^{+0.074}$ with a period of 3.36 ± 0.03 yr. We find model-independent masses of 11.7 ± 1.4 and 8.5 ± 0.5 M Jup for β Pic b and c, respectively, assuming coplanarity. The mass of β Pic b is consistent with the hottest start evolutionary models, at an age of 25 ± 3 Myr. A direct detection of β Pic c would provide a second calibration measurement in a coeval system.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Global energy spectrum of the general oceanic circulation

Abstract Advent of satellite altimetry brought into focus the pervasiveness of mesoscale eddies $${{{{{{{\bf{{{{{{{{\mathcal{O}}}}}}}}}}}}}}}}({100})$$ O ( 100 ) km in size, which are the ocean’s analogue of weather systems and are often regarded as the spectral peak of kinetic energy (KE). Yet, understanding of the ocean’s spatial scales has been derived mostly from Fourier analysis in small "representative” regions that cannot capture the vast dynamic range at planetary scales. Here, we use a coarse-graining method to analyze scales much larger than what had been possible before. Spectra spanning over three decades of length-scales reveal the Antarctic Circumpolar Current as the spectral peak of the global extra-tropical circulation, at ≈ 10 4 km, and a previously unobserved power-law scaling over scales larger than 10 3 km. A smaller spectral peak exists at ≈ 300 km associated with mesoscales, which, due to their wider spread in wavenumber space, account for more than 50% of resolved surface KE globally. Seasonal cycles of length-scales exhibit a characteristic lag-time of ≈ 40 days per octave of length-scales such that in both hemispheres, KE at 10 2 km peaks in spring while KE at 10 3 km peaks in late summer. These results provide a new window for understanding the multiscale oceanic circulation within Earth’s climate system, including the largest planetary scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Thermodynamic properties of ε -Fe with thermal electronic excitation effects on vibrational spectra

The thermodynamic properties of hexagonal-close-packed iron (ε–Fe) are essential for investigating the internal structure and dynamic properties of planetary cores. Despite their importance to planetary sciences, experimental investigations of ε–Fe at relevant conditions are still challenging. Therefore, ab initio calculations are crucial to elucidating the thermodynamic properties of this system. Here, we use a free energy calculation scheme based on the phonon gas model compatible with temperature-dependent phonon frequencies. We investigate the effects of electronic thermal excitations, which introduces a temperature dependence on phonon frequencies, and the implication for the thermodynamic properties of ε–Fe at extreme pressure (P) and temperature (T) conditions. We disregard phonon-phonon interactions, i.e., anharmonicity and their effect on phonon frequencies. Nevertheless, the current scheme is also applicable to T -dependent anharmonic frequencies. We conclude that the impact of thermal electronic excitations on vibrational properties is not significant up to ~4000 K at 200 GPa but should not be ignored at higher temperatures or pressures. However, the static free energy F st must always include the effect of thermal excitation fully in a continuum of T. Furthermore, our results for isentropic equations of state show good agreement with data from recent ramp compression experiments up to 1400 GPa conducted at the National Ignition Facility.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

A Shallow-Water Model for Convective Self-Aggregation

Randomly distributed convective storms can self-aggregate in the absence of large-scale forcings. Here we present a 1D shallow-water model to study the convective self-aggregation. This model simulates the dynamics of the planetary boundary layer and represents convection as a triggered process. Once triggered, convection lasts for finite time and occupies finite length. Here we show that the model can successfully simulate self-aggregation, and that the results are robust to a wide range of parameter values. In the simulations, convection excites gravity waves. The gravity waves then form a standing wave pattern, separating the domain into convectively active and inactive regions. We analyze the available potential energy (APE) budget and show that convection generates APE, providing energy for self-aggregation. By performing dimensional analysis, we develop a scaling theory for the size of convective aggregation, which is set by the gravity wave speed, damping time scale, and number density of convective storms. This paper provides a simple modeling framework to further study convective self-aggregation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Planets Across Space and Time (PAST). II. Catalog and Analyses of the LAMOST–Gaia–Kepler Stellar Kinematic Properties

The Kepler telescope has discovered over 4000 planets (candidates) by searching ~200,000 stars over a wide range of distance (order of kpc) in our Galaxy. Characterizing the kinematic properties (e.g., Galactic component membership and kinematic age) of these Kepler targets (including the planet candidate hosts) is the first step toward studying Kepler planets in the Galactic context, which will reveal fresh insights into planet formation and evolution. In this paper, the second part of the Planets Across the Space and Time (PAST) series, by combining the data from the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) and Gaia and then applying the revised kinematic methods from PAST I, we present a catalog of kinematic properties (i.e., Galactic positions, velocities, and the relative membership probabilities among the thin disk, thick disk, Hercules stream, and the halo) as well as other basic stellar parameters for 35,835 Kepler stars. Further analyses of the LAMOST–Gaia–Kepler catalog demonstrate that our derived kinematic age reveals the expected stellar activity-age trend. Furthermore, we find that the fraction of thin (thick) disk stars increases (decreases) with the transiting planet multiplicity (N {sub p} = 0, 1, 2 and 3+) and the kinematic age decreases with N {sub p}, which could be a consequence of the dynamical evolution of planetary architecture with age. The LAMOST–Gaia–Kepler catalog will be useful for future studies on the correlations between the exoplanet distributions and the stellar Galactic environments as well as ages.

47 OTHER INSTRUMENTATION↗