Search NASA⌕ Search

SEARCH · Search NASA

Results for “Approximate computing”

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 919 records · Page 51

SGD-Net: Efficient Model-Based Deep Learning with Theoretical Guarantees

Deep unfolding networks have recently gained popularity for solving imaging inverse problems. However, the computational and memory complexity of data-consistency layers within traditional deep unfolding networks scales with the number of measurements, limiting their applicability to large-scale imaging inverse problems. We propose SGD-Net as a new methodology for improving the efficiency of deep unfolding through stochastic approximations of the data-consistency layers. Our theoretical analysis shows that SGD-Net can be trained to approximate batch deep unfolding networks to an arbitrary precision. Our simulations on intensity diffraction tomography and sparse-view computed tomography show that SGD-Net can match the performance of the traditional batch network at a fraction of training and testing complexity.Deep unfolding networks have recently gained popularity for solving imaging inverse problems. However, the computational and memory complexity of data-consistency layers within traditional deep unfolding networks scales with the number of measurements, limiting their applicability to large-scale imaging inverse problems. We propose SGD-Net as a new methodology for improving the efficiency of deep unfolding through stochastic approximations of the data-consistency layers. Our theoretical analysis shows that SGD-Net can be trained to approximate batch deep unfolding networks to an arbitrary precision. Our simulations on intensity diffraction tomography and sparse-view computed tomography show that SGD-Net can match the performance of the traditional batch network at a fraction of training and testing complexity.

97 MATHEMATICS AND COMPUTING↗

Mapping Thermoelectric Transport in a Multicomponent Alloy Space

Abstract Interest in high entropy alloy thermoelectric materials is predicated on achieving ultralow lattice thermal conductivity κ L through large compositional disorder. However, here it is shown that for a given mechanism, such as mass contrast phonon scattering, κ L will be minimized along the binary alloy with highest mass contrast, such that adding an intermediate mass atom to increase atomic disorder can increase thermal conductivity. Only when each component adds an independent scattering mechanism (such as adding strain fluctuation to an existing mass fluctuation) is there a benefit. In addition, both charge carriers and heat‐carrying phonons are known to experience scattering due to alloying effects, leading to a trade‐off in thermoelectric performance. Analytic transport models are applied, based on perturbation and effective medium theories, to predict how alloy scattering will affect the thermal and electronic transport across the full compositional range of several pseudo‐ternary and pseudo‐quaternary alloy systems. To do so, a multicomponent extension is demonstrated to both thermal and electronic binary alloy scattering models based on the virtual crystal approximation. Finally, it is shown that common functional forms used in computational thermodynamics can be applied to this problem to further generalize the scattering behavior that is modeled.

36 MATERIALS SCIENCE↗

Sm 2 Ru 3 Sn 5 : A Noncentrosymmetric Cubic Member of the Ln 2 M 3 X 5 Family

An optimized synthetic method is presented for Sm 2 Ru 3 Sn 5 and investigate its physical properties and electronic structure. Sm 2 Ru 3 Sn 5 is prepared by arc-melting stoichiometric ratios of the elements and is confirmed by single crystal and powder X-ray diffraction. An antiferromagnetic transition is observed at T N = 3.8 K. A modified Curie-Weiss fit to the data in the range 50–150 K yields a Curie-Weiss temperature: θ CW = −36.6 K and an effective magnetic moment: μ eff = 0.83 μ B , in agreement with a Sm 3+ oxidation state. Field-dependent magnetization up to H = 7 T at 2 K shows a maximum response of 0.06 μ B , which is significantly lower than the expected Sm 3+ saturation moment (0.71 μ B ). Resistivity measurements indicate metallic behavior, and analysis of the magnetic entropy from the heat capacity reveals a doublet ground state due to crystal electric field splitting. The electronic structure and density of states are calculated with density function theory and further supported by the local density approximation with dynamical mean-field theory. Finally, the experimental and computational results highlight localized Sm 3+ moments and suggest a possible interplay between Ruddelman–Kitel–Kasuya–Yosida and Kondo interactions, positioning Sm 2 Ru 3 Sn 5 as a promising material for studying topology and complex physical phenomena.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Subleading corrections to the S 3 free energy of necklace quiver theories dual to massive IIA

We investigate the S 3 free energy of N = 3 Chern-Simons-matter quiver gauge theories with gauge group U(N) r (r ≥ 2) where the sum of Chern-Simons levels does not vanish, beyond the leading order in the large-N limit. We take two different approaches to explore the sub-leading structures of the free energy. First we evaluate the matrix integral for the partition function in the ’t Hooft limit using a saddle point approximation. Second we use an ideal Fermi-gas model to compute the same partition function, but in the limit of fixed Chern-Simons levels. The resulting expressions for the free energy F = – log Z are then compared in the overlapping parameter regime. The Fermi-gas approach also hints at a universal 1/6 log N correction to the free energy. Since the quiver gauge theories we consider are dual to massive Type IIA theory, we expect the sub-leading correction of the planar free energy in the large ’t Hooft parameter limit to match higher-derivative corrections to the tree-level holographic dual free energy, which have not yet been fully investigated.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Half Metallic Ferromagnetic Character in ZnXP{sub 2} (X = Ge, Si) Chalcopyrites Doped with Mn

The electronic and the magnetic properties of ZnXP{sub 2} (X = Ge, Si) were studied using the Korring-Kohn-Rostoker (KKR) method combined with the coherent potential approximation (CPA). The total and partial density of state (DOS) are computed for different Mn concentrations. The total magnetic moment, FM and DLM energies, and their variation as well as Curie temperature are also given. It is shown that the substitution of X cations (X = Ge, Si) by Mn atoms in ZnXP{sub 2} chalcopyrites leads to half metallic ferromagnetic character with double exchange mechanism. Thus, the critical temperature can be controlled by varying the concentration of manganese impurity.

36 MATERIALS SCIENCE↗

Estimation of stress and stress-induced permeability change in a geological nuclear waste repository in a thermo-hydrologically coupled simulation

Geologic disposal is a promising solution for a safe permanent isolation of accumulated high-level nuclear waste from nuclear powerplants. The behavior of host rock is highly coupled thermally, hydromechanically, and chemically. Numerical simulations of such coupled phenomena for the extremely long term (>100,000 years) and large length scale (>1 km) of geologic disposal remain to be computationally challenging. In this study, a methodology has been developed to approximate stress and stress-induced permeability change in host rock using only thermo-hydrological (TH) variables. A coupled thermo-hydromechanical (THM) simulation is carried out using TOUGH-FLAC simulator to model THM behaviors of a generic nuclear waste repository, in order to evaluate the performance of the developed methodology, which is implemented in a coupled TH simulation. Results show that stress and permeability change estimated by the developed methodology in the TH-coupled simulation match those calculated in the THM-coupled simulation over the simulated timeframe of over 10,000 years. Details about errors in stress and permeability estimates accrued by the developed methodology are discussed here. The developed methodology will help incorporate stress-induced permeability change into existing TH simulators for the long-term radionuclide transport in geologic disposal.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

A higher-order finite-element implementation of the nonlinear Fokker–Planck collision operator for charged particle collisions in a low density plasma

Collisions between particles in a low density plasma are described by the Fokker–Planck collision operator. In applications, this nonlinear integro-differential operator is often approximated by linearised or ad-hoc model operators due to computational cost and complexity. In this work, we present an implementation of the nonlinear Fokker–Planck collision operator written in terms of Rosenbluth potentials in the Rosenbluth–MacDonald–Judd (RMJ) form. The Rosenbluth potentials may be obtained either by direct integration or by solving partial differential equations (PDEs) similar to Poisson's equation: we optimise for performance and scalability by using sparse matrices to solve the relevant PDEs. We represent the distribution function using a tensor-product continuous-Galerkin finite-element representation and we derive and describe the implementation of the weak form of the collision operator. We present tests demonstrating a successful implementation using an explicit time integrator and we comment on the speed and accuracy of the operator. Finally, we speculate on the potential for applications in the current and next generation of kinetic plasma models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Spectrally accelerated edge and scrape-off layer gyrokinetic turbulence simulations

This paper presents the first gyrokinetic (GK) simulations of edge and scrape-off layer (SOL) turbulence accelerated by a velocity-space spectral approach in the full-f GK code GENE-X. Building upon the original grid velocity-space discretization, we derive and implement a new spectral formulation and verify the numerical implementation using the method of manufactured solution. We conduct a series of spectral turbulence simulations focusing on the TCV-X21 reference case (Oliveira et al., 2022 [26]) and compare these results with previously validated grid simulations (Ulbl et al., 2023 [25]). The spectral approach reproduces the outboard midplane (OMP) profiles (density, temperature, and radial electric field), dominated by trapped electron mode (TEM) turbulence, with excellent agreement and significantly lower velocity-space resolution. As a consequence, the spectral approach reduces the computational cost (CPUh) by at least an order of magnitude, of approximately 50 for the TCV-X21 case. This enables high-fidelity GK simulations to be performed within a few days on modern CPU-based supercomputers for medium-sized devices and establishes GENE-X as a powerful tool for studying edge and SOL turbulence, moving towards reactor-relevant devices like ITER.

Gyrokinetic↗

A hybrid deep neural operator/finite element method for ice-sheet modeling

One of the most challenging and consequential problems in climate modeling is to provide probabilistic projections of sea level rise. A large part of the uncertainty of sea level projections is due to uncertainty in ice sheet dynamics. At the moment, accurate quantification of the uncertainty is hindered by the cost of ice sheet computational models. In this work we develop a hybrid approach to approximate existing ice sheet models at a fraction of their cost. Our approach consists of replacing the finite element model for the momentum equations for the ice velocity, the most expensive part of an ice sheet model, with a Deep Operator Network, while we retain a classic finite element discretization for the evolution of the ice thickness. We show that the resulting hybrid model is very accurate and it is an order of magnitude faster than the traditional finite element model. Further, a distinctive feature of the proposed model, compared to other neural network approaches, is that it can handle high-dimensional parameter spaces (parameter fields) such as the basal friction at the bed of the glacier and can therefore be used for generating samples for uncertainty quantification. Further, we study the impact of hyper-parameters, number of unknowns and correlation length of the parameter distribution on the training and accuracy of the Deep Operator Network on a synthetic ice sheet model. We then target the evolution of the Humboldt glacier in Greenland and show that our hybrid model can provide accurate statistics of the glacier mass loss and can be effectively used to accelerate the quantification of uncertainty.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Probing the consistency of cosmological contours for supernova cosmology

As the scale of cosmological surveys increases, so does the complexity in the analyses. This complexity can often make it difficult to derive the underlying principles, necessitating statistically rigorous testing to ensure the results of an analysis are consistent and reasonable. This is particularly important in multi-probe cosmological analyses like those used in the Dark Energy Survey (DES) and the upcoming Legacy Survey of Space and Time, where accurate uncertainties are vital. In this paper, we present a statistically rigorous method to test the consistency of contours produced in these analyses and apply this method to the Pippin cosmological pipeline used for type Ia supernova cosmology with the DES. We make use of the Neyman construction, a frequentist methodology that leverages extensive simulations to calculate confidence intervals, to perform this consistency check. A true Neyman construction is too computationally expensive for supernova cosmology, so we develop a method for approximating a Neyman construction with far fewer simulations. We find that for a simulated dataset, the 68% contour reported by the Pippin pipeline and the 68% confidence region produced by our approximate Neyman construction differ by less than a percent near the input cosmology; however, they show more significant differences far from the input cosmology, with a maximal difference of 0.05 in $Ω$ M and 0.07 in w. In conclusion, this divergence is most impactful for analyses of cosmological tensions, but its impact is mitigated when combining supernovae with other cross-cutting cosmological probes, such as the cosmic microwave background.

79 ASTRONOMY AND ASTROPHYSICS↗

Strategies to Obtain Reliable Energy Landscapes from Embedded Multireference Correlated Wavefunction Methods for Surface Reactions

Embedded correlated wavefunction (ECW) theory is a powerful tool for studying ground- and excited-state reaction mechanisms and associated energetics in heterogeneous catalysis. Several factors are important to obtaining reliable ECW energies, critically the construction of consistent active spaces (ASs) along reaction pathways when using a multireference correlated wavefunction (CW) method that relies on a subset of orbital spaces in the configuration interaction expansion to account for static electron correlation, e.g., complete AS self-consistent field theory, in addition to the adequate partitioning of the system into a cluster and environment, as well as the choice of a suitable basis set and number of states included in excited-state simulations. Here, in this work, we conducted a series of systematic studies to develop best-practice guidelines for ground- and excited-state ECW theory simulations, utilizing the decomposition of NH 3 on Pd(111) as an example. We determine that ECW theory results are relatively insensitive to cluster size, the aug-cc-pVDZ basis set provides an adequate compromise between computational complexity and accuracy, and that a fixed-clean-surface approximation holds well for the derivation of the embedding potential. Additionally, we demonstrate that a merging approach, which involves generating ASs from the molecular fragments at each configuration, is preferable to a creeping approach, which utilizes ASs from adjacent structures as an initial guess, for the generation of consistent potential energy curves involving open-d-shell metal surfaces, and, finally, we show that it is essential to include bands of excited states in their entirety when simulating excited-state reaction pathways.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Efficient Implementation of NOCI-MP2 Using the Resolution of the Identity Approximation with Application to Charged Dimers and Long C–C Bonds in Ethane Derivatives

An efficient implementation of the perturb-then-diagonalize nonorthogonal configuration interaction method with second-order Møller-Plesset perturbation theory (NOCI-MP2) is introduced. Relative to other low scaling multireference perturbation theories, NOCI-MP2 often requires a much smaller active space because of the use of nonorthogonal reference configurations. Reworking the NOCI-MP2 equations with the resolution of the identity (RI) approximation enables the method to have the same memory requirements and computational scaling as single reference RI-MP2. The working equations are extended to include single substitutions as required when the reference determinants do not satisfy the Hartree-Fock equations. A detailed computational algorithm is presented along with timings to establish the performance of the implementation. NOCI-MP2 is applied to the binding energy and charge resonance energy in dication and monocation π dimers, as well as didiamantane ethane, and hexaphenylethane. A well-defined set of nonorthogonal determinants are obtained using absolutely localized molecular orbitals (ALMOs), as solutions to the self-consistent field for molecular interactions (SCF-MI) equations corresponding to covalent and ionic determinants. Agreement with experimental information where available, and other multireference methods, is satisfactory, with the use of an 0.3 au level shift to guard against large MP2 amplitudes. For didiamantane ethane and hexaphenylethane, large dispersion forces help stabilize the molecules despite the steric repulsion. However, in the case of hexaphenylethane, the energy penalty from the geometric distortion of the fragments critically weakens the bond.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Conical Intersection Influences the Ground State Rearrangement of Fulvene to Benzene

The rearrangement of fulvene to benzene is believed to play an important role in the formation of soot during hydrocarbon combustion. Previous work has identified two possible mechanisms for the rearrangement-a unimolecular path and a hydrogen-atom-assisted, bimolecular path. Computational results to date have suggested that the unimolecular mechanism faces a barrier of about 74 kcal/mol, which makes it unable to compete with the bimolecular mechanism under typical combustion conditions. This computed barrier is about 10 kcal/mol higher than the experimental value, which is an unusually large discrepancy for modern electronic structure theory. In the present work, we have reinvestigated the unimolecular mechanism computationally, and we have found a second transition state that is approximately 10 kcal/mol lower in energy than the previously identified one and, therefore, in excellent agreement with the experimental value. The existence of two transition states for the same rearrangement arises because there is a conical intersection between the two lowest singlet states which occurs in the vicinity of the reaction coordinates. The two possible paths around the cone on the lower adiabatic surface give rise to the two distinct saddle points. The lower barrier for the unimolecular mechanism now makes it competitive with the bimolecular one, according to our calculations. In support of this conclusion, we have reanalyzed some previous experimental results on anisole pyrolysis, which leads to benzene as a significant product and have shown that the unimolecular and bimolecular mechanisms for fulvene ? benzene must be occurring competitively in that system. Finally, we have identified that similar conical intersections arise during the isomerizations of benzofulvene and isobenzofulvene to naphthalene.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

New Dimension in Ab Initio Electronic Structure Theory: Temperature, Pressure, and Chemical Potential

Ab initio electronic structure theory has transformed gas-phase molecular science with its predictive ability. In the attempt to bring such predictive ability to macroscopic systems and condensed matter, the theory must integrate quantum mechanics with statistical thermodynamics, so that thermodynamic functions such as free energy, internal energy, entropy, and chemical potentials are computed as functions of temperature in a systematically converging series of approximations. Here, a general, versatile strategy of elevating ab initio electronic structure theory to nonzero temperatures is introduced and discussed.

74 ATOMIC AND MOLECULAR PHYSICS↗

Advanced data analysis in inertial confinement fusion and high energy density physics

Bayesian analysis enables flexible and rigorous definition of statistical model assumptions with well-characterized propagation of uncertainties and resulting inferences for single-shot, repeated, or even cross-platform data. This approach has a strong history of application to a variety of problems in physical sciences ranging from inference of particle mass from multi-source high-energy particle data to analysis of black-hole characteristics from gravitational wave observations. The recent adoption of Bayesian statistics for analysis and design of high-energy density physics (HEDP) and inertial confinement fusion (ICF) experiments has provided invaluable gains in expert understanding and experiment performance. In this Review, we discuss the basic theory and practical application of the Bayesian statistics framework. We highlight a variety of studies from the HEDP and ICF literature, demonstrating the power of this technique. Due to the computational complexity of multi-physics models needed to analyze HEDP and ICF experiments, Bayesian inference is often not computationally tractable. Two sections are devoted to a review of statistical approximations, efficient inference algorithms, and data-driven methods, such as deep-learning and dimensionality reduction, which play a significant role in enabling use of the Bayesian framework. We provide additional discussion of various applications of Bayesian and machine learning methods that appear to be sparse in the HEDP and ICF literature constituting possible next steps for the community. We conclude by highlighting community needs, the resolution of which will improve trust in data-driven methods that have proven critical for accelerating the design and discovery cycle in many application areas.

47 OTHER INSTRUMENTATION↗

Artificial-intelligence-driven shot reduction in quantum measurement

Variational Quantum Eigensolver (VQE) provides a powerful solution for approximating molecular ground state energies by combining quantum circuits and classical computers. However, estimating probabilistic outcomes on quantum hardware requires repeated measurements (shots), incurring significant costs as accuracy increases. Optimizing shot allocation is thus critical for improving the efficiency of VQE. Current strategies rely heavily on hand-crafted heuristics requiring extensive expert knowledge. This paper proposes a reinforcement learning (RL)-based approach that automatically learns shot assignment policies to minimize total measurement shots while achieving convergence to the minimum of the energy expectation in VQE. The RL agent assigns measurement shots across VQE optimization iterations based on the progress of the optimization. This approach reduces VQE's dependence on static heuristics and human expertise. When the RL-enabled VQE is applied to a small molecule, a shot reduction policy is learned. The policy demonstrates transferability across systems and compatibility with other wavefunction Ansätze. In addition to these specific findings, this work highlights the potential of RL for automatically discovering efficient and scalable quantum optimization strategies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Loss-of-Coolant Accident Analysis of a High-Burnup Pressurized Water Reactor Core Design Using Gadolinia-Doped UO 2

High-burnup and extended enrichment fuels are of interest for extending the cycle lengths of pressurized water reactors from 18 to 24 months. Changes to the fuel design and core loading scheme have potentially significant safety implications due to power distribution effects, reduced thermal conductivity, and/or increased plenum pressures due to additional burnable poison loadings. Additionally, higher burnups result in increased material degradation and risk of fuel fragmentation, relocation, and dispersal (FFRD). A representative 24-month core design using gadolinia-doped UO 2 was analyzed for performance under large-break (LB) loss-of-coolant accident (LOCA) conditions using PARCS, RELAP5-3D, and BISON. Furthermore, all considered acceptance criteria were met, with no cases exceeding the 1477 K maximum cladding temperature or the post-quench ductility oxidation limit of 17% equivalent cladding reacted. Full-core FFRD susceptibility was estimated to be approximately 455 kg, though high uncertainties exist with current approaches for computing susceptibility. Undoped fuel rods are more likely to be limiting due to higher linear heat rates. Relatively high burnup and linear heat rate rods located in second batch assemblies are of greatest safety significance during LB LOCA for this high-burnup core design.

High burnup↗