Search NASA⌕ Search

SEARCH · Search NASA

Results for “planetary 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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Constraints on Chondrule Generation, Disk Dynamics, and Asteroid Accretion from the Compositions of Carbonaceous Meteorites

The elemental and isotopic compositions of meteorites are expected to reflect several key processes that occurred in the early solar system, including the migration of gas and dust throughout the protoplanetary disk, the formation of chondrules, and the accretion of the first planetary bodies. However, the specific origins of the various compositions measured among these rocks are currently poorly constrained, limiting our understanding of these processes. Here, we use previously measured elemental and isotopic compositions of chondrites and iron meteorites to identify that carbonaceous (CC) meteorites are mixtures of noncarbonaceous (NC) material, calcium–aluminum-rich inclusion (CAI) material, and CI (Ivuna-like) material, in varying proportions. These trends indicate that chondrules in CO (Ornans-like), CM (Mighei-like), CV (Vigarano-like), and TL (Tagish Lake) chondrites share near-identical average proportions of CI material, arguing that they were generated through the remelting of preexisting NC chondrules all in the same disk environment. Because this proportion likely evolved over space and time throughout the disk, this similarity argues that these chondrules originate from a restricted spatial region and time interval, favoring their generation through a localized event. Moreover, the compositions of CR (Renazzo-like) chondrites indicate that their constituents formed through mechanisms different from those in CO, CM, CV, and TL chondrites. Furthermore, the recovered proportions of CI material in CC iron meteorites and chondrites together also argue for evolution in either the predominant direction of dust and gas motion in the first ~10 au of the disk or the radial distance of planetesimal accretion throughout the CC reservoir.

58 GEOSCIENCES↗
DOE OSTI2025

Improving the Capabilities and Computational Efficiency of the RTE+RRTMGP Radiation Code (Final Report)

This report details progress on the RTE+RRTMGP radiation codes made during the period of performance. RTE+RRTMGP is a set of codes for computing radiative fluxes in planetary atmospheres. RRTMGP uses a k-distribution to provide an optical description (absorption and possibly Rayleigh optical depth) of the gaseous atmosphere, along with the relevant source functions, on a pre-determined spectral grid given temperatures, pressures, and gas concentration. RTE computes fluxes given spectrally-resolved optical descriptions and source functions. Spectrally-resolved fluxes are summarized (“reduced”) via a user extensible class. The initial release of the code and the design choices are described in Pincus et al. 2019; the codes are available on Github. Although RRTMGP was based on current (at the time) empirical spectroscopic data, RTE and RRTMGP were developed in large part to modernize software practices. The design focused on flexibility broadly interpreted: by separating code from data and allowing data to drive computation; in coupling to the host model (e.g. the coupling of clouds to radiative fluxes is user-controlled); with respect to programming languages (computational tasks are accessed via widely-compatible C interfaces); and with respect to hardware (the codes run on a range of CPU and GPU architectures). The code also puts an emphasis on modularity and clarity. RTE+RRTMGP v1.0 was released in September 20219. This award supported the evolution of the RTE+RRTMGP code base to support greater flexibility, accuracy, and efficiency.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Planet Occurrence Rate Correlated to Stellar Dynamical History: Evidence from Kepler and Gaia

The dynamical history of stars influences the formation and evolution of planets significantly. To explore the influence of dynamical history on the planet formation and evolution using observations, we assume stars that experienced significantly different dynamical histories tend to have different relative velocities. Utilizing the accurate Gaia–Kepler Stellar Properties Catalog, we select single main-sequence stars and divide these stars into three groups according to their relative velocities, i.e., high-V, medium-V, and low-V stars. After considering the known biases from Kepler data and adopting prior and posterior correction to minimize the influence of stellar properties on planet occurrence rate, we find that high-V stars have a lower occurrence rate of super-Earths and sub-Neptunes (1–4 R {sub ⊕}, P < 100 days) and a higher occurrence rate of sub-Earth (0.5–1 R {sub ⊕}, P < 30 days) than low-V stars. Additionally, high-V stars have a lower occurrence rate of hot Jupiter-sized planets (4–20 R {sub ⊕}, P < 10 days) and a slightly higher occurrence rate of warm or cold Jupiter-sized planets (4–20 R {sub ⊕}, 10 < P < 400 days). After investigating multiplicity and eccentricity, we find that high-V planet hosts prefer a higher fraction of multiplanet systems and lower average eccentricity, which are consistent with the eccentricity–multiplicity dichotomy of Kepler planetary systems. All of these statistical results favor the scenario that high-V stars with large relative velocity may experience fewer gravitational events, while low-V stars may be influenced by stellar clustering significantly.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

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↗
DOE OSTI2022

“Gray Zone” Simulations Using a Three-Dimensional Planetary Boundary Layer Parameterization in the Weather Research and Forecasting Model

Abstract Generating accurate weather forecasts of planetary boundary layer (PBL) properties is challenging in many geographical regions, oftentimes due to complex topography or horizontal variability in, for example, land characteristics. While recent advances in high-performance computing platforms have led to an increase in the spatial resolution of numerical weather prediction (NWP) models, the horizontal gridcell spacing (Δ x ) of many regional-scale NWP models currently fall within or are beginning to approach the gray zone (i.e., Δ x ≈ 100–1000 m). At these gridcell spacings, three-dimensional (3D) effects are important, as the most energetic turbulent eddies are neither fully parameterized (as in traditional mesoscale simulations) nor fully resolved [as in traditional large-eddy simulations (LES)]. In light of this modeling challenge, we have implemented a 3D PBL parameterization for high-resolution mesoscale simulations using the Weather Research and Forecasting Model. The PBL scheme, which is based on the algebraic model developed by Mellor and Yamada, accounts for the 3D effects of turbulence by calculating explicitly the momentum, heat, and moisture flux divergences in addition to the turbulent kinetic energy. In this study, we present results from idealized simulations in the gray zone that illustrate the benefit of using a fully consistent turbulence closure framework under convective conditions. While the 3D PBL scheme reproduces the evolution of convective features more appropriately than the traditional 1D PBL scheme, we highlight the need to improve the turbulent length scale formulation. Significance Statement The spatial resolution of weather models continues to increase at a rapid rate in accordance with the enhancement of computing power. As a result, smaller-scale atmospheric features become more explicitly resolved. However, most numerical models still ignore the impact of horizontal weather variations on boundary layer flows, which becomes more important at these smaller spatial scales. To address this issue, we have implemented a new modeling approach, using fundamental principles, which accounts for horizontal variability. Our results show that including three-dimensional effects of turbulence is necessary to achieve realistic boundary layer characteristics. This novel technique may be useful for many applications including complex terrain flows, pollutant dispersion, and surface–atmosphere interaction studies.

Juliano, Timothy W.↗
DOE OSTI2025

Best estimate of the planetary boundary layer height from multiple remote sensing measurements

Remote sensing measurements have been widely used to estimate the planetary boundary layer height (PBLHT). Each remote sensing approach offers unique strengths and faces different limitations. In this study, we use machine learning (ML) methods to produce a best-estimate PBLHT (PBLHT-BE-ML) by integrating four PBLHT estimates derived from remote sensing measurements at the Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) observatory. Three ML models – random forest (RF) classifier, RF regressor, and light gradient-boosting machine (LightGBM) – were trained on a dataset from 2017 to 2023 that included radiosonde, various remote sensing PBLHT estimates, and atmospheric meteorological conditions. Evaluations indicated that PBLHT-BE-ML from all three models improved alignment with the PBLHT derived from radiosonde data (PBLHT-SONDE), with LightGBM demonstrating the highest accuracy under both stable and unstable boundary layer conditions. Feature analysis revealed that the most influential input features at the SGP site were the PBLHT estimates derived from (a) potential temperature profiles retrieved using Raman lidar (RL) and atmospheric emitted radiance interferometer (AERI) measurements (PBLHT-THERMO), (b) vertical velocity variance profiles from Doppler lidar (PBLHT-DL), and (c) aerosol backscatter profiles from micropulse lidar (PBLHT-MPL). The trained models were then used to predict PBLHT-BE-ML at a temporal resolution of 10 min, effectively capturing the diurnal evolution of PBLHT and its significant seasonal variations, with the largest diurnal variation observed over summer at the SGP site. We applied these trained models to data from the ARM Eastern Pacific Cloud Aerosol Precipitation Experiment (EPCAPE) field campaign (EPC), where the PBLHT-BE-ML, particularly with the LightGBM model, demonstrated improved accuracy against PBLHT-SONDE. Analyses of model performance at both the SGP and EPC sites suggest that expanding the training dataset to include various surface types, such as ocean and ice-covered areas, could further enhance ML model performance for PBLHT estimation across varied geographic regions.

Zhang, Damao [Pacific Northwest National Laborator↗
DOE OSTI2020

KELT-25 b and KELT-26 b: A Hot Jupiter and a Substellar Companion Transiting Young A Stars Observed by TESS

We present the discoveries of KELT-25 b (TIC 65412605, TOI-626.01) and KELT-26 b (TIC 160708862, TOI-1337.01), two transiting companions orbiting relatively bright, early A stars. The transit signals were initially detected by the KELT survey and subsequently confirmed by Transiting Exoplanet Survey Satellite (TESS) photometry. KELT-25 b is on a 4.40 day orbit around the V = 9.66 star CD-24 5016 (T{sub eff}=8280{sub -180}{sup +440} K, M&lowast; = 2.18{sub -0.11}{sup +0.12} M{sub ⊙}), while KELT-26 b is on a 3.34 day orbit around the V = 9.95 star HD 134004 (T{sub eff} = 8640{sub -240}{sup +500} K, M&lowast; = 1.93{sub -0.16}{sup +0.14} M{sub ⊙}), which is likely an Am star. We have confirmed the substellar nature of both companions through detailed characterization of each system using ground-based and TESS photometry, radial velocity measurements, Doppler tomography, and high-resolution imaging. For KELT-25, we determine a companion radius of R{sub P} = 1.64{sub -0.043}{sup +0.039} R{sub J} and a 3σ upper limit on the companion’s mass of ~64 M{sub J}. For KELT-26 b, we infer a planetary mass and radius of M{sub P} = 1.41{sub -0.51}{sup +0.43} M{sub J} and R{sub P} = 1.94{sub -0.058}{sup +0.060} R{sub J}. From Doppler tomographic observations, we find KELT-26 b to reside in a highly misaligned orbit. This conclusion is weakly corroborated by a subtle asymmetry in the transit light curve from the TESS data. KELT-25 b appears to be in a well-aligned, prograde orbit, and the system is likely a member of the cluster Theia 449.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Biologically Cooling the Planet: A BRISK Approach that Directly Connects the Ground to Deep Space

Climate change is a grave threat to our ability to thrive and prosper on Earth. Accumulation of greenhouse gases in the atmosphere stifle the Earth’s ability to cool itself, which leads to a rise of the average temperature – Global warming. The planetary cooling is exclusively radiative – thermal radiation escapes into cold outer space via a flux of photons at the low-infrared range. The penetration of these photons through Earth’s atmosphere is dictated by atmospheric transparency, which, serendipitously possesses a window centered around the peak of the Plank spectrum for a black body at about room temperature – the average temperature of Earth. This window spanning about 8 --13 μm wavelengths allows most of the photons in this range to escape into the Cosmos carrying our planet’s heat away. The efficiency of this outer heat flux – and thus planet’s cooling ability – can be controlled via emissivity – a function of wavelength which multiplies the fundamental Plank spectrum of a black body at a given temperature. In the presence of non-trivial atmospheric transparency it can be shown that radiative cooling can be enhanced with properly “designed” emissivity function. For the last 100 years, or so, people have been looking into ways to engineer materials to take advantage of such enhancement for various cooling applications, from making ice to cool clothes, and substantial progress has been achieved, including cooling in direct sunlight. The beauty if this approach is that it is entirely passive – no energy input is required. The objective of this project was to investigate whether biological systems possess such advanced cooling capacity, developed via millennia of evolution, to enhance plant’s adaptation to a hot environment, to increase cooling capability, and possibly modulate spectral emissivities on short or long time scales. We wanted to know if BRISK (Biological Radiative Sky Cooling) occurs and to what extent. It is unknown whether biological systems can perform this "spectral redistribution." It is reasonable to expect, however, that they might: in particular, desert plants under high heat stress and arid conditions might evolve this mechanism in order to maintain a tolerable internal temperature. Our intended approach was twofold: extensive literature search and supporting experimental investigations on select plants. We obtained the answers from the former, while the latter proved to be more complicated to perform with required accuracy in the timespan/budget of a small seedling project. After one year of research, we concluded that among the plant species studied there seem to be no evidence of BRISK, but it’s not the end of the story. Our planet’s biosphere is immense, and the evolution is continuing. With our preliminary results and a better understanding of what to look for further discoveries await, considering the accelerated pace of infrared imaging and spectroscopy technology development, from microscopic to satellite-based studies of plants and canopies. In the end, Earth’s energy balance is ultimately tied to radiative cooling so that further investigation into this problem is directly tied to solving the Climate Change challenge.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Lyα Observations of Comet C/2013 A1 (Siding Spring) Using MAVEN IUVS Echelle

The close approach of comet C/2013 A1 (Siding Spring) to Mars in 2014 October provided a unique opportunity to observe a dynamically new Oort cloud comet with potential for interaction with a planet’s atmosphere. The water-originating hydrogen coma of the comet extended to over 20 million km from the nucleus. Determining the properties of this coma contributes to characterizing the comet’s water content and production rate. The present study analyzes a unique data set of high spectral resolution UV observations of comet C/2013 A1 Siding Spring measured by the Mars Atmosphere and Volatile Evolution spacecraft. The Siding Spring observations capture Lyα emissions from the Martian corona, the interplanetary medium, as well as the cometary H and D reservoirs. The isolated cometary spectra are analyzed to reveal a velocity distribution of H atoms that are consistent with model estimates of H{sub 2}O photodissociated H emissions and of OH photodissociated H emissions, Doppler shifted from the main comet H emission line center by 18 km s{sup −1} and 8 km s{sup −1}, respectively. The variations in comet H brightness with distance from the nucleus are used to constrain cometary water production to a rate of 0.5 × 10{sup 28} molecules s{sup −1} at a time when Siding Spring was at 1.5 au, pre-perihelion.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Mesoscale Organization in Cumulus-Coupled Stratocumulus

Marine cloud systems cover a substantial portion of the world’s oceans. Most of these clouds form relatively close to the ocean surface, typically within one to two kilometers, a region referred to by meteorologists as the marine boundary layer. They are composed predominantly of liquid water, although ice particles can occur in mid- and high-latitude marine clouds during winter. In satellite imagery, these clouds appear bright against the darker ocean surface below, reflecting a large fraction of incoming sunlight back into space that would otherwise warm the ocean. Because marine boundary layer clouds cover such an extensive area of the ocean, they exert a significant influence on Earth’s overall transfer of solar energy absorbed by the surface and thermal energy emitted to space, a balance known as the planetary radiation budget. Marine boundary layer clouds are typically thin, and their formation and dissipation depend on a delicate balance between processes acting at the ocean surface below and the warm, dry air above. They are notoriously difficult to simulate accurately in weather forecast models, which often produce too few marine low clouds in midlatitudes and clouds in tropical regions that are excessively bright, meaning they reflect too much solar radiation. The marine boundary layer is frequently characterized by widespread overcast cloud cover that often transitions from a continuous, single-layer deck to more broken cloud fields toward the tropics. These transitions typically proceed through an intermediate stage in which shallow, broken clouds form beneath the overlying stratiform cloud deck. Once broken clouds develop below the overcast, they frequently self-organize into cloud clusters known as marine boundary layer convective complexes (MBLCCs), although the mechanisms governing the formation and organization of MBLCCs remain poorly understood. Accurately representing these transitions in long-range weather forecast models is essential because they influence the properties of air masses advected over the continental United States and Europe, and they become increasingly important for forecasts on seasonal and longer timescales. We employed two complementary approaches to investigate the processes controlling MBLCCs and their impact on marine cloud cover. Long-term observations from the U.S. Department of Energy’s Eastern North Atlantic (ENA) Observatory provided a unique dataset that allowed us to characterize fundamental properties of MBLCCs, including their typical size and frequency of occurrence. These observations were combined with high-resolution numerical simulations performed on supercomputers to examine the evolution of MBLCCs during cold-air outbreaks over the ENA region.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

New Candidates for Planetary-mass Brown Dwarfs in IC 348

We have used infrared images obtained with the Wide Field Camera 3 on board the Hubble Space Telescope to search for planetary-mass brown dwarfs in the star-forming cluster IC 348. In those images, we have identified 12 objects that have colors indicative of spectral types later than M8, corresponding to masses of ≲30 M {sub Jup} at the age of IC 348. The four brightest candidates have been observed with spectroscopy, all of which are confirmed to have late types. Two of those candidates appear to be young, and thus are likely members of the cluster, while the ages and membership of the other two candidates are uncertain. One of the former candidates is the faintest known member of IC 348 in extinction-corrected K {sub s} and is expected to have a mass of 4–5 M {sub Jup} based on evolutionary models and an assumed age of 3 Myr. Four of the remaining eight candidates have ground-based photometry that further supports their candidacy as brown dwarfs, some of which are fainter (and potentially less massive) than the known members.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

High repetition rate mapping of the interaction between a laser plasma and magnetized background plasma via laser induced fluorescence

In this work, the laminar coupling of energy between a laser-produced plasma and a background magnetized plasma was investigated via planar laser induced fluorescence diagnostic and magnetic flux probes. Experiments performed on the Large Plasma Device at the University of California, Los Angeles, mapped out the two-dimensional spatiotemporal evolution of the laser-plasma (debris) ion velocity distribution function (VDF) to assess debris-background coupling in a sub-Alfvénic regime. The acquisition of these data necessitates high repetition rate (1 Hz) as each dataset is the accumulation of thousands of laser shots, which would not be feasible in single-shot experiments. Fully kinetic, three-dimensional particle-in-cell simulations are compared to the measured VDFs to provide a framework in which we can understand the coupling of a sub-Alfvénic plasma flow through a preformed, magnetized plasma. The simulations display the same departure from the expected gyromotion of the debris plasma as observed in the experimental data, and in conjunction with the measured magnetic field traces, have led to the direct observation of the collisionless coupling via laminar fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Oxidation of the Interiors of Carbide Exoplanets

Astrophysical measurements have shown that some stars have sufficiently high carbon-to-oxygen ratios such that the planets they host would be mainly composed of carbides instead of silicates. We studied the behavior of silicon carbide in the presence of water under the high pressure–temperature conditions relevant to planetary interiors in the laser-heated diamond-anvil cell. When reacting with water, silicon carbide converts to silica (stishovite) and diamond at pressures up to 50 GPa and temperatures up to 2500 K: $\mathrm{SiC}+2{{\rm{H}}}_{2}{\rm{O}}\to {\mathrm{SiO}}_{2}+{\rm{C}}+2{{\rm{H}}}_{2}$. Therefore, if water can be incorporated into carbide planets during their formation or through later delivery, they could be oxidized and have mineralogy dominated by silicates and diamond in their interiors. The reaction could produce CH 4 at shallower depths and H 2 at greater depths that could be degassed from the interior, causing the atmospheres of the converted carbon planets to be rich in reducing gases. Excess water after the reaction can be stored in dense silica polymorphs in the interiors of the converted carbon planets. Such conversion of mineralogy to diamond and silicates would decrease the density of carbon-rich planet, making the converted planets distinct from silicate planets in mass–radius relations for the 2–8 Earth mass range.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Polarization and variability of compact sources measured in Planck time-ordered data

This paper introduces a new Planck Catalog of Polarized and Variable Compact Sources (PCCS-PV) comprising 153 sources, the majority of which are extragalactic. The data include both the total flux density and linear polarization measured by Planck with frequency coverage from 30 to 353 GHz, and temporal spacing ranging from days to years. We classify most sources as beamed, extragalactic radio sources; the catalog also includes several radio galaxies, Seyfert galaxies, and Galactic and Magellanic Cloud sources, including H II regions and planetary nebulae. An advanced extraction method applied directly to the multifrequency Planck time-ordered data, rather than the mission sky maps, was developed to allow an assessment of the variability of polarized sources. Our analysis of the time-ordered data from the Planck mission, tod2flux, allowed us to catalog the time-varying emission and polarization properties for these sources at the full range of polarized frequencies employed by Planck, 30-353 GHz. PCCS-PV provides the time-and frequency-dependent, polarized flux densities for all 153 sources. To illustrate some potential applications of the PCCS-PV, we conducted preliminary comparisons of our measurements of selected sources with published data from other astronomical instruments. In summary, we find general agreement between the Planck and the Institut de Radioastronomie Millimétrique (IRAM) polarization measurements as well as with the Metsähovi 37 GHz values at closely similar epochs. Here, these combined measurements also show the value of PCCS-PV results and the PCCS2 catalog for filling in missing spectral (or temporal) coverage and helping to define the spectral energy distributions of extragalactic sources. In turn, these results provide useful clues as to the physical properties of the sources.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Mars Upper Atmospheric Temperature and Atomic Oxygen Density Derived from the O i 130.4 nm Emission Observed by NASA’s MAVEN Mission

The abundance of atomic oxygen is a major unknown for the Mars upper atmosphere, and systematic measurements of this key species is a primary objective of NASA’s Mars Atmosphere and Volatile Evolution (MAVEN) mission. Moreover, the Martian dayside temperature structure and its variability has been a subject of considerable debate in the past. In this study, an inversion analysis of the periapsis limb scans of the O i 130.4 nm emission observed by MAVEN is performed to quantify the neutral temperature and atomic oxygen abundance in the Mars upper atmosphere. Our analysis yields unprecedented resolution of the spatial variation of Mars atomic oxygen, showing that its exobase density decreases with the solar zenith angle (SZA) from ∼0° to 60° by about an order of magnitude, from ∼10{sup 8} to ∼10{sup 7} cm{sup −3} near perihelion (2015 April) and from ∼10{sup 7} to ∼10{sup 6} cm{sup −3} near aphelion (2017 November). For a given SZA the exobase densities near perihelion are about an order of magnitude higher than those near aphelion. The [O]/[CO{sub 2}] mixing ratio at 130 km altitude can vary from ∼0.1% to 4.6%, which also decreases with increasing SZA. Moreover, it is shown that the Mars exobase temperature does not vary significantly with SZA from ∼0° to ∼60°, which is estimated to be ∼239.2 ± 27.6 K near perihelion and ∼162.4 ± 19.9 K near aphelion, with a difference of 76.8 K between these two seasonal extremes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A Laboratory-driven Multiscale Investigation of X-Ray Induced Mass Loss and Photochemical Evolution in Cosmic Carbon and Silicate Dust

We present the results of an integrated laboratory and modeling investigation into the impact of stellar X-rays on cosmic dust. Carbonaceous grains were prepared in a cooled (< 200 K) supersonic expansion from aromatic molecular precursors, and were later irradiated with 970 eV X-rays. Silicate (enstatite) grains were prepared via laser ablation, thermally annealed, and later irradiated with 500 eV X-rays. Infrared spectra of the 3.4 μm band of the carbon sample prepared with benzene revealed 84% ± 5% band area loss for an X-ray dose of 5.2 x 10 23 eV.cm -2 . Infrared spectra of the 8 -12 μm Si-O band of the silicate sample revealed band area loss up to 63% ± 5% for doses of 2.3 x 10 23 eV.cm -2 . A hybrid Monte Carlo particle trajectory approach was used to model the impact of X-rays and ensuing photoelectrons, Auger and collisionally ionized electrons through the bulk. As a result of x-ray ionization and ensuing Coulomb explosions on surface molecules, the calculated mass loss is 60% for the carbonaceous sample and 46% for the silicate sample, within a factor of 2 of the IR band loss, supporting an X-ray induced mass-loss mechanism. We apply the laboratory X-ray destruction rates to estimate the lifetimes of dust grains in protoplanetary disks surrounding 1 M ⊙ and 0.1 M ⊙ G and M stars. In both cases, X-ray destruction timescales are short (a few million years) at the disk surface, but are found to be much longer than typical disk lifetimes (≳10 Myr) over the disk bulk.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Simulations of axion minihalos

The axion, motivated as a solution to the strong CP problem, is also a viable dark matter candidate. The axion field takes random values in causally disconnected regions if the symmetry breaking that establishes the particle occurs after inflation, leading to white-noise density fluctuations at low wave numbers and forming dense minihalos with sub-planetary masses subsequently. There have been two recent proposals that appear capable of testing this scenario, namely using pulsar timing arrays and studying cosmological microlensing caustics. Motivated by these proposals, we use N-body simulations to study the formation of substructures from white-noise density fluctuations. The density profiles of our relaxed axion minihalos can be described by the Navarro–Frenk–White profile, and the minihalos’ concentration number agrees well with a simple, physically-motivated model. We develop a semi-analytic formula to fit the mass function from our simulation, which agrees broadly at different redshifts and only differs at factor of two level from classic halo mass functions. This analytic mass function allows us to consider uncertainties in the post-inflation axion scenario, as well as extrapolate our high-redshift simulations results to the present. Our work estimates the present-day abundance of axion substructures, as is necessary for predicting their effect on cosmological microlensing caustics and pulsar timing. Our calculations suggest that if pulsar timing and microlensing probes can reach recent sensitivity forecasts, they may be sensitive to the post-inflation axion dark matter scenario, even when accounting for uncertainties pertaining to axion strings. For pulsar timing, the most significant caveat is whether axion minihalos are disrupted by stars, which our estimates show is mildly important at the most relevant masses. In conclusion, as our gravitational simulations are scale invariant, the results can be extended to models where the dark matter is comprised of other axion-like particles and even clusters of primordial black holes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Subsurface fluvial sediments beneath InSight on Mars from geophysical constraints

Subsurface structure investigation on Mars is crucial for understanding its geological evolution and past hydrological conditions. Elysium Planitia (EP), located near the hypothesized ancient ocean shorelines, could contain clues for past water activity and paleoclimate. Here we present better-constrained subsurface models beneath InSight extending to ~800 m depth, obtained from joint inversion of seismic and seismoacoustic coupling data, and use the well-resolved subsurface structure to explore the lithological profile through rock physics models. The derived subsurface lithology agrees well with local geological context and exhibits a shallow 60-m-thick low-rigidity layer consistent with hydrated sedimentary materials. Despite possible contributions of aeolian and volcanic deposits, we favor the interpretation that the low-rigidity layer originated from fluvial activity in EP during the Hesperian or Hesperian-to-Amazonian epoch, as supported by adjacent paleo-shoreline morphology observations. These results hint at a period of warmer paleoclimate at low latitudes, possibly during high-obliquity phases of Mars’ rotational axis.

58 GEOSCIENCES↗