Search NASA⌕ Search

SEARCH · Search NASA

Results for “Functional expansion”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Machine learned potential for high-throughput phonon calculations of metal—organic frameworks

Metal–organic frameworks (MOFs) are highly porous and versatile materials studied extensively for applications such as carbon capture and water harvesting. However, computing phonon-mediated properties in MOFs, like thermal expansion and mechanical stability, remains challenging due to the large number of atoms per unit cell, making traditional Density Functional Theory (DFT) methods impractical for high-throughput screening. Recent advances in machine learning potentials have led to foundation atomistic models, such as MACE-MP-0, that accurately predict equilibrium structures but struggle with phonon properties of MOFs. In this work, we developed a workflow for computing phonons in MOFs within the quasi-harmonic approximation with a fine-tuned MACE model, MACE-MP-MOF0. The model was trained on a curated dataset of 127 representative and diverse MOFs. The fine-tuned MACE-MP-MOF0 improves the accuracy of phonon density of states and corrects the imaginary phonon modes of MACE-MP-0, enabling high-throughput phonon calculations with state-of-the-art precision. The model successfully predicts thermal expansion and bulk moduli in agreement with DFT and experimental data for several well-known MOFs. These results highlight the potential of MACE-MP-MOF0 in guiding MOF design for applications in energy storage and thermoelectrics.

Elena, Alin Marin↗
DOE OSTI2022

Symmetry-correct bonding in density functional theory calculations for delta phase Pu

The long-held conclusion that magnetic order in delta phase Pu makes the structure mechanically unstable in density functional theory calculations is incorrect if the magnetic order is allowed to be non-collinear. The non-collinear 3Q spin structure is intrinsically cubic and makes all structurally equivalent bonds equivalent in their bonding character. Applied to density functional theory calculations on delta phase Pu, the 3Q spin structure results in elastic constants and phonons with the correct symmetry. Finally, the calculated phonon dispersion agrees better with experiment than previous calculations, and the calculated thermal expansion shows the unique behavior seen experimentally: it is negative.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Developing Machine Learning Interatomic Potential for Fe-Cr-Ni Alloys

Accurate prediction of creep and fatigue behavior of stainless steel at elevated temperatures in hydrogen environment requires fundamental understanding of alloy-hydrogen interaction at cross-scale including bulk lattice and key defects such as vacancies, grain boundaries, surfaces, stacking faults, dislocations, and precipitates. This project aims to predict creep behavior of 347H stainless steel with H using machine learning interatomic potentials based on first-principles density functional theory simulations. The Moment Tensor Potentials platform is adopted for this work since it demonstrates a fine balance between model accuracy and computational efficiency. The potential is well trained based on large amount of high-fidelity density functional theory calculations. The validation is carried out by comparing various important properties including short range order, coefficient of thermal expansion, elastic properties, stacking fault energy, grain boundary energy, and surface energy. This work lays the foundation for reliable atomistic simulation of high temperature hydrogen attack of stainless steel.

density functional theory (DFT)↗
DOE OSTI2023

A prediction of the thermodynamic, thermophysical, and mechanical properties of CrTaO 4 from first principles

It is reported that the self-forming CrTaO 4 oxide scale can protect refractory high-entropy alloys from oxidation, superior to Cr 2 O 3 . In this paper, the phase stability, mechanical, and thermal properties of three polymorphous phases of CrTaO 4 are systematically investigated from first-principles density functional theory calculations. The mechanical properties predicted using the strain–energy methods indicated that all three phases are mechanically stable. The temperature dependence of elastic constants and polycrystalline moduli of three phases demonstrated the thermal softening as temperature increase. The Helmholtz-free energies as a function of volume and temperature are derived from phonon dispersions within the quasi-harmonic approximation at six strained volumes. The calculated apparent bulk coefficients of thermal expansion of these three phases are evaluated, the highest value approximately 13.4× 10 –6 K –1 within a temperature range of 500–2000 K for the rutile I4 1 md phase. The lattice thermal conductivity calculated by the Debye–Callaway model suggested that the rutile type I4 1 md phase has the lowest value of approximately 2.1 W/m/K at 1800 K. The other two phases, C2/m and P2/c, exhibit higher values due to relatively lower Grüneisen parameters and larger phonon velocities. The melting point of CrTaO 4 is predicted to be between 1975 and 2449 K using ab initio molecular dynamics simulations. Finally, this work provides a comprehensive theoretical understanding of the thermodynamic, mechanical, and thermal properties for the new material CrTaO 4 and serves as an example of a viable computational design strategy for improved oxidation resistance of refractory alloys at high temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Constraints on standard model effective field theory for a Higgs boson produced in association with W or Z bosons in the H → $\textrm{b}\overline{\textrm{b}}$ decay channel in proton-proton collisions at $\sqrt{s}$ = 13 TeV

A standard model effective field theory (SMEFT) analysis with dimension-six operators probing nonresonant new physics effects is performed in the Higgs-strahlung process, where the Higgs boson is produced in association with a W or Z boson, in proton-proton collisions at a center-of-mass energy of 13 TeV. The final states in which the W or Z boson decays leptonically and the Higgs boson decays to a pair of bottom quarks are considered. The analyzed data were collected by the CMS experiment between 2016 and 2018 and correspond to an integrated luminosity of 138 fb −1 . An approach designed to simultaneously optimize the sensitivity to Wilson coefficients of multiple SMEFT operators is employed. Likelihood scans as functions of the Wilson coefficients that carry SMEFT sensitivity in this final state are performed for different expansions in SMEFT. The results are consistent with the predictions of the standard model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Heptagon functions and seven-gluon amplitudes in multi-Regge kinematics

We compute all 2 → 5 gluon scattering amplitudes in planar $ \mathcal{N} $ N = 4 super-Yang-Mills theory in the multi-Regge limit that is sensitive to the non-trivial (“long”) Regge cut. We provide the amplitudes through four loops and to all logarithmic accuracy at leading power, in terms of single-valued multiple polylogarithms of two variables. To obtain these results, we leverage the function-level results for the amplitudes in the Steinmann cluster bootstrap. To high powers in the series expansion in the two variables, our results agree with the recently conjectured all-order central emission vertex used in the Fourier-Mellin representation of amplitudes in multi-Regge kinematics. Our results therefore provide a resummation of the Fourier-Mellin residues into single-valued polylogarithms, and constitute an important cross-check between the bootstrap approach and the all-orders multi-Regge proposal.

1/N expansion↗
DOE OSTI2021

A unified system for Coleman–De Luccia transitions

Highlights: • Monotonicity controls whether radius of CDL bounce is negative mode direction. • Decays from de Sitter and Anti-de Sitter analyzed in single geometric framework. • de Sitter to Anti-de Sitter instantons come in two types, similar to type-A vs type-B. We show that the presence or absence of negative modes associated with the expansion or contraction of Coleman–De Luccia instantons is controlled by the monotonicity properties of a single function parameterized by the tension and two vacuum energies. This approach simplifies and unifies certain aspects of phase transitions between two de Sitter vacua, between two Anti-de Sitter, and between one of each. This strategy may serve as a guide for identifying negative modes in more complicated gravitational theories, for instance supergravity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

The cold high-pressure approach to hydrogen delivery

In view of the very expensive and wasteful nature of today's approaches to H 2 delivery, in this work we explore the possibility of transporting cold (200 K) high pressure (875 bar) H 2 in thermally insulated trailers and dispensing H 2 directly from the trailer, with the potential to eliminate station compressor, cascade, and refrigerator, leading to major reductions in station complexity, maintenance, electricity consumption, and cost, while improving functionality by enabling essentially unlimited back to back refuels, and improving safety due to reduced H 2 expansion energy at low temperature. Detailed techno-economic analysis shows promise for substantial delivery cost reductions through cold high pressure H 2 dispensed directly from the trailer. Results indicate that: (1) Terminal operations for cold high pressure H 2 delivery are $\$$0.32/kg H 2 more expensive than for 350 bar compressed gas delivery (today's lowest cost H 2 delivery technology) due to higher level of pressurization (to 1000 bar) and chilling needs (to 165 K). (2) Trailer cost drops slightly ($\$$0.73 vs. $\$$0.81/kg H 2 for a 350 bar trailer) due to increased capacity (1035 kg H 2 delivered vs. 700 kg) compensating for increased capital cost ($\$$906,900 for cold high pressure H 2 vs. $\$$634,000 for 350 bar trailer). (3) Cold hydrogen delivery presents major advantages in fueling station cost, reduced from $\$$1.27 to $\$$0.46/kg H 2 due to elimination of major system components: compressor, cascade, and chiller. (4) Total compression cost (terminal + station) drops from $\$$0.92/kg H 2 ($\$$0.32 terminal and $\$$0.60 station) for 350 bar trailers to $\$$0.55/kg H 2 (all at the terminal) for cold high pressure H 2 . (5) Elimination of small-scale station compressors is the main contributor to reduced delivery cost due to their inefficiency, capital expense, and maintenance needs. In summary, total delivery cost reduction vs. 350 bar trailer equals $\$$0.58/kg H 2 (from $\$$2.96 to 2.38/kg H 2 ), equivalent to 24% of the total delivery cost. This large cost advantage will improve the economics of H 2 vehicles facilitating the transition to a future of zero emission transportation.

08 HYDROGEN↗
DOE OSTI2021

Using Computationally-Determined Properties for Machine Learning Prediction of Self-Diffusion Coefficients in Pure Liquids

The ability to predict transport properties of liquids quickly and accurately will greatly improve our understanding of fluid properties both in bulk and complex mixtures, as well as in confined environments. Such information could then be used in the design of materials and processes for applications ranging from energy production and storage to manufacturing processes. As a first step, we consider the use of machine learning (ML) methods to predict the diffusion properties of pure liquids. Recent results have shown that Artificial Neural Networks (ANNs) can effectively predict the diffusion of pure compounds based on the use of experimental properties as the model inputs. In the current study, a similar ANN approach is applied to modeling diffusion of pure liquids using fluid properties obtained exclusively from molecular simulations. A diverse set of 102 pure liquids is considered, ranging from small polar molecules (e.g., water) to large nonpolar molecules (e.g., octane). Self-diffusion coefficients were obtained from classical molecular dynamics (MD) simulations. Since nearly all the molecules are organic compounds, a general set of force field parameters for organic molecules was used. The MD methods are validated by comparing physical and thermodynamic properties with experiment. Computational input features for the ANN include physical properties obtained from the MD simulations as well as molecular properties from quantum calculations of individual molecules. Furthermore, fluid properties describing the local liquid structure were obtained from center of mass radial distribution functions (COM-RDFs). Feature sensitivity analysis revealed that isothermal compressibility, heat of vaporization, and the thermal expansion coefficient were the most impactful properties used as input for the ANN model to predict the MD simulated self-diffusion coefficients. The MD-based ANN successfully predicts the MD self-diffusion coefficients with only a subset (2 to 3) of the available computationally determined input features required. A separate ANN model was developed using literature experimental self-diffusion coefficients as model targets. Although this second ML model was not as successful due to a limited number of data points, a good correlation is still observed between experimental and ML predicted self-diffusion coefficients.

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

A multi-ancestry GWAS of Fuchs corneal dystrophy highlights the contributions of laminins, collagen, and endothelial cell regulation

Fuchs endothelial corneal dystrophy (FECD) is a leading indication for corneal transplantation, but its molecular etiology remains poorly understood. We performed genome-wide association studies (GWAS) of FECD in the Million Veteran Program followed by multi-ancestry meta-analysis with the previous largest FECD GWAS, for a total of 3970 cases and 333,794 controls. We confirm the previous four loci, and identify eight novel loci: SSBP3, THSD7A, LAMB1, PIDD1, RORA, HS3ST3B1, LAMA5, and COL18A1. We further confirm the TCF4 locus in GWAS for admixed African and Hispanic/Latino ancestries and show an enrichment of European-ancestry haplotypes at TCF4 in FECD cases. Among the novel associations are low frequency missense variants in laminin genes LAMA5 and LAMB1 which, together with previously reported LAMC1, form laminin-511 (LM511). AlphaFold 2 protein modeling, validated through homology, suggests that mutations at LAMA5 and LAMB1 may destabilize LM511 by altering inter-domain interactions or extracellular matrix binding. Finally, phenome-wide association scans and colocalization analyses suggest that the TCF4 CTG18.1 trinucleotide repeat expansion leads to dysregulation of ion transport in the corneal endothelium and has pleiotropic effects on renal function.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Fermions, quantum gravity, and holography in two dimensions

We study a model comprising N flavors of Kähler Dirac fermion propagating on a triangulated two-dimensional disk which is constrained to have a negative average bulk curvature. Dirichlet boundary conditions are chosen for the fermions. Quantum fluctuations of the geometry are included by summing over all possible triangulations consistent with these constraints. We show in the limit N → ∞ that the partition function is dominated by a regular triangulation of two-dimensional hyperbolic space. We use strong coupling expansions and Monte Carlo simulation to show that in this limit boundary correlators of the fermions have a power law dependence on boundary separation as one expects from holography. However, we argue that this behavior breaks down for any finite number of massive fields in the thermodynamic limit and quantum fluctuations of the bulk geometry drive the theory into a nonholographic phase. In contrast, for massless fermions, we find evidence that the boundary is conformal even for finite N . This is consistent with theoretical results in quantum Liouville theory. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2024

Structure Factors for Hot Neutron Matter from Ab Initio Lattice Simulations with High-Fidelity Chiral Interactions

We present the first ab initio lattice calculations of spin and density correlations in hot neutron matter using high-fidelity interactions at next-to-next-to-next-to-leading order in chiral effective field theory. These correlations have a large impact on neutrino heating and shock revival in core-collapse supernovae and are encapsulated in functions called structure factors. Unfortunately, calculations of structure factors using high-fidelity chiral interactions were well out of reach using existing computational methods. In this Letter, we solve the problem using a computational approach called the rank-one operator (RO) method. The RO method is a general technique with broad applications to simulations of fermionic many-body systems. It solves the problem of exponential scaling of computational effort when using perturbation theory for higher-body operators and higher-order corrections. Using the RO method, we compute the vector and axial static structure factors for hot neutron matter as a function of temperature and density. Here, the ab initio lattice results are in good agreement with virial expansion calculations at low densities but are more reliable at higher densities. Random phase approximation codes used to estimate neutrino opacity in core-collapse supernovae simulations can now be calibrated with ab initio lattice calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Performance of Microreactor Test Article with Embedded Sensors During Testing in The Single Primary Heat Extraction and Removal Emulator

The nuclear industry is pursuing microreactors that can be factory assembled and deployed to remote regions for reliable power generation. One class of microreactors uses a monolithic metal core block coupled to heat pipes that use passive heat flow, increasing the surface area for heat transfer without requiring active coolant flow through the reactor core. Additional experimental testing is required to understand the heat rejection limitations of heat pipes and thermal stresses in the monolithic core block due to significant temperature gradients. This report describes the initial characterization and testing of a stainless-steel test article that was fabricated with embedded sensors to measure heat pipe performance limits, as well as spatially distributed temperatures and strains during electrically heated thermal testing to simulate nuclear heating. The electrically heated testing was performed in the Single Primary Heat Extraction and Removal Emulator facility located at Idaho National Laboratory. The ultimate goals of this work are to (1) accurately monitor temperature and strain distributions that result from differential thermal expansion in the test articles and (2) quantify the heat rejection limits of heat pipes as a function of operating temperature and working fluid during steady-state and transient operations. Initial tests focused on the feasibility of using advanced fiber-optic sensors and other sensor technologies to improve the understanding of temperature and strain distributions within the electrically heated experiments. However, these sensing capabilities could benefit the broader microreactor community if the sensors could be used to monitor component and system health during nuclear operations to inform a limited number of microreactor operators, ultimately reducing operation and maintenance costs and moving toward semiautonomous operation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Design, Selection, and Evaluation of Reinforcement Learning Single Agents for Ground Target Tracking

Previous approaches for small fixed-wing unmanned air systems that carry strapdown rather than gimbaled cameras achieved satisfactory ground target tracking performance using both standard and deep reinforcement learning algorithms. However, these approaches have significant restrictions and abstractions to the dynamics of the vehicle, such as constant airspeed and constant altitude, because the number of states and actions was necessarily limited. Thus, extensive tuning was required to obtain good tracking performance. The expansion from 4 state–action degrees of freedom to 15 enabled the agent to exploit previous reward functions that produced novel yet undesirable emergent behavior. This paper investigates the causes of and various potential solutions to undesirable emergent behavior in the ground target tracking problem. A combination of changes to the environment, reward structure, action space simplification, command rate, and controller implementation provides insight into obtaining stable tracking results. Consideration is given to reward structure selection and refinement to mitigate undesirable emergent behavior. Results presented in the paper for a simulated environment of a single unmanned air system tracking a randomly moving single ground target show that a soft actor–critic algorithm can produce feasible tracking trajectories without limiting the state space and action space, provided that the environment is properly posed.

Engineering↗
DOE OSTI2020

Implementing Chiral Three-Body Forces in Terms of Medium-Dependent Two-Body Forces

Three-nucleon (3N) forces are an indispensable ingredient for accurate few-body and many- body nuclear structure and reaction theory calculations. While the direct implementation of chiral 3N forces can be technically very challenging, a simpler approach is given by employing instead a medium-dependent NN interaction V med that reflects the physics of three-body forces at the two-body normal-ordered approximation. We review the derivation and construction of Vmed from the chiral 3N interaction at next-to-next-to-leading order (N2LO), consisting of a long-range 2π-exchange term, a mid-range 1π-exchange component and a short-range contact-term. Several applications of V med to the equation of state of cold nuclear and neutron matter, the nucleon single-particle potential in nuclear matter, and the nuclear quasiparticle interaction are discussed. We also explore differences in using local vs. nonlocal regulating functions on 3N forces and make direct comparisons to exact results at low order in perturbation theory expansions for the equation of state and single-particle potential. We end with a discussion and numerical calculation of the in-medium NN potential V med from the next-to-next-to-next-to-leading order (N3LO) chiral 3N force, which consists of a series of long-range and short-range terms.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Influence of cation species on thermal expansion of Y 2 Si 2 O 7 –Gd 2 Si 2 O 7 solid solutions

Mixtures of Y 2 Si 2 O 7 and Gd 2 Si 2 O 7 were synthesized by solid-state reaction at 1600°C and characterized via in situ x-ray diffraction (XRD) to determine their coefficients of thermal expansion (CTE). All solid solutions within the system exhibited the orthorhombic δ-RE 2 Si 2 O 7 (Pna2 1 ) structure. Thermal expansion measurements of Y 2 Si 2 O 7 and Gd 2 Si 2 O 7 correlated well with reported values in literature, and all synthesized solid solutions exhibited CTEs between Y 2 Si 2 O 7 and Gd 2 Si 2 O 7 . Generally, there was a slight decrease in CTE exhibited by the materials with increasing Gd 2 Si 2 O 7 content, with Gd 2 Si 2 O 7 having the lowest CTEs and Y 2 Si 2 O 7 the highest CTEs. Here, the decrease in CTE was attributed to stronger bonds of Gd-O over Y-O, as determined by calculated crystal orbital Hamilton populations using density functional theory. However, such differences were very small and crystal structure was the dominating factor in CTE trends.

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

Variational quantum solver employing the PDS energy functional

In our previous work (J. Chem. Phys. 2020, 153, 201102),we reported a new class of quantum algorithms that are based on the quantum computation of the connected moment expansion to find the ground and excited state energies. In particular, the Peeters-Devreese-Soldatov (PDS) formulation is found variational and bearing the potential for further combining with the existing variational quantum infrastructure. Following this direction, here we propose a variational quantum solver employing the PDS energy gradient. In comparison with the usual variational quantum eigensolver (VQE) and the original static PDS approach, the proposed variational quantum solver offers an effective approach to achieve high ac-curacy at finding the ground state and its energy through the rotation of the trial wave function of modest quality guided by the low order PDS energy gradients, thus improves the ac-curacy and efficiency of the quantum simulation. We demonstrate the performance of the proposed variational quantum solver for toy models, H2molecule, and strongly correlated planar H4system in some challenging situations. In all the case studies, the proposed variational quantum approach out-performs the usual VQE and static PDS calculations even at the lowest order.

Peng, Bo↗
DOE OSTI2026

Energy-Screened Many-Body Expansion for Protein–Ligand Interactions: Examining Convergence for Metalloenzymes Through Seven–Body Interactions

Fragment-based quantum chemistry is a powerful strategy for calculating protein−ligand interaction energies using quantum chemistry methods. Rigorous convergence often requires hundreds of atoms in the protein binding-site model, especially if that model is constructed using distance-based criteria to select amino acid residues, while three- and four-body calculations exhibit instability related to combinatorial proliferation in the number of subsystem calculations. Here, we report an energy-based screening protocol for the many-body expansion applied to protein−ligand interactions, implemented in the open-source FRAGME∩T code. Using a combination of aggressive screening based on semiempirical quantum chemistry, with an improved graph-theoretical algorithm to eliminate unimportant subsystems, we are able to perform n-body calculations up to n = 7 using density functional theory in triple-ζ basis sets. Distance cutoffs further reduce the cost without compromising accuracy. Rapid and stable convergence of the many-body expansion is obtained by n = 4, for a pair of metalloenzymes in which a divalent ion coordinates directly to the ligand. As compared to previous results that relied solely on distance cutoffs, oscillations in the n-body corrections are reduced or eliminated, although residual errors remain in one case. This work demonstrates that benchmark-quality protein−ligand interaction energies can be systematically converged using a method with excellent parallel efficiency and scalability.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗