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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Electronic origin of reorganization energy in interfacial electron transfer

Electron transfer (ET) reactions underpin energy conversion and chemical transformations in both biological and a biological systems. The efficiency of any ET process relies on achieving a desired ET rate within an optimal driving force range. Marcus theory provides a microscopic framework for understanding the activation free energy—and therefore the rate—of ET in terms of a key parameter: the reorganization energy. For electrified solid–liquid interfaces, it has long been conventionally understood that only factors in the electrolyte phase are responsible for determining the reorganization energy and that the electronic density of states (DOS) of the electrode only serves to dictate the number of thermally accessible channels for ET. Here we show instead that the electrode DOS plays a central role in governing the reorganization energy, far outweighing its conventionally assumed role. Using atomically layered heterostructures, we tune the DOS of graphene and measure outer-sphere ET kinetics. We find the ensuing variation in ET rate arises from strong modulation in a reorganization energy associated with image potential localization in the electrode. Here we redefine the traditional paradigm of heterogeneous ET kinetics, revealing a deeper role of the electrode electronic structure in interfacial reactivity.

Electrochemistry↗
DOE OSTI2026

Living on the edge: a non-perturbative resolution to the negativity of bulk entropies

Lin, Maldacena, Rozenberg, and Shan (LMRS) presented a new information paradox in black hole physics by noticing that the entanglement and Rényi entropies in a two-sided black hole can become negative when the geometry contains a very large number of matter excitations behind the black hole horizon. While originally this puzzle was presented in the context of BPS two-sided black holes in two-dimensional supergravity, the negativity in fact persists for more general two-sided black holes in the presence of a large number of matter excitations. Since the entanglement and Rényi entropies in ordinary quantum systems cannot be negative, resolving this puzzle is a necessary step towards understanding the quantum mechanical description of black holes. In this paper, we explain how to address the entanglement negativity puzzle, both in the original setting discussed by LMRS and in more general non-supersymmetric settings, by summing over all non-perturbative contributions to the gravitational path integral. We then interpret this result from the point of view of a dual matrix integral, which we use to extend our analysis beyond the regime of validity of the genus re-summation performed in the gravitational path integral. In this regime, positivity is rescued by new saddles of the matrix integral, a one-eigenvalue instanton and a two-eigenvalue instanton. Finally, we formulate a similar puzzle and its resolution using random tensor network techniques.

2D Gravity↗
DOE OSTI2026

Time-dependent-bases with local CUR decomposition method for accelerating turbulent combustion simulations

Here, this study presents a novel reduced-order modeling framework, Time-Dependent Bases with Local CUR decomposition (TDB-L-CUR), designed to efficiently and accurately approximate the species transport equations in reacting flow simulations. The method extends the existing TDB-CUR approach for chemically reacting flows (Jung et al. Comput. Methods Appl. Mech. Engrg. 437 (2025) 117758), which leverages matrix decomposition techniques to form a global-in-space, time-dependent low-dimensional manifold. While TDB-CUR performs well in homogeneous systems, it may be less well-suited to spatially heterogeneous systems such as turbulent flames, where higher-rank approximations are typically required. The proposed TDB-L-CUR framework introduces two methodological extensions to the baseline approach. First, it applies unsupervised clustering to partition the physical domain into distinct regions, enabling spatially localized manifold construction, thereby reducing the rank required for the reduced-order representation. Second, it incorporates a computational singular perturbation (CSP)-based scheme for identifying and penalizing fast species, allowing for spatio-temporally adaptive mitigation of chemical stiffness. The proposed framework is validated on a hierarchy of test cases, including a one-dimensional premixed flame, a two-dimensional nonpremixed ignition case with vortex interaction, and a three-dimensional turbulent premixed flame. TDB-L-CUR significantly improves accuracy over TDB-CUR while further reducing computational cost. The fully on-the-fly formulation of TDB-L-CUR (i.e., requiring no offline training or prior knowledge) makes it a robust and scalable tool for reduced-order modeling of reactive flows.

Local manifold↗
DOE OSTI2026

Enabling probabilistic learning on manifolds through double diffusion maps

Here, we present a generative learning framework for probabilistic sampling that extends Probabilistic Learning on Manifolds (PLoM), which is designed to generate statistically consistent realizations of a random vector in a finite-dimensional Euclidean space, informed by a (representative) set of observations. In its original form, PLoM constructs a reduced-order probabilistic model by combining three main components: (a) kernel density estimation to approximate the underlying probability measure, (b) Diffusion Maps to characterize the manifold of the data, and (c) a reduced-order Itô Stochastic Differential Equation (ISDE) to sample from the learned distribution. However, its sampling dynamics are posed in the ambient space and the retained number of reduced coordinates is chosen by projection-reconstruction error. In practice, this often (i) requires more coordinates than the data’s intrinsic dimension to achieve stable sampling and (ii) lacks a smooth, basis-independent lifting back to the data domain; moreover, standard Diffusion Maps emphasize harmonic eigenfunctions and can miss non-harmonic latent structure. We address these limitations by decoupling geometry learning from sampling: a first Diffusion Maps pass identifies non-harmonic coordinates on which we formulate a full-order ISDE directly in the latent space, while Double Diffusion Maps captures multiscale geometric features and Geometric Harmonics (GH) learns a smooth lifting map to the ambient variables that is independent of the particular diffusion basis. This hybrid design preserves the system’s dynamical richness with a compact geometric representation and enables principled out-of-sample inference. The effectiveness and robustness of the proposed method are illustrated through two numerical studies: one based on data generated from two-dimensional Hermite polynomial functions and another based on high-fidelity simulations of a detonation wave in a reactive flow.

Double diffusion maps↗
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

Pseudocapacitive Titanium Oxynitride Nanowires for Ultrahigh Capacitance Supercapacitors

High-quality, multifunctional two-dimensional (2D) titanium oxynitide (TiNO) thin films and one-dimensional (1D) TiNO nanowires have been synthesized using a pulsed laser deposition, a simple, fast, and congruent evaporation method. First-principles calculations as a function of surface orientation and termination indicate that surface oxidation of TiNO nanowires can stabilize the (110) orientation observed experimentally. The specific capacitance value for the TiNO nanowire samples (2725 mF/cm 2 ) has been found to be nearly six times more than that of the TiNO thin film samples (400 mF/cm 2 ), which is attributed to the high packing density of TiNO nanowires over a given area. The nanowire samples have also been found to exhibit a significantly higher energy density (1.35 μWh/cm 2 ) than the TiNO thin-film samples (0.33 μWh/cm 2 ). Thus, the TiNO material system in thin-film and nanowire forms has been demonstrated to be a promising candidate for use as an electrode material in supercapacitors and other charge-storage applications.

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

Diagonal isometric form for tensor network states in two dimensions

Isometric tensor network states (isoTNS) generalize the isometric form of the one-dimensional matrix product states to tensor networks in two and higher dimensions. Here, we introduce an alternative isometric form for isoTNS by incorporating auxiliary tensors to represent the orthogonality hypersurface. We implement the time evolving block decimation algorithm on this new isometric form and benchmark the method by computing ground states and the real time evolution of the transverse field Ising model in two dimensions on large square lattices of up to 1250 sites. Our results demonstrate that isoTNS can efficiently capture the entanglement structure of two-dimensional area law states. The short-time dynamics is also accurately reproduced even at the critical point. Our isoTNS formulation further allows for a natural extension to different lattice geometries, such as the honeycomb or kagome lattice.

2-dimensional systems↗
DOE OSTI2026

High-performance nanodevices based on WGe 2⁢ N 4 monolayer

Two-dimensional (2D) 𝑀⁢𝐴 2 ⁢𝑍 4 -family monolayers (MLs) have emerged as promising semiconductors due to their element tunability and rich electronic and optoelectronic properties. In this work, using first-principles calculations, we investigate the electronic, mechanical, transport, and optoelectronic properties of WGe 2⁢ N 4 ML with a small indirect bandgap. Various nanodevices based on WGe 2 ⁢N 4 ML are studied, including pn-junction diodes, pin-junction field-effect transistors (FETs), and phototransistors. The present results reveal that the WGe 2⁢ N 4 ML exhibits high rigidity, thermal stability, and remarkable light absorption. These nanodevices demonstrate excellent performance: (1) the pn-junction diode shows a high rectification ratio and a near-Shockley-limit ideality factor, (2) the pin-junction FET exhibits significant gate voltage modulation capability with subthreshold swing as low as 71 mV/dec (close to the theoretical limit calculated based on the Boltzmann distribution), and (3) the phototransistor displays strong optoelectronic responses in the visible and ultraviolet regions. Furthermore, these findings establish WGe 2 ⁢N 4 ML as a versatile platform for developing high-performance, multifunctional nanoelectronic and optoelectronic devices, significantly expanding the application potential of 𝑀⁢𝐴 2⁢ 𝑍 4 -family materials.

2-dimensional systems↗
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↗
DOE OSTI2026

Dynamics of ballistic photocurrents driven by Coulomb scattering in a two-dimensional material

First-principles real-time time-dependent density-functional theory (rt-TDDFT) calculations reveal the existence of ballistic photocurrents generated by Coulomb scattering, which has not previously been considered as a mechanism for the bulk photovoltaic effect. With monolayer GeS as an example, it is predicted that ballistic currents can be comparable to shift currents under experimentally accessible conditions.

2-dimensional systems↗
DOE OSTI2026

Site-decorated model for unconventional frustrated magnets: Ultranarrow phase crossover and two-dimensional spin reversal transition

Here, the site-decorated Ising model is introduced to advance the understanding and experimental realization of the recently discovered one-dimensional (1D) finite-temperature ultranarrow phase crossover in an external magnetic field, while mitigating the geometric complexities of traditional bond-decorated models. The unconventional frustration and physics are clarified by exactly mapping the 1D site-decorated Ising model in a magnetic field onto a zero-field bond-decorated 𝐽 1 −𝐽 2 Ising model with conventional geometrical frustration. Furthermore, although higher-dimensional Ising models in an external field remain unsolved exactly, an exact solution for a spin-reversal transition—driven by an exotic, hidden half-ice, half-fire state induced by site decoration—is derived. This transition, triggered by a slight variation in temperature or magnetic field—without changing its direction—even in the weak-field limit, offers a promising route toward energy-efficient applications such as data storage and processing. The results suggest that site decoration offers an avenue for materials and device design, particularly in systems such as mixed 𝑑−𝑓 compounds, optical lattices, and neural networks, calling for further studies with site-decorated Heisenberg models. In addition, the site-decorated model offers a rigorous test ground for artificial intelligence (AI) in science, as the analytic derivation of the present results was not only validated but also improved by a general-purpose large language model, inspiring the use of AI as scientific discoverer.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

2D reactive transport model of shale chemical weathering and biogeochemical fluxes along a mountainous hillslope, East River Watershed, Colorado: Input files and simulation results

This data package contains input files and simulation results for a two-dimensional (2D) reactive transport model used to quantitatively analyze the coupled hydrological and biogeochemical processes governing shale weathering and associated biogeochemical fluxes under realistic environmental conditions in the high-elevation East River Watershed. These data support the conclusions presented in Stolze et al. (Water Resources Research, under review), "Model-based interpretation of solute exports and carbon partitioning during shale weathering in a mountainous hillslope". The model simulates atmospheric-subsurface gas exchange, subsurface water flow, and shale weathering processes under dynamic, year-scale conditions along a shale-underlain hillslope located in the East River watershed. The simulations were performed using the PFLOTRAN flow and reactive transport code and executed on the Perlmutter supercomputer to leverage its large-scale parallel computing capabilities. The data package contains two zipped folders, "model_input_files" and "simulation_results", and one readme.txt file. "model_input_files" contains the necessary input files to run the calibrated base-base model presented in Stolze et al. (Water Resources Research, under review). "simulation_results" contains a single hdf5 file ("Output_2D_hillslope_model.h5") which includes the results of simulation performed using the base-case model. This file can be opened with HDFView 3.1.4, Python, or MATLAB. "readme.txt" contains relevant information about the base-case model and provides guidelines on how to run the associated input files provided in the folder "model_input_files". Furthermore, readme.txt provides information regarding the model results provided in "Output_2D_hillslope_model.h5" such as matrix dimensionality and output units. Field datasets used to evaluate model performance were collected at three monitoring wells located along a hillslope transect (PLM1, PLM2, and PLM3). Dissolved ion concentration data were collected from November 2016 to October 2021 for Ca, Mg, DIC, Na, K, SO4 (Dong et al., 2025 - dic_npoc_data_2014_2024.zip - DOI:10.15485/1660459; Williams et al., 2025 - anion_data_2014_2024.zip - DOI:10.15485/1668054; Dong et al., 2025 - cation_data_2014_2024.zip - DOI:10.15485/1668055). Note that we used the files named er_PLM1_xx_yy, er_PLM2_xx_yy, and er_PLM3_xx_yy where xx stands for the name of the aqueous species and yy stands for the depth where the measurements were performed. Soil water content ([0 - 1] m) and water table depth were collected from November 2016 to October 2021 (Wan et al., 2024 - Dynamic_water_table__depthsFig2b.csv and Soil_water_content_Fig4e.csv - DOI:10.15485/2322567). Gaseous CO2 concentration were collected from October 2020 to December 2021(Wan et al., 2024 - Soil_CO2_concentrations_Fig4h.csv - DOI:10.15485/2322567) Gaseous CO2 flux from the subsurface to the atmosphere were collected in the vicinity of PLM2 from October 2019 to May 2022 (Wu et al., 2025). Soil microbial biomass concentration was measured from August 2016 to June 2017 (Sorensen et al., 2019 - 2017_East_River_Pumphouse_Microbial_Biomass__1_.csv - DOI:10.15485/1577267) All field data are published as CSV files compatible with Microsoft Excel, MATLAB, and Python, or as text files. The coordinates of the monitoring wells and the CO2(g) flux sensor in the coordinate system WGS84 are: -PLM1: [38.9197710 ; -106.9492750] -PLM2: [38.9201580 ; -106.9487170] -PLM3: [38.9207843 ; -106.9483668] -PLM4: 38.9210060 ; -106.9479528] -CO2(g) flux sensor: [38.9199180 ; -106.9489906] ------------------------------------------------------------------------------------------- This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231. This research used resources of the National Energy Research Scientific Computing Center (NERSC), a Department of Energy User Facility using NERSC award BER-ERCAP 23980, BER-ERCAP 28550, and BER-ERCAP 33789.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Casimir-Polder potential on an excited atom near an atomic array

We develop a microscopic description of the fluctuation-mediated Casimir-Polder (CP) shifts on a 'test' two-level atom placed near a two-dimensional atomic array of two-level atoms. We derive the resonant and off-resonant CP potentials experienced by the excited test atom using fourth-order perturbation theory, under the assumption that the test atom resonance is far detuned from those of the array atoms. The total potential on the test atom can be described as the sum of the pairwise resonant and off-resonant potentials resulting from its interaction with the individual atoms of the array. We analyze the asymptotic scaling of CP shifts as a function of the test atom-array separation, and its dependence on various system parameters: array spacing and size, and dipole orientation of the array atoms. Our results bridge the description of CP potential across two distinct regimes: (i) from a single-atom limit where we recover the well-known two-atom Van der Waals potential, (ii) to a macroscopic boundary limit, where we demonstrate new asymptotic scaling laws. We demonstrate that these scaling laws can be tuned via the microscopic parameters of the atomic array, establishing atomically-controlled arrays as a versatile platform for tailoring fluctuation-induced QED phenomena.

FOS: Physical sciences↗
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↗