Search NASASearch

SEARCH · Search NASA

Results for “planetary core”

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

Thermoelastic Anomaly of Iron Carbonitride Across the Spin Transition and Implications for Planetary Cores

Carbon and nitrogen are considered as candidate light elements present in planetary cores. However, there is limited understanding regarding the structure and physical properties of Fe-C-N alloys under extreme conditions. Here diamond anvil cell experiments were conducted, revealing the stability of hexagonal-structured Fe 7 (N 0.75 C 0.25 ) 3 up to 120 GPa and 2100 K, without undergoing any structural transformation or dissociation. Notably, the thermal expansion coefficient and Grüneisen parameter of the alloy exhibit a collapse at 55–70 GPa. First-principles calculations suggest that such anomaly is associated with the spin transition of iron within Fe 7 (N 0.75 C 0.25 ) 3 . Our modeling indicates that the presence of ∼1.0 wt% carbon and nitrogen in liquid iron contributes to 9–12% of the density deficit of the Earth's outer core. The thermoelastic anomaly of the Fe-C-N alloy across the spin transition is likely to affect the density and seismic velocity profiles of (C,N)-rich planetary cores, thereby influencing the dynamics of such cores.

Geology

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

Iron Bonding with Light Elements: Implications for Planetary Cores Beyond the Binary System

Light element alloying in iron is required to explain density deficit and seismic wave velocities in Earth’s core. However, the light element composition of the Earth’s core seems hard to constrain as nearly all light element alloying would reduce the density and sound velocity (elastic moduli). The alloying light elements include oxidizing elements like oxygen and sulfur and reducing elements like hydrogen and carbon, yet their chemical effects in the alloy system are less discussed. Moreover, Fe-X-ray Absorption Near Edge Structure (Fe-XANES) fingerprints have been studied for silicate materials with ferrous and ferric ions, while not many X-ray absorption spectroscopy (XAS) studies have focused on iron alloys, especially at high pressures. To investigate the bonding nature of iron alloys in planetary interiors, we presented X-ray absorption spectroscopy of iron–nitrogen and iron–carbon alloys at high pressures up to 50 GPa. Together with existing literature on iron–carbon, –hydrogen alloys, we analyzed their edge positions and found no significant difference in the degree of oxidation among these alloys. Pressure effects on edge positions were also found negligible. Our theoretical simulation of the valence state of iron, alloyed with S, C, O, N, and P also showed nearly unchanged behavior under pressures up to 300 GPa. This finding indicates that the high pressure bonding of iron alloyed with light elements closely resembles bonding at the ambient conditions. We suggest that the chemical properties of light elements constrain which ones can coexist within iron alloys.

Advanced Photon Source

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

Dynamic strength of iron under pressure-temperature conditions of Earth’s inner core

Iron (Fe) is a primary constituent of terrestrial planetary cores, yet its rheological properties under extreme conditions remain uncertain. Here we present direct measurements of Fe strength at 310-430 GPa pressures and 3700-5800 K temperatures, obtained using Rayleigh-Taylor (RT) instability experiments at the National Ignition Facility. Single-crystal α-Fe samples with [001] and [111] orientations are shock-ramp compressed past the α-ε transition along paths approaching Earth’s inner core conditions. We find that ε-Fe derived from [001] α-Fe is consistently stronger (11-20 GPa) than that from [111] α-Fe (8-18 GPa), contrary to the trend at ambient conditions. Large-scale molecular dynamics simulations reproduce this atypical strength anisotropy and attribute it to microstructural evolution during the phase transition and subsequent ε-phase plasticity. Ripple growth analysis further constrains viscosities of 100-170 Pa·s under the driven conditions. Our results provide experimental benchmarks for Fe rheology at inner-core conditions, with implications for seismic anisotropy and the geodynamo.

Condensed-matter physics

Liquid Structure of Iron and Iron–Nitrogen–Carbon Alloys Within the Cores of Small Terrestrial Bodies

Abstract Nitrogen has been proposed to be stored within planetary cores, but its effects on the structure and density of molten Fe–alloys have not been explored experimentally. Using energy‐dispersive X‐ray diffraction, we determined the structure of Fe–N(–C) liquids at core conditions (1–7 GPa and 1700–1900°C) within a Paris‐Edinburgh press. Variation of N up to 7 wt.% and C up to 1.5 wt.% results in near‐linear changes in Fe–Fe atom distances and structure factor with increasing light element content. We did not observe a significant pressure‐driven structural transition in Fe–N(–C) liquids. We model the expansion of the Fe–Fe bonds using a modified Birch‐Murnaghan equation of state. With this model, we demonstrate that N or C contamination could lead to an overestimation of the Fe–Fe distances of pure Fe. We observe that the incorporation of 1 wt.% N or C into Fe results in a change in Fe–Fe distances that is twice as significant as the effect of 1 GPa. By approximating the change in volume, we infer that N and C incorporated in liquid iron could contribute to the density deficit observed in the cores of terrestrial bodies.

Pease, Allison [Department of Earth and Environmen

Sound velocities and thermal equation of state of fcc -iron-nickel alloys at high pressure and high temperature: Implications for the cores of Moon and several planets

Fcc-Fe-Ni alloy is believed to be the most dominant solid constitute of moderate-sized terrestrial planetary cores. Investigating the physical properties, especially the density and sound velocity of Fe-Ni alloys and comparing them with seismic observations is an indispensable approach to constructing compositional models for planetary interiors. In this study, we conducted sound velocity measurements on Fe-Ni alloys with 10 wt.% and 20 wt.% Ni up to ∼13.5 GPa and 1073 K, using the ultrasonic interferometry technique in a multi-anvil apparatus in conjunction with synchrotron radiation. By fitting the experimental data to finite strain equations, the bulk and shear moduli and their pressure and temperature derivatives are derived, yielding K S0 =145.8(14) GPa, G 0 = 73.2(7) GPa, K S0 ’ = 5.89(24), G 0 ’ = 2.89(8), (∂K S /∂T) P = -0.0181(12) GPa/K and (∂G/∂T) P = -0.0393(10) GPa/K for fcc-Fe 80 Ni 20 . An examination of the density-velocity relationship shows that compressional wave velocity is insensitive to temperature within the current pressure and temperature range, while shear wave velocity exhibits a large reduction with increasing temperature. Here, extrapolation of the sound velocities following the finite strain theories suggests that much slower Vs should be expected at pressure and temperature conditions corresponding to those of the lunar core. Possible core density and velocity profiles for other moderate planets and satellites, such as Mars, Mercury, and Ganymede are also calculated.

Equation of state

Thermovelocimetric characterization of liquid metal convection in a rotating slender cylinder

Rotating turbulent convection occurs ubiquitously in natural convective systems encompassing planetary cores, oceans and atmospheres, as well as in many industrial applications. While the global heat and mass transfer of water-like rotating Rayleigh–Bénard convection is well-documented, the dynamics in low-Prandtl-number liquid metals remain less understood. In this study, we experimentally investigate rotating Rayleigh-Bénard convection in liquid gallium (Prandtl number Pr ≈ 0.027) within a slender cylinder (diameter-to-height aspect ratio Γ = D/H =1/2) using novel thermovelocimetric diagnostic techniques that integrate simultaneous multi-point thermometry and ultrasonic Doppler velocity measurements. Our results reveal the formation of a stable, global-scale azimuthal wavenumber m = 2 quadrupolar vortex at low supercriticality. We propose that enhanced wall modes facilitated by the slender cylinder geometry interact with the bulk flow to create these large-scale axialized vortices. Furthermore, our findings imply a distinct scaling behavior for the wall-mode precession frequency in liquid metals, extending previous results obtained for moderate-Pr fluids. In conclusion, this provides new insights into wall-bulk coupling mechanisms of low-Pr rotating convective turbulence.

Magnetohydrodynamics

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

Effect of pressure cycling and compression rate on the bcc-hcp transition in an FeNi alloy

Here, we investigate the body-centered cubic (bcc) to hexagonal close-packed (hcp) phase transition in Fe-10wt. %Ni alloy, combining pressure cycling and fast compression with time-resolved synchrotron x-ray diffraction in a dynamic diamond anvil cell. Three pressure cycles were conducted with compression rates ranging from 0.1 to nearly 10 3 GPa/s. During the first cycle with the slowest compression, the observed orientations in the bcc and hcp phases are consistent with the Burgers mechanism, followed by c-axis rotation of the hcp phase consistent with {$10$ $\bar{1}2$} twinning. During the following cycles with fast compression at 10 2 –10 3 GPa/s, the hcp phase exhibits negligible c-axis rotation with a nearly constant c/a ratio of ∼1.61 up to ∼30 GPa, indicating suppression of plastic deformation (especially, twinning) due to sample confinement from the gasket. Notably, the onset pressure of the transition decreases with additional pressure cycling and faster compression, which normally leads to over-pressurization. This suggests that defects or shear induced from the pressure cycling reduces the transition kinetics even during fast compression. These insights into the deformation and transition behavior in an FeNi alloy under multiple dynamic loading cycles can offer guidance for future design of advanced structural alloys and improve our understanding of planetary core processes.

FeNi alloy

Dynamical structure factors of warm dense matter from time-dependent orbital-free and mixed-stochastic-deterministic density functional theory

Abstract We present the first calculations of the inelastic part of the dynamical structure factor (DSF) for warm dense matter (WDM) using time-dependent orbital-free density functional theory (TD-OF-DFT) and mixed-stochastic-deterministic (mixed) Kohn Sham TD-DFT (KS TD-DFT). WDM is an intermediate phase of matter found in planetary cores and laser-driven experiments, where the accurate calculation of the DSF is critical for interpreting x-ray Thomson scattering measurements. Traditional TD-DFT methods, while highly accurate, are computationally expensive, motivating the exploration of TD-OF-DFT and mixed TD-KS-DFT as more efficient alternatives. We applied these methods to experimentally measured WDM systems, including solid-density aluminum and beryllium, compressed beryllium, and carbon–hydrogen mixtures. Our results show that TD-OF-DFT requires a dynamical kinetic energy potential in order to qualitatively capture the plasmon response. Additionally, it struggles with capturing bound electron contributions. In contrast, mixed TD-KS-DFT offers greater accuracy in distinguishing bound and free electron effects, aligning well with experimental data, though at a higher computational cost. This study highlights the trade-offs between computational efficiency and accuracy, demonstrating that TD-OF-DFT remains a valuable tool for rapid scans of parameter space, while mixed TD-KS-DFT should be preferred for high-fidelity simulations. Our findings provide insight into the future development of DFT methods for WDM and suggest potential improvements for TD-OF-DFT.

36 MATERIALS SCIENCE

Electrical conductivities of (Mg,Fe)O at extreme pressures and implications for planetary magma oceans

During planet formation, planets undergo many impacts that can generate magma oceans. When these crystallize, part of the magma densifies via iron enrichment and migrates to the core–mantle boundary, forming an iron-rich basal magma ocean (BMO). The BMO could generate a dynamo in early Earth and super-Earths if the electrical conductivity of the BMO, which is thought to be sensitive to its Fe content, is sufficiently high. To test this hypothesis, here we conduct laser-driven shock experiments on ferropericlase (Mg x ,Fe 1−x )O (0.95 ≤ x ≤ 1) as an Fe-rich BMO analogue, perform density functional theory molecular dynamics simulations on MgO and calculate the long-term evolution of super-Earths. We find that the d.c. conductivities of MgO and (Mg,Fe)O are indistinguishable between 467 GPa and 1,400 GPa, despite previous predictions. Here, we predict that super-Earths larger than 3–6 Earth masses can produce BMO-driven dynamos that are almost one order of magnitude stronger than core-driven dynamos for several billion years.

Exoplanets

Revealing the Moon's Taurus‐Littrow Landslide via Integrated Analysis of Pristine Apollo 17 Soil Core 73001/2

The “light mantle” deposit at the base of South Massif in the Moon's Taurus‐Littrow Valley was a primary science target for the Apollo 17 exploration. The possibility that it was a landslide triggered by ejecta from Tycho Crater is critical for establishing the age of Tycho and constraining recent lunar impact chronology; however, the mechanism of emplacement of the deposit has recently been questioned. The newly opened 73001/73002 double drive tube from Station 3 sampled 70.6 cm deep into the regolith and represents the first stratigraphic section of an extraterrestrial landslide deposit returned to Earth. Here we provide an overview of the stratigraphy of the 73001/73002 core based on top to bottom variations revealed by coordinated laboratory analyses and explore constraints on the emplacement of the light mantle deposit. Briefly, the upper ∼10 cm of 73002 contains a disturbed zone from space weathering and emplacement of ejecta from a nearby crater that excavated and ejected basaltic material. Below 10 cm is a nearly uniform unit of immature regolith. These data support a single event for the emplacement of the deposit at this location, followed by weathering and mixing of materials from nearby crater ejecta in the upper 10 cm. Slight variations in chemistry and clast components may reflect the relative stratigraphy of the South Massif slope, with material toward the bottom of 73001 originating from lower slopes and material from higher up in the core representing regolith from higher up the South Massif slopes.

58 GEOSCIENCES

Shock compression of FeOOH and implications for iron-water interactions in super-earth magma oceans

Iron(Fe)-water reactions in a magma ocean can influence water storage and density of planets. These reactions can form Fe-O-H phases, whose density, melting, and electronic properties at planetary interior conditions are important for informing planetary models. Here, we study natural goethite (α-FeOOH) that is shock-compressed along its principal Hugoniot. Analysis of our velocity interferometer system for any reflector (VISAR) results extends the equation of state to over 800 GPa. X-ray diffraction and VISAR reflectivity results indicate the onset of melting occurs at ~95 GPa with complete melting by 166 GPa, which may be relevant to low seismic velocity anomalies observed above the core-mantle boundary. Analysis of X-ray emission spectroscopy results up to 285 GPa shows the spin crossover of Fe, with dominantly low spin Fe above ~265 GPa in the melt, supporting formation of dense basal magma oceans in terrestrial planets. Using our measured FeOOH densities, we model planetary interiors up to 10 Earth masses. Assuming FeOOH forms via iron-water reactions, the radius decreases by up to 28%, while the density increases by up to 165% compared to the unreacted cases, providing an avenue to investigate water storage and evolution in super-Earths and sub-Neptunes.

Zhang, Yanyao [Stanford Univ., CA (United States)]

Building wet planets through high-pressure magma–hydrogen reactions

Close-in transiting sub-Neptunes are abundant in our Galaxy. Planetary interior models based on their observed radius–mass relationship suggest that sub-Neptunes contain a discernible amount of either hydrogen (dry planets) or water (wet planets) blanketing a core composed of rocks and metal. Water-rich sub-Neptunes have been believed to form farther from the star and then migrate inwards to their present orbits. Here we report experimental evidence of reactions between warm, dense hydrogen fluid and silicate melt that release silicon from the magma to form alloys and hydrides at high pressures. We found that oxygen liberated from the silicate melt reacts with hydrogen, producing an appreciable amount of water up to a few tens of weight per cent, which is much greater than previously predicted based on low-pressure ideal gas extrapolation. Consequently, these reactions can generate a spectrum of water contents in hydrogen-rich planets, with the potential to reach water-rich compositions for some sub-Neptunes, implying an evolutionary relationship between hydrogen-rich and water-rich planets. Therefore, detection of a large amount of water in exoplanet atmospheres may not be the optimal evidence for planet migration in the protoplanetary disk, calling into question the assumed link between composition and planet formation location.

Horn, H. W. [Arizona State University, Tempe, AZ (

Thermally induced structural competitiveness and metastability of body-centered-cubic iron under nonequilibrium conditions

The structure and stability of iron near melting at multi-megabar pressures are of significant interest in high pressure physics and earth and planetary sciences. While the body-centered cubic (BCC) phase is generally recognized as unstable at lower temperatures, its stability relative to the hexagonal close-packed (HCP) phase at high temperatures (approximately 0.5 eV) in the Earth’s inner core (IC) remains a topic of ongoing theoretical and experimental debate. Motivated by a series of studies attempting to understand the state of the IC and resulting in contradictory findings, we aim to explore the physical factors that may contribute to this confusion in a broader context. Our ab initio calculations show a significant drop in energy, the emergence of a plateau and a local minimum in the potential energy surface, and stabilization of all phonon modes of BCC iron at elevated electron temperatures (>1–1.5 eV). These effects increase the competition among the BCC, HCP, and the face-centered cubic (FCC) phases and lead to the metastability of the BCC structure. Furthermore, at pressure-temperature conditions relevant to the IC, we can expect substantial coexistence of FCC and HCP iron due to their similar free energies and low transition barrier (∼0.2 eV/atom), and the thermodynamic stability of BCC iron is enhanced by its substantial lattice vibration entropy. In conclusion, these findings provide a clear theoretical framework for understanding iron phase relations and solidification processes, both experimentally and in the IC.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND

A small core in Vesta inferred from Dawn’s observations

Vesta’s large-scale interior structure had previously been constrained primarily using the gravity and shape data from the Dawn mission. However, these data alone still allow a wide range of possibilities for the differentiation state of the body. The moment of inertia is arguably the most diagnostic parameter related to the radial density distribution of a planetary body, making it crucial for assessing the body’s state of internal differentiation. Determining the moment of inertia requires additional measurements of the amplitudes of small rotational motions, such as precession and nutation. Here we report an updated estimate of the moment of inertia of Vesta inferred from Dawn’s Doppler tracking via the Deep Space Network and onboard imaging data. The recovered value for Vesta’s normalized polar moment of inertia is $\overline{C}$/MR 2 = 0.4208 ± 0.0047 (where M is the mass of Vesta and R is the reference radius), which is only 6.6% lower than the homogeneous value of 0.4505. This value, combined with the gravity field and global shape, suggests that Vesta’s interior has limited density stratification beneath its howardite–eucrite–diogenite-dominated crust. We propose two possible origin scenarios that are consistent with the observed constraints. In the first scenario, Vesta’s interior did not undergo full differentiation due to late accretion. In the second scenario, Vesta originated as an impact remnant of a larger differentiated body re-accreted with non-chondritic bulk composition produced from a catastrophic impact. Vesta did not experience complete differentiation in either scenario, suggesting that its current state reflects a complex interplay between its accretion timing, thermal evolution, redistribution of 26 Al bearing melt and/or impact processes.

CNEOS 2014-01-08 bolide

The Psyche Gamma-Ray and Neutron Spectrometer

A Gamma-Ray and Neutron Spectrometer (GRNS) instrument has been developed as part of the science payload for NASA’s Discovery Program Psyche mission to the M-class asteroid (16) Psyche. The GRNS instrument is designed to measure the elemental composition of Psyche with the goal to understand the origin of this mysterious, potentially metal-rich planetary body. The GRNS will measure the near-surface abundances for the elements Ni, Fe, Si, K, S, Al, and Ca, as well as the spatial distribution of Psyche’s metal-to-silicate fraction (or metal fraction). These measurements address three of the five Psyche mission science objectives: determine if Psyche is a core; determine whether small metal bodies incorporate light elements into the metal phase; and determine whether Psyche was formed under reducing conditions. The Gamma-Ray Spectrometer (GRS) uses a cryocooled, high-purity Ge (HPGe) sensor to detect cosmic-ray generated gamma rays in the 60 to 9000-keV energy range. The HPGe sensor is surrounded by a borated plastic anticoincidence shield that provides three functions: active background rejection from charged particle interactions in the HPGe sensor; fast neutron measurements; and direct measurements of the incident galactic cosmic ray flux. The Neutron Spectrometer (NS) uses three 3 He gas proportional sensors, each with different material wraps to measure thermal (<0.4 eV), low-energy epithermal (0.4 eV to 1 keV), and high-energy epithermal (up to 100 keV) neutrons. This paper provides an overview of the Psyche GRNS, including: its science and measurement objectives; the design of the instrument hardware, software, and operation; pre-launch performance measurements and its initial performance in space; and an overview of its data products and expected operation for different Psyche mission phases.

Engineering - Instrumentation related to nuclear s