SEARCH · Search NASA
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.
The catalytic mechanism and reaction network of oxidative coupling of methane over La 2 O 3 are thoroughly investigated by density functional theory calculations.
We report accurate identification of descriptors for catalytic activities has long been essential to the in-depth understanding of catalysis and recently to set the basis for catalyst screening. However, commonly used methods suffer from low accuracy in predictability. This study reports an enhanced approach to accurately identify the descriptors from a kinetic dataset using a machine learning (ML) surrogate model. CO hydrogenation to methanol over Cu-based catalysts was taken as a case study. Our model captures not only the contribution from individual elementary steps but also the interaction between relevant steps within a reaction network, which was found to be essential for high accuracy. As a result, six effective descriptors are identified, which are accurate enough to ensure the trained gradient boosted regression (GBR) model for good prediction of the methanol turnover frequency (TOF) over metal (M)-doped Cu(111) model surfaces (M = Au, Cu, Pd, Pt, Ni). More importantly, going beyond the purely mathematical ML model, the catalytic role of each identified descriptor can be revealed by using model-agnostic interpretation tools, which enhances the insight into the promoting effect of alloying. The trained GBR model outperforms the conventional derivative-based methods in terms of both the predictability and the mechanism understanding. It opens alternative possibilities toward accurate descriptor-based rational catalyst optimization.
To understand the gas–surface chemistry above the thermal protection system of a hypersonic vehicle, it is necessary to map out the kinetics of key elementary reaction steps. In this work, extensive periodic density functional theory (DFT) calculations are performed to elucidate the interaction of atomic oxygen and nitrogen with both the basal plane and edge sites of highly oriented pyrolytic graphite (HOPG). Reaction energies and barriers are determined for adsorption, desorption, diffusion, recombination, and several reactions. These DFT results are compared with the most recent finite-rate model for air–carbon ablation. Our DFT results corroborated some of the parameters used in the model but suggest that further refinement may be necessary for others. The calculations reported here will help to establish a predictive kinetic model for the complex reaction network present under hypersonic flight conditions.
Electrochemical CO reduction can serve as a sequential step in the transformation of CO 2 into multicarbon fuels and chemicals. In this study, we provide insights on how to steer selectivity for CO reduction almost exclusively toward a single multicarbon oxygenate by carefully controlling the catalyst composition and its surrounding reaction conditions. Under alkaline reaction conditions, we demonstrate that planar CuAg electrodes can reduce CO to acetaldehyde with over 50% Faradaic efficiency and over 90% selectivity on a carbon basis at a modest electrode potential of −0.536 V vs. the reversible hydrogen electrode. The Faradaic efficiency to acetaldehyde was further enhanced to 70% by increasing the roughness factor of the CuAg electrode. Density functional theory calculations indicate that Ag ad-atoms on Cu weaken the binding energy of the reduced acetaldehyde intermediate and inhibit its further reduction to ethanol, demonstrating that the improved selectivity to acetaldehyde is due to the electronic effect from Ag incorporation. These findings will aid in the design of catalysts that are able to guide complex reaction networks and achieve high selectivity for the desired product.
Hemostasis relies on a reaction network of serine proteases and their cofactors to form a blood clot. Coagulation factor IXa (protease) plays an essential role in hemostasis as evident from the bleeding disease associated with its absence. RNA aptamers specifically targeting individual coagulation factors have potential as anticoagulants and as probes of the relationship between structure and function. Here, we report X-ray structures of human factor IXa without a ligand bound to the active site either in the apo-form or in complex with an inhibitory aptamer specific for factor IXa. The aptamer binds to an exosite in the catalytic domain and allosterically distorts the active site. Our studies reveal a conformational ensemble of IXa states, wherein large movements of Trp 215 near the active site drive functional transitions between the closed (aptamer-bound), latent (apo), and open (substrate-bound) states. The latent state of the apo-enzyme may bear on the uniquely poor catalytic activity of IXa compared to other coagulation proteases. The exosite, to which the aptamer binds, has been implicated in binding VIIIa and heparin, both of which regulate IXa function. Our findings reveal the importance of exosite-driven allosteric modulation of IXa function and new strategies to rebalance hemostasis for therapeutic gain.
The design and implementation of future flowing liquid-lithium plasma-facing components (LLPFCs) will be dependent on several factors. Of course, one of the most important is the need to be able to deal with high heat fluxes incident on the surface of the LLPFCs, but there are also several other important liquid-metal behaviors that have been identified for their critical impact on the feasibility of a LLPFC. One of these is the ability to constantly wet 100% of the plasma-facing component area and the best way to achieve that. Another key point is knowing and understanding the erosion and corrosion of the surfaces subject to a flowing liquid-lithium system and the ability for hydrogen and helium uptake by the system. The Center for Plasma Material Interactions (CPMI) has been tasked with looking at these various issues. The Mock-up Entry module for EAST device was used to investigate wetting and erosion effects and to design a suitable distribution and collection system with a liquid-lithium loop. The vapor shielding effects of lithium on the surface were also modeled and studied. A model coupling CRANE, an open-source global reaction network solver, and Zapdos, a plasma transport solver, is being developed to better understand the dynamics of the vapor cloud. Experiments on the Magnum-PSI at the Dutch Institute for Fundamental Energy Research have been carried out to study the vapor shielding effect and obtain experimental benchmarks to verify the model. Also, initial experiments using the Hybrid Illinois Device for Research and Applications have been performed to understand the pumping effects of lithium on helium. Experiments with a drop of liquid lithium (~100 mg) into a helium plasma have shown the ability of lithium to take out the cold recycling helium gas as well as hydrogen and oxygen impurity gases. The improvement in plasma performance was significant, and further understanding of this effect will have impacts on how future LLPFCs will be designed. Further investigation into the exact mechanism for helium pumping by lithium needs to be performed in the future. Finally, this paper presents a summary of the results obtained at the CPMI.
Abstract Compared to other computational physics areas such as codes for general computational fluid dynamics, the documentation of verification methods for plasma fluid codes remains under developed. Current analytical solutions for plasma are often highly limited in terms of testing highly coupled physics, due to the harsh assumptions needed to derive even simple plasma equations. This work highlights these limitations, suggesting the method of manufactured solutions (MMSs) as a potential option for future verification efforts. To demonstrate the flexibility of MMS in verifying these highly coupled systems, the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework was utilized. Thanks to the MOOSE framework’s robustness and modularity, as well as to its physics module capabilities and ecosystem applications (i.e. Zapdos and the chemical reaction network) developed for plasma physics modeling and simulation, this report lays the groundwork for a structured method of conducting plasma fluid code verification.
At temperatures below the QCD phase transition, any substantial lepton number in the Universe can only be present within the neutrino sector. In this work, we systematically explore the impact of a non-vanishing lepton number on Big Bang Nucleosynthesis (BBN) and the Cosmic Microwave Background (CMB). Relying on our recently developed framework based on momentum averaged quantum kinetic equations for the neutrino density matrix, we solve the full BBN reaction network to obtain the abundances of primordial elements. We find that the maximal primordial total lepton number L allowed by BBN and the CMB is -0.12 (-0.10) ≤ L ≤ 0.13 (0.12) for NH (IH), while specific flavor directions can be even more constrained. This bound is complementary to the limits obtained from avoiding baryon overproduction through sphaleron processes at the electroweak phase transition since, although numerically weaker, it applies at lower temperatures and is obtained completely independently. We publicly release the C++ code COFLASY-C on GitHub (https://github.com/mariofnavarro/COFLASY/tree/COFLASY-C) which solves for the evolution of the neutrino quantum kinetic equations numerically.
Abstract The 2017 detection of a kilonova coincident with gravitational-wave emission has identified neutron star mergers as the major source of the heaviest elements and dramatically constrained alternative theories of gravity. Observing a population of such sources has the potential to transform cosmology, nuclear physics, and astrophysics. However, with only one confident multi-messenger detection currently available, modelling the diversity of signals expected from such a population requires improved theoretical understanding. In particular, models that are quick to evaluate and are calibrated with more detailed multi-physics simulations are needed to design observational strategies for kilonovae detection and to obtain rapid-response interpretations of new observations. We use grey-opacity models to construct populations of kilonovae, spanning ejecta parameters predicted by numerical simulations. Our modelling focuses on wavelengths relevant for upcoming optical surveys, such as the Rubin Observatory Legacy Survey of Space and Time (LSST). In these simulations, we implement heating rates that are based on nuclear reaction network calculations. We create a Gaussian-process emulator for kilonova grey opacities, calibrated with detailed radiative transfer simulations. Using recent fits to numerical relativity simulations, we predict how the ejecta parameters from binary neutron star (BNS) mergers shape the population of kilonovae, accounting for the viewing-angle dependence. Our simulated population of BNS mergers produce peak i-band absolute magnitudes of −20 ≤ Mi ≤ −11. A comparison with detailed radiative transfer calculations indicates that further improvements are needed to accurately reproduce spectral shapes over the full light curve evolution.
The presence of a strong shallow heat source of unknown origin in accreting neutron star crusts has been inferred by analyzing x-ray observations of their cooling in quiescence. Here, we model the cooling of KS 1731–260 using realistic crust compositions and nuclear heating and cooling sources from detailed nuclear reaction network calculations. We find that the required strength of the shallow heat source in KS 1731–260 is reduced by more than a factor of 3 compared to previous analysis, a 5𝜎 difference that alleviates the need for exotic solutions. Our analysis also suggests the existence of an impure nuclear pasta layer in the inner crust of KS 1731–260, though future observations will provide more stringent constraints. In addition, we obtain constraints on the dominant surface burning modes of KS 1731–260 over its history.
Fission data are essential inputs to reaction networks involved in nucleosynthesis simulations and nuclear forensics. In such applications as well as in the description of multichance fission, the characteristics of fission for odd-mass nuclei are just as important as those for even-even nuclei. The fission theories that aim at explicitly describing fission dynamics are typically based on some variant of the nuclear mean-field theory. In such cases, the treatment of systems with an odd number of particles is markedly more involved, both formally and computationally. In this article, we use the blocking prescription of the Hartree-Fock-Bogoliubov theory with Skyrme energy functionals to compute the deformation properties of odd-mass uranium isotopes. We show that the resulting fission fragment distributions depend quite significantly on the spin of the odd neutron. By direct calculation of the spin distribution of the fissioning nucleus, we propose a methodology to rigorously predict the charge and mass distributions in odd-mass nuclei.
The nuclear reaction network within the interior of the Sun is an efficient MeV physics factory and can produce long-lived particles generic to dark sector models. In this work we consider the sensitivity of satellite instruments, primarily the RHESSI spectrometer, that observe the quiet Sun in the MeV regime where backgrounds are low. We find that quiet Sun observations offer a powerful and complementary probe in regions of parameter space, where the long-lived particle decay length is longer than the radius of the Sun and shorter than the distance between the Sun and Earth. We comment on connections to recent model-building work on heavy neutral leptons coupled to neutrinos and high-quality axions from mirror symmetries.
The Rare-RI Ring (R3) is a recently commissioned cyclotronlike storage ring mass spectrometer dedicated to mass measurements of exotic nuclei far from stability at Radioactive Isotope Beam Factory (RIBF) in RIKEN. The first application of mass measurement using the R3 mass spectrometer at RIBF is reported. Rare isotopes produced at RIBF— 127 Sn, 126 In, 125 Cd, 124 Ag, 123 Pd—were injected in R3. Masses of 126 In, 125 Cd, and 123 Pd were measured whereby the mass uncertainty of 123 Pd was improved. This is the first reported measurement with a new storage ring mass spectrometry technique realized at a heavy-ion cyclotron and employing individual injection of the preidentified rare nuclei. The latter is essential for the future mass measurements of the rarest isotopes produced at RIBF. The impact of the new 123 Pd result on the solar r-process abundances in a neutron star merger event is investigated by performing reaction network calculations of 20 trajectories with varying electron fraction Y e . It is found that the neutron capture cross section on 123 Pd increases by a factor of 2.2 and β-delayed neutron emission probability, P 1n , of 123 Rh increases by 14%. The neutron capture cross section on 122 Pd decreases by a factor of 2.6 leading to pileup of material at A = 122, thus reproducing the trend of the solar r-process abundances. The trend of the two-neutron separation energies (S 2n ) was investigated for the Pd isotopic chain. The new mass measurement with improved uncertainty excludes large changes of the S 2n value at N = 77. Finally, such large increase of the S 2n values before N = 82 was proposed as an alternative to the quenching of the N = 82 shell gap to reproduce r-process abundances in the mass region of A = 112 – 124.
Here, we present a joint-likelihood analysis of big bang nucleosynthesis (BBN) and cosmic microwave background (CMB) data, consistently combining likelihoods and taking into account uncertainties in nuclear reaction rates for the first time. Bayesian inference is performed on the baryon abundance and the effective number of neutrino species, 𝑁 eff , using a CMB Boltzmann solver in combination with LINX , a new flexible and efficient BBN code. We marginalize over Planck nuisance parameters and nuclear rates to find 𝑁 eff =3.08$^{+0.15}_{−0.14}$, 2.94$^{+0.16}_{−0.15}$, or 2.96$^{+0.13}_{−0.14}$, for three separate reaction networks. This framework enables robust testing of the lambda cold dark matter paradigm and its variants with CMB and BBN data.
Contemporary computation is expensive, with large language models and artificial intelligence becoming more common in daily life. However, high-performance computing is reaching the limits in speed and energy expenditure, and domain science requires ever-increasing computational capacity, with simulations and data analysis pipelines ever-growing in complexity. As we progress towards post-exascale computation, with the associated high energy costs, new methods of energy-conscious computation are required. Novel analog and hybrid digital-analog systems can overcome these challenges, and chemical reactions offer a promising avenue. Computers based on chemistry can provide compact desktop devices with immense computational power. These devices are readily scalable by considering greater reaction systems or vessels, meeting the high-performance requirements for scientific workflows. In this article, we present ChemComp, a compilation pipeline for the conversion of ordinary differential equations into implementable chemical reactions. We then demonstrate the solving capabilities of ChemComp by emulating a potential chemical reservoir device. We leverage the multi-layer intermediate representation (MLIR) compiler framework to implement an expressive chemical reaction abstraction and propose a path for chemical reaction networks (CRNs) to represent mathematical problems effectively. Combined, we demonstrate a potential workflow that can harness chemistry’s computing power to create energy-efficient, high-performance computation systems for contemporary computing needs.
Organic materials are an attractive choice for structural components due to their light weight and versatility. However, because they decompose at low temperatures relative to tradiational materials they pose a safety risk due to fire and loss of structural integrity. To quantify this risk, analysts use chemical kinetics models to describe the material pyrolysis and oxidation using thermogravimetric analysis. This process requires the calibration of many model parameters to closely match experimental data. Previous efforts in this field have largely been limited to finding a single best-fit set of parameters even though the experimental data may be very noisy. Furthermore the chemical kinetics models are often simplified representations of the true de- composition process. The simplification induces model-form errors that the fitting process cannot capture. In this work we propose a methodology for calibrating decomposition models to thermogravimetric analysis data that accounts for uncertainty in the model-form and experimental data simultaneously. The methodology is applied to the decomposition of a carbon fiber epoxy composite with a three-stage reaction network and Arrhenius kinetics. The results show a good overlap between the model predictions and thermogravimetric analysis data. Uncertainty bounds capture devia- tions of the model from the data. The calibrated parameter distributions are also presented. In conclusion, the distributions may be used in forward propagation of uncertainty in models that leverage this material.
Despite the explosion of soil metagenomic data, we lack a synthesized understanding of patterns in the distribution and functions of soil microorganisms. These patterns are critical to predictions of soil microbiome responses to climate change and resulting feedbacks that regulate greenhouse gas release from soils. To address this gap, we assay 1,512 manually curated soil metagenomes using complementary annotation databases, read-based taxonomy, and machine learning to extract multidimensional genomic fingerprints of global soil microbiomes. Our objective is to uncover novel biogeographical patterns of soil microbiomes across environmental factors and ecological biomes with high molecular resolution. We reveal shifts in the potential for (i) microbial nutrient acquisition across pH gradients; (ii) stress-, transport-, and redox-based processes across changes in soil bulk density; and (iii) greenhouse gas emissions across biomes. We also use an unsupervised approach to reveal a collection of soils with distinct genomic signatures, characterized by coordinated changes in soil organic carbon, nitrogen, and cation exchange capacity and in bulk density and clay content that may ultimately reflect soil environments with high microbial activity. Genomic fingerprints for these soils highlight the importance of resource scavenging, plant-microbe interactions, fungi, and heterotrophic metabolisms. Across all analyses, we observed phylogenetic coherence in soil microbiomes—more closely related microorganisms tended to move congruently in response to soil factors. Collectively, the genomic fingerprints uncovered here present a basis for global patterns in the microbial mechanisms underlying soil biogeochemistry and help beget tractable microbial reaction networks for incorporation into process-based models of soil carbon and nutrient cycling.
The use of differential scanning calorimetry (DSC) to measure the thermal behavior of individual components and electrolyte/electrode combinations is common. However, here we focus on DSC tests on an anode, cathode, and electrolyte (ACE) component combination over a temperature range that includes many of the phase transitions and key reactions (i.e., to 500 °C) that contribute to thermal runaway. This method can help quantify the complex reaction network in a full cell, thereby informing potential safety issues. Here, we used DSC heat flow data from a solid-state Li 0.43 CoO 2 +C+PVDF | LLZO | Li metal ACE sample and its components to quantify key factors affecting results. We focused on three areas: (1) ACE sample preparation and assembly in DSC pans, (2) DSC measurement parameters, and (3) heat flow analysis. Key points include the choice of component ratios (e.g., commercially relevant N:P capacity ratio), the importance of conductive carbon and binder, type of pan used, DSC ramp rate, and integration method used when dealing with broad and overlapping exothermic peaks. This work deepens the scientific basis and best practices for obtaining heat flow data from ACE samples for early-stage evaluation of solid-state and beyond battery safety.