Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “TWO-DIMENSIONAL SYSTEMS”

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

Geometric Delocalization in Two Dimensions

We demonstrate the existence of transient two-dimensional surfaces where a random-walking particle escapes to infinity in contrast to localization in standard flat two-dimensional space. We first prove that any rotationally symmetric two-dimensional membrane embedded in flat three-dimensional space cannot be transient. Then we formulate a criterion for the transience of a general asymmetric two-dimensional membrane. We use it to explicitly construct a class of transient two-dimensional manifolds with a nontrivial metric and height function but “zero average curvature,” which we dub “tablecloth manifolds.” The absence of the logarithmic infrared divergence of the Laplace-Beltrami operator in turn implies the absence of weak localization, nonexistence of bound states in shallow potentials, and breakdown of the Mermin-Wagner theorem and Kosterlitz-Thouless transition on the tablecloth manifolds, which may be realizable in both quantum simulators and corrugated two-dimensional materials.

Anderson localization↗
DOE OSTI2026

Integration of Online Cross-Section Generation Capability with Depletion and Transient Solvers in Griffin

Griffin is a Multiphysics Object-Oriented Simulation Environment (MOOSE)-based reactor multiphysics analysis application jointly developed by Argonne and Idaho National Laboratories under the DOENE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. In FY25, an online crosssection generation capability based on the Self-Shielding Application Programming Interface (SSAPI) was demonstrated for TRISO-fueled reactor problems under steady-state conditions. This fiscal year, that capability was extended to support depletion and transient multiphysics calculations, enabling high-fidelity analyses that generate self-shielded cross sections on the fly from the actual evolving composition and temperature states rather than from pre-tabulated libraries. For depletion, a two-way coupling was established in which SSAPI computes compact-averaged self-shielded cross sections that the depletion solver then uses to advance the Bateman equations, with the updated compositions returned to SSAPI at each step; the depletion module was refactored to support both library-based and SSAPI-based cross sections, and additional logic was added to track daughter isotopes and to exclude minor isotopes for efficiency. For transient analysis, the SSAPI multigroup library was extended with the kinetics data required for time-dependent calculations, the Improved Quasi-Static (IQS) scheme was coupled with SSAPI, and several supporting capabilities were implemented, including a self-shielding treatment that lets control rods and drums move within a self-shielded model, which had previously been impossible and had ruled out rod- and drum-movement transients with on-the-fly cross sections altogether, a new mixing scheme for delayed-neutron precursor decay constants, a checkpoint-based restart workflow, and performance improvements such as pointwise cross-section interpolation and the bypassing of unnecessary Dancoff factor calculations. The implemented capabilities were verified against Serpent Monte Carlo solutions. For depletion, a prismatic pin-cell problem based on a Next Generation Nuclear Plant (NGNP) Very High Temperature Reactor benchmark showed excellent agreement, with eigenvalue differences within 200 pcm over the entire burnup range (up to 140 MWD/kgU) and fission-product and actinide inventories agreeing to within 0.8% and 2.5%, respectively; a heat-pipe microreactor assembly problem with a much higher fuel loading confirmed the same behavior and quantified the bias introduced when the multigroup equivalence effect is neglected. For transient analysis, a pin-cell problem with a step reactivity insertion and temperature feedback reproduced the analytically expected asymptotic power and showed close agreement between the direct and IQS solutions, and a two-dimensional microreactor core problem with control-drum rotation exercised the new moving-drum self-shielding treatment and demonstrated successful coupling of the online crosssection generation with both the direct and IQS transient methods. The capability was further exercised on a full-core pebble-bed problem, in which Griffin was coupled with the System Analysis Module (SAM) to simulate load-following operation of the gPBR with the Doppler feedback resolved at the TRISO fuel kernel temperature. These developments in Griffin provide a convenient, high-fidelity approach to cross-section generation for advanced thermal reactors with geometrically complex and highly heterogeneous configurations, including TRISO-fueled prismatic and pebble-bed systems, and support steady-state, depletion, and transient multiphysics calculations. They also enable self-shielded cross sections to be evaluated directly at the actual coupled state of the system, thereby establishing a foundation for high-fidelity, fully coupled multiphysics analysis of advanced reactors

Park, H.↗
NASA NTRSDate not supplied

Toward Large Field-of-view X-ray Spectrometers with Magnetic Microcalorimeters

A future X-ray telescope mission like Lynx, a X-ray flagship mission concept, requires a large field of view, high energy resolution, and high spatial resolution. The Magnetic Microcalorimeter (MMC) detectors we are developing for such a mission feature a 5-arcminute field of view, approximately 100,000 pixels, each spanning 0.5–1 arcsecond. These detectors require the readout of 8,000 sensors. Three different array types have been designed to meet scientific goals, two of which utilize thermal multiplexing to reduce the number of readout channels in a scheme called “Hydra.” We are currently fabricating detectors on a recently completed buried wiring wafer with connections to every sensor. The sensors are read out using a microwave multiplexer (μMUX) with flux ramp modulation. A single μMUX can read out approximately 1,000 sensors with only two feed lines. We will report on the test of the first two-dimensional μMUX chip, consisting of 78 resonators. μMUX chips that pass screening will be directly bonded to MMC array chips using indium bump bonding. In parallel, we have tested a kinetic inductance traveling wave parametric amplifier (KI-TWPA) to amplify signals at low temperatures. The KI-TWPA exhibits extremely low internal noise, approaching the standard quantum limit. Since the energy resolution of an MMC is primarily limited by the readout system, the KI-TWPA is an ideal amplifier for MMCs with μMUX readout. We recently measured the readout noise of a μMUX with a KI-TWPA, demonstrating its capability for reading out MMC arrays. We present the design and fabrication of the latest MMC detectors, along with its readout performance.

Magntic Microcalorimeters↗
DOE OSTI2026

Magnetism and magnetoelastic effect in the two-dimensional van der Waals multiferroic CuCrP 2 ⁢S 6

Here, we report a magnetic and neutron diffraction study on the ground state magnetism and field evolution of single-crystal van der Waals multiferroic CuCrP 2 ⁢S 6 . The ordered moments align along the 𝑏 axis in the "𝐴-type" antiferromagnetic configuration with a spin-flop transition along the same direction. Field application along 𝑎 introduces a smooth transition to a fully polarized ferromagnetic state via in-plane spin rotation. These findings resolve the ambiguity of the ground state magnetization direction in CuCrP 2 ⁢S 6 and uncover its field responses, providing a firm basis for future magnetoelectric study. A magnetoelastic coupling effect connecting the interlayer spacing and the magnetic order was further revealed, highlighting the out-of-plane strain as an effective control knob for tuning magnetism both in this system and in related van der Waals magnets.

Guo, Jiasen [Oak Ridge National Laboratory (ORNL),↗
DOE OSTI2026

Avoidance of disruptions on KSTAR due to vertical displacement events via novel real-time stability assessment

Disruption avoidance via the DECAF approach has been achieved on KSTAR using a novel real-time vertical stability assessment and a multiactuator feedback control strategy. The development of disruption avoidance strategies with reactor-relevant reliability is an urgent activity, enabling future fusion power plants. The stability metric employed is based on a new formulation of a vertical force gradient balance metric evaluated across the poloidal cross section of the plasma, with parameters tuned using historical data. Evaluation of this metric on a validation set of 400 recent KSTAR shots indicates >82% of Vertical displacement events can be avoided via feedback control. Essential to its calculation is the two-dimensional toroidal current density distribution in the plasma. Measurement of this profile faster than fully-converged equilibrium reconstructions can deliver is found to improve forecaster performance and is achieved with a surrogate model that takes as input magnetic diagnostic measurements and outputs the current profile on a basis comprising the top principal components of historical current profiles (from past equilibrium reconstructions). This method solves the non-uniqueness problem typically faced when reconstructing current profiles directly from diagnostics, while improving computational time and accuracy. On average, profiles produced by this model reach coefficients of determination of >0.99 with respect to those from equilibrium reconstructions. The avoidance actuators employed include poloidal field coils and an electron cyclotron current drive system. The multiactuator approach, as shown in this first demonstration, allows disruption avoidance while minimizing impact to operational performance. This ability, along with its flexibility and speed, makes this new approach an attractive option for avoiding these types of disruptions in reactors.

Tobin, Matthew [Columbia Univ., New York, NY (Unit↗
DOE OSTI2026

Synthesis and Characterization of Layered Actinide (U, Np, Pu) Oxide and Hydroxide Phases

Systematic structural comparisons across the early actinides remain limited by the scarcity of well-defined transuranic layered oxide and oxyhydroxide phases. Here, we report the synthesis and single-crystal characterization of new layered actinide compounds spanning U, Np, and Pu obtained under mild hydrothermal conditions in concentrated alkali hydroxide media. These include hydrated oxides α-Cs 2 U 2 O 7 •0.5H 2 O and Rb 2 An 2 O 7 •0.5H 2 O (An = Np, Pu), oxy-hydroxides Rb 6 [(AnO 2 ) 6 O 8 (OH) 2 ]•xH 2 O x = 0, 0.5 (An = U, Np) and Rb 4 [(UO 2 ) 5 O 6 (OH) 2 ]•2H 2 O, as well as nitrate-intercalated compounds Cs 6 [(AnO 2 ) 3 O 4 (OH) 2 ](NO 3 ) 2 (An = Np, Pu). Single-crystal X-ray diffraction studies reveal extended two-dimensional architectures constructed from edge- and vertex-sharing actinyl polyhedra, with systematic evolution in equatorial coordination, hydration, and anionic sheet topology across the U–Np–Pu series. Incorporation of nitrate anions within the interlayer region of Cs 6 [(AnO 2 ) 3 O 4 (OH) 2 ](NO 3 ) 2 establishes a previously unobserved structural motif in layered transuranic oxyhydroxides, demonstrating an additional pathway for anion-mediated framework stabilization. Correlation of crystallographic metrics with single-crystal Raman spectroscopy provides new vibrational benchmarks linking differences in An═O yl bond lengths to the equatorial coordination and the interstitial cations. These findings expand the structural hierarchy of layered actinide materials and address clarifying periodic trends governing topology, bonding, and vibrational signatures in high-valent 5ƒ oxide systems.

actinides↗
DOE OSTI2026

Molecularly engineered ZnO–carbon nanosheets from fumaric acid precursors for efficient photocatalytic water purification

The photocatalytic breakdown of organic contaminants is crucial for the development of water purification technology. Zinc oxide (ZnO) is an extensively researched photocatalyst; however, its efficacy is hindered by rapid charge recombination and limited utilization of UV irradiation. Resolving these issues necessitates integrating ZnO with conductive carbon phases via scalable, low-temperature synthesis. We provide a molecularly designed sol–gel method that converts zinc–fumarate coordination networks into two-dimensional ZnO–carbon nanosheets utilizing solely aqueous precursors and mild annealing temperatures (400–600 °C). This method utilizes fumaric acid as a dicarboxylate linker and polyvinylpyrrolidone (PVP) as a structural carbon source to produce ultrathin wurtzite ZnO nanosheets embedded inside an amorphous carbon matrix. The resultant ZnO-C hybrid achieves nearly complete methylene blue degradation within 10 min under UV-A illumination, demonstrating first-order kinetics and outstanding recyclability. Compared with commercial ZnO, the ZnO–C nanosheets exhibit comparable rapid photocatalytic degradation, enhanced adsorption behavior, a porous nanosheet morphology, and an integrated ZnO–carbon interfacial structure. These findings provide a viable molecular-templating approach to fabricating various metal oxide–carbon photocatalysts and underscore substantial enhancements in semiconductor efficacy in eco-friendly water treatment systems.

Ozcan, Muca [ORNL] (ORCID:0000000320020474)↗
DOE OSTI2026

Quantum Filtering and Analysis of Multiplicities in Eigenvalue Spectra

Fine-grained spectral properties of quantum Hamiltonians, including both eigenvalues and their multiplicities, provide useful information for characterizing many-body quantum systems as well as for understanding phenomena such as topological order. Extracting such information with small additive error is #BQP-complete in the worst case. In this work, we introduce QFAMES (quantum filtering and analysis of multiplicities in eigenvalue spectra), a quantum algorithm that efficiently identifies clusters of closely spaced dominant eigenvalues and determines their multiplicities under physically motivated assumptions, which allows us to bypass worst-case complexity barriers. QFAMES also enables the estimation of observable expectation values within targeted energy clusters, providing a powerful tool for studying quantum phase transitions and other physical properties. We validate the effectiveness of QFAMES through numerical demonstrations, including its applications to characterizing quantum phases in the transverse-field Ising model and estimating the ground-state degeneracy of a topologically ordered phase in the two-dimensional toric code model. We also generalize QFAMES to the setting of mixed initial states. Our approach offers rigorous theoretical guarantees and significant advantages over existing subspace-based quantum spectral analysis methods, particularly in terms of the sample complexity and the ability to resolve degeneracies.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

i -Wave Symmetry Altermagnetism and Anomalous Hall Conductivity in Monolayer FeCl 3

In this work, the FeCl 3 monolayer is proposed to host a rare and exotic i -wave symmetry altermagnetism. Using first-principles calculations, we performed a detailed investigation of the structural stability and electronic, magnetic, and spintronic properties of monolayer FeCl 3 . Our results suggest the system is semiconducting in nature and exhibits i-wave symmetry nonrelativistic spin splitting (NRSS), originating from its hexagonal crystal structure with roto-inversion symmetries and anisotropic crystal field. The momentum-resolved spin-splitting further supports the presence of spin-dependent topological features. Further, the monolayer exhibits sizable conventional spin-Hall and anomalous Hall conductivities, highlighting its promise as a versatile platform for next-generation flexo-spintronic applications.

altermagnetism↗
DOE OSTI2026

Creation of depth-confined, shallow nitrogen-vacancy centers in diamond with tunable density

Engineering shallow nitrogen-vacancy (NV) centers in diamond holds the key to unlocking new advances in nanoscale quantum sensing. We find that the creation of near-surface NVs through delta doping during diamond growth allows for tunable control over both NV depth confinement (with a twofold improvement relative to low-energy ion implantation) and NV density, ultimately resulting in highly sensitive single defects and ensembles with coherence limited by NV–NV interactions. Additionally, we demonstrate the utility of our shallow delta-doped NVs by imaging magnetism in few-layer CrSBr, a two-dimensional magnet. Here, we anticipate that the control afforded by near-surface delta doping will enable new developments in NV quantum sensing from nanoscale nuclear magnetic resonance to entanglement-enhanced metrology.

Chemical vapor deposition↗
DOE OSTI2024

Analytical Model for Atomic Relaxation in Twisted Moiré Materials

By virtue of being atomically thin, the electronic properties of heterostructures built from two-dimensional materials are strongly influenced by atomic relaxation. The atomic layers behave as flexible membranes rather than rigid crystals. Here we develop an analytical theory of lattice relaxation in twisted moiré materials. We obtain analytical results for the lattice displacements and corresponding pseudo gauge fields, as a function of twist angle. We benchmark our results for twisted bilayer graphene and twisted WSe 2 bilayers using large-scale molecular dynamics simulations. Our single-parameter theory is valid in graphene bilayers for twist angles 𝜃 ≳ 0.7°, and in twisted WSe 2 for 𝜃 ≳ 1.6°. Furthermore, we also investigate how relaxation alters the electronic structure in twisted bilayer graphene, providing a simple extension to the continuum model to account for lattice relaxation.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Ion Transport in Charged Membranes: Linking Electric-Field-Driven Mechanisms to Pore Size via Perturbation Analysis

Ion-exchange membranes are a critical component in electrochemical systems. Nevertheless, the understanding and modeling of ion transport within these porous structures have been limited by particular complexity reductions, either ignoring the dimensionality of their porous network architecture or imposing geometric assumptions (i.e., overlapping double layers). Before addressing this morphology-transport gap, a framework that relates the driving forces of transport to the geometry of a single pore is required. In this work, our modeling domain consists of a two-dimensional single pore with charged walls, connecting two identical electrolyte reservoirs. Using the Poisson-Nernst-Planck equations and regular perturbation theory, we decouple the electric fields and analyze the driving forces of ion transport, specifically electromigration and induced electroosmosis within the pore. These processes are described as analytical functions of the interaction aspect ratio, ?, defined as the ratio of the pore radius to the Debye length. Using this parameter, our study (i) describes the interplay between electromigrative and electroosmotic mechanisms that set ionic conductivity, (ii) identifies a dimensionless group of intrinsic electrolyte properties that indicates the predominant driving force, and (iii) provides a qualitative, confinement-dependent perspective on selectivity in ion-conducting membranes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Imaging a terahertz superfluid plasmon in a two-dimensional superconductor

The superconducting gap defines the fundamental energy scale for the emergence of dissipationless transport and collective phenomena in a superconductor. In layered high-temperature cuprate superconductors, in which the Cooper pairs are confined to weakly coupled two-dimensional (2D) copper–oxygen (CuO 2 ) planes, terahertz (THz) spectroscopy at subgap millielectronvolt (meV) energies has provided crucial insights into the collective superfluid response perpendicular to the superconducting layers. However, within the CuO 2 planes, the collective superfluid response manifests as plasmonic charge oscillations at energies far exceeding the superconducting gap, obscured by strong dissipation. Here, in this study, we present spectroscopic evidence of a below-gap, 2D superfluid plasmon in few-layer Bi 2 Sr 2 CaCu 2 O 8+x and spatially resolve its deeply subdiffractive THz electrodynamics. By placing the superconductor in the near field of a spintronic THz emitter, we reveal this distinct resonance—absent in bulk samples and observed only in the superconducting phase—and determine its plasmonic nature by mapping the geometric anisotropy and dispersion. Crucially, these measurements offer a direct view of the momentum-dependent and frequency-dependent superconducting transition in two dimensions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Real-time observation of toroidal current redistributions induced by three-dimensional MHD phenomena triggering vertical displacement events in tokamak plasmas

Three-dimensional MHD instabilities, including edge-localized modes (ELMs) and internal reconnection events (IREs), have been observed to precipitate loss of vertical stability in tokamak plasmas, resulting in vertical displacement events (VDEs). This vertical destabilization can occur due to toroidal current redistributions and/or shape changes resulting from these phenomena. Using a recently introduced method for rapidly reconstructing the two-dimensional toroidal plasma current density profile in real-time, results are presented that demonstrate the specific current distribution changes that occur during ELMs (on KSTAR) and IREs (on MAST-U) that lead to loss of vertical control. The method most efficiently reconstructs the toroidal current density profile by doing so on a basis of principal components of historical profiles. These principal components isolate dominant current profile dynamics, improving interpretability, increasing speed, and reducing dimensionality of the profile computation. On KSTAR, this computation is executed in the real-time plasma control system at a rate of 10 kHz (limited by available CPU cycle times), allowing the current profile evolution to be assessed at several times over the course of each ELM event. Further, by incorporating the reconstructions into a novel vertical stability metric, the contribution of specific current profile dynamics to the loss of vertical stability can be assessed in real-time for VDE avoidance and improved understanding of the causal relationship between three-dimensional MHD phenomena and VDEs. The success of this method in approximating toroidal current density profiles from kinetic equilibrium reconstructions is also presented ($R^2=0.990$), along with its capability to produce other equilibrium quantities of interest in real-time at high time resolution.

edge-localized modes↗
DOE OSTI2026

MATEY: multiscale adaptive transformer models for spatiotemporal physical systems

Accurate representation of the multiscale features in spatiotemporal physical systems using vision transformer architectures requires extremely long, computationally prohibitive token sequences. To address this issue, we propose two novel adaptive tokenization schemes that dynamically adjust patch sizes based on local features: one ensures convergent behavior to uniform patch refinement, while the other offers better computational efficiency. Moreover, we present a set of spatiotemporal attention schemes, where the temporal or axial spatial dimensions are decoupled, to evaluate their baseline computational and data efficiencies and to determine whether adaptive tokenization can improve this performance. We assess the performance of the proposed multiscale adaptive model, MATEY, in a sequence of experiments. Compared to a full spatiotemporal attention scheme or a scheme that decouples only the temporal dimension, we find that fully decoupled axial attention is less efficient and expressive, requiring more training time and model parameters to achieve the same accuracy. The experiments on the adaptive tokenization schemes show that, compared to a uniformly refined model, the proposed schemes achieve comparable or improved accuracy at a much lower cost in the tested two-dimensional settings. While the asymptotic analysis suggests the potential for favorable scaling, empirical validation at substantially longer sequence lengths remains to be performed in future work. Finally, we demonstrate in two fine-tuning tasks featuring different physics that models pretrained on PDEBench data outperform the ones trained from scratch, especially in the low data regime with frozen attention.

adaptive tokenization↗
DOE OSTI2026

Josephson dynamics in two-dimensional ring-shaped condensates

We investigate Josephson transport in a fully closed, two-dimensional superfluid circuit formed by a ring-shaped 87 Rb Bose-Einstein condensate that contains two optical barriers acting as movable weak links. Translating these barriers at controlled speeds imposes a steady bias current, enabling direct mapping of the current-chemical-potential (𝐼−Δ⁢𝜇) characteristics. For narrow junctions (𝑤 ≈ 1µ⁢m), the circuit exhibits a pronounced dc branch that terminates at a critical current 𝐼 𝑐 = 9⁢(1) × 10 3 s −1 ; above this threshold, the system switches to an ac, resistive regime. Classical-field simulations that include the moving barriers quantitatively reproduce both the nonlinear 𝐼−Δ⁢𝜇 curve and the measured 𝐼 𝑐 , validating the underlying microscopic picture. Analysis of the ensuing phase dynamics shows that dissipation is mediated by the nucleation and traversal of vortex-antivortex pairs through the junctions, while the bulk condensate remains globally phase locked—direct evidence of the ring's topological constraint enforcing quantized circulation. These results establish a cold-atom analog of a superconducting quantum interference device in which Josephson dynamics can be resolved at the single-vortex level, providing a versatile platform for atomtronic circuit elements, nonreciprocal Josephson devices, and on-chip Sagnac interferometers for multiaxis rotation sensing.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

JIMWLK on a quantum computer

We propose a method for solving the Jalilian-Marian-Iancu-McLerran-Weigert-Leonidov-Kovner (JIMWLK) evolution equation on quantum computers. Our approach exploits the reformulation of the JIMWLK equation as a Lindblad master equation governing the rapidity evolution of the hadronic density matrix, as established in prior work. To render the problem tractable for quantum simulation, we introduce several approximations: the two-dimensional transverse plane is reduced to a one-dimensional radial lattice by assuming azimuthal symmetry of the jump operators; the gauge group is restricted to SU(2); and the infinite Wilson lines of the JIMWLK equation are replaced by finite Wilson links along the light-cone direction. The resulting bosonic Hilbert space is truncated using the electric field basis familiar from Hamiltonian lattice gauge theory, with states restricted to angular momenta 𝑗 ≤ 𝑗 max . We derive the matrix elements of the JIMWLK Lindblad jump operators in this basis. As a benchmark, we demonstrate rapid convergence of the fundamental dipole expectation value with 𝑗 max for both pure and mixed Gaussian initial density matrices. For the simplest truncation, 𝑗 max =1/2, we implement the Lindblad evolution using a quantum simulation algorithm verified with the Qiskit statevector simulator by decomposing the non-unitary evolution operator into a linear combination of unitaries. This work establishes a concrete pathway toward quantum simulation of high-energy quantum chromodynamics evolution equations, with direct relevance to the physics program of the Electron-Ion Collider.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2026

Molecular Inventory and Comparative Organic Profiling of A Homogenized Aggregate Sample From Asteroid (101955) Bennu

Evaluating the distribution and diversity of solvent-soluble organic molecules in extraterrestrial materials provide information about the physicochemical environments where prebiotic chemistry occurred and is essential for assessing the inventory of prebiotic compounds available to the early Earth. In this work, we conducted a comprehensive molecular characterization of solvent-soluble organic compounds in a homogenized aggregate sample (ID: OREX-800107-128) of unsorted regolith from asteroid (101955) Bennu returned by NASA’s OSIRIS-REx mission. Using two-dimensional gas chromatography coupled with high-resolution time-of-flight mass spectrometry (GC×GC-HRMS), we identified 87 organic compounds across several molecular classes, including alkanes, polycyclic aromatic hydrocarbons, and nitrogen-, oxygen-, and sulfur-bearing species. Comparison with the organic inventory of CI (Orgueil), C2-ungrouped (Tagish Lake and Tarda), and CM2 (Murchison) carbonaceous chondrites revealed that Bennu sample OREX-800107-128 is most similar to Tagish Lake and Tarda yet compositionally distinct from all of these meteorites. This Bennu sample exhibits a relatively high abundance of low-molecular-weight aromatics (e.g., toluene), pyridine derivatives, and sulfur species such as dimethyl sulfite, consistent with episodic low-temperature aqueous alteration. The presence of these organic compounds suggests that this material experienced chemically diverse processes in the parent body conducive to complex organic synthesis, including aldol condensation reactions responsible for N- and O-bearing compounds. We also compared the molecular distributions of this Bennu sample and previously analyzed samples from (162173) Ryugu to assess potential common origins between these two carbonaceous asteroids. Bennu shows more chemically diverse and nitrogen-rich organic materials than Ryugu. While both asteroids share broadly similar primitive compositions, the observed chemical differences indicate distinct evolutionary histories within a likely common formation environment though whether Bennu and Ryugu derive from a single compositionally heterogeneous parent body or from multiple closely related progenitors remains unresolved. Further detailed comparative analyses of material from Bennu, Ryugu, and specific carbonaceous meteorites (e.g., Tagish Lake and Tarda) will provide new insights into how the N-rich molecular profile observed in Bennu samples was formed and evolved in the outer Solar System.

sample return↗