Search NASA⌕ Search

SEARCH · Search NASA

Results for “Formal methods”

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

A thermodynamically consistent discretization of 1D thermal-fluid models using their metriplectic 4-bracket structure

Thermodynamically consistent models in continuum physics, i.e. models which satisfy the first and second laws of thermodynamics, may be expressed using the metriplectic formalism. In this work, we leverage the structures underlying this modeling formalism to preserve thermodynamic consistency in discretizations of a fluid model. The procedure relies (1) on ensuring that the spatial semi-discretization retains certain symmetries and degeneracies of the Poisson and metriplectic 4-brackets, and (2) on the use of an appropriate energy conserving time-stepping method. Here, the minimally simple yet nontrivial example of a one-dimensional thermal-fluid model is treated. It is found that preservation of the requisite symmetries and degeneracies of the 4-bracket is relatively simple to ensure in Galerkin spatial discretizations, suggesting a path forward for thermodynamically consistent discretizations of more complex fluid models using more specialized Galerkin methods.

Hamiltonian structure↗
DOE OSTI2024

tLaSDI: Thermodynamics-informed latent space dynamics identification

Here we propose a latent space dynamics identification method, namely tLaSDI, that embeds the first and second principles of thermodynamics. The latent variables are learned through an autoencoder as a nonlinear dimension reduction model. The latent dynamics are constructed by a neural network-based model that precisely preserves certain structures for the thermodynamic laws through the GENERIC formalism. An abstract error estimate is established, which provides a new loss formulation involving the Jacobian computation of autoencoder. The autoencoder and the latent dynamics are simultaneously trained to minimize the new loss. Computational examples demonstrate the effectiveness of tLaSDI, which exhibits robust generalization ability, even in extrapolation. In addition, an intriguing correlation is empirically observed between a quantity from tLaSDI in the latent space and the behaviors of the full-state solution.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Stabilization of [(N 5 ) 2 BX] 2– and [(N 5 ) 2 B 2 X 2 ] 2– (X = H, F, Cl, Br) by Conjugation and Hyperconjugation Effects

The isolation of nucleophilic boron bases has led to a paradigm shift in boron chemistry. Previous studies of the bis(carbene) borylene complexes revealed that these compounds possess strong donor abilities, and their reaction inertness is due to the large steric hindrance between boron reagents and reactant. In the present study, we have theoretically studied the [(N 5 ) 2 BX] 2– and [(N 5 ) 2 B 2 X 2 ] 2– compounds (X = H, F, Cl, Br). Their electronic structures and properties are discussed by using the NBO, LOL, and ELF methods. We found that both π-conjugation and hyperconjugation effects can effectively stabilize the substituted nucleophilic anionic boron compounds [(N 5 ) 2 BX] 2– and [(N 5 ) 2 B 2 X 2 ] 2– . Substituents, especially X = H, stabilize the boron center through highly delocalized π-bonding, involving the formally “empty” in-plane p orbitals of the boron atom. While the halogen substituents have high electron withdrawal ability, leading to systems being less stable, we suggest the borinium anions [(N 5 ) 2 BH] 2– and [(N 5 ) 2 B 2 H 2 ] 2– as possible synthetic targets of novel environmentally friendly catalysts.

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

Photoluminescence line shapes of nanocrystals: Contributions from first- and second-order vibronic couplings

Here, we present a microscopic, parameter-free approach for computing the photoluminescence spectra of a single semiconductor nanocrystal. The method derives exciton–phonon coupling directly from the semi-empirical pseudopotential framework and systematically incorporates both diagonal and off-diagonal exciton-phonon interactions, expanded to second-order in the phonon coordinates. The dipole–dipole correlation function was calculated using a Dyson expansion within the Kubo–Toyozawa formalism, enabling a consistent description of the role of pure dephasing and population transfer on the photoluminescence spectral features. Applied to CdSe/CdS core–shell nanocrystals, the approach quantitatively reproduces experimental photoluminescence spectra over a wide temperature range, revealing that quadratic phonon couplings account for nearly half of the homogeneous linewidth above ≈ 100−150 K, while off-diagonal couplings leading to exciton thermalization play only a minor role and only as T → 300 K.

Dephasing Rate↗
DOE OSTI2025

A Parallel Alternative for Energy-Efficient Neural Network Training and Inferencing

Energy efficiency of training and inferencing with large neural network models is a critical challenge facing the future of sustainable large-scale machine learning workloads. This paper introduces an alternative strategy, called phantom parallelism, to minimize the net energy consumption of traditional tensor (model) parallelism, the most energy-inefficient component of large neural network training. The approach is presented in the context of feed-forward network architectures as a preliminary, but comprehensive, proof-of-principle study of the proposed methodology. We derive new forward and backward propagation operators for phantom parallelism, implement them as custom autograd operations within an end-to-end phantom parallel training pipeline and compare its parallel performance and energy-efficiency against those of conventional tensor parallel training pipelines. Formal analyses that predict lower bandwidth and FLOP counts are presented with supporting empirical results on up to 256 GPUs that corroborate these gains. Experiments are shown to deliver ∼50% reduction in the energy consumed to train FFNs using the proposed phantom parallel approach when compared with conventional tensor parallel methods. Additionally, the proposed approach is shown to train smaller phantom models to the same model loss on smaller GPU counts as larger tensor parallel models on larger GPU counts offering the possibility for even greater energy savings.

Seal, Sudip [ORNL] (ORCID:0000000332330656)↗
DOE OSTI2024

Thermodynamics-informed latent space dynamics identification

This software showcases a latent space dynamics identification method, namely tLaSDI, that embeds the first and second principles of thermodynamics. The latent variables are learned through an autoencoder as a nonlinear dimension reduction model. The latent dynamics are constructed by a neural network-based model that precisely preserves certain structures for the thermodynamic laws through the GENERIC formalism. An abstract error estimate is established, which provides a new loss formulation involving the Jacobian computation of autoencoder. The autoencoder and the latent dynamics are simultaneously trained to minimize the new loss. Computational examples demonstrate the effectiveness of tLaSDI, which exhibits robust generalization ability, even in extrapolation. In addition, an intriguing correlation is empirically observed between a quantity from tLaSDI in the latent space and the behaviors of the full-state solution.

Cheung, Siu Wun↗
DOE OSTI2026

Hydride and Seek: Comparing Crystallographic Hydride Placement Techniques with an Open-Shell Cobalt Complex

Locating hydrides is crucial in organometallic chemistry but difficult to do accurately using X-ray diffraction. Electron diffraction has been proposed as a way to overcome this problem but has not been systematically compared to neutron diffraction and to quantum crystallography (Hirshfeld atom refinement, HAR) to test this hypothesis. Here, we present a comparative analysis of methods for a terminal cobalt hydride complex by comparing a single-crystal neutron diffraction reference structure to results from single-crystal X-ray diffraction with and without Hirshfeld atom refinement (HAR, NoSpherA2), density functional theory (DFT), and electron diffraction (3D-ED/MicroED) refined under kinematical and dynamical formalisms. Conventional X-ray diffraction gives lower precision than neutron diffraction as expected. Despite expected improvements, HAR gives systematic deviation from the neutron benchmark. Interestingly, optimized DFT equilibrium geometries are closer to the neutron value than the value from HAR. On the other hand, electron diffraction with a high-quality data set coupled with dynamical refinement localizes the hydride in difference maps and gives excellent agreement with the neutron data. Dynamical refinement is crucial, as kinematical refinement does not allow assignment of a hydride peak. This cross-modal comparison defines the conditions under which 3D-ED/MicroED delivers high-precision metal–hydride distances for this open-shell cobalt hydride.

anions↗
DOE OSTI2025

Theory of tunneling between two-dimensional electron layers driven by spin pumping: Adiabatic regime and beyond

Tunneling spectroscopy between parallel two-dimensional (2D) electronic systems provides a powerful method to probe the underlying electronic properties by measuring tunneling conductance. Here, in this work, we present a theoretical framework for spin transport in 2D-to-2D tunneling systems, driven by spin pumping. This theory applies to a vertical heterostructure where two layers of metallic 2D electron systems are separated by an insulating barrier, with one layer exchange coupled to a magnetic layer driven at resonance. Utilizing a nonperturbative Floquet-Keldysh formalism, we derive general expressions for the tunneling spin and charge currents across a broad range of driving frequencies, extending beyond the traditional adiabatic pumping regime. At low frequencies, we obtain analytical results that recover the known behaviors in the adiabatic regime. However, at higher frequencies, our numerical findings reveal significant deviations in the dependence of spin and charge currents on both frequency and precession angle. This work offers fresh insights into the role of magnetization dynamics in tunneling transport, opening up new avenues for exploring nonadiabatic spin pumping phenomena.

Green's function methods↗
DOE OSTI2026

Linear-Scaling Asymmetric Triples Correction through the Solution of the DLPNO–CCSD Lambda Equations: DLPNO–CCSD(T) Λ

In this research, we derive equations for solving for the stationary points of the DLPNO–CCSD Lagrangian, in the t 1 -transformed formalism introduced earlier and as currently implemented in the P SI 4 quantum chemistry software package. These lambda equations in the local pair natural orbital basis allow for the evaluation of CCSD(T) Λ energetics with linear-scaling computational effort, also known as the asymmetric triples correction. This DLPNO–CCSD(T) Λ method allows for accurate triples contributions to be computed for larger molecules, especially in cases that CCSD(T) is known to be insufficient, such as with multireference systems and bond-breaking systems. We showcase the accuracy of our code on reaction energies, barrier heights, and noncovalent interaction energies. Also showcased are the capabilities of our code by evaluating DLPNO–CCSD(T) Λ energetics on large noncovalent dimers up to 112 atoms, as well as a rhodium catalyst complex containing 66 atoms.

Cluster chemistry↗
DOE OSTI2025

$\gamma$-ray angular correlations, distributions and linear polarization in tracking arrays

With $\gamma $-ray tracking arrays, the angular correlation and, especially, the linear polarization of $\gamma $ rays can be measured with high precision. However, because of the complexity of these spectrometers, measuring angular distributions is more involved. The techniques and associated formalisms were developed and established in the 1950s. The experimental setups of that era were relatively simple when compared to the advanced $\gamma $-ray tracking arrays available today. This paper aims to review and refine these techniques, unify the nomenclatures and, specifically, adapt them for use with the newest tracking arrays. In addition, new methods for producing the required response functions for these arrays are introduced.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

SF6 / Acetone Separation & Purification

Sulfur hexafluoride (SF₆) is recognized as the most potent greenhouse gas used in the power transmission and semiconductor sectors. In the past decade, its global emissions have increased sharply, largely due to the absence of effective disposal and recovery pathways. This concern was formally addressed under the Kyoto Protocol, which urged participating nations to adopt measures to limit SF₆ emissions and mitigate its role in climate change. In recent years, various abatement strategies have been explored, including non-thermal plasma (NTP) technologies such as radio frequency, microwave, dielectric barrier discharge, and electron beam systems. While these methods can decompose SF₆, they also produce hazardous by-products like sulfur oxyfluorides, sulfur dioxide, hydrofluoric acid, and fluorine gas. Because these compounds pose both environmental and health risks, replacing traditional disposal practices with modern, more effective treatment methods has become an urgent priority.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Multipole responses in fissioning nuclei and their uncertainties

Electromagnetic multipole responses are key inputs to model the structure, decay, and reactions of atomic nuclei. With the introduction of the finite amplitude method (FAM), large-scale calculations of the nuclear linear response in heavy deformed nuclei have become possible. This work provides a detailed study of multipole responses in actinide nuclei with Skyrme energy density functionals. We quantify both systematic and statistical uncertainties induced by the functional parametrization in FAM calculations. We also extend the FAM formalism to perform blocking calculations with the equal filling approximation for odd-mass and odd-odd nuclei, and analyze the impact of blocking configurations on the response. In conclusion, by examining the entire plutonium isotopic chain from the proton to the neutron drip line, we find a large variability of the response with the neutron number and study how it correlates with the deformation of the nuclear ground state.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

SYCL for Performance Portability: Application Experience with Coupled Cluster Formalism in Quantum Chemistry on Exascale Systems

The exascale computing has brought unprecedented heterogeneity in node architectures, with systems such as Frontier and Aurora featuring diverse GPU accelerators, network connectivity among others. Ensuring performance portability across these platforms is a key challenge. To address this, we employ the SYCL programming model to develop portable, high-performance quantum chemistry workloads. As a representative application, we focus on the non-iterative Triples component of the coupled-cluster CCSD(T) method, a key driver in quantum chemistry. In this work, we report on our experience deploying SYCL-based implementations using both DPC++ and AdaptiveCPP across two flagship exascale platforms: OLCF Frontier with AMD MI250X GPUs and ALCF Aurora with Intel GPUs. Our results demonstrate that SYCL enables efficient, single-source implementations that scale to thousands of nodes, delivering performance on par with vendor-optimized HIP solutions. We highlight key insights into runtime behavior, kernel portability, and scaling characteristics, showing that SYCL offers a viable path for performance-portable computing.

Bagusetty, Abhishek [Argonne National Laboratory (↗
DOE OSTI2025

Metal–Metal Bonding Influences Hydride Reactivity in [Sn–Rh] 3+ and [Sn–Ni] 2+ Bimetallics

Heavier group 14 metal hydrides serve as key intermediates in catalytic transformations, such as hydroboration. Regenerating such intermediates via a clean hydride source like dihydrogen could provide catalytic processes with a more economical alternative to, e.g., silanes and hydridoborane reagents. Herein, we report our efforts toward this goal using a [Sn–Rh] 3+ bimetallic system with the formal Rh I center acting as a potential dihydrogen activator. Targeting the introduction of a hydride ligand to the bimetallic core, our reactivity studies have revealed a preference for the hydride ligand to be bound to the Rh center, instead of the Sn center. Finally, examination of the electronic structures of these complexes via theoretical and experimental methods has revealed the electron acceptor nature of the Rh center within the bimetallic core and offers an explanation for the localization of the hydride between metal centers.

anions↗
DOE OSTI2025

Impossibility of obtaining time-independent, three-dimensional, spherically symmetric densities of confined systems of relativistically moving constituents

The quantum-mechanical definition of probability, the uncertainty principle, and Poincaré invariance provide strong basic restrictions on the ability to define spatial densities associated with form factors describing the properties of confined systems of relativistically moving constituents. Despite this, many papers ignore one or more of these restrictions. Here I show how to obtain time-independent , two-dimensional densities that are consistent with the stated restrictions. This is done using the light-front, infinite momentum frame formalism. Two-dimensional density interpretations of the axial-vector form factor and all three gravitational form factors is obtained. The resulting mass radius is smaller than the charge radius. Additionally, an expression of a two-dimensional mass density related to the trace of the energy momentum tensor is obtained. I also show that all known methods for finding three-dimensional densities—using the Breit frame, Abel transformations, Wigner distributions, and spherically symmetric wave packets with vanishing spatial extent—violate the basic restrictions in different ways. Furthermore, the use of the latter leads to densities that vanish almost everywhere in space as time increases from an initial value.

form factors↗
DOE OSTI2025

Measuring the Burgers vector of dislocations with dark-field X-ray microscopy

The subsurface dynamics of dislocations are essential to many properties of bulk crystalline materials. However, it is challenging to characterize a bulk crystal by conventional transmission electron microscopy (TEM) due to the limited penetration depth of electrons. A novel X-ray imaging technique – dark-field X-ray microscopy (DFXM) – was developed to image hierarchical dislocation structures in bulk crystals. While today's DFXM can effectively map the line structures of dislocations, it is still challenging to quantify the Burgers vectors, the key characterization governing the dislocation behaviors. Here, we extend the 'invisibility criterion' formalism from the TEM theory to the geometrical-optics model of DFXM and demonstrate the consistency between DFXM and dark-field TEM using multi-diffraction-peak imaging for a single edge dislocation. Due to the practical difficulty of multi-peak DFXM experiments, we further study how the Burgers vector effect is encoded for a single-peak DFXM experiment. Using the geometrical-optics DFXM simulation, we explore the asymmetry of rocking tilt scans at different rolling tilts and develop a new method to characterize the Burgers vector. The conclusions of this study advance our understanding of the use of DFXM in characterizing individual dislocations, enabling the connection from bulk DFXM imaging to dislocation mechanics.

36 MATERIALS SCIENCE↗
DOE OSTI2026

AthenaK: A Performance-portable Version of the Athena++ Adaptive Mesh Refinement Framework

We describe AthenaK: a new implementation of the Athena++ block-based adaptive mesh refinement framework using the Kokkos programming model. Finite volume methods for Newtonian, special relativistic, and general relativistic (GR) hydrodynamics and magnetohydrodynamics (MHD), and GR-radiation hydrodynamics and MHD, as well as a module for evolving Lagrangian tracer or charged test particles (e.g., cosmic rays) are implemented using the framework. In two companion papers, we describe (1) a new solver for the Einstein equations based on the Z4c formalism, and (2) a GRMHD solver in dynamical spacetimes also implemented using the framework, enabling new applications in numerical relativity. By adopting Kokkos, the code can be run on virtually any hardware, including CPUs, GPUs from multiple vendors, and emerging Advanced RISC Machine processors. AthenaK shows excellent performance and weak scaling, achieving over 1 billion cell updates per second for hydrodynamics in three dimensions on a single NVIDIA Grace Hopper processor. It does this with a typical parallel efficiency of 80% on 65,536 AMD GPUs on the OLCF Frontier system. Such performance portability enables AthenaK to leverage modern exascale computing systems for challenging applications in astrophysical fluid dynamics, numerical relativity, and multimessenger astrophysics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Gaunt and Breit two-electron contributions to mean-field transformations and fine structure splitting

Materials utilized by novel energy systems are often studied using weakly correlated mean-field theories. However, if these systems incorporate heavy elements, relativistic effects must be included. Therefore, a Kramers unrestricted coupled cluster with singles and doubles excitation formalism within a molecular mean-field exact two-component framework (X2C mmf ) using a four-component Dirac–Hartree–Fock (DHF) reference state is presented. The exact X2C mmf transformed normal-order Hamiltonian incorporates all one-electron and two-electron (2e) contributions from the Coulomb, Gaunt, and Breit operators and is used with the equation of motion method to calculate the excitation energies of the alkali group of elements. Using this framework, the effects of 2e Gaunt and Breit integrals are studied. Results demonstrate growing contributions from these integrals to the generated X2C mmf mean-fields and electronic fine structure calculations with increasing atomic number. Overall, this paper outlines the method, its effect within the X2C mmf approach, and lays the foundation for future theoretical development of relativistic calculations within this framework.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗