Search NASA⌕ Search

SEARCH · Search NASA

Results for “dependable 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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Toward Quantitatively Accurate Calculation of the Redox-Associated Acid–Base and Ligand Binding Equilibria of Aquacobalamin

In redox processes in complex transition metal-containing species are often intimately associated with changes in ligand protonation states and metal coordination number. Moreover, a major challenge is therefore to develop consistent computational approaches for computing pH-dependent redox and ligand dissociation properties of organometallic species. Reduction of the Co center in the vitamin B12 derivative aquacobalamin can be accompanied by ligand dissociation, protonation, or both, making these properties difficult to compute accurately. We examine this challenge here by using density functional theory and continuum solvation to compute Co ligand binding equilibrium constants (Kon/off), pKas and reduction potentials for models of aquacobalamin in aqueous solution. We consider two models for cobalamin ligand coordination: the first follows the hexa, penta, tetra coordination scheme for Co III , Co II , and Co I species, respectively, and the second model features saturation of each vacant axial coordination site on Co II and Co I species with a single, explicit water molecule to maintain six directly interacting ligands or water molecules in each oxidation state. Comparing these two coordination schemes in combination with five dispersion-corrected density functionals, we find that the accuracy of the computed properties is largely independent of the scheme used, but including only a continuum representation of the solvent yields marginally better results than saturating the first solvation shell around Co throughout. PBE performs best, displaying balanced accuracy and superior performance overall, with RMS errors of 80 mV for seven reduction potentials, 2.0 log units for five pK a s and 2.3 log units for two log K on/off values for the aquacobalamin system. Furthermore, we find that the BP86 functional commonly used in corrinoid studies suffers from erratic behavior and inaccurate descriptions of Co axial ligand binding, leading to substantial errors in predicted pK a s and K on/off values. Finally, these findings demonstrate the effectiveness of the present approach for computing electrochemical and thermodynamic properties of a complex transition metal-containing cofactor.

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

Real-time data processing for serial crystallography experiments

We report the use of streaming data interfaces to perform fully online data processing for serial crystallography experiments, without storing intermediate data on disk. The system produces Bragg reflection intensity measurements suitable for scaling and merging, with a latency of less than 1 s per frame. Our system uses the CrystFEL software in combination with the ASAP::O data framework. In a series of user experiments at PETRA III, frames from a 16 megapixel Dectris EIGER2 X detector were searched for peaks, indexed and integrated at the maximum full-frame readout speed of 133 frames per second. The computational resources required depend on various factors, most significantly the fraction of non-blank frames ('hits'). The average single-thread processing time per frame was 242 ms for blank frames and 455 ms for hits, meaning that a single 96-core computing node was sufficient to keep up with the data, with ample headroom for unexpected throughput reductions. Further significant improvements are expected, for example by binning pixel intensities together to reduce the pixel count. We discuss the implications of real-time data processing on the `data deluge' problem from recent and future photon-science experiments, in particular on calibration requirements, computing access patterns and the need for the preservation of raw data.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Dependence of the kinetic energy absorption capacity of bistable mechanical metamaterials on impactor mass and velocity

Using an alternative mechanism to dissipation or scattering, bistable structures and mechanical metamaterials have shown promise for mitigating the detrimental effects of impact by reversibly locking energy into strained material. Herein, we extend prior works on impact absorption via bistable metamaterials to computationally explore the dependence of kinetic energy transmission on the velocity and mass of the impactor, with strain rates exceeding 10 2 s –1 . We observe a large dependence on both impactor parameters, ranging from significantly better to worse performance than a comparative linear material. We then correlate the variability in performance to solitary wave formation in the system and give analytical estimates of idealized energy absorption capacity under dynamic loading. In addition, we find a significant dependence on damping accompanied by a qualitative difference in solitary wave propagation within the system. The complex dynamics revealed in this study offer potential future guidance for the application of bistable metamaterials to applications including human and engineered system shock and impact protection devices.

42 ENGINEERING↗
DOE OSTI2017

CompNanoTox2015: novel perspectives from a European conference on computational nanotoxicology on predictive nanotoxicology

A first European Conference on Computational Nanotoxicology, CompNanoTox, was held in November 2015 in Benahavís, Spain with the objectives to disseminate and integrate results from the European modeling and database projects (NanoPUZZLES, ModENPTox, PreNanoTox, MembraneNanoPart, MODERN, eNanoMapper and EU COST TD1204 MODENA) as well as to create synergies within the European NanoSafety Cluster. This conference was supported by the COST Action TD1204 MODENA on developing computational methods for toxicological risk assessment of engineered nanoparticles and provided a unique opportunity for crossfertilization among complementary disciplines. The efforts to develop and validate computational models crucially depend on high quality experimental data and relevant assays which will be the basis to identify relevant descriptors. The ambitious overarching goal of this conference was to promote predictive nanotoxicology, which can only be achieved by a close collaboration between the computational scientists (e.g. database experts, modeling experts for structure, (eco) toxicological effects, performance and interaction of nanomaterials) and experimentalists from different areas (in particular toxicologists, biologists, chemists and material scientists, among others). The main outcome and new perspectives of this conference are summarized here.

nanoqsar↗
DOE OSTI2026

Machine learning approach for vibronically renormalized electronic band structures

Here, we present a machine learning (ML) method for efficient computation of vibrational thermal expectation values of physical properties from first principles. Our approach is based on the nonperturbative frozen phonon formulation in which stochastic Monte Carlo algorithm is employed to sample configurations of nuclei in a supercell at finite temperatures based on a first-principles phonon model. A deep-learning neural network is trained to accurately predict physical properties associated with sampled phonon configurations, thus bypassing the time-consuming ab initio calculations. To incorporate the point-group symmetry of the electronic system into the ML model, group-theoretical methods are used to develop a symmetry-invariant descriptor for phonon configurations in the supercell. We apply our ML approach to compute the temperature dependent electronic energy gap of silicon based on density functional theory (DFT). We show that, with less than a hundred DFT calculations for training the neural network model, an order of magnitude larger number of sampling can be achieved for the computation of the vibrational thermal expectation values. Our work highlights the promising potential of ML techniques for finite temperature first-principles electronic structure methods.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

A MultiGPU Performance-Portable Solution for Array Programming Based on Kokkos

Today, multiGPU nodes are widely used in high-performance computing and data centers. However, current programming models do not provide simple, transparent, and portable support for automatically targeting multiple GPUs within a node on application areas of array programming. In this paper, we describe a new application programming interface based on the Kokkos programming model to enable array computation on multiple GPUs in a transparent and portable way across both NVIDIA and AMD GPUs. We implement different variations of this technique to accommodate the exchange of stencils (array boundaries) among different GPU memory spaces, and we provide autotuning to select the proper number of GPUs, depending on the computational cost of the operations to be computed on arrays, that is completely transparent to the programmer. We evaluate our multiGPU extension on Summit (#5 TOP500), with six NVIDIA V100 Volta GPUs per node, and Crusher that contains identical hardware/software as Frontier (#1 TOP500), with four AMD MI250X GPUs, each with 2 Graphics Compute Dies (GCDs)for a total of 8 GCDs per node. We also compare the performance of this solution against the use of MPI + Kokkos, which is the cur-rent de facto solution for multiple GPUs in Kokkos. Our evaluation shows that the new Kokkos solution provides good scalability for many GPUs and a faster and simpler solution (from a programming productivity perspective) than MPI + Kokkos.

Valero Lara, Pedro↗
DOE OSTI2016

Ab initio calculations of the lattice parameter and elastic stiffness coefficients of bcc Fe with solutes

Here, we present an efficient methodology for computing solute-induced changes in lattice parameters and elastic stiffness coefficients Cij of single crystals using density functional theory. We also introduce a solute strain misfit tensor that quantifies how solutes change lattice parameters due to the stress they induce in the host crystal. Solutes modify the elastic stiffness coefficients through volumetric changes and by altering chemical bonds. We compute each of these contributions to the elastic stiffness coefficients separately, and verify that their sum agrees with changes in the elastic stiffness coefficients computed directly using fully optimized supercells containing solutes. Computing the two elastic stiffness contributions separately is more computationally efficient and provides more information on solute effects than the direct calculations. We compute the solute dependence of polycrystalline averaged shear and Young's moduli from the solute dependence of the single-crystal Cij. We then apply this methodology to substitutional Al, B, Cu, Mn, Si solutes and octahedral interstitial C and N solutes in bcc Fe. Comparison with experimental data indicates that our approach accurately predicts solute-induced changes in the lattice parameter and elastic coefficients. The computed data can be used to quantify solute-induced changes in mechanical properties such as strength and ductility, and can be incorporated into mesoscale models to improve their predictive capabilities.

36 MATERIALS SCIENCE↗
DOE OSTI2015

The Multi-Step CADIS method for shutdown dose rate calculations and uncertainty propagation

Shutdown dose rate (SDDR) analysis requires (a) a neutron transport calculation to estimate neutron flux fields, (b) an activation calculation to compute radionuclide inventories and associated photon sources, and (c) a photon transport calculation to estimate final SDDR. In some applications, accurate full-scale Monte Carlo (MC) SDDR simulations are needed for very large systems with massive amounts of shielding materials. However, these simulations are impractical because calculation of space- and energy-dependent neutron fluxes throughout the structural materials is needed to estimate distribution of radioisotopes causing the SDDR. Biasing the neutron MC calculation using an importance function is not simple because it is difficult to explicitly express the response function, which depends on subsequent computational steps. Furthermore, the typical SDDR calculations do not consider how uncertainties in MC neutron calculation impact SDDR uncertainty, even though MC neutron calculation uncertainties usually dominate SDDR uncertainty.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Initial-State Dependent Optimization of Controlled Gate Operations with Quantum Computer

There is no unique way to encode a quantum algorithm into a quantum circuit. With limited qubit counts, connectivity, and coherence times, a quantum circuit optimization is essential to make the best use of near-term quantum devices. We introduce a new circuit optimizer called AQCEL, which aims to remove redundant controlled operations from controlled gates, depending on initial states of the circuit. Especially, the AQCEL can remove unnecessary qubit controls from multi-controlled gates in polynomial computational resources, even when all the relevant qubits are entangled, by identifying zero-amplitude computational basis states using a quantum computer. As a benchmark, the AQCEL is deployed on a quantum algorithm designed to model final state radiation in high energy physics. For this benchmark, we have demonstrated that the AQCEL-optimized circuit can produce equivalent final states with much smaller number of gates. Moreover, when deploying AQCEL with a noisy intermediate scale quantum computer, it efficiently produces a quantum circuit that approximates the original circuit with high fidelity by truncating low-amplitude computational basis states below certain thresholds. Our technique is useful for a wide variety of quantum algorithms, opening up new possibilities to further simplify quantum circuits to be more effective for real devices.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Extending SEER for Extreme Heterogeneity

Heterogeneous and multi-device nodes are increasingly common in high-performance computing and data centers, yet existing programming models often lack simple, transparent, and portable support for these diverse architectures. The main contribution of this work is the development of novel SEER capabilities to address this challenge by providing a descriptive programming model that allows applications to seamlessly leverage heterogeneous nodes across various device types. SEER uses efficient memory management and can select the proper device[s] depending on the computational cost of the applications. This is completely transparent to the programmer, thereby providing a highly productive programming environment. Integrating extreme heterogeneity into the SEER library as shown with the use of NVIDIA and AMD GPUs simultaneously allows it to expand and exploit the performance possibilities. Our analysis based on the well-known Conjugate Gradient algorithm reports accelerations above 1.5 × on computationally demanding steps of such an algorithm by using both architectures simultaneously.

Teranishi, Keita [ORNL] (ORCID:0000000166472690)↗
DOE OSTI2015

Electronic spectra from TDDFT and machine learning in chemical space

Due to its favorable computational efficiency, time-dependent (TD) density functional theory (DFT) enables the prediction of electronic spectra in a high-throughput manner across chemical space. Its predictions, however, can be quite inaccurate. We resolve this issue with machine learning models trained on deviations of reference second-order approximate coupled-cluster (CC2) singles and doubles spectra from TDDFT counterparts, or even from DFT gap. We applied this approach to low-lying singlet-singlet vertical electronic spectra of over 20 000 synthetically feasible small organic molecules with up to eight CONF atoms. The prediction errors decay monotonously as a function of training set size. For a training set of 10 000 molecules, CC2 excitation energies can be reproduced to within ±0.1 eV for the remaining molecules. Analysis of our spectral database via chromophore counting suggests that even higher accuracies can be achieved. Based on the evidence collected, we discuss open challenges associated with data-driven modeling of high-lying spectra and transition intensities.

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

Real-Time Time-Dependent Nuclear–Electronic Orbital Approach: Dynamics beyond the Born–Oppenheimer Approximation

The quantum mechanical treatment of both electrons and nuclei is crucial in nona- diabatic dynamical processes such as proton-coupled electron transfer. The nuclear electronic orbital (NEO) method provides an elegant framework for including nuclear quantum effects beyond the Born-Oppenheimer approximation. To enable the study of nonequilibrium properties, we derive and implement a real-time NEO (RT-NEO) approach based on time-dependent Hatree-Fock or density functional theory, in which both the electronic and nuclear degrees of freedom are propagated in a time-dependent variational framework. Nuclear and electronic spectral features can be resolved from the time-dependent dipole moment computed using the RT-NEO method. The test cases demonstrate the dynamical interplay between the quantum nuclei and the electrons through vibronic coupling. Moreover, vibrational excitation in the RT-NEO approach is demonstrated by applying a resonant driving eld, and electronic excitation is demonstrated by simulating excited state intramolecular proton transfer. This work shows that the RT-NEO approach is a promising tool to study nonadiabatic quantum dynamical processes within a time-dependent variational description for the coupled electronic and nuclear degrees of freedom.

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

Modeling molecule-plasmon interactions using quantized radiation fields within time-dependent electronic structure theory

We present a combined cavity quantum electrodynamics/ab initio electronic structure approach for simulating plasmon-molecule interactions in the time domain. The simple Jaynes-Cummings-type model Hamiltonian typically utilized in such simulations is replaced with one in which the molecular component of the coupled system is treated in a fully ab initio way, resulting in a computationally efficient description of general plasmon-molecule interactions. Mutual polarization effects are easily incorporated within a standard ground-state Hartree-Fock computation, and time-dependent simulations carry the same formal computational scaling as real-time time-dependent Hartree-Fock theory. As a proof of principle, we apply this generalized method to the emergence of a Fano-like resonance in coupled molecule-plasmon systems; this feature is quite sensitive to the nanoparticle-molecule separation and the orientation of the molecule relative to the polarization of the external electric field.

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

A Decomposition-Based Learn-To-Optimize Approach with Feasibility Layer Assistance for Sub-Hourly Unit Commitment

Sub-hourly unit commitment (UC) with 15-min intervals is gaining significant attention as a way to respond rapidly to the fluctuations in electricity supply and demand introduced by renewable resources. However, the increased temporal resolution and complex inter-temporal dependencies pose substantial computational challenges for traditional optimization methods. To this end, this paper explores a decomposition-based learn-to-optimize approach. Building on recent advances in machine learning, our method revisits the long- overlooked Lagrangian relaxation framework, which is a classical decomposition technique that enables tractable subproblem solving. These smaller subproblems are inherently well-suited for machine learning, as their reduced dimensionality and structural regularity allow predictive models to efficiently learn and generalize solution patterns. We thus propose a generic predictive model, which embeds Gated Recurrent Units (GRUs) and Attention in the encoder-decoder structure, and integrate a rule-based feasibility layer to capture temporal dependencies, reduce training effort, and improve feasibility w.r.t. unit-level constraints. Our method has been validated on the IEEE 118-bus system, demonstrating promising performance in solving sub-hourly UC problems efficiently and feasibly.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Dipole Moment Calculations Using Multiconfiguration Pair-Density Functional Theory and Hybrid Multiconfiguration Pair-Density Functional Theory

The dipole moment is the molecular property that most directly indicates molecular polarity. The accuracy of computed dipole moments depends strongly on the quality of the calculated electron density, and the breakdown of single-reference methods for strongly correlated systems can lead to poor predictions of the dipole moments in those cases. Here, we derive the analytical expression for obtaining the electric dipole moment by multiconfiguration pair density functional theory (MCPDFT), and we assess the accuracy of MC-PDFT for predicting dipole moments at equilibrium and nonequilibrium geometries. We show that MC-PDFT dipole moment curves have reasonable behavior even for stretched geometries, and they significantly improve upon the CASSCF results by capturing more electron correlation. Herein, the analysis of a dataset consisting of 18 first-row transition metal diatomics and 6 main-group polyatomic molecules with multireference character suggests that MC-PDFT and its hybrid extension (HMC-PDFT) perform comparably to CASPT2 and MRCISD+Q methods and have a mean unsigned deviation of 0.2–0.3 D with respect to the best available dipole moment reference values. We explored the dependence of the predicted dipole moments upon the choice of the on-top density functional and active space, and we recommend the tPBE and hybrid tPBE0 on-top choices for the functionals combined with the moderate correlated participating orbital scheme for selecting the active space. With these choices, the mean unsigned deviations (in debyes) of the calculated equilibrium dipole moments from the best estimates are 0.77 for CASSCF, 0.29 for MC-PDFT, 0.24 for HMCPDFT, 0.28 for CASPT2, and 0.25 for MRCISD+Q. These results are encouraging because the computational cost of MC-PDFT or HMC-PDFT is largely reduced compared to the CASPT2 and MRCISD+Q methods.

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

Large change of interlayer vibrational coupling with stacking in Mo 1-x W x Te 2

Stacking variations in quasi-two-dimensional materials can have an important influence on material properties, such as changing the topology of the band structure. Unfortunately, the weakness of van der Waals (vdW) interactions makes it difficult to compute the stacking dependence of properties, and even in a material as simple as graphite the stacking energetics remain unclear. Mo 1 - x W x Te 2 is a material in which three differently stacked phases are conveniently accessible by temperature changes: 1 T ' , T d * , and the reported Weyl semimetal phase T d . Additionally, the transitions proceed via layer sliding, and the corresponding interlayer shear mode (ISM) is relevant not just for the stacking energetics but also for understanding the relationship between Weyl physics and structural changes. However, the interlayer interactions of Mo 1 - x W x Te 2 are not well understood, with wide variation in computed properties. We report inelastic neutron scattering of the ISM in a Mo 0.91 W 0.09 Te 2 crystal. The ISM energies are generally consistent with the linear chain model, as expected given the weak interlayer interaction, though there are some discrepancies from predicted intensities. However, the interlayer force constants K x in the T d * and 1 T ' phases are substantially weaker than that of T d at 75(3) and 83(3)%, respectively. Considering that the relative positioning of atoms in neighboring layers is approximately the same regardless of overall stacking, our results suggest that longer-range influences, such as stacking-induced electronic band-structure changes, may be responsible for the substantial change in the interlayer vibrational coupling and thus the C 55 elastic constant. These findings should elucidate the stacking energetics of Mo 1 - x W x Te 2 and other vdW layered materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗