Search NASA⌕ Search

SEARCH · Search NASA

Results for “primitive bodies”

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.

Quantum many-body linear algebra, Hamiltonian moments, and a coupled-cluster inspired framework

Here, we propose a general strategy to develop quantum many-body approximations of primitives in linear algebra algorithms. As a practical example, we introduce a coupled-cluster inspired framework to produce approximate Hamiltonian moments and demonstrate its application in various linear algebra algorithms for ground state estimation. Through numerical examples, we illustrate the difference between the ground-state energies arising from quantum many-body linear algebra and those from the analogous many-body perturbation theory. Our results support the general idea of designing quantum many-body approximations outside of perturbation theory, providing a route to new algorithms and approximations.

Algorithms and data structure↗

The Non-carbonaceous–Carbonaceous Meteorite Dichotomy

The isotopic dichotomy between non-carbonaceous (NC) and carbonaceous (CC) meteorites indicates that meteorite parent bodies derive from two genetically distinct reservoirs, which presumably were located inside (NC) and outside (CC) the orbit of Jupiter and remained isolated from each other for the first few million years of the solar system. Here we review the discovery of the NC–CC dichotomy and its implications for understanding the early history of the solar system, including the formation of Jupiter, the dynamics of terrestrial planet formation, and the origin and nature of Earth’s building blocks. The isotopic difference between the NC and CC reservoirs is probably inherited from the solar system’s parental molecular cloud and has been maintained through the rapid formation of Jupiter that prevented significant exchange of material from inside (NC) and outside (CC) its orbit. The growth and/or migration of Jupiter resulted in inward scattering of CC bodies, which accounts for the co-occurrence of NC and CC bodies in the present-day asteroid belt and the delivery of presumably volatile-rich CC bodies to the growing terrestrial planets. Earth’s primitive mantle, at least for siderophile elements like Mo, has a mixed NC–CC composition, indicating that Earth accreted CC bodies during the final stages of its growth, perhaps through the Moon-forming giant impactor. The late-stage accretion of CC bodies to Earth is sufficient to account for the entire budget of Earth’s water and highly volatile species.

58 GEOSCIENCES↗

ExaTN: Scalable GPU-Accelerated High-Performance Processing of General Tensor Networks at Exascale

We present ExaTN (Exascale Tensor Networks), a scalable GPU-accelerated C++ library which can express and process tensor networks on shared- as well as distributed-memory high-performance computing platforms, including those equipped with GPU accelerators. Specifically, ExaTN provides the ability to build, transform, and numerically evaluate tensor networks with arbitrary graph structures and complexity. It also provides algorithmic primitives for the optimization of tensor factors inside a given tensor network in order to find an extremum of a chosen tensor network functional, which is one of the key numerical procedures in quantum many-body theory and quantum-inspired machine learning. Numerical primitives exposed by ExaTN provide the foundation for composing rather complex tensor network algorithms. We enumerate multiple application domains which can benefit from the capabilities of our library, including condensed matter physics, quantum chemistry, quantum circuit simulations, as well as quantum and classical machine learning, for some of which we provide preliminary demonstrations and performance benchmarks just to emphasize a broad utility of our library.

97 MATHEMATICS AND COMPUTING↗

Discontinuous Galerkin discretization for quantum simulation of chemistry

All-electron electronic structure methods based on the linear combination of atomic orbitals method with Gaussian basis set discretization offer a well established, compact representation that forms much of the foundation of modern correlated quantum chemistry calculations—on both classical and quantum computers. Despite their ability to describe essential physics with relatively few basis functions, these representations can suffer from a quartic growth of the number of integrals. Recent results have shown that, for some quantum and classical algorithms, moving to representations with diagonal two-body operators can result in dramatically lower asymptotic costs, even if the number of functions required increases significantly. We introduce a way to interpolate between the two regimes in a systematic and controllable manner, such that the number of functions is minimized while maintaining a block-diagonal structure of the two-body operator and desirable properties of an original, primitive basis. Techniques are analyzed for leveraging the structure of this new representation on quantum computers. Empirical results for hydrogen chains suggest a scaling improvement from O(N 4.5 ) in molecular orbital representations to O(N 2.6 ) in our representation for quantum evolution in a fault-tolerant setting, and exhibit a constant factor crossover at 15 to 20 atoms. Moreover, we test these methods using modern density matrix renormalization group methods classically, and achieve excellent accuracy with respect to the complete basis set limit with a speedup of 1–2 orders of magnitude with respect to using the primitive or Gaussian basis sets alone. These results suggest our representation provides significant cost reductions while maintaining accuracy relative to molecular orbital or strictly diagonal approaches for modest-sized systems in both classical and quantum computation for correlated systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Organic matter in carbonaceous chondrite lithologies of Almahata Sitta: Incorporation of previously unsampled carbonaceous chondrite lithologies into ureilitic regolith

The Almahata Sitta (AhS) meteorite is a unique polymict ureilite. Recently, carbonaceous chondritic lithologies were identified in AhS. Organic matter (OM) is ubiquitously found in primitive carbonaceous chondrites. The molecular and isotopic characteristics of this OM reflect its origin and parent body processes, and are particularly sensitive to heating. The C1 lithologies AhS 671 and AhS 91A were investigated, focusing mainly on the OM. We found that the OM in these lithologies is unique and contains primitive isotopic signatures, but experienced slight heating possibly by short-term heating event(s). These characteristics support the idea that one or more carbonaceous chondritic bodies were incorporated into the ureilitic parent body. Here, the uniqueness of the OM in the AhS samples implies that there were large variations in primitive carbonaceous chondritic materials in the solar system other than known primitive carbonaceous chondrite groups such as CI, CM, and CR chondrites.

58 GEOSCIENCES↗

Incorporation of 16 O-rich anhydrous silicates in the protolith of highly hydrated asteroid Ryugu

The abundant phyllosilicate and carbonate minerals characterizing most of the returned particles from asteroid Ryugu suggest a history of extensive aqueous alteration on its parent body, similar to the rare mineralogically altered, but chemically primitive, CI (Ivuna-type) chondrite meteorites. Particle C0009 differs mineralogically from other Ryugu particles examined so far by containing anhydrous silicates at a level of ~0.5 vol%, and thus can help shed light on the unaltered original materials that constituted Ryugu’s protolith. In situ oxygen isotope measurements of the most Mg-rich olivine and pyroxene in C0009 reveal two populations of Δ 17 O: –25‰ to –15‰ and –8‰ to –3‰. The former and the latter populations correlate well with silicate morphologies similar to those seen in amoeboid olivine aggregates and chondrule phenocrysts, respectively, both of which are abundant in less aqueously altered carbonaceous chondrites. Here, this result also highlights the presence of olivine with Δ 17 O close to the solar value in either a CI chondrite or an asteroid with CI-chondrite characteristics, and provides strong evidence that amoeboid olivine aggregates and Mg-rich chondrules accreted into Ryugu’s protolith. Our data also raise the possibility that the protoliths of CI and other carbonaceous chondrites incorporated similar anhydrous silicates.

58 GEOSCIENCES↗

Nitrogen- and oxygen-rich organic material indicative of polymerization in pre-aqueous cryochemistry on Bennu’s parent body

Nitrogen-containing organic compounds play key biological roles, and their identification in primitive astromaterials such as meteorites can shed light on the origin of life. However, meteorites are typically contaminated by uncontrolled exposure to Earth. Here we show that pristine samples returned from asteroid Bennu contain polymeric organics exceptionally rich in nitrogen and oxygen. These polymers contain a variety of functional groups including amines, amides, N-heterocycles, and aliphatic and aromatic hydrocarbons, among others. They are seen in a carbonaceous vein with mineral inclusions and in multilayered organic sheets. Their morphology and composition indicate formation from pre-aqueous N-rich precursors and later modification during aqueous alteration. These findings demonstrate that asteroids like Bennu contain complex nitrogen-rich organic phases formed by pre-aqueous and aqueous processes, and they expand the known inventory of potential prebiotic extraterrestrial compounds.

Sandford, Scott A↗

Bayesian constraints on the origin and geology of exoplanetary material using a population of externally polluted white dwarfs

ABSTRACT White dwarfs that have accreted planetary bodies are a powerful probe of the bulk composition of exoplanetary material. In this paper, we present a Bayesian model to explain the abundances observed in the atmospheres of 202 DZ white dwarfs by considering the heating, geochemical differentiation, and collisional processes experienced by the planetary bodies accreted, as well as gravitational sinking. The majority (>60 per cent) of systems are consistent with the accretion of primitive material. We attribute the small spread in refractory abundances observed to a similar spread in the initial planet-forming material, as seen in the compositions of nearby stars. A range in Na abundances in the pollutant material is attributed to a range in formation temperatures from below 1000 K to higher than 1400 K, suggesting that pollutant material arrives in white dwarf atmospheres from a variety of radial locations. We also find that Solar System-like differentiation is common place in exoplanetary systems. Extreme siderophile (Fe, Ni, or Cr) abundances in eight systems require the accretion of a core-rich fragment of a larger differentiated body to at least a 3σ significance, whilst one system shows evidence that it accreted a crust-rich fragment. In systems where the abundances suggest that accretion has finished (13/202), the total mass accreted can be calculated. The 13 systems are estimated to have accreted masses ranging from the mass of the Moon to half that of Vesta. Our analysis suggests that accretion continues for 11 Myrs on average.

79 ASTRONOMY AND ASTROPHYSICS↗

Physics-Driven Construction of Compact Primitive Gaussian Density Fitting Basis Sets

We present a model-assisted density fitting (MADF) basis set generator, an algorithm for generating primitive atomic Gaussian density fitting (DF) basis sets (DFBSs) from a contracted Gaussian orbital basis set (OBS). The MADF algorithm produces DFBSs suitable for accurate robust DF approximation of 2-particle interactions in mean-field and correlated electronic structures. The algorithm is designed to (a) saturate the OBS product space by a large regularized set of primitive solid-harmonic Gaussian shells with nonuniform distribution of exponents, followed by (b) pruning of the shells according to their contributions to the 2- body energy of a correlated atomic ensemble. Building the DFBS generator model almost exclusively on mathematical and physical principles allows one to limit the number of parameters that control the density fitting error to three, with a single set of parameters sufficient for computations with all basis cardinal numbers, with and without correlation of core electrons, with and without scalar and spin-dependent relativistic effects, spanning almost all of the Periodic Table. Performance assessment included basis sets up to quadruple-ζ quality from several major basis set families, using molecules composed of main-group, d-block, and f-block elements. The resulting DF errors in Hartree−Fock and second-order MP2 energies (with relativistic all-electron treatments, when appropriate) were on the order of 20 and 10 μE h per electron, respectively.

Approximation↗

Point containment algorithms for constructive solid geometry with unbounded primitives

Here, we present several algorithms for evaluating point containment in constructive solid geometry (CSG) trees with unbounded primitives. Three algorithms are presented based on postfix, prefix, and infix notations of the CSG binary expression tree. We show that prefix and infix notations enable short-circuiting logic, which reduces the number of primitives that must be checked during point containment. To evaluate the performance of the algorithms, each algorithm was implemented in the OpenMC Monte Carlo particle transport code, which relies on CSG to represent solid bodies through which subatomic particles travel. Two sets of tests were carried out. First, the execution time to generate a rasterized image of a 2D slice of three CSG models of varying complexity was measured. Use of both prefix and infix notations offered significant speedup over the postfix notation that has traditionally been used in particle transport codes, with infix resulting in a 6 x reduction in execution time relative to postfix for a model of a tokamak fusion device. We then measured the execution time of neutron transport simulations of the same three models using each of the algorithms. The results and performance improvements reveal the same trends as for the rasterization test, with a 5.52 x overall speedup using the infix notation relative to the original postfix notation in OpenMC for the tokamak model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Fast inversion, preconditioned quantum linear system solvers, fast Green's-function computation, and fast evaluation of matrix functions

Preconditioning is the most widely used and effective way for treating ill-conditioned linear systems in the context of classical iterative linear system solvers. We introduce a quantum primitive called fast inversion, which can be used as a preconditioner for solving quantum linear systems. The key idea of fast inversion is to directly block encode a matrix inverse through a quantum circuit implementing the inversion of eigenvalues via classical arithmetics. We demonstrate the application of preconditioned linear system solvers for computing single-particle Green's functions of quantum many-body systems, which are widely used in quantum physics, chemistry, and materials science. We analyze the complexities in three scenarios: the Hubbard model, the quantum many-body Hamiltonian in the plane-wave-dual basis, and the Schwinger model. We also provide a method for performing Green's function calculation in second quantization within a fixed-particle manifold and note that this approach may be valuable for simulation more broadly. Aside from solving linear systems, fast inversion also allows us to develop fast algorithms for computing matrix functions, such as the efficient preparation of Gibbs states. Furthermore, we introduce two efficient approaches for such a task, based on the contour-integral formulation and the inverse transform, respectively.

97 MATHEMATICS AND COMPUTING↗

Reducing Memory Consumption in Calico with Shared Memory

This document details the work to reduce memory consumption in Calico. Calico is SimTools’ Constructive Solid Geometry (CSG) and geometry painting library. It is primarily used to paint material volume fractions in the Eulerian meshes of the physics codes. Calico provides point-in-body checks for the geometry supplied by an Oso model, which are then aggregated by the host codes. In addition, Calico can be used to build Oso models and is used by Ingen for that purpose. Oso models, and thus Calico, provide support for various CSG primitives such as spheres, cylinders, surfaces generated by rotating tabular curve data, and STL files as well as binary combinations of those primitives. Prior to refactoring Calico will run out of memory on CTS-1 machines when 36 MPI ranks are used per node when reading STL models on the order of 1.5 GB. This limitation is a bottleneck in designer workflow. This problem has been alleviated through the use of data structures to both reduce memory consumption and to leverage MPI-3 shared memory. This report details the data structures targeted for refactoring in Calico, the methods and implementation details for reducing memory consumption and leveraging shared memory, and results for one test problem. Results show a memory reduction when loading a 1 GB STL file by a factor of 27.5, from 93.4 to 3.4 GB.

97 MATHEMATICS AND COMPUTING↗

Meteor over New York City: Brines in a primitive CM asteroid

The CI (Ivuna-type) carbonaceous material returned from asteroids Ryugu and Bennu contain mobilized sodium from the evaporation or freezing of liquid water into brines, shedding light on the internal structure of ice-rich CI-type worlds and the formation of prebiotic organic compounds. The formation of brines has not been demonstrated in CM (Mighei-type) carbonaceous chondrites, which also supplied organic matter to the early Earth. Here, we announce the fall of a primitive meteorite from a daytime fireball over the New York metropolitan area in July 2024. It is a CM2 breccia that contains unique CM1 clasts rich in water and sodium. The meteorite contains abundant amino acids and other products of organic chemistry in brines that reveal subsurface processes on CM-type asteroid parent bodies.

Geosciences↗

Crystallization histories of the group IIF iron meteorites and Eagle Station pallasites

The group IIF iron meteorites and Eagle Station pallasites (PES) have highly siderophile element abundances (HSE; Re, Os, Ir, Ru, Pt, and Pd) of metal that are consistent with formation in planetesimal cores by fractional crystallization with minor to major solid metal-liquid metal mixing. Modeling of HSE abundances of the IIF irons indicates a complex formation history that included the mixing of primitive and evolved solid and liquid metals. By contrast, modeling of HSE abundances of PES metal suggests these meteorites formed mainly as equilibrium solids from a common liquid. Abundances of some of the siderophile elements in the IIF irons and PES are permissive of a common core origin, however, the abundances of W and Ni indicate the PES ultimately formed on a more oxidized body. The PES most likely formed by the injection of olivine present at the core-mantle boundary into a metallic core liquid as a result of impact. The core then crystallized inward, trapping the olivine.

58 GEOSCIENCES↗

Compositional evidence for chondrule origins of low‐Ca pyroxenes in comet Wild 2 and a giant cluster IDP

Abstract A literature compilation of 1136 low‐Ca pyroxene compositions from chondrules from 12 primitive type 2–3 carbonaceous, ordinary and enstatite chondrite groups define unique regions on an Al 2 O 3 and Cr 2 O 3 diagram when compared to low‐Ca pyroxenes from equilibrated type 4‐6 chondrites. Measured compositions of 100 low‐Ca pyroxenes from comet Wild 2 and a giant cluster IDP of probable cometary origin are similar to each other and fall in the type 2–3 chondrite chondrule region suggesting that most of the pyroxenes likely formed in the solar nebula like conventional chondrules. The data imply that most low Ca‐pyroxenes from comet Wild 2 and the giant cluster IDP formed from igneous crystallization processes and did not experience significant thermal metamorphism, indicating that the low‐Ca pyroxenes were unlikely incorporated into large parent bodies prior to accretion in their respective comet bodies. An intriguing group of nine low‐Ca pyroxenes from comet Wild 2 with low Cr and Al that fall where type 4–6 chondrites are located are interpreted as products of condensation. The compositional data combined with previously measured oxygen isotopes on 17 low‐Ca pyroxenes support earlier conclusions that comet samples have links with carbonaceous, ordinary, and possibly enstatite chondrite groups. Our results provide additional evidence that comets accreted materials from multiple chondrule reservoirs throughout the solar nebula.

Geochemistry & Geophysics↗

Triple oxygen and hydrogen stable isotope composition of hydroxyl water in Orgueil (CI-type) and Tagish Lake (C2-type) carbonaceous chondrites

The primitive carbonaceous chondrites are of interest to cosmochemical science because they contain relatively large amounts of ‘water’ (H 2 O and/or OH – ) within phyllosilicate minerals. This water is evidence for the accretion of ices by their parent planetesimals, and thus represents an archive of the isotopic compositions of H 2 O in the protoplanetary environments. Here, in this study, we used thermogravimetry-enabled laser spectroscopy (TGA-IRIS) analyses of the Orgueil and Tagish Lake meteorites to make δ 2 H, δ 18 O, and Δ′ 17 O measurements of the H 2 O and OH – contained in the different hydrous minerals that comprise each meteorite. In Orgueil, we measured mass-weighted averages of OH – in the saponite and serpentine phyllosilicate matrix to be δ 2 H = 192 ‰, δ 18 O = 1.5 ‰, which are unquestionably of extraterrestrial origin with Δ′ 17 O = 1.0 ‰. For Tagish Lake, analogous values of OH– in the saponite and serpentine phyllosilicate matrix are δ 2 H = 704 ‰, δ 18 O = 11.3 ‰, and are similarly unambiguously extraterrestrial with Δ′ 17 O = 0.82 ‰. We estimate that the parent H 2 O involved in aqueous alteration of Orgueil had δ 18 O value ≥ +23 ‰. In Orgueil, we interpret the phyllosilicate petrographic relationships, and the δ 18 O values of OH – in saponite and serpentine to indicate that saponite formed first, at a lower temperature by 35 to 53 °C than serpentine. This suggests that the Orgueil parent body experienced increasing temperature during the phase of active aqueous alteration (prograde) which set the δ 18 O OH values of the serpentine and saponite. In the case of Tagish Lake, serpentine formed at a lower temperature by 32 to 60 °C than saponite, for the simplest case with constant δ 18 O H2O values. If serpentine formed first, followed by saponite formation at 32 to 60 °C °C higher temperature, this suggests that Tagish Lake sample TL1 underwent prograde aqueous alteration as the parent body heated up. We find evidence that the parent H 2 O for Orgueil, Tagish Lake sample TL1, and Murchison (based on data from a previous study) may have had Δ′ 17 O values of >0.64 ‰.

Aqueous alteration↗

Activity coefficients of siderophile elements in Fe-Si liquids at high pressure

Metallic core formation in differentiated bodies in the inner solar system can take place between low pressures (near 1 bar) to much higher pressures (up to 100 GPa). Most thermodynamic models of metal-silicate equilibria utilise activity coefficients for metallic tracers in Fe liquids, nearly all of which have been carried out at low pressures. This study focuses on the effect of pressure on activity coefficients for Au, P, V, Mn, Ga, Zn, Cd, Sn, W, Pb, and Nb in liquid Fe-Si alloys. From a series of experiments at 10 GPa, 2373 K containing variable Si content in a metallic liquid we have derived epsilon interaction parameters in Fe-Si liquids ($ε^{Si}_{M}$). Comparison of 1 GPa and 10 GPa data shows no difference except for Nb. Epsilon parameters derived from low pressure experiments can thus be used to calculate activity coefficients for application to higher pressure processes (at least to 10 GPa).

58 GEOSCIENCES↗

The iron oxidation state of Ryugu samples

The Hayabusa2 mission sampled Ryugu, an asteroid that did not suffer extensive thermal metamorphism, and returned rocks to the Earth with no significant air exposure. It therefore offers a unique opportunity to study the redox state of carbonaceous Cb-type asteroids and evaluate the overall redox state of the most primitive rocks of the solar system. An analytical framework was developed to investigate the iron mineralogy and valence state in extraterrestrial material at the micron scale by combining x-ray diffraction, conventional Mössbauer (MS), and nuclear forward scattering (NFS) spectroscopies. An array of standard minerals was analyzed and cross-calibrated between MS and NFS. Then, MS and NFS spectra on three Ryugu grains were collected at the bulk and the micron scales. In Ryugu samples, iron is essentially accommodated in magnetite, clay minerals (serpentine–smectite), and sulfides. Only a single set of Mössbauer parameters was necessary to account for the entire variability observed in MS and NFS spectra, at all spatial scales investigated. These parameters therefore make up a fully consistent iron mineralogical model for the Ryugu samples. As far as MS and NFS spectroscopies are concerned, Ryugu grains are overall similar to each other and share most of their mineralogical features with CI-type chondrites. In detail however, no ferrihydrite is found in Ryugu particles even at the very sensitive scale of Mössbauer spectroscopy. The typical Fe 3+ /Fe tot of clay minerals is much lower than typical redox ratios measured in CI chondrites (Fe 3+ /Fe tot = 85%–90%). Furthermore, magnetite from Ryugu is stoichiometric with no significant maghemite component, whereas up to 12% of maghemite was previously identified in the Orgueil's so-called magnetite. These differences suggest that most CI meteorites suffered terrestrial alteration and that the preterrestrial composition of these carbon-rich samples was less oxidized than previously measured. However, it is not clear yet whether or not the parent bodies of CI chondrites were as reduced as Ryugu. Finally, the high spatial resolution of NFS allows to disentangle the redox state and the crystal chemistry of iron accommodated in serpentine and smectite. The most likely polytype of serpentine is lizardite, containing <35% of Fe 3+ , a fraction of which being tetrahedrally coordinated. Smectite is more oxidized (Fe 3+ /Fe tot > 65%) and mainly contains octahedral ferric iron. In conclusion, this finding implies that these clays formed from highly alkaline fluids and the spatial variability highlighted here may suggest a temporal evolution or a spatial variability of the nature of this fluid.

58 GEOSCIENCES↗