Search NASA⌕ Search

SEARCH · Search NASA

Results for “function approximation”

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 505 records · Page 28

Active learning for SNAP interatomic potentials via Bayesian predictive uncertainty

Bayesian inference with a simple Gaussian error model is used to efficiently compute prediction variances for energies, forces, and stresses in the linear SNAP interatomic potential. Here, the prediction variance is shown to have a strong correlation with the absolute error over approximately 24 orders of magnitude. Using this prediction variance, an active learning algorithm is constructed to iteratively train a potential by selecting the structures with the most uncertain properties from a pool of candidate structures. The relative importance of the energy, force, and stress errors in the objective function is shown to have a strong impact upon the trajectory of their respective net error metrics when running the active learning algorithm. Batched training of different batch sizes is also tested against singular structure updates, and it is found that batches can be used to significantly reduce the number of retraining steps required with only minor impact on the active learning trajectory.

97 MATHEMATICS AND COMPUTING↗

The separatrix electron density in JET, ASDEX upgrade and alcator C-Mod H-mode plasmas: A common evaluation procedure and correlation with engineering parameters

The separatrix electron density is an important parameter for core-edge scenario integration in tokamak devices, as it influences plasma confinement, divertor detachment and disruption avoidance. This quantity has been measured in H-mode discharges on JET, ASDEX Upgrade and Alcator C-Mod by applying the same fitting function to Thomson scattering measurements, and by employing the same analysis technique based on scrape-off layer power balance. To estimate the power crossing the separatrix, the inter-ELM time derivative of the plasma energy d W /d t has been experimentally evaluated and found to be approximately a constant fraction of the absorbed heating power. Correlations between n e,sep and engineering parameters have been investigated, revealing that n e,sep scales with the divertor neutral pressure p 0,div in a similar manner across all devices. Additionally, when n e,sep is normalized to the obtained n e,sep dependency, no clear correlation with the plasma current is found. These observations are in agreement with the 2-point model, which suggests that the upstream separatrix density is mainly set by the recycling at the divertor target.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Crystal structure of valbenazine, C 24 H 38 N 2 O 4

The crystal structure of valbenazine has been solved and refined using synchrotron X-ray powder diffraction data and optimized using density functional theory techniques. Valbenazine crystallizes in space groupP2 1 2 1 2 1 (#19) witha= 5.260267(17),b= 17.77028(7),c= 26.16427(9) Å,V= 2445.742(11) Å 3 , andZ= 4 at 295 K. The crystal structure consists of discrete molecules and the mean plane of the molecules is approximately (8,−2,15). There are no obvious strong intermolecular interactions. There is only one weak classical hydrogen bond in the structure, from the amino group to the ether oxygen atom. Two intramolecular and one intermolecular C–H⋯O hydrogen bonds also contribute to the lattice energy. The powder pattern has been submitted to ICDD for inclusion in the Powder Diffraction File™ (PDF®)

Materials Science↗

Thermodynamics of Tritium Trapping by Point Defects in Intermetallic Al 12 (TM) 2.35 Aluminide Coating Phase

Density functional theory simulations have been carried out to investigate the potential for tritium trapping by metal vacancies in intermetallic Al 12 (TM) 2.35 phase (TM = Fe, Cr, and Ni) as function of temperature and tritium partial pressure. It was found that tritium could be favorably trapped by Fe and Ni vacancies and not favorably trapped by Al and Cr vacancies. However, due to the presence of partially occupied Al sites in bulk Al 12 (TM) 2.35 , leading to the approximate number of ~255 Al atoms in the unit cell, 86 sites were found energetically favorable to the creation of an Al vacancy. While adding a tritium atom in an Al vacancy is not energetically favorable, the tritiated defect still has a negative Gibbs free energy because the energy gain for creating an Al vacancy overcome the energy cost of adding the tritium species. Based on the calculated Gibbs free energy, the first tritiation of a metal vacancy, at conditions relevant to in-reactor operations, should be more favorable for Al, followed Fe, Ni, and Cr vacancies. By comparing the behavior of tritium in Al 12 (TM) 2.35 with previously studied Fe-Al coating phases (i.e., FeNiAl 5 , Fe 4 Al 13 , and Fe 2 Al 5.6 ), we found that there is a correlation between interstitial tritium solubility and the potential for vacancy trapping. The current trend suggests that if the insertion of an interstitial tritium cost more than 0.3 eV, then trapping by metal vacancies should be preferred. By combining the simulations results obtained to date, we noticed different trapping mechanisms of tritium in the Al coating. Tritium is mostly trapped by Fe and Ni vacancies in the outer Fe-Al coating phase Al 12 (TM) 2.35 while tritium should be preferentially trapped by Al and Fe vacancies for the inner Fe-Al coating phases (FeNiAl 5 , Fe 4 Al 13 , Fe 2 Al 5.6 ). Altogether, these studies show that tritium interacts differently with the various Fe-Al aluminide phases, they also suggest that tritium trapping and retention could be more efficient if metal defects are present and if the solubility of interstitial tritium in the different phases is low.

36 MATERIALS SCIENCE↗

Effect of dilute Rh on oxygen dissociation, spillover, and the oxidation of Cu across many orders of magnitude pressure

Knowledge of how trace amounts of more reactive metals influence the oxidation rate and mechanism of Cu surfaces is essential for developing strategies to optimize the performance of Cu-based catalysts. We find that the addition of 1% Rh to Cu(111) increases the initial O 2 dissociation rate by approximately 9-fold. CO poisoning experiments reveal that single Rh atoms activate O 2 and facilitate the spillover of atomic oxygen to Cu sites. Scanning tunneling microscopy (STM) and in situ X-ray photoelectron spectroscopy (XPS) support this mechanism, showing enhanced surface oxygen near Rh atoms. Here, a density functional theory (DFT)-based model demonstrates that Rh binds the O 2 precursor 0.15 eV more strongly than Cu(111) and lowers the O 2 dissociation barrier by 0.02 eV. Both single-crystal and nanoparticle experiments show that at low oxygen pressures, Rh enhances Cu oxidation, whereas at higher pressures, it inhibits deeper oxidation, as evidenced by in situ ultraviolet-visible (UV-vis) spectra.

36 MATERIALS SCIENCE↗

Pentaquarks made of light quarks and their admixture to baryons

This paper is a continuation of our studies of multiquark hadrons. The antisymmetrization of their wave functions required by Fermi statistics is nontrivial, as it mixes orbital, color, spin, and flavor structures. In our previous papers we developed a method to find them based on the representations of the permutation group, and derived the explicit wave functions for baryons excited to the first and second shells (L = 1, 2), tetraquarks $qq$$\overline{q}$$\overline{q}$ and hexaquarks (6q). Now we apply it to light pentaquarks ($qqq$$\overline{q}$), in the S- and P-shells (L = 0, 1). Using Jacobi coordinates, one can use the hyperdistance approximation in 12-dimensional space. We further address the issue of “unquenching” of baryons, by considering their mixing with pentaquarks, via two channels, through the addition of σ-like or π-like $\overline{q}$$q$ pairs. This mixing is central for understanding of the observed flavor asymmetry of the antiquark sea, the amount of orbital motion issue as well as other nucleon properties.

Baryons↗

Quasiperiodic potassium adlayer on decagonal Al–Ni–Co quasicrystal

Quasiperiodicity in free-electron-like metals is a subject of significant interest within the scientific community. In this work, utilizing scanning tunneling microscopy (STM), low energy electron diffraction (LEED), and density functional theory (DFT), we demonstrate the formation of a quasiperiodic potassium monolayer on the tenfold surface of decagonal Al–Ni–Co quasicrystal. A dispersed growth comprising of isolated K adatoms is observed at sub-monolayer coverage, which coalesce with increasing coverage and forms pentagonal and decagonal quasiperiodic motifs. LEED demonstrates distinct rings of spots displaying decagonal symmetry. Furthermore, our DFT calculation using the W-approximant surface that closely resembles d-Al–Ni–Co shows that sizable adsorbate-substrate interaction makes the potassium adatoms bind to the quasiperiodically dispersed favorable adsorption sites resulting in the formation of quasiperiodic potassium monolayer. Notably, the experimental motifs obtained from STM measurements align remarkably well with the DFT predictions, underscoring the intricate relationship between the electronic structure of the substrate that drives the quasiperiodic growth of the potassium adlayer.

alkali metal↗

An efficient explicit implementation of a near-optimal quantum algorithm for simulating linear dissipative differential equations

We propose an efficient block-encoding technique for the implementation of the Linear Combination of Hamiltonian Simulations (LCHS) for simulating dissipative initial-value problems. This algorithm approximates a target nonunitary operator as a weighted sum of Hamiltonian evolutions, thereby emulating a dissipative problem by mixing various time scales. We introduce an efficient encoding of the LCHS into a quantum circuit based on a simple coordinate transformation that turns the dependence on the summation index into a trigonometric function. Classically, this method is equivalent to the use of a highly accurate Fejér-Clenshaw-Curtis quadrature formula. Quantumly, this significantly simplifies block-encoding of a dissipative problem and allows one to perform an exponential number of Hamiltonian simulations by a single Quantum Signal Processing (QSP) circuit. The resulting LCHS circuit has high success probability and the selector scales logarithmically with the number of terms in the LCHS sum and linearly with time. Careful analysis of error convergence proves that this method is more efficient than other LCHS circuits that have recently appeared in the literature. We verify the quantum circuit and its scaling by simulating it on a digital emulator of fault-tolerant quantum computers and, as a test problem, solve the advection-diffusion equation. The proposed algorithm can be used for simulating a wide class of nonunitary initial-value problems including the Liouville equation with added dissipation and linear embeddings of nonlinear systems, such as the Koopman-von Neumann and Carleman embeddings.

Novikau, I [Lawrence Livermore National Laboratory↗

nf-core/proteinfamilies: a scalable pipeline for the generation of protein families

The growth of metagenomics-derived amino acid sequence data has transformed our understanding of protein function, microbial diversity, and evolutionary relationships. However, the vast majority of these proteins remain functionally uncharacterized. Grouping the millions of such uncharacterized sequences with the few experimentally characterized ones allows the transfer of annotations, while the inspection of conserved residues with multiple sequence alignments can provide clues to function, even in the absence of existing functional information. To address the challenges associated with this data surge and the need to group sequences, we present a scalable, open-source, parametrizable Nextflow pipeline (nf-core/proteinfamilies) that generates nascent protein families or assigns new proteins to existing families. The computational benchmarks demonstrated that resource usage scales approximately linearly with input size, and the biological benchmarks showed that the generated protein families closely resemble manually curated families in widely used databases.

Nextflow↗

Final Report for FE0032098: Improving the cost-effectiveness of algal CO2 utilization by synergistic integration with power plant and wastewater treatment operations

The overall goal of this project was to demonstrate an engineering-scale open raceway pond algae cultivation system (approximately 180 m²), including the integration of technologies that utilized carbon dioxide (CO₂) from a coal-fired power plant and wastewater-derived nutrient inputs for cost-effective and environmentally friendly biomass production. The key advantages associated with the innovative algae cultivation system and its integration with wastewater treatment functions, as described herein, had been demonstrated in previous bench- and pilot-scale work by the project team partners. This project combined those approaches to maximize practical benefits and available synergies, resulting in a significant reduction in the net cost of producing algal biomass products. The primary target algal species for the project was Spirulina, which served as a high-protein content ingredient for food and animal feed. Spirulina was selected because it had already been approved by the FDA, was in use as a food ingredient, and commanded prices of up to $30/kg. It had a typical protein content of 50–75%, comparable to other high-protein concentrates, featured high digestibility without requiring pretreatment, and offered a high conversion ratio in animal feed applications. In addition, Spirulina had a relatively high content of the blue pigment phycocyanin, which could be extracted as a high-value co-product prior to using the remaining biomass for nutritional purposes.

Schideman, Lance [University of Illinois]↗

Data-driven reduced-order models for port-Hamiltonian systems with operator inference

Hamiltonian operator inference has been developed in Sharma et al. (2022) to learn structure-preserving reduced-order models (ROMs) for Hamiltonian systems. The method constructs a low-dimensional model using only data and knowledge of the functional form of the Hamiltonian. The resulting ROMs preserve the intrinsic structure of the system, ensuring that the mechanical and physical properties of the system are maintained. In this work, we extend this approach to port-Hamiltonian systems, which generalize Hamiltonian systems by including energy dissipation, external input, and output. Based on snapshots of the system’s state and output, together with the information about the functional form of the Hamiltonian, reduced operators are inferred through optimization and are then used to construct data-driven ROMs. To further alleviate the complexity of evaluating nonlinear terms in the ROMs, a hyper-reduction method via discrete empirical interpolation is applied. Accordingly, we derive error estimates for the ROM approximations of the state and output. Lastly, we demonstrate the structure preservation, as well as the accuracy of the proposed port-Hamiltonian operator inference framework, through numerical experiments on a linear mass–spring-damper problem and a nonlinear Toda lattice problem.

97 MATHEMATICS AND COMPUTING↗

Completely Multipolar Model for Many-Body Water–Ion and Ion–Ion Interactions

This work constructs an advanced force field, the Completely Multipolar Model (CMM), to quantitatively reproduce each term of an energy decomposition analysis (EDA) for aqueous solvated alkali metal cations and halide anions and their ion pairings. We find that all individual EDA terms remain well-approximated in the CMM for ion-water and ion-ion interactions, except for polarization, which shows errors due to the partial covalency of ion interactions near their equilibrium. We quantify the onset of the dative bonding regime by examining the change in molecular polarizability and Mayer bond indices as a function of distance, showing that partial covalency manifests by breaking the symmetry of atomic polarizabilities while strongly damping them at short-range. This motivates an environment-dependent atomic polarizability parameter that depends on the strength of the local electric field experienced by the ions to account for strong damping, with anisotropy introduced by atomic multipoles. The resulting CMM model for ions provides accurate dimer surfaces and three-body polarization and charge transfer compared to EDA, and shows excellent performance on various ion benchmarks including vibrational frequencies and cluster geometries.

Heindel, Joseph P↗

Estimating Electron Temperature and Density Using Van Allen Probe Data: Typical Behavior of Energetic Electrons in the Inner Magnetosphere

Abstract The Earth's inner magnetosphere contains multiple electron populations influenced by different factors. The cold electrons of the plasmasphere, warm plasma that contributes to the ring current, and the relativistic plasma of the radiation belts often seem to behave independently. Using omni‐directional flux and energy measurements from the HOPE and Magnetic Electron Ion Spectrometer instruments aboard the Van Allen Probes, we provide a detailed density and temperature description of the inner magnetosphere, offering a comprehensive statistical analysis of the entire Van Allen Probe era. While number density and temperature data at geosynchronous orbit are available, this study focuses on the warm plasma in the inner magnetosphere . Values of density and temperature are extracted by fitting energy and phase space density to obtain the distribution function. The fitted distributions are related to the zeroth and second moments to estimate the number density and temperature. Analysis has indicated that a two Maxwellian fit is sufficient over a wide range of and that there are two independent plasma populations. The more energetic population has a median number density of approximately and a temperature of around 130 keV, with a temperature peak observed between L * = 4 and L * = 4.5. This population is relatively uniform in magnetic local time (MLT). In contrast, the less energetic warm electron population has a median number density of about and a temperature of 7.4 keV. Strong statistical trends in density and temperature across both L * and MLT are presented, along with potential sources driving these variations.

58 GEOSCIENCES↗

Unveiling Hidden Lyman Alpha Emitters in the DESI DR1 Data

We present an automatic method based on machine-learning convolutional neural network (CNN) architecture to detect Lyman alpha emitters (LAE) hidden in the Data Release 1 spectroscopic dataset of the Dark Energy Spectroscopic Instrument (DESI). Those LAEs mostly have incorrect redshift estimations because the current DESI pipeline is not designed to detect and measure the redshifts of galaxies at $z>2$. To uncover those sources, we first visually inspect thousands of DESI spectra and construct a sample, consisting of both LAEs and non-LAEs, for training and testing the CNN-based model to (1) detect LAEs in DESI spectra and (2) determine their Ly$α$ redshifts. The final model yields $95.2\%$ purity and $95.9\%$ completeness for detecting LAEs. We apply this model to approximately $2\times10^{6}$ spectra of sources targeted as emission-line galaxies and detect 19,685 LAEs from $z\sim2$ to $3.5$ within 12 minutes with a single GPU, illustrating the high efficiency of this model for identifying LAEs. The detected LAEs are mostly at the bright end of the luminosity function with Ly$α$ luminosity $L_{\rm Lyα} \gtrsim 10^{43}$ erg/s. The high signal-to-noise composite spectrum of the detected LAEs further shows various spectral features, including P-Cygni profiles of metal lines and MgII emission lines, possible indicators of Lyman continuum escape fraction, revealing the rich astrophysical information in this LAE sample. Finally, this sample can be used to train and validate the pipelines for redshift determination of LAEs for the preparation of the DESI-II survey.

Chan, Jui-Kuan [Taiwan, Natl. Taiwan U.] (ORCID:00↗

Beyond Single-Reference Fixed-Node Approximation in Ab Initio Diffusion Monte Carlo Using Antisymmetrized Geminal Power Applied to Systems with Hundreds of Electrons

Diffusion Monte Carlo (DMC) is an exact technique to project out the ground state (GS) of a Hamiltonian. Since the GS is always bosonic, in Fermionic systems, the projection needs to be carried out while imposing antisymmetric constraints, which is a nondeterministic polynomial hard problem. In practice, therefore, the application of DMC on electronic structure problems is made by employing the fixed-node (FN) approximation, consisting of performing DMC with the constraint of having a fixed, predefined nodal surface. How do we get the nodal surface? The typical approach, applied in systems having up to hundreds or even thousands of electrons, is to obtain the nodal surface from a preliminary mean-field approach (typically, a density functional theory calculation) used to obtain a single Slater determinant. This is known as single reference. In this paper, we propose a new approach, applicable to systems as large as the C 60 fullerene, which improves the nodes by going beyond the single reference. In practice, we employ an implicitly multireference ansatz (antisymmetrized geminal power wave function constraint with molecular orbitals), initialized on the preliminary mean-field approach, which is relaxed by optimizing a few parameters of the wave function determining the nodal surface by minimizing the FN-DMC energy. We highlight the improvements of the proposed approach over the standard single-reference method on several examples and, where feasible, the computational gain over the standard multireference ansatz, which makes the methods applicable to large systems. We also show that physical properties relying on relative energies, such as binding energies, are affordable and reliable within the proposed scheme.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Thermophysical Properties of NaCl–UCl 3 –PuCl 3 Molten Salts: A Combined Computational and Experimental Study

Actinide-bearing molten salts for use as fuels are an essential part of next generation molten salt reactors. Yet, numerous multicomponent salt mixtures are underdeveloped or have not been investigated. Here, this study, based on a combination of experimental and modeling techniques, is dedicated to determining and understanding a variety of properties of the ternary system of NaCl–UCl 3 –PuCl 3 , which represents a scenario for burnup of NaCl–UCl 3 fuel, at two compositions (∼10 and 5 mol % PuCl 3 in eutectic NaCl–UCl 3 pseudobinary) and a range of temperatures. Evaluation of the heat flow and mass loss data showed the 0.61NaCl–0.30UCl 3 –0.09PuCl 3 salt had a melting temperature of 551 ± 5 °C. Two additional thermal effects were observed occurring at approximately 410 and 494 °C. The transition occurring at 410 °C may be due to the presence of oxide in the salt. Extrapolation of thermodynamic data indicates the transition occurring at 494 °C is due to the formation of a liquid phase. Experimental testing determined the density of this system is a linear function of temperature and can be represented by the equation ρ = 4.014–0.0010T(°C), R 2 = 0.992. Additionally, by using atomistic modeling, we found that increasing the PuCl 3 content from 5 to 10 mol % led to the formation of larger Pu 3+ clusters and slower transport of ions.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Breaking the Million-Electron and 1 EFLOP/s Barriers: Biomolecular-Scale Ab Initio Molecular Dynamics Using MP2 Potentials

The accurate simulation of complex biochemical phenomena has historically been hampered by the computational requirements of high-fidelity molecular-modeling techniques. Quantum mechanical methods, such as ab initio wave-function (WF) theory, deliver the desired accuracy, but have impractical scaling for modeling biosystems with thousands of atoms. Combining molecular fragmentation with MP2 perturbation theory, this study presents an innovative approach that enables biomolecular-scale ab initio molecular dynamics (AIMD) simulations at WF theory level. Leveraging the resolution-of-the-identity approximation for Hartree-Fock and MP2 gradients, our approach eliminates computationally intensive four-center integrals and their gradients, while achieving near-peak performance on modern GPU architectures. The introduction of asynchronous time steps minimizes time step latency, overlapping computational phases and effectively mitigating load imbalances. Utilizing up to 9,400 nodes of Frontier and achieving 59% (1006.7 PFLOP/s) of its double-precision floating-point peak, our method enables us to break the million-electron and 1EFLOP/s barriers for AIMD simulations with quantum accuracy.

Kurzak, Jakub↗

Neptunium (V) Solubility in WIPP Brine: Presence and Absence of WIPP-Relevant Organic and Inorganic Ligands

The oxidation state and solubility of neptunium (Np) in 5 M NaCl and high ionic strength synthetic WIPP (Waste Isolation Pilot Plant) brines were systematically investigated as a function of pCH+ (7-11) in the presence and absence of WIPP-relevant organic ligands, including EDTA (Ethylenediaminetetraacetic acid), oxalate, citrate, and acetate, as well as inorganic ligands such as borate and carbonate at a temperature of 23 ± 2°C. This study was conducted through long-term batch solubility experiments spanning approximately 950 days, using an undersaturation approach. A comprehensive set of experimental and spectroscopic techniques, including UV-VIS-NIR spectroscopy and Np L III -edge X-ray Absorption Spectroscopy (XAS), was employed to identify the solubility-controlling Np solid phases and the predominant aqueous Np species present in the samples.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗