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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9

All sources NASA NTRS DOE OSTI

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↗

DOE OSTI2024

Fine structure of current noise spectra in nanoelectromechanical resonators

Here, we study the frequency-dependent noise of a suspended carbon nanotube quantum dot nanoelectromechanical resonator induced by electron-vibration coupling. Using a rigorous Keldysh diagrammatic technique, we establish a formal framework connecting the vibrational properties to electrical measurements. We find that the noise power spectrum exhibits a narrow resonant peak at the frequency of the vibrational modes. However, this fine structure tends to disappear due to a coherent cancellation effect when the tunneling barriers are tuned to a symmetric point. Notably, measuring the electrical current noise spectra provides a sensitive alternative method for detecting the damping and dephasing of quantum vibrational modes.

Electric measurements↗

DOE OSTI2025

Orthogonality broadcasting and quantum position verification

The no-cloning theorem leads to information-theoretic security in various quantum cryptographic protocols. However, this security typically derives from a possibly weaker property that classical information encoded in certain quantum states cannot be broadcast. To formally capture this property, we introduce the study of ‘orthogonality broadcasting.’ When attempting to broadcast the orthogonality of two different qubit bases, we establish that the power of classical and quantum communication is equivalent. However, quantum communication is shown to be strictly more powerful for broadcasting orthogonality in higher dimensions. We then relate orthogonality broadcasting to quantum position verification and provide a new method for establishing error bounds in the no pre-shared entanglement model that can address protocols previous methods could not. Our key technical contribution is an uncertainty relation that uses the geometric relation of the states that undergo broadcasting rather than the non-commutative aspect of the final measurements.

quantum cryptography↗

DOE OSTI2024

Sparse-Stochastic Fragmented Exchange for Large-Scale Hybrid Time-Dependent Density Functional Theory Calculations

Here we extend our recently developed sparse-stochastic fragmented exchange formalism for ground-state near-gap hybrid DFT to calculate absorption spectra within linear-response time-dependent generalized Kohn-Sham DFT (LR-GKS-TDDFT) for systems consisting of thousands of valence electrons within a grid-based/plane-wave representation. A mixed deterministic/fragmented-stochastic compression of the exchange kernel, here using long-range explicit exchange functionals, provides an efficient method for accurate optical spectra. Both real-time propagation as well as frequency-resolved Casida-equation-type approaches for spectra are presented, and the method is applied to large molecular dyes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Porting Classical Approaches for Quantum Simulations to Quantum Computers

Simulating quantum many-body systems is one of the most promising problems in which we might anticipate that quantum computers should show quantum advantage. Unfortunately, there is still a gap between this promise and actual practice. New quantum algorithms need to be developed and the current quantum algorithms have various difficulties - e.g efficient state preparation - which must be overcome and improved upon. In many cases, classical approaches need to be ported over to quantum devices. In this project we have developed a suite of new quantum algorithms which makes progress in this regard. We developed a new optimization scheme for variational quantum eigensolvers, UBOS, which mitigates problems with local minimas and barren plateaus while improving convergence to the ground state by an order of magnitude. We developed a new way to utilize qubitization to find ground states of nearly frustration-free Hamiltonians faster than all previous methods. We developed a series of state preparation techniques which helps initialize parameterized quantum circuits into reasonable starting points on which quantum algorithms are then applied. In addition to the development of novel algorithms, it is critical to have classical simulation techniques for approximately simulating quantum circuits which can be used to benchmark and understand quantum algorithms. Toward that end, we developed a novel POVM formalism to simulate quantum circuits as well as exemplify the massive parallelization of tensor network methodologies. Finally, we developed physical understanding of entanglement phase transitions such as many-body localization and random tensor networks.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Interpretable, extensible linear and symbolic regression models for charge density prediction using a hierarchy of many-body correlation descriptors

Here, density functional theory (DFT) is routinely used to make electronic structure predictions for high-throughput screening of materials and molecules for technologically relevant areas, like the identification of better catalysts, electronic materials, and drug discovery. However, the DFT formalism is limited by (a) its poor (quadratic-to-quartic) scaling, and (b) the need to perform repeated eigenvalue computations of the electronic Hamiltonian as part of its self-consistent field (SCF) iteration procedure to obtain the converged ground state electron density, ρ (r). Approaches that directly predict ρ (r) of a structure with high accuracy can accelerate conventional SCF calculations and can also be used in linearly scaling methods such as orbital-free DFT. To this end, we present a procedure to predict the ground state electron density of molecular and periodic three-dimensional systems directly from the atomic structure with a particular emphasis on physical interpretability. In our framework, ρ (r) is modeled using many-body correlation descriptors that accurately capture the effects of local atomic arrangements in the neighborhood of a grid point. Our use of a linear regression scheme to fit to charge density data enables transparent analysis of the relative contributions of various types of local atomic correlations. By systematically including increasingly complex correlations, our model is shown to accurately predict ρ (r) for a variety of chemically and electronically diverse systems — amorphous Ge, Al(001) slab, crystalline Ga 2 O 3 , molecular benzene, and polyethylene. We then demonstrate a symbolic regression-based protocol to construct easily computable, interpretable features from lower-order correlations that significantly improves our electron density predictions with effectively no increase in the computational cost.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

Actinide 5f Occupations: The Case of PuO 2

In actinide chemistry, the formal number of open-shell 5f electrons, n open , is a well-defined quantity with an integer value. The effective 5f occupation, n f , additionally takes donation and back-donation into account, generally has a non-integer value, and has varying numerical definitions. The present study explores the important distinction between n f and n open in actinide chemistry with the example of PuO2, by using electronic structure methods with a relativistic Hamiltonian in combination with experimental Pu M 5 -edge high-energy-resolution X-ray absorption and emission spectroscopic data. The total donation to the metal in PuO 2 is between 3.1 and 2.4 electrons, depending on the type of calculation, most of which is to the Pu 6d and 5f shells. The donation into 5f is sensitive to the approximations in the electronic structure model but likely amounts to 1.6/0.8 electrons when the diffuse regions of the 5f shell are included/excluded. Valence band resonant inelastic X-ray scattering experiments demonstrate that Pu 5f electron density is present in the valence band; thus, there is a clear experimental signature of covalent bonding in PuO 2 . Pu M 5 -edge and M 3 -edge high-energy-resolution X-ray absorption near-edge structures for Pu 3+ and Pu 4+ in an aqueous solution are compared to PuO 2 , showing that Pu in PuO 2 has an nf closer to Pu 4+ (aq).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Accelerated data-driven materials science with the Materials Project

The Materials Project was launched formally in 2011 to drive materials discovery forwards through high-throughput computation and open data. More than a decade later, the Materials Project has become an indispensable tool used by more than 600,000 materials researchers around the world. This Perspective describes how the Materials Project, as a data platform and a software ecosystem, has helped to shape research in data-driven materials science. We cover how sustainable software and computational methods have accelerated materials design while becoming more open source and collaborative in nature. Next, we present cases where the Materials Project was used to understand and discover functional materials. We then describe our efforts to meet the needs of an expanding user base, through technical infrastructure updates ranging from data architecture and cloud resources to interactive web applications. Finally, we discuss opportunities to better aid the research community, with the vision that more accessible and easy-to-understand materials data will result in democratized materials knowledge and an increasingly collaborative community.

Horton, Matthew K↗

DOE OSTI2026

Correlation function metrology for warm dense matter: Recent developments and practical guidelines

X-ray Thomson scattering (XRTS) has emerged as a valuable diagnostic for matter under extreme conditions, as it captures the intricate many-body physics of the probed sample. Recent advances, such as the model-free temperature diagnostic of Dornheim et al. [Nat. Commun. 13 , 7911 (2022)], have demonstrated how much information can be extracted directly within the imaginary-time formalism. However, since the imaginary-time formalism is a concept often difficult to grasp, we provide here a systematic overview of its theoretical foundations and explicitly demonstrate its practical applications to temperature inference, including relevant subtleties. Furthermore, we present recent developments that enable the determination of the absolute normalization, Rayleigh weight, and density from XRTS measurements without reliance on uncontrolled model assumptions. Finally, we outline a unified workflow that guides the extraction of these key observables, offering a practical framework for applying the method to interpret experimental measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Data-driven particle dynamics: Structure-preserving coarse-graining for emergent behavior in non-equilibrium systems

Multiscale systems are ubiquitous in science and technology, but are notoriously challenging to simulate as short spatiotemporal scales must be appropriately linked to emergent bulk physics. When expensive high-dimensional dynamical systems are coarse-grained into low-dimensional models, the entropic loss of information leads to emergent physics which are dissipative, history-dependent, and stochastic. To machine learn coarse-grained dynamics from time-series observations of particle trajectories, we propose a framework using the metriplectic bracket formalism that preserves these properties by construction; most notably, the framework guarantees discrete notions of the first and second laws of thermodynamics, conservation of momentum, and a discrete fluctuation-dissipation balance crucial for capturing non-equilibrium statistics. We introduce the mathematical framework abstractly before specializing to a particle discretization. As labels are generally unavailable for entropic state variables, we introduce a novel self-supervised learning strategy to identify emergent structural variables. We validate the method on benchmark systems and demonstrate its utility on two challenging examples: (1) coarse-graining star polymers at challenging levels of coarse-graining while preserving non-equilibrium statistics, and (2) learning models from high-speed video of colloidal suspensions that capture coupling between local rearrangement events and emergent stochastic dynamics. We provide open-source implementations in both PyTorch and LAMMPS, enabling large-scale inference and extensibility to diverse particle-based systems.

Computational Engineering, Finance, and Science (c↗

DOE OSTI2024

Assessing spin-density wave formation in La 3 Ni 2 O 7 from electronic structure calculations

Here, we employ correlated density-functional theory methods (DFT + Hubbard U) to investigate the spin-density wave state of the bilayer Ruddlesden-Popper (RP) nickelate La 3 Ni 2 O 7 which becomes superconducting under pressure. We predict that the ground state of this bilayer RP material has traits of both the double spin-stripe and the single spin-charge stripe phases proposed in the literature as it corresponds to in-plane up/up′/down/down′ diagonal stripes with up/down being high spin (formally Ni 2+ : d 8 ), and up′ /down′ being low spin (formally Ni 3+ : d 7 ). The main feature of this solution (that is insulating even at U = 0) is the dominant role of the d x 2 −y 2 bands around the Fermi level, which would become doped with the introduction of electrons via oxygen vacancies. In spite of the similarity with cuprates in terms of the dominant role of d x 2 −y 2 bands, some differences are apparent in the magnetic ground state of La 3 Ni 2 O 7 : the antiferromagnetic out-of-plane coupling within the bilayer (linked to the d z 2 orbitals forming a spin-singlet-like configuration) is found to be the dominant one while in-plane interactions are reduced due to the stripe order of the ground state. With pressure, this striped magnetic ground state remains similar in nature but the increase in bandwidth quickly transitions La 3 Ni 2 O 7 into a metallic state with all the activity close to the Fermi level involving, to a large extent, d x 2 −y 2 orbitals. This is reminiscent of the cuprates and may provide key insights into how superconductivity arises in this material under pressure.

LaBollita, Harrison↗

DOE OSTI2025

The Zintl pnictides Yb10CdSb9 and Yb14CdSb11: New candidate thermoelectric materials

The synthesis of new materials is the lifeline of solid-state science, and it continues to offer us unique opportunities for testing various theoretical formulations and models on a practical material. Such an avenue, therefore, provides a breeding ground for technological innovations and advancements that can completely revolutionize our world. Here, we report the results of our exploratory syntheses in the Yb–Cd–Sb compositional space that lead to the identification of two new Zintl antimonides, namely, Yb10CdSb9 and Yb14CdSb11. Their crystal structures were established via single-crystal X-ray diffraction methods; the basic electronic and transport properties of the new materials were also characterized. Yb10CdSb9 crystallizes in a disordered variant of the tetragonal Ca10LiMgSb9 structure type with unit cell parameters a = 11.8473(8) Å and c = 17.1302(12) Å (space group P42/mnm). Yb14CdSb11 crystallizes in the tetragonal Ca14AlSb11 structure type with unit cell parameters: a = 16.605(3) Å and c = 12.144(7) Å (space group I41/acd). Although the structures of both compounds can be rationalized within the framework of the Zintl formalism, based on the partitioning of the valence electrons in the much disordered Yb10CdSb9 phase, the charge is indicative of a slightly electron-rich composition. Electronic structure calculations in both cases support the notion of intrinsic semiconductor behavior, as expected for a Zintl phase. The temperature dependence of the electrical resistivity of a single crystal of Yb10CdSb9 is in line with that, and the evolution of the Seebeck coefficient indicates an electron-dominated transport mechanism, and a respectable power factor of 0.71 μW/cm K2 at 460 K can be calculated for Yb10CdSb9. The electrical resistivity of Yb14CdSb11, however, evolves in a semimetallic manner, which could suggest an overdoped sample or degenerate semiconducting behavior.

Ogunbunmi, Michael O. (ORCID:0000000253409198)↗

DOE OSTI2025

𝐴𝑏 initio density-matrix approach to exciton coherence: Phonon scattering, Coulomb interactions, and radiative recombination

Relaxation processes following light excitation in semiconductors are key in materials-based quantum technology applications. These processes are broadly studied in atomically thin transition-metal dichalcogenides, quasi-two-dimensional excitonic semiconductors in which atomistic design allows for tunable excited-state properties, such as relaxation lifetimes and photoinduced coherence. In this work, we present a density-matrix-based approach to compute exciton relaxation within a many-body ab initio perspective. We expand our previously developed Lindblad density-matrix formalism to capture multichannel electron-hole pair relaxation processes, including phonon and Coulomb scattering as well as radiative recombination, and we study their effect on the time-resolved excited-state propagation. Using monolayer MoSe 2 as a prototypical example, we examine many-body effects on the time-dependent dynamics of photoactive excitations, exploring how the electron-hole pair interactions are reflected in variations of the excitation energy, spectral signature, and state coherence. In conclusion, our method supplies a detailed understanding of exciton relaxation mechanisms in realistic materials, offering a previously unexplored pathway to study excited-state dynamics in semiconductors from first principles.

Band structure methods↗