Search NASA⌕ Search

SEARCH · Search NASA

Results for “molecular decomposition”

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 OSTI2020

Realization of an Asymmetric Non‐Aqueous Redox Flow Battery through Molecular Design to Minimize Active Species Crossover and Decomposition

Abstract This communication presents a mechanism‐based approach to identify organic electrolytes for non‐aqueous redox flow batteries (RFBs). Symmetrical flow cell cycling of a pyridinium anolyte and a cyclopropenium catholyte resulted in extensive capacity fade due to competing decomposition of the pyridinium species. Characterization of this decomposition pathway enabled the rational design of next‐generation anolyte/catholyte pairs with dramatically enhanced cycling performance. Three factors were identified as critical for slowing capacity fade: (1) separating the anolyte–catholyte in an asymmetric flow cell using an anion exchange membrane (AEM); (2) moving from monomeric to oligomeric electrolytes to limit crossover through the AEM; and (3) removing the basic carbonyl moiety from the anolyte to slow the protonation‐induced decomposition pathway. Ultimately, these modifications led to a novel anolyte–catholyte pair that can be cycled in an AEM‐separated asymmetric RFB for 96 h with >95 % capacity retention at an open circuit voltage of 1.57 V.

Shrestha, Anuska↗
DOE OSTI2023

The Energetic Origins of Pi–Pi Contacts in Proteins

Accurate potential energy models of proteins must describe the many different types of noncovalent interactions that contribute to a protein’s stability and structure. Pi–pi contacts are ubiquitous structural motifs in all proteins, occurring between aromatic and nonaromatic residues and play a nontrivial role in protein folding and in the formation of biomolecular condensates. Guided by a geometric criterion for isolating pi–pi contacts from classical molecular dynamics simulations of proteins, we use quantum mechanical energy decomposition analysis to determine the molecular interactions that stabilize different pi–pi contact motifs. We find that neutral pi–pi interactions in proteins are dominated by Pauli repulsion and London dispersion rather than repulsive quadrupole electrostatics, which is central to the textbook Hunter–Sanders model. This results in a notable lack of variability in the interaction profiles of neutral pi–pi contacts even with extreme changes in the dielectric medium, explaining the prevalence of pi-stacked arrangements in and between proteins. We also find interactions involving pi-containing anions and cations to be extremely malleable, interacting like neutral pi–pi contacts in polar media and like typical ion–pi interactions in nonpolar environments. Like-charged pairs such as arginine–arginine contacts are particularly sensitive to the polarity of their immediate surroundings and exhibit canonical pi–pi stacking behavior only if the interaction is mediated by environmental effects, such as aqueous solvation.

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

Decomposition characteristics of C4F7N-based SF6-alternative gas mixtures

C4F7N [2,3,3,3-tetrafluoro-2-(trifluoromethyl)propanenitrile]/CO2 gas mixtures are being developed as an eco-friendly electrical insulator to replace SF6, the most potent greenhouse industrial gaseous dielectric. However, recent studies have reported complicated and often conflicting decomposition pathways for C4F7N/CO2 gas mixtures, which has raised concerns. In this work, the decomposition characteristics of C4F7N/CO2 gas mixtures were studied comprehensively by both designed computations and experiments. Computations were performed starting from fundamental propositions of C4F7N/CO2 decompositions, which were further experimentally verified by pyrolysis, long-term thermal aging with/without catalytic materials (industrial-grade molecular sieves 4A), and electrical decomposition by spark discharge. The results of both computations and experiments suggest that in an ideal thermal decomposition, C4F7N is likely to decompose into C2F6 and small fluoronitriles first at high temperatures. The generation of C3F6 and C2N2 from C4F7N thermal decomposition at lower temperatures appears because of the catalytic effect of incompatible materials, for example, the industrial-grade molecular sieves 4A that we tested. The electron impact dissociation of C4F7N plays an important role in C4F7N electrical decomposition, leading to additional formation of distinctive small molecules of CF4 and C2N2 of low concentrations. It was pointed out based on a real arcing test in a load disconnector that the decomposition of C4F7N gas mixtures in real applications will be at a much moderate and manageable rate than what was obtained from the highly accelerated laboratory tests presented in this work. The signatures of decomposition products extracted in this study provide invaluable guidance for developing decomposition-based diagnosis and fixation of decomposition byproducts toward SF6-free power grids.

Physics↗
DOE OSTI2025

Temperature-based reactive flow model for triaminotrinitrobenzene (TATB) plastic bonded explosives

A new reactive flow model is presented for triaminotrinitrobenzene (TATB)-based plastic bonded explosives, applicable to shock initiation and steady detonation problems of differing initial temperature. Temperature disequilibrium is assumed between unreacted explosive, material in the vicinity of compressed defects (called hot spots), and reaction products. The model incorporates temperature-dependent decomposition reaction rates. Particularly, Arrhenius model parameters were derived from quantum-based molecular dynamics simulations of TATB decomposition. Further, a model of detonation carbon aggregation is incorporated, describing the slow release of energy inherent to detonation in TATB-based materials. Model parameters were calibrated against gas gun shock initiation experiments and steady detonation rate stick tests. The predictive ability of the model in the shock initiation regime is tested against recent thin pulse experiments. The model is found to perform equally well in predicting the size-effect curve of ambient, cold, and hot rate sticks. The present work demonstrates the viability of incorporating results from subscale simulations into a continuum-scale reactive flow model.

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

Rapid loss of complex polymers and pyrogenic carbon in subsoils under whole-soil warming

Subsoils contain more than half of soil organic carbon (SOC) and are expected to experience rapid warming in the coming decades. Yet our understanding of the stability of this vast carbon pool under global warming is uncertain. In particular, the fate of complex molecular structures (polymers) remains debated. Here we show that 4.5 years of whole-soil warming (+4 °C) resulted in less polymeric SOC (sum of specific polymers contributing to SOC) in the warmed subsoil (20–90 cm) relative to control, with no detectable change in topsoil. Warming stimulated the subsoil loss of lignin phenols (–17 ± 0%) derived from woody plant biomass, hydrolysable lipids cutin and suberin, derived from leaf and woody plant biomass (–28 ± 3%), and pyrogenic carbon (–37 ± 8%) produced during incomplete combustion. Given that these compounds have been proposed for long-term carbon sequestration, it is notable that they were rapidly lost in warmed soils. We conclude that complex polymeric carbon in subsoil is vulnerable to decomposition and propose that molecular structure alone may not protect compounds from degradation under future warming.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Revealing the etching process of water-soluble Au 25 nanoclusters at the molecular level

Etching (often considered as decomposition) is one of the key considerations in the synthesis, storage, and application of metal nanoparticles. However, the underlying chemistry of their etching process still remains elusive. Here, we use real-time electrospray ionization mass spectrometry to study the reaction dynamics and size/structure evolution of all the stable intermediates during the etching of water-soluble thiolate-protected gold nanoclusters (Au NCs), which reveal an unusual “recombination” process in the oxidative reaction environment after the initial decomposition process. Interestingly, the sizes of NC species grow larger and their ligand-to-metal ratios become higher during this recombination process, which are distinctly different from that observed in the reductive growth of Au NCs (e.g., lower ligand-to-metal ratios with increasing sizes). The etching chemistry revealed in this study provides molecular-level understandings on how metal nanoparticles transform under the oxidative reaction environment, providing efficient synthetic strategies for new NC species through the etching reactions.

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

Hybrid biological-chemical strategy for converting polyethylene into a recyclable plastic monomer using engineered Corynebacterium glutamicum

Converting polyethylene (PE) into valuable materials, particularly ones that are better for the environment than the incumbent plastics, not only helps mitigate environmental issues caused by plastic waste but also alleviates the long-standing problem of microbial fermentation competing with food supplies. However, the inherent robustness of PE due to its strong carbon-carbon bonds and high molecular weight necessitates harsh decomposition conditions, resulting in diverse decomposition outcomes that present significant challenges for downstream applications, especially for bioconversion. In this study, we demonstrate a hybrid biological-chemical conversion process for PE, converting its decomposition products, namely short-chain diacids, into a monomer, β-keto-δ-lactone (BKDL), for highly recyclable polydiketoenimine plastics using engineered Corynebacterium glutamicum. Since BKDL synthesis requires a substantial supply of malonyl-CoA, we employed an alternative biosynthesis pathway that leverages C. glutamicum's natural proficiency in amino acid production. We optimized this pathway in vivo by minimizing carbon loss to CO2 and byproducts, improving the transporter system, and maximizing co-factor regeneration. Furthermore, we co-optimized the PE deconstruction process to produce predominantly C4 to C6 diacids and integrated three catabolic pathways into the engineered strain to enhance diacid utilization, maximizing the carbon conversion from PE. Finally, an engineered polyketide synthase was introduced into C. glutamicum to enable BKDL synthesis. This work demonstrates the potential of a chemo-biological hybrid strategy for recycling plastic waste, highlighting its promise in addressing environmental challenges and promoting sustainable materials.

Zhan, Chunjun↗
DOE OSTI2024

Observation of a super-tetrahedral cluster of acetonitrile-solvated dodecaborate dianion via dihydrogen bonding

We launched a combined negative ion photoelectron spectroscopy and multiscale theoretical investigation on the geometric and electronic structures of a series of acetonitrile-solvated dodecaborate clusters, i.e., B 12 H 12 2- ·nCH 3 CN (n = 1–4). The electron binding energies of B 12 H 12 2- ·nCH 3 CN are observed to increase with cluster size, suggesting their enhanced electronic stability. B3LYP-D3(BJ)/ma-def2-TZVP geometry optimizations indicate each acetonitrile molecule binds to B 12 H 12 2- via a threefold dihydrogen bond (DHB) B3–H3 ⋮⋮⋮ H3C–CN unit, in which three adjacent nucleophilic H atoms in B 12 H 12 2- interact with the three methyl hydrogens of acetonitrile. The structural evolution from n = 1 to 4 can be rationalized by the surface charge redistributions through the restrained electrostatic potential analysis. Notably, a super-tetrahedral cluster of B 12 H 12 2- solvated by four acetonitrile molecules with 12 DHBs is observed. The post-Hartree–Fock domain-based local pair natural orbital- coupled cluster singles, doubles, and perturbative triples [DLPNO-CCSD(T)] calculated vertical detachment energies agree well with the experimental measurements, confirming the identified isomers as the most stable ones. Furthermore, the nature and strength of the intermolecular interactions between B 12 H 12 2- and CH 3 CN are revealed by the quantum theory of atoms-in-molecules and the energy decomposition analysis. Ab initio molecular dynamics simulations are conducted at various temperatures to reveal the great kinetic and thermodynamic stabilities of the selected B 12 H 12 2- ·CH 3 CN cluster. The binding motif in B 12 H 12 2- ·CH 3 CN is largely retained for the whole halogenated series B 12 X 12 2- ·CH 3 CN (X = F–I). This study provides a molecular-level understanding of structural evolution for acetonitrile-solvated dodecaborate clusters and a fresh view by examining acetonitrile as a real hydrogen bond (HB) donor to form strong HB interactions.

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

Leveraging High-resolution Molecular Composition of Soil Organic Matter to Enhance Carbon Cycling Modeling

Soils store more carbon than the atmosphere and vegetation combined, yet Earth system models still struggle to predict how this vast reservoir will respond to environmental change. A central limitation is that most soil biogeochemical models represent organic matter using bulk conceptual pools or chemically homogeneous fractions, preventing direct use of rapidly expanding molecular-scale datasets. Here we develop and test a new soil decomposition framework that explicitly integrates high-resolution information on organic matter composition. First, we construct a molecularly informed litter decomposition module in which plant inputs are partitioned into five functional compound classes—carbohydrates, proteins, lignin-like aromatics, lipids, and carbonyls—using a molecular mixing model calibrated to solid-state 13 C Nuclear Magnetic Resonance (NMR) spectra. Class-specific kinetics, lignin-dependent physical protection, and substrate-driven microbial carbon use efficiency allow the module to capture metabolic tradeoffs associated with enzyme production and nutrient limitation. We then embed this litter module within a microbially explicit whole-soil model that tracks the transformation of these compound classes through particulate organic matter, dissolved organic matter, mineral-associated organic matter, and microbial biomass. High-resolution Fourier Transform Ion Cyclotron Resonance mass spectrometry (FTICR-MS) data are used to link internal pools to measurable soil organic matter fractions and to constrain key process parameters. Applications at soil-core and ecosystem scales demonstrate that the new model reproduces observed soil respiration dynamics while providing mechanistic attribution of CO 2 fluxes to specific chemical classes and pools. Compared to existing frameworks such as the Community Land Model soil biogeochemistry module and the Millennial model, our approach maintains competitive predictive skill while substantially improving interpretability and opportunities for data–model integration. This work illustrates a viable pathway for leveraging molecular-scale observations to reduce structural uncertainty in soil carbon–climate feedback projections.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Pyrolysis of bio-derived dioxolane fuels: A ReaxFF molecular dynamics study

Alkyl-substituted 1,3-dioxolanes, including 4,5-dimethyl-2-pentan-3-yl-1,3-dioxolane (Fuel 1), 4,5-dimethyl-2-pentyl-1,3-dioxolane (Fuel 2), and 2-(heptan-3-yl)-4,5-dimethyl-1,3-dioxolane (Fuel 3), have been recently suggested as potential biodiesels. In this paper, we investigate the initial pyrolysis of the alkyl-substituted 1,3-dioxolanes at high temperatures using ReaxFF molecular dynamics (MD) simulations. We analyze the decomposition rate, reaction mechanism, and product distribution in the pyrolysis of the three alkyl 1,3-dioxlanes. The three fuels primarily decompose to 4,5-dimethyl-1,3-dioxolane radical and hydrocarbons derived from the alkyl side-chains. The further decomposition of 4,5-dimethyl-1,3-dioxolane radical primarily leads to 2-C 4 H 8 and CO 2 within a few decomposition steps. The hydrocarbon product distribution is significantly affected by the molecular structure of the alkyl side-chain, which would have a strong influence on the sooting tendency of these fuels. The ReaxFF simulations predict that the order of sooting tendency would be Fuel 3 > Fuel 1 > Fuel 2, which agrees with the measured sooting tendency trend. Based on the pyrolysis mechanism identified by ReaxFF, we propose a new alkyl dioxolane, 4-hexyl-5-methyl-1,3-dioxolane (Fuel 4), which might produce even less soot, by modifying the molecular structure of Fuel 2. Our ReaxFF simulation shows that Fuel 4 produce much less C 4 H 8 , an effective non-aromatic soot precursor, than Fuel 2. Moreover, more carbon atoms are bonded to each oxygen atom in Fuel 4 than Fuel 2, which would help reduce soot yield by removing more carbon atoms from the soot-producing pool of species. The major decomposition pathways identified in this work can be used to develop chemical kinetic models for 1,3-dioxolane based compounds, as biodiesel components, applicable to combustion engine simulations. We also demonstrate that the chemical kinetic insight offered by ReaxFF simulations can be used to design new fuel molecules with more desired properties.

09 BIOMASS FUELS↗
DOE OSTI2026

SPRUCE FT-ICR MS, Bulk Chemistry, and Mass Loss from Litter Decomposition Study in Experimental Plots, Marcell Experimental Forest, Minnesota, 2015-2017

This dataset contains molecular, bulk chemical, and mass loss measurements from a litter decomposition study at the Spruce and Peatland Responses Under Changing Environments (SPRUCE) experimental site within the Marcell Experimental Forest in northern Minnesota, USA. This site is in a Sphagnum spp. ombrotrophic bog forest. Litterbags were deployed into the peat in September 2015 across three warming levels (+0, +4.5, and +9°C) under ambient and elevated carbon dioxide (CO₂ - +500 ppm) and retrieved after roughly 0.5, 1, and 2 years of field incubation (2015-09-23 to 2017-08-02). Litterbags containing six peatland litter types: black spruce needles (Picea mariana - SPL), spruce fine roots (SPR), Sphagnum angustifolium (ANG), Sphagnum magellanicum (MAG), Labrador tea leaves (Rhododendron groenlandicum - LTL), and Labrador tea roots (LTR). Molecular composition of water-soluble organic matter extracts was characterized using Fourier Transform Ion Cyclotron Resonance Mass Spectrometry (FT-ICR MS) at 9.4 Tesla, operated in negative ion mode with electrospray ionization, providing molecular formula assignments and compound-class distributions across the decomposition time series. Bulk chemical characterization included elemental analysis (percent carbon, nitrogen, and phosphorus) and Fourier Transform Infrared Spectroscopy (FTIR) to quantify functional group composition. Litter mass loss was tracked gravimetrically at each retrieval interval, expressed as percent mass remaining relative to initial dry mass for each litter type and treatment combination. These data are valuable for understanding how vegetation shifts driven by increased atmospheric CO2 and temperature in peatlands alter litter inputs and organic matter stabilization trajectories, with implications for projecting and modeling peatland carbon cycling. This dataset contains two data files in comma-separated value (.csv) format. Additional metadata are provided: two data dictionaries and a file-level metadata file in comma separate (.csv) format and a user guide in PDF (*.pdf) format.

decomposition↗
DOE OSTI2025

Deciphering decomposition pathways of high explosives with cryogenic X-ray Raman spectroscopy

We employed cryogenic X-ray Raman spectroscopy to investigate the early-stage decomposition of the high explosive molecule hexanitrohexaazaisowurtzitane (CL-20). By systematically varying the radiation dose under cryogenic conditions, we induced the decomposition of the molecule using ionizing radiation and observed the evolution of spectral features at the carbon, nitrogen, and oxygen K edges. Through extensive first-principles calculations, we identified key intermediates in the early stages of the decomposition process, resulting from C–C and C–N bond cleavage which leads to the opening of the internal cage structure. A detailed analysis of spectral trends and fingerprints provided evidence supporting N–NO 2 homolytic cleavage as the primary initial decomposition pathway. The combination of advanced core-level spectroscopy methods and state-of-the-art theoretical calculations enabled a comprehensive characterization of the molecular changes induced by controlled radiation dose exposures. In conclusion, our findings establish a benchmark for understanding the decomposition chemistry of high-explosive materials, offering important insights into their stability and reactivity under extreme conditions.

X-ray Raman spectroscopy↗
DOE OSTI2025

Toward an Understanding of Linear Scaling Relations through Energy Decomposition Analysis

The discovery of linear scaling relations has fundamentally changed the field of heterogeneous catalysis. The scaling relations have been rationalized based on the d-band theory, specifically a separation of sp and d electron contributions to adsorption energies. Within the framework of energy decomposition analysis, a full understanding of such a separation would require one to further break down the adsorption energy into distinct energy components such as electrostatics, polarization, charge transfer, and van der Waals interactions, and to examine the sp and d contributions to each of them. As a step in this direction, we analyzed the interaction energy between CH x (x = 1–4) adsorbates and fcc(100) transition metal surfaces (M = Cu, Ag, Au, Rh, and Pt), with the surfaces represented both as slabs in plane-wave density functional theory (pw-DFT) calculations and as atomic clusters in atomic-orbital basis density functional theory (ao-DFT) calculations. Through an absolutely localized molecular orbital (ALMO) based energy decomposition analysis of the ao-DFT adsorption energy, each of the interaction energy components (electrostatics, polarization, van der Waals, and charge transfer) was found to follow its own scaling relations, with an intricate interplay among these energy components yielding the overall scaling relations for the total adsorption energies. Using the recently introduced ALMO-based polarization and charge-transfer analysis schemes, we further dissected polarization into metal surface and adsorbate contributions, and charge transfer into metal → adsorbate and adsorbate → metal contributions. The contributions from the sp and d electrons of the metal to these terms were further quantified, and the dominant role of the metal d electrons was reaffirmed. These results shed light on how CHx adsorbates interact with metal surfaces and further reveal the physical origin of the scaling relations.

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

Detecting Thermally-Induced Spinodal Decomposition with Picosecond Ultrasonics in Cast Austenitic Stainless Steels

Destructive techniques to monitor nuclear reactor component health may not always be available during service, as they are time-consuming and often require pre-installed inspection coupons. Non-destructive evaluation (NDE) techniques can bridge this gap by rapidly identifying the state of mission-critical reactor components, via inference between NDE-measurable material properties and those of ultimate interest, such as ductility and toughness. Here, we demonstrate one such inference about the health of thermally aged cast austenitic stainless steels. Observations of surface acoustic wave peak (SAW) splitting correlate with spinodal decomposition-induced embrittlement as destructively measured by Charpy impact energy. Elastodynamic calculations and molecular dynamics simulations of the effects of spinodal decomposition on elastic moduli support that the new acoustic modes present are due to stiffening in the δ-ferrite domains. Finally, this discovery enables one to probe structure-property relationships in materials in a greatly accelerated manner, suggesting that similar inference methods can be used to determine material fitness-for-service, or to quickly uncover new structure-property relationships.

304-type stainless steel↗
DOE OSTI2023

Temperature-based reactive flow model for TATB plastic bonded explosives

A new reactive flow model is presented for TATB-based plastic bonded explosives, applicable to shock initiation and detonation problems of differing initial temperature. Temperature disequilibrium is assumed between unreacted explosive, material in the vicinity of compressed defects (called hot spots), and detonation products. The model incorporates temperature-dependent decomposition reaction rates. Arrhenius model parameters were derived from subscale quantum-based molecular dynamics simulations of shock-induced TATB decomposition. Further, a model of detonation carbon aggregation is incorporated, describing the late-time slow release of energy inherent to detonation in TATB-based materials. The predictive ability of the model in the shock initiation regime is tested against recent thin pulse experiments on PBX 9502. Here, the model is found to perform equally well in predicting the size effect curve of ambient, cold, and hot PBX 9502 rate sticks.

Chemistry - Chemical explosives↗
DOE OSTI2025

Similarity Metric for Data Optimization and Efficient Training of Reactive Machine Learning Force Fields for Hydrocarbon Radiolysis

Radiolysis is a common approach to sterilize polymers, chemically modify them for upcycling, and accelerate their decomposition for recycling purposes. Reactive molecular dynamics (MD) simulations provide a powerful tool to generate atomic-level trajectories of the reactive processes and quantify radiolytic chemical degradation pathways. For this, machine learning (ML) surrogate models for reactive force fields with quantum mechanical accuracy are now widely used, which require ML training data sets that can provide information on atomic environments for target chemical systems. However, radiolysis chemistry can be highly complex and diverse, which poses significant challenges for generating training data to parametrize ML models. In this regard, we developed a method for optimizing the training data set using a cosine similarity metric to help guide training set selection for radiolysis of polyethylene, a model hydrocarbon polymer, as well as to enhance the transferability of our reactive ML force field (MLFF) to a variety of molecular and polymeric systems. Our approach performs atom-by-atom comparisons between local atomic environments to pinpoint important data points associated with rare and localized events, such as radiolysis damage within structures. We apply this approach to train the Chebyshev Interaction Model for Efficient Simulation (ChIMES) MLFF model, which expresses the atomic interaction potentials in terms of linear combinations of many-body Chebyshev polynomials. We first show that our method can reduce our training set size by ∼70% while improving overall accuracy compared to more standard MD model fitting approaches. We then validate our optimum model against diverse hydrocarbon simulation data, including simple alkanes and systems with unsaturated carbon bonds, over a wide range of thermodynamic conditions. Finally, we use our ChIMES model to perform MD simulations of radiolytic damage with large-scale systems that help avoid system size effects. Overall, our approach yields an MD force field that retains most of the accuracy of the underlying quantum method while yielding many orders of improvement in computational efficiency. In conclusion, our efforts will have impact on future hydrocarbon polymer radiolysis studies, where the chemical details of the polymer–radiation interactions can have a strong effect on the resulting products observed in experiments.

Hydrocarbons↗
DOE OSTI2026

Scalable Quantum Monte Carlo Method for Polariton Chemistry via Mixed Block Sparsity and Tensor Hypercontraction Method

We present a reduced-scaling auxiliary-field quantum Monte Carlo (AFQMC) framework designed for large molecular systems and ensembles, with or without coupling to optical cavities. Our approach leverages the natural block sparsity of the Cholesky decomposition (CD) of electron repulsion integrals in molecular ensembles and employs tensor hypercontraction (THC) to efficiently compress low-rank Cholesky blocks. By representing the Cholesky vectors in a mixed format, keeping high-rank blocks in block-sparse form and compressing low-rank blocks with THC, we reduce the scaling of exchange-energy evaluation from quartic to robust cubic in the number of molecular orbitals N, while lowering memory from cubic toward quadratic. Benchmark analyses on one-, two-, and three-dimensional molecular ensembles (up to ∼1,200 orbitals) show that (a) the number of nonzeros in Cholesky tensors grows linearly with system size across dimensions; (b) the average numerical rank increases sublinearly and does not saturate at these sizes; and (c) rank heterogeneity─some blocks nearly full rank and many low rank, naturally motivates the proposed mixed block sparsity and THC scheme for efficient calculation of exchange energy. In conclusion, we demonstrate that the mixed scheme yields cubic wall-time scaling with favorable prefactors and preserves AFQMC accuracy.

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

Tracing priming effects in palsa peat carbon dynamics using a stable isotope-assisted metabolomics approach

Introduction: Peatlands store up to a third of global soil carbon, and in high latitudes their litter inputs are increasing and changing in composition under climate change. Although litter significantly influences peatland carbon and nutrient dynamics by changing the overall lability of peatland organic matter, the physicochemical mechanisms of this impact—and thus its full scope—remain poorly understood. Methods: We applied multimodal metabolomics (UPLC-HRMS, 1 H NMR) paired with 13 C Stable Isotope-Assisted Metabolomics (SIAM) to track litter carbon and its potential priming effects on both existing soil organic matter and carbon gas emissions. Through this approach, we achieved molecule-specific tracking of carbon transformations at unprecedented detail. Results: Our analysis revealed several key findings about carbon dynamics in palsa peat. Microbes responded rapidly to litter addition, producing a short-term increase in CO 2 emissions, fueled nearly exclusively by transformations of litter carbon. Litter inputs significantly contributed to the organic nitrogen pool through amino acids and peptide derivatives, which served as readily accessible nutrient sources for microbial communities. We traced the fate of plant-derived polyphenols including flavonoids like rutin, finding evidence of their degradation through heterocyclic C-ring fission, while accumulation of some polyphenols suggested their role in limiting overall decomposition. The SIAM approach detected subtle molecular changes indicating minimal and transient priming activity that was undetectable through conventional gas measurements alone. This transient response was characterized by brief microbial stimulation followed by rapid return to baseline metabolism. Pre-existing peat organic matter remained relatively stable; significant priming of its consumption was not observed, nor was its structural alteration. Discussion: This suggests that while litter inputs temporarily increase CO 2 emissions, they don’t sustain long-term acceleration of stored carbon decomposition or substantially decrease peat’s carbon store capacity. Our findings demonstrate how technological advancements in analytical tools can provide a more detailed view of carbon cycling processes in complex soil systems.

54 ENVIRONMENTAL SCIENCES↗