SEARCH · Engineering Papers
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.
Machine learning (ML)-based feature analysis reveals universal design rules regardless of density functional choices. Using the consensus among multiple functionals, we identify robust lead complexes in ML-accelerated chemical discovery.
NOTE: The manuscript associated with this data package is currently in review. The data may be revised based on reviewer feedback. Upon manuscript acceptance, this data package will be updated with the final dataset and additional metadata. This data package is associated with the manuscript “Fire Severity-Dependent Shifts in Vegetation Parameter Sensitivity: A Pre- and Post-Fire Analysis Using ELM-FATES and Explainable AI” submitted to Journal of Advances in Modeling Earth Systems (Zahura et al. 2026). The study examines vegetation physiological parameters controlling pre-fire and post-fire vegetation dynamics. To support this analysis, 73 vegetation parameters in Functionally Assembled Terrestrial Ecosystem Simulator (FATES) (Fisher et al., 2018) , which is coupled with E3SM (Energy Exascale Earth System Model) land model (ELM, ELM-FATES), were perturbed using a Sobol sequence to generate 1,024 ensemble members for two plant functional types: needleleaf evergreen extratropical trees (NEET) and C3 grass. Simulations were conducted for the pre-fire period (2016) and post-fire period (2018–2023). Burn severity was represented by modifying the Nesterov index in FATES to 75,000, 150,000, and 300,000 for low, moderate, and high severity, respectively. A no-fire scenario was also included. Simulations were performed for 16 grid cells in the American River Watershed across different burn severities and plant functional types. XGBoost (eXtreme Gradient Boosting) models were trained using the parameter ensembles and ELM-FATES-simulated outputs, including leaf area index (LAI), gross primary productivity (GPP), aboveground biomass, vegetation evaporation, transpiration, and soil evaporation. Models were trained separately for each year and burn severity, followed by SHAP (SHapley Additive exPlanations) analysis to identify changes in dominant parameters after fire disturbance. For details on how to navigate data packages generated by this project, see https://data.ess-dive.lbl.gov/portals/PNNLRiverCorridorSFA/About. The data package contains the ELM-FATES simulation data. The scripts and data related to the analysis will be added later. The inputs and outputs from ELM-FATES are inside the “FATES” folder. “FATES_domain_surface” contains the domain and surface netcdfs that were used to run ELM-FATES in the study area. “FATES_parameters” contains the 1024 ensembles that were generated using Sobol sequence. “FATES_outputs” folder contains ELM-FATES simulated variables. All files are .csv and .nc (NetCDF).
Van der Waals (vdW) magnets both allow exploration of fundamental 2D physics and offer a route toward exploiting magnetism in next generation information technology, but vdW magnets with complex, noncollinear spin textures are currently rare. We report here the syntheses, crystal structures, magnetic properties and magnetic ground states of four bulk vdW metal–organic magnets (MOMs): FeCl 2 (pym), FeCl 2 (btd), NiCl 2 (pym), and NiCl 2 (btd), pym = pyrimidine and btd = 2,1,3-benzothiadiazole. Using a combination of neutron diffraction and bulk magnetometry we show that these materials are noncollinear magnets. Although only NiCl 2 (btd) has a ferromagnetic ground state, we demonstrate that low-field hysteretic metamagnetic transitions produce states with net magnetization in zero-field and high coercivities for FeCl 2 (pym) and NiCl 2 (pym). By combining our bulk magnetic data with diffuse scattering analysis and broken-symmetry density-functional calculations, we probe the magnetic superexchange interactions, which when combined with symmetry analysis allow us to suggest design principles for future noncollinear vdW MOMs. These materials, if delaminated, would prove an interesting new family of 2D magnets.
Finite-temperature quantum field theories are formulated in terms of Green's functions and self-energies on the Matsubara axis. In multiorbital systems, these quantities are related to positive semidefinite matrix-valued functions of the Carathéodory and Schur class. Analysis, interpretation, and evaluation of derived quantities such as real-frequency response functions requires analytic continuation of the off-diagonal elements to the real axis. Here we derive the criteria under which such functions exist for given Matsubara data and present an interpolation algorithm that intrinsically respects their mathematical properties. For small systems with precise Matsubara data, we find that the continuation exactly recovers all off-diagonal and diagonal elements. In real-materials systems, we show that the precision of the continuation is sufficient for the analytic continuation to commute with the Dyson equation, and we show that the truncation of the off-diagonal self-energy elements leads to considerable approximation artifacts.
ABSTRACT In this work, we present the galaxy clustering measurements of the two DES lens galaxy samples: a magnitude-limited sample optimized for the measurement of cosmological parameters, maglim, and a sample of luminous red galaxies selected with the redmagic algorithm. maglim/redmagic sample contains over 10 million/2.5 million galaxies and is divided into six/five photometric redshift bins spanning the range z ∈ [0.20, 1.05]/z ∈ [0.15, 0.90]. Both samples cover 4143 $\deg ^2$ over which we perform our analysis blind, measuring the angular correlation function with an S/N ∼ 63 for both samples. In a companion paper, these measurements of galaxy clustering are combined with the correlation functions of cosmic shear and galaxy–galaxy lensing of each sample to place cosmological constraints with a 3 × 2pt analysis. We conduct a thorough study of the mitigation of systematic effects caused by the spatially varying survey properties and we correct the measurements to remove artificial clustering signals. We employ several decontamination methods with different configurations to ensure the robustness of our corrections and to determine the systematic uncertainty that needs to be considered for the final cosmology analyses. We validate our fiducial methodology using lognormal mocks, showing that our decontamination procedure induces biases no greater than 0.5σ in the (Ωm, b) plane, where b is the galaxy bias.
Localized bases play an important role in understanding electronic structure. In periodic insulators, a natural choice of localized basis is given by the Wannier functions which depend on a choice of unitary transform known as a gauge transformation. Over the past few decades, there have been many works that have focused on optimizing the choice of the gauge so that the corresponding Wannier functions are maximally localized or reflect some symmetry of the underlying system. In this work, we consider fully nonperiodic materials where the usual Wannier functions are not well defined and gauge optimization is impractical. To tackle the problem of calculating exponentially localized generalized Wannier functions in both periodic and nonperiodic systems, we discuss the ‘iterated projected position (IPP)” algorithm. Here, the IPP algorithm is based on matrix diagonalization and therefore unlike optimization-based approaches, it does not require initialization and cannot get stuck at a local minimum. Furthermore, the IPP algorithm is guaranteed by a rigorous analysis to produce exponentially localized functions under certain mild assumptions. We numerically demonstrate that the IPP algorithm can be used to calculate exponentially localized bases for the Haldane model, the Kane-Mele model (in both Z 2 invariant even and Z 2 invariant odd phases), and the p x + ip y model on a quasicrystal lattice.
The application of spatially resolved mass spectrometry (MS) and MS imaging approaches for studying biomolecular processes in the kidney is rapidly growing. These powerful methods, which enable label-free and multiplexed detection of many molecular classes across omics domains (including metabolites, drugs, proteins and protein post-translational modifications), are beginning to reveal new molecular insights related to kidney health and disease. Further, the complexity of the kidney often necessitates multiple scales of analysis for interrogating biofluids, whole organs, functional tissue units, single cells and subcellular compartments. Various MS methods can generate omics data across these spatial domains and facilitate both basic science and pathological assessment of the kidney. Optimal processes related to sample preparation and handling for different MS applications are rapidly evolving. Emerging technology and methods, improvement of spatial resolution, broader molecular characterization, multimodal and multiomics approaches and the use of machine learning and artificial intelligence approaches promise to make these applications even more valuable in the field of nephology. Overall, spatially resolved MS and MS imaging methods have the potential to fill much of the omics gap in systems biology analysis of the kidney and provide functional outputs that cannot be obtained using genomics and transcriptomic methods.
Single-atom catalysts (SACs) have aroused tremendous interest over the past decade, particularly in the community of energy and environment-related electrocatalysis. A rapidly growing number of recent publications have recognized it as a promising candidate with maximum atomic utilization, distinct activity, and selectivity in comparison to bulk catalysts and nanocatalysts. However, the complexity of localized coordination environments and the dispersion of isolated sites lead to significant difficulties when it comes to gaining insight into the intrinsic behavior of electrocatalytic reactions. Furthermore, the low metal loadings of most SACs make conventional ensemble measurements less likely to be accurate on the subnanoscale. Thus, it remains challenging to probe the activity and properties of individual atomic sites by available commercial instruments and analytical methods. In spite of this, continuing efforts have lately focused on the development of advanced measurement methodologies, which are very useful to the fundamental understanding of SACs. There have recently been a number of in situ/operando techniques applied to SACs, such as electron microscopy, spectroscopy, and other analysis methods, which support relevant functions to identify the active sites and reaction intermediates and to investigate the dynamic behavior of localized structures of the catalytic sites. This Account aims to present recent electrochemical probing techniques which can be used to identify single-atomic catalytic sites within solid supports. First, we describe the basic principles of molecular probe methods for the study and analysis of electrocatalytic site behavior. In particular, the in situ probing technique enabled by surface interrogation scanning electrochemical microscopy (SISECM) can measure the active site density and kinetic rate with high resolution. An alternative electrochemical probing technique is further demonstrated on the basis of single-entity electrochemistry, which allows the unique electrochemical imaging of the size and catalytic rate of single atoms, molecules, and clusters. Next, the merits and limitations of different electrochemical techniques are then discussed, along with perspectives for future prospects. Apart from this, we further showcase the powerful capability of emerging electrochemical probing techniques for determining significant effects and properties of SACs for various electrocatalytic reactions, including oxygen reduction and evolution, hydrogen evolution, and nitrate reduction. Overall, electrochemical techniques with atomic resolution have greatly increased opportunities for observing, measuring, and understanding the surface and interface chemistry during energy conversion. In the future, it is anticipated that the development of electrochemical probing techniques will be advanced with innovative perspectives on the behavior and features of SACs. We hope that this Account can contribute in several ways to promoting the fundamental knowledge and technical progress of emerging electrochemical measurements for studying SACs.
Pseudomonas aeruginosa CtpA is a carboxyl-terminal processing protease that partners with the outer membrane lipoprotein LbcA to degrade at least five cell wall-associated proteins, four of which are cell wall hydrolases. This activity plays an important role in supporting P. aeruginosa virulence in a mouse model of acute pneumonia. However, almost nothing is known about the molecular mechanisms underlying CtpA and LbcA function. Here, we used structural analysis to show that CtpA alone assembles into an inactive hexamer comprising a trimer of dimers, which limits its substrate access and prevents nonspecific degradation. The adaptor protein LbcA is a right-handed open spiral with 11 tetratricopeptide repeats, which might wrap around a substrate to deliver it to CtpA for degradation. By structure-guided mutagenesis and functional assays, we also showed that the interfaces of the CtpA trimer of dimers and an N-terminal helix of LbcA are important for LbcA-mediated substrate degradation by CtpA both in vitro and in vivo. This work improves our understanding of the molecular mechanism of the LbcA-CtpA proteolytic system and reveals some striking differences from the arrangements found in some other bacterial CTPs.
We perform a global QCD analysis of unpolarized parton distribution functions (PDFs) in the proton, including new 𝑊+ charm production data from 𝑝𝑝 collisions at the LHC and semi-inclusive pion and kaon production data in lepton-nucleon deep-inelastic scattering, both of which have been suggested for constraining the strange quark PDF. Compared with a baseline global fit that does not include these datasets, the new analysis reduces the uncertainty on the strange quark distribution over the range 0.01 < 𝑥 < 0.3, and provides a consistent description of processes sensitive to strangeness in the proton. Including the new datasets, the ratio of strange to nonstrange sea quark distributions is $R_s = (s + \bar{s})/(\bar{u} +\bar{d})$ $=$ {$0.72^{+0.52}_{−0.34}, 0.46^{+0.30}_{−0.20}, 0.32^{+0.23}_{−0.15}$} for 𝑥 ={$0.01, 0.04, 0.1$} at 𝑄 2 $=$ 4 GeV 2 . The data place more stringent constraints on the strange asymmetry $(s - \bar{s})$, which is found to be consistent with zero in this range.
We perform a simulation-based analysis on the potential of emerging ferroelectric tunnel junctions (FTJs) as a memory device for crossbar arrays. Though FTJs are promising due to their low power switching characteristics compared to other emerging technologies, the greatest challenge for FTJs is the tradeoff between integration density and read performance. Our analysis highlights the need to co-optimize the ferroelectric thickness of the FTJ and read/write voltages to achieve proper functionality at large array sizes. Our analysis shows that FTJ-based crossbar achieves 93% higher sense margin at isoread power of 116 nW (per bit), but this FTJ design comes at a cost of 9.28× higher write power at isowrite time of 250 ns. In response, we study the potential tradeoffs of design points outside the feasible region to understand what device characteristics are desired to overcome such challenges.
Germanium selenide (GeSe) is a van der Waals–bonded layered material with promising optoelectronic properties, which has been experimentally synthesized for two-dimensional (2D) semiconductor applications. In the monolayer, due to reduced dimensionality and, thus, screening environment, perturbations such as the presence of defects have a significant impact on its properties. We apply density functional theory and many-body perturbation theory to understand the electronic and optical properties of GeSe containing a single selenium vacancy in the –2 charge state. We predict that the vacancy results in midgap “trap states” that strongly localize the electron and hole density and lead to sharp, low-energy optical absorption peaks below the predicted pristine optical gap. Analysis of the exciton wave function reveals that the 2D Wannier-Mott exciton of the pristine material monolayer is strongly modified by the presence of the defect. The lowest-energy exciton wave function is highly localized around the defect, with the Bohr radius reduced by a factor of 4 when compared to the pristine material and a dipole moment along the out-of-plane axis due to the defect-induced symmetry breaking. Altogether, these results suggest that the vacancy significantly modifies the excitonic properties of the system, demonstrating the importance of considering defects in the context of material design.
The tensile failure stress of a brittle-material sphere is important to the survival of the sphere itself but is particularly relevant to systems that employ such spheres for applications in which the sphere exterior is subjected to tension during system fabrication or operation. That maximum hoop or tangential stress can be measured by diametrally compressing a C-sphere specimen directly machined from such spheres. This combination of specimen geometry and test method was conceived at Oak Ridge National Laboratory around 20 years ago. Although the C-sphere geometry has been tested with larger spheres, the ease and efficacy of its testing are unknown for submillimeter sphere diameters. To consider this, C-sphere specimens were machined from yttria-stabilized zirconia (YSZ) spheres of ∼500 µm diameter. The maximum (hoop or tangential) tensile stress was determined as a function of compressive failure force and specimen dimensions using finite element analysis. Companion effective area and volume as a function of Weibull modulus were computed. Hoop tensile failure stresses were then measured and fractography conducted to examine flaw type and location and to estimate the mirror constant of this YSZ grade. This novel study demonstrates that hoop tensile failure stress can be determined with sub-millimeter brittle-material spheres.
Whittle index policy is a heuristic to the intractable restless multi-armed bandits (RMAB) problem. Although it is provably asymptotically optimal, finding Whittle indices remains difficult. In this paper, we present Neural-Q-Whittle, a Whittle index based Q-learning algorithm for RMAB with neural network function approximation, which is an example of nonlinear two-timescale stochastic approximation with Q-function values updated on a faster timescale and Whittle indices on a slower timescale. Despite the empirical success of deep Q-learning, the non-asymptotic convergence rate of Neural-Q-Whittle, which couples neural networks with two-timescale Q-learning largely remains unclear. This paper provides a finite-time analysis of Neural-Q-Whittle, where data are generated from a Markov chain, and Q-function is approximated by a ReLU neural network. Our analysis leverages a Lyapunov drift approach to capture the evolution of two coupled parameters, and the nonlinearity in value function approximation further requires us to characterize the approximation error. Combing these provide Neural-Q-Whittle with convergence rate, where is the number of iterations.
As a first step in preparing for the return of samples from the Moon by the Artemis Program, NASA initiated the Apollo Next Generation Sample Analysis Program (ANGSA). ANGSA was designed to function as a low-cost sample return mission and involved the curation and analysis of samples previously returned by the Apollo 17 mission that remained unopened or stored under unique conditions for 50 years. These samples include the lower portion of a double drive tube previously sealed on the lunar surface, the upper portion of that drive tube that had remained unopened, and a variety of Apollo 17 samples that had remained stored at -27 °C for approximately 50 years. ANGSA constitutes the first preliminary examination phase of a lunar “sample return mission” in over 50 years. It also mimics that same phase of an Artemis surface exploration mission, its design included placing samples within the context of local and regional geology through new orbital observations collected since Apollo and additional new “boots-on-the-ground” observations, data synthesis, and interpretations provided by Apollo 17 astronaut Harrison Schmitt. ANGSA used new curation techniques to prepare, document, and allocate these new lunar samples, developed new tools to open and extract gases from their containers, and applied new analytical instrumentation previously unavailable during the Apollo Program to reveal new information about these samples. Most of the 90 scientists, engineers, and curators involved in this mission were not alive during the Apollo Program, and it had been 30 years since the last Apollo core sample was processed in the Apollo curation facility at NASA JSC. There are many firsts associated with ANGSA that have direct relevance to Artemis. ANGSA is the first to open a core sample previously sealed on the surface of the Moon, the first to extract and analyze lunar gases collected in situ, the first to examine a core that penetrated a lunar landslide deposit, and the first to process pristine Apollo samples in a glovebox at -20 °C. All the ANGSA activities have helped to prepare the Artemis generation for what is to come. The timing of this program, the composition of the team, and the preservation of unopened Apollo samples facilitated this generational handoff from Apollo to Artemis that sets up Artemis and the lunar sample science community for additional successes.
We present a kinetic stability analysis of the solar wind electron distribution function consisting of the Maxwellian core and the magnetic-field aligned strahl, a superthermal electron beam propagating away from the sun. We use an electron strahl distribution function obtained as a solution of a weakly collisional drift-kinetic equation, representative of a strahl affected by Coulomb collisions but unadulterated by possible broadening from turbulence. This distribution function is essentially non-Maxwellian and varies with the heliospheric distance. The stability analysis is performed with the Vlasov–Maxwell linear solver leopard. We find that depending on the heliospheric distance, the core-strahl electron distribution becomes unstable with respect to sunward-propagating kinetic-Alfvén, magnetosonic, and whistler modes, in a broad range of propagation angles. The wavenumbers of the unstable modes are close to the ion inertial scales, and the radial distances at which the instabilities first appear are on the order of 1 au. However, we have not detected any instabilities driven by resonant wave interactions with the superthermal strahl electrons. Instead, the observed instabilities are triggered by a relative drift between the electron and ion cores necessary to maintain zero electric current in the solar wind frame (ion frame). Contrary to strahl distributions modelled by shifted Maxwellians, the electron strahl obtained as a solution of the kinetic equation is stable. Our results are consistent with the previous studies based on a more restricted solution for the electron strahl.
Heavy-fermion materials are compounds in which localized f orbitals hybridize with delocalized d ones, leading to quasiparticles with large renormalized masses. The presence of strongly correlated f electrons at the Fermi level may also lead to long-range order, such as magnetism, or unconventional superconductivity. From a theoretical point of view, the “standard model” for heavy-fermion compounds is the periodic Anderson model (PAM). Despite being extensively scrutinized, its thermodynamic properties in three-dimensional (3D) lattices have not been carefully addressed by unbiased methodologies. Here we investigate the 3D PAM employing state-of-the-art finite-temperature auxiliary-field quantum Monte Carlo simulations. Furthermore, we present the behavior of the kinetic energy, the entropy, the specific heat, and the double occupancy as functions of the temperature and the hybridization strength. From these quantities, and by the analysis of the spin-spin correlation functions, we investigate the occurrence of magnetic phase transitions at finite temperatures, and we determine the phase diagram of the model, including the behavior of the Néel temperature as a function of the external parameters.
In previous years, SGTech funded enhancements to the isotopics routine in the software called Gamma Detector Response and Analysis Software-Detector Response Function (GADRAS-DRF), including the addition of peak fit customization capabilities. A project was also funded that focused on implementing a peak-based model fitting routine, allowing model fitting to be performed without dependence on export-controlled cross-sections. In FY25 significant improvements were made to the custom peak fitting interface, accompanied by several validation studies within GADRAS-DRF. These studies encompassed IsotopeID performance, distributed source analysis, isotopics validation, and activity estimation. Additionally, the peak-only model fitting option was validated using an HPGe measurement of a rotating drum with line sources.