Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Planetary”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Dimorphos’s Material Properties and Estimates of Crater Size from the DART Impact

On 2022 September 26, the Double Asteroid Redirection Test (DART) spacecraft intentionally collided with Dimorphos, the moon of the binary asteroid system 65803 Didymos. This collision provided the first full-scale test of a kinetic impactor for planetary defense. Images from DART’s DRACO camera revealed Dimorphos to be an oblate spheroid covered in boulders of varying sizes and shapes. Very little was known about Dimorphos prior to DART’s impact, including its shape, structure, and material properties. Approach observations and those following the DART impact have provided crucial knowledge that narrows the parameter space relevant to modeling the impact into Dimorphos. Here we present the results of a suite of hydrocode simulations of the DART impact on Dimorphos. Despite remaining uncertainties, initial models of DART’s kinetic impact provide important information about the results of DART (e.g., potential crater size and morphology, ejecta mass) and the properties of Dimorphos. Simulations here suggest that Dimorphos has near-surface strength ranging from a few Pascals to tens of kPa, which corresponds to crater sizes of ~40–60 m. Simulated crater sizes provide a crucial comparison metric for the European Space Agency Hera mission when it arrives at the Didymos system. Hera’s measurement of crater size in combination with measurement of Dimorphos’s mass will allow us to assess our simulations and provide the information needed to make the DART impact experiment both the first test of a planetary defense mitigation mission and the first full-scale planetary defense simulation validation exercise.

36 MATERIALS SCIENCE↗
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↗
DOE OSTI2024

Air-coupled tsunamis generated from impacts and airbursts: Our understanding before Hunga-Tonga Hunga-Ha'apai

The effort to prevent or mitigate the effects of an impact on Earth is known as planetary defense. A significant component of planetary defense research involves risk assessment. Much of our understanding of the risk from near-Earth objects comes from the geologic record in the form of impact craters, but not all asteroid impacts are crater-forming events. Small asteroids explode before reaching the surface, generating an airburst, and most impacts into the ocean do not penetrate the water to form a crater in the sea floor. The risk from these non-crater-forming ocean impacts and airbursts is difficult to quantify and represents a significant uncertainty in our assessment of the overall threat. We are currently working to better understand impact scenarios that can generate dangerous tsunamis. One of the suggested mechanisms for the production of asteroid–generated tsunamis is by direct coupling of the pressure wave to the water, analogous to the means by which a moving weather front can generate a meteotsunami. To test this hypothesis, we ran a series of airburst simulations and provided time-resolved pressure and wind profiles to use as source functions for tsunami models. We used the CTH hydrocode to model the various airburst scenarios to compare to the results of other simulations and provide time dependent boundary conditions as input to shallow-water wave propagation codes. The strongest and most destructive meteotsunamis are generated by atmospheric pressure oscillations with amplitudes of only a few hPa1 (mbar), corresponding to changes in sea level of a few cm. The resulting wave is strongest when there is a resonance between the ocean and the atmospheric forcing. A Proudman resonance takes place when the atmospheric disturbance’s translational speed (U) equals the longwave phase speed $\sqrt{gh}$ of shallow water wave. Coupling is strongest when the Froude number (Fr=U/c) is unity. A weather front propagates much slower than the speed of sound, so meteotsunamis are most common and dangerous in shallow bodies of water such as the Mediterranean Sea or Lake Michigan. By contrast, the blast wave from an airburst or crater-forming impact propagates at a speed faster than a tsunami in the deepest ocean, and a Proudman resonance cannot be achieved even though the overpressures are orders of magnitude greater. However, blast wave profiles are N-waves in which a sharp shock wave leading to overpressure is followed by a more gradual rarefaction to a much longer-duration underpressure phase. Even though the blast outruns the water wave it is forcing, the tsunami should continue to be driven by the out-of-resonance gradient associated with the suction phase, which may depend strongly on the details of the airburst or impact scenario. The open question is whether there are any conditions under which such an airburst-driven tsunami can be dangerous enough to contribute to the overall impact risk. We have also identified other potential mechanisms for airburst-generated tsunamis: 1) reaction force at the surface from the plume ejected into space, which carries significant momentum, 2) expanding toroidal vortices at the surface, which travel more slowly than the shock wave and can generate a Proudman resonance in relatively shallow ocean (such as continental shelf), and 3) steam explosion from seawater ablation by a “Type II” (Libyan Desert Glass-type) airburst in which the hot vapor jet descends to the surface. On January 15, 2022, the Hunga-Tonga Hunga-Ha’apai volcano, located approximately 60 km north of Tongatapu, the main island of Tonga, violently erupted with a powerful explosion, culminating the period of volcanic activity that started in December of 2021. This event and resulting tsunamis provided an existence proof for the air pressure wave coupling mechanism we proposed. It also suggests that it can be stronger and more significant over much greater distances than we contemplated, leading to global tsunamis associated with impact events on land as well as in the water. Large atmospheric explosions generate global Lamb waves with larger amplitudes, longer periods, and slower speeds than the local and regional blast waves we modeled prior to that event. This paper reviews our analysis and modeling of airburst-driven tsunamis prior to the 2022 Hunga-Tonga Hunga-Ha’apai tsunami, which was the subject of two presentations at the 2023 Planetary Defense Conference and is the subject of another paper currently in preparation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

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)]↗
DOE OSTI2026

Final Technical Report for DoE award DE‐SC0023367 “Energetic Electron Transport in Magnetized Plasma with Magnetic Islands”

This project investigated how plasmas interact with energetic particles and solid materials under extreme conditions relevant to fusion energy, space plasmas, and planetary environments. Using experiments on the DIII-D National Fusion Facility, the research first examined how high-energy electrons move, become trapped, and are released in plasmas containing magnetic islands—structures commonly found in fusion reactors and Earth’s magnetosphere—providing new insight into particle transport and acceleration processes. The project also explored plasma-driven chemical reactions that can occur during meteoroid entry into planetary atmospheres, demonstrating that simple molecules such as ammonia can be produced and survive in high-temperature plasma conditions. Together, these results improve understanding of plasma behavior across laboratory, space, and planetary systems while informing fusion plasma control and plasma–material interaction studies. The project additionally contributed to workforce development by training graduate students, undergraduates, and early-career researchers and by disseminating results through peer-reviewed publications and international scientific conferences.

Orlov, Dmitri Mikhailovich [UC San Diego] (ORCID:0↗
DOE OSTI2024

Air-coupled tsunamis generated from impacts and airbursts: Our understanding before Hunga-Tonga Hunga-Ha'apai

The effort to prevent or mitigate the effects of an impact on Earth is known as planetary defense. A significant component of planetary defense research involves risk assessment. Much of our understanding of the risk from near-Earth objects comes from the geologic record in the form of impact craters, but not all asteroid impacts are crater-forming events. Small asteroids explode before reaching the surface, generating an airburst, and most impacts into the ocean do not penetrate the water to form a crater in the sea floor. The risk from these non-crater-forming ocean impacts and airbursts is difficult to quantify and represents a significant uncertainty in our assessment of the overall threat. We are currently working to better understand impact scenarios that can generate dangerous tsunamis.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

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↗
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 OSTI2025

Advancing Aerosol Chemical Characterization and Vertical Profiling over the Southern Great Plains Using Uncrewed Aerial Sampling and Offline Aerosol Mass Spectrometry

Recent advancements in uncrewed aerial systems (UASs) and particulate matter (PM) analytical techniques have provided opportunities for atmospheric research. In this study, we deployed the Department of Energy’s fixed-wing ArcticShark UAS to examine PM 2.5 composition at varying altitudes─within and above the planetary boundary layer (PBL)─over the Southern Great Plains atmospheric observatory (SGP). A total of 22 flights were conducted across March, June, and August 2023. Composite filter samples were collected during each flight and analyzed with offline aerosol mass spectrometry (AMS), complemented by on-board real-time sensors and ground-based instrumentation, to provide a comprehensive view of regional aerosol characteristics. Results show clear vertical and seasonal differences in the aerosol composition. Relative to ground-level measurements, aloft samples exhibited shifts in the distribution of organic and inorganic PM, with the organic composition varying distinctly across seasons. Particulate organic nitrogen (ON) was elevated, with bulk compositions similar in March and June but strongly altered in August, likely driven by biomass burning and enhanced photochemical activity. Combined AMS and chemical ionization mass spectrometry analyses detected amines, amides, and amino acids. PM above the planetary boundary layer was enriched in oxidized organic aerosols, while ground-level PM contained higher nitrate and sulfate. Seasonal differences in aqueous-phase processing were also observed, which were strongest in March during persistent cloud cover and weaker in the drier August period, suggesting a shift from aqueous- to gas-phase SOA formation. In conclusion, these findings highlight the value of UAS in advancing PM measurements and vertical profiling of aerosol composition.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

A Scale‐Dependent Analysis of the Barotropic Vorticity Budget in a Global Ocean Simulation

Abstract The climatological mean barotropic vorticity budget is analyzed to investigate the relative importance of surface wind stress, topography, planetary vorticity advection, and nonlinear advection in dynamical balances in a global ocean simulation. In addition to a pronounced regional variability in vorticity balances, the relative magnitudes of vorticity budget terms strongly depend on the length‐scale of interest. To carry out a length‐scale dependent vorticity analysis in different ocean basins, vorticity budget terms are spatially coarse‐grained. At length‐scales greater than 1,000 km, the dynamics closely follow the Topographic‐Sverdrup balance in which bottom pressure torque, surface wind stress curl and planetary vorticity advection terms are in balance. In contrast, when including all length‐scales resolved by the model, bottom pressure torque and nonlinear advection terms dominate the vorticity budget (Topographic‐Nonlinear balance), which suggests a prominent role of oceanic eddies, which are of km in size, and the associated bottom pressure anomalies in local vorticity balances at length‐scales smaller than 1,000 km. Overall, there is a transition from the Topographic‐Nonlinear regime at scales smaller than 1,000 km to the Topographic‐Sverdrup regime at length‐scales greater than 1,000 km. These dynamical balances hold across all ocean basins; however, interpretations of the dominant vorticity balances depend on the level of spatial filtering or the effective model resolution. On the other hand, the contribution of bottom and lateral friction terms in the barotropic vorticity budget remains small and is significant only near sea‐land boundaries, where bottom stress and horizontal viscous friction generally peak.

Khatri, Hemant↗
DOE OSTI2025

Mechanisms in Regulating the Quasi‐Biennial Oscillation in Exascale Earth System Model Version 2

Abstract The update to the deep convection parameterization in the energy exascale Earth system model version 2 (E3SMv2) makes the simulated quasi‐biennial oscillation (QBO) have a shorter period than without the update when other tunable parameters the same since the update makes convection more intense but less frequent while leaving the time‐mean convective heating tendency almost unchanged. Amplitudes of momentum fluxes of parameterized gravity waves (GWs) are intensified since they are determined by the square of convective heating tendency. In addition, stronger planetary waves are simulated in E3SMv2 with the convection scheme update, partially due to the tropospheric precipitation change. Furthermore, there is evidence that planetary waves are intensified in the stratosphere due to the dissipation of enhanced parameterized GWs. These factors are found to be responsible for modulating the QBO simulated in E3SMv2.

Li, Yuanpu [Climate and Global Dynamics Division U↗
DOE OSTI2024

Structural evolution of liquid silicates under conditions in Super-Earth interiors

Molten silicates at depth are crucial for planetary evolution, yet their local structure and physical properties under extreme conditions remain elusive due to experimental challenges. In this study, we utilize in situ X-ray diffraction (XRD) at the Matter in Extreme Conditions (MEC) end-station of the Linear Coherent Linac Source (LCLS) at SLAC National Accelerator Laboratory to investigate liquid silicates. Using an ultrabright X-ray source and a high-power optical laser, we probed the local atomic arrangement of shock-compressed liquid (Mg,Fe)SiO 3 with varying Fe content, at pressures from 81(9) to 385(40) GPa. We compared these findings to ab initio molecular dynamics simulations under similar conditions. Results indicate continuous densification of theO-O and Mg-Si networks beyond Earth’s interior pressure range, potentially altering melt properties at extreme conditions. This could have significant implications for early planetary evolution, leading to notable differences in differentiation processes between smaller rocky planets, such as Earth and Venus, and super-Earths, which are exoplanets withmasses nearly three times that of Earth.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Contrasting biological production trends over land and ocean

Terrestrial and marine ecosystems constitute the primary components of the Earth’s biosphere, yet their photosynthetic productions are typically studied separately, which limits understanding of planetary carbon uptake and biosphere health. Here, in this study, using multiple satellite-derived products, we identify contrasting net primary production (NPP) trends between land and ocean, probably reflecting their differential sensitivity to climate warming, especially in tropical regions. Planetary NPP shows an overall increase of 0.11 ± 0.13 PgC yr −1 (P = 0.05) from 2003 to 2021, driven by a significant terrestrial enhancement of 0.20 ± 0.07 PgC yr −1 (P < 0.001) and partially offset by an oceanic decline of −0.12 ± 0.12 PgC yr −1 (P = 0.07). While land contributes to the strong upwards NPP trend, the interannual variability in global NPP is predominantly driven by the ocean, especially during strong El Niño–Southern Oscillation events. Our findings highlight the resilience and potential vulnerability of biosphere primary productivity in a warming climate, calling for integrated land–ocean monitoring and assessment to support climate mitigation initiatives.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Exploring toroidal anvil profiles for larger sample volumes above 4 Mbar

Abstract With the advent of toroidal and double-stage diamond anvil cells (DACs), pressures between 4 and 10 Mbar can be achieved under static compression, however, the ability to explore diverse sample assemblies is limited on these micron-scale anvils. Adapting the toroidal DAC to support larger sample volumes offers expanded capabilities in physics, chemistry, and planetary science: including, characterizing materials in soft pressure media to multi-megabar pressures, synthesizing novel phases, and probing planetary assemblages at the interior pressures and temperatures of super-Earths and sub-Neptunes. Here we have continued the exploration of larger toroidal DAC profiles by iteratively testing various torus and shoulder depths with central culet diameters in the 30–50 µm range. We present a 30 µm culet profile that reached a maximum pressure of 414(1) GPa based on a Pt scale. The 300 K equations of state fit to our P–V data collected on gold and rhenium are compatible with extrapolated hydrostatic equations of state within 1% up to 4 Mbar. This work validates the performance of these large-culet toroidal anvils to > 4 Mbar and provides a promising foundation to develop toroidal DACs for diverse sample loading and laser heating.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Magnetohydrodynamic simulations preliminarily predict the habitability and radio emission of TRAPPIST-1e

TRAPPIST-1e, an Earth-sized exoplanet in the habitable zone of the nearby M dwarf TRAPPIST-1, may experience magnetospheric responses that vary with stellar space weather, which could potentially influence both its habitability and radio emissions. Our objective is to investigate how different Earth-like magnetospheric configurations of TRAPPIST-1e – specifically variations in dipolar magnetic field strength and axial tilt – respond to diverse stellar space weather conditions, including events analogous to coronal mass ejections (CMEs), and to assess their implications for potential habitability and expected radio emissions. We conducted 3D magnetohydrodynamic simulations of the TRAPPIST-1e system using the PLUTO code in spherical coordinates. The planetary magnetic field was modelled as dipolar, with equatorial strengths from Earth-like to several times stronger. The dipole axis spans a representative range of axial tilts. We investigate four stellar wind environments, from sub-Alfvénic flow to CME-like disturbances. Planetary shielding was quantified based on the magnetopause standoff distance, and radio powers were estimated via empirical scaling laws. Our simulations show that both shielding and radio power depend strongly on the magnetic configuration. Stronger fields increase protection, while larger tilts reduce it. Radio power increases with both field strength and tilt across all wind regimes. An Earth-like magnetic field can provide effective shielding even under intense CMEs, whereas high tilts require stronger fields. Predicted radio powers reach ~10 20 erg s −1 during CMEs, making bursts from close-in, magnetised planets more detectable. However, for TRAPPIST-1e, the maximum cyclotron frequency lies below the Earth’s ionospheric cutoff (~10 MHz), making ground-based detection currently infeasible.

Sun: magnetic fields↗
DOE OSTI2025

An accurate measurement of the spectral resolution of the JWST Near Infrared Spectrograph

The spectral resolution (R ≡ λ/Δλ) of spectroscopic data is crucial information for accurate kinematic measurements. In this letter we present a robust measurement of the spectral resolution of the JWST Near Infrared Spectrograph (NIRSpec) in fixed slit (FS) and integral field spectroscopy (IFS) modes. Due to the similarity of the utilized slit dimension in the FS mode to that of the shutters in the multi-object spectroscopy (MOS) mode, our resolution measurements in the FS mode can also be used for the MOS mode in principle. We modeled H and He lines of the planetary nebula SMP LMC 58 using a Gaussian line spread function (LSF) to estimate the wavelength-dependent resolution for multiple disperser and filter combinations. We corrected for the intrinsic width of the planetary nebula’s H and He lines due to its expansion velocity by measuring it from a higher-resolution X-shooter spectrum. We find that NIRSpec’s in-flight spectral resolutions exceed the pre-launch estimates provided in the JWST User Documentation by 11–53% in the FS mode and by 1–24% in the IFS mode across the covered wavelengths. We recover the expected trend that the resolution increases with the wavelength within a configuration. The robust and accurate LSFs presented in this letter will enable high-accuracy kinematic measurements using NIRSpec for applications in cosmology and galaxy evolution.

methods: data analysis↗
DOE OSTI2025

Completion of lunar magma ocean solidification at 4.43 Ga

Crystallization of the lunar magma ocean yielded a chemically unique liquid residuum named KREEP. This component is expressed as a large patch on the near side of the Moon and a possible smaller patch in the northwest portion of the Moon’s South Pole-Aitken basin on the far side. Thermal models estimate that the crystallization of the lunar magma ocean (LMO) could have spanned from 10 and 200 My, while studies of radioactive decay systems have yielded inconsistent ages for the completion of LMO crystallization covering over 160 My. Here, we show that the Moon achieved >99% crystallization at 4,429 ± 76 Ma, indicating a lunar formation age of ~4,450 Ma or possibly older. Using the 176 Lu– 176 Hf decay system (t 1/2 = 37 Gy), we found that the initial 176 Hf/ 177 Hf ratios of lunar zircons with varied U–Pb ages are consistent with their crystallization from a KREEP-rich reservoir with a consistently low 176 Lu/ 177 Hf ratio of 0.0167 that emerged ~140 My after solar system formation. The previously proposed younger model age of ~4.33 Ga for the source of mare basalts (240 My after solar system formation) might reflect the timing of a large impact. Our results demonstrate that lunar magma ocean crystallization took place while the Moon was still battered by planetary embryos and planetesimals leftover from the main stage of planetary accretion. The study of Lu–Hf model ages for samples brought back from the South Pole-Aitken basin will help to assess the lateral continuity of KREEP and further understand its significance in the early history of the Moon.

Science & Technology - Other Topics↗
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)↗