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 163 records · Page 9

Vibronic coupling in energy transfer dynamics and two-dimensional electronic–vibrational spectra

We introduce a heterodimer model in which multiple mechanisms of vibronic coupling and their impact on energy transfer can be explicitly studied. We consider vibronic coupling that arises through either Franck-Condon activity in which each site in the heterodimer has a local electron-phonon coupling or Herzberg-Teller activity in which the transition dipole moment coupling the sites has an explicit vibrational mode-dependence. We have computed two-dimensional electronic-vibrational (2DEV) spectra for this model while varying the magnitude of these two effects and find that 2DEV spectra contain static and dynamic signatures of both types of vibronic coupling. Franck-Condon activity emerges through a change in the observed excitonic structure, while Herzberg-Teller activity is evident in the appearance of significant side-band transitions that mimic the lower-energy excitonic structure. A comparison of quantum beating patterns obtained from analysis of the simulated 2DEV spectra shows that this technique can report on the mechanism of energy transfer, elucidating a means of experimentally determining the role of specific vibronic coupling mechanisms in such processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

FlameNEST: explicit profile likelihoods with the Noble Element Simulation Technique

We present FlameNEST, a framework providing explicit likelihood evaluations in noble element particle detectors using data-driven models from the Noble Element Simulation Technique. FlameNEST provides a way to perform statistical analyses on real data with no dependence on large, computationally expensive Monte Carlo simulations by evaluating the likelihood on an event-by-event basis using analytic probability elements convolved together in a single TensorFlow multiplication. Furthermore, this robust framework creates opportunities for simple inter-collaboration analyses which will be fundamental for the future of experimental dark matter physics.

47 OTHER INSTRUMENTATION↗

Reinforcement learning pulses for transmon qubit entangling gates

The utility of a quantum computer is highly dependent on the ability to reliably perform accurate quantum logic operations. For finding optimal control solutions, it is of particular interest to explore model-free approaches, since their quality is not constrained by the limited accuracy of theoretical models for the quantum processor—in contrast to many established gate implementation strategies. In this work, we utilize a continuous control reinforcement learning algorithm to design entangling two-qubit gates for superconducting qubits; specifically, our agent constructs cross-resonance and CNOT gates without any prior information about the physical system. Using a simulated environment of fixed-frequency fixed-coupling transmon qubits, we demonstrate the capability to generate novel pulse sequences that outperform the standard cross-resonance gates in both fidelity and gate duration, while maintaining a comparable susceptibility to stochastic unitary noise. We further showcase an augmentation in training and input information that allows our agent to adapt its pulse design abilities to drifting hardware characteristics, importantly, with little to no additional optimization. Our results exhibit clearly the advantages of unbiased adaptive-feedback learning-based optimization methods for transmon gate design.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Dynamical correlations and nonequilibrium sum rules in photodoped Hubbard ladders

Using matrix product state techniques we study the nonequilibrium dynamical response of the half-filled Hubbard ladder when subject to an optical pump. Optical pumping offers a way of producing and manipulating new strongly correlated phenomena by suppressing existing magnetic correlations. The ladder allows the effects of pump directionality to be investigated, and compared to a single chain it has strong spin-charge coupling and a fully gapped excitation spectrum, promising different nonequilibrium physics. We compute time-dependent correlations, including the nonequilibrium dynamical structure factors for spin and charge. By deriving a combined spin-charge sum rule that applies both in and out-of-equilibrium, we show that spectral weight is pumped directly from the antiferromagnetic spin response into a low energy 𝜔 ∼ 0 charge response below the Mott gap. The transfer of weight is pump direction dependent: pumping directed along the legs disrupts magnetic correlations more than pumping in the rung-direction, even if the post-pump energy density is similar. The charge correlation length is dramatically enhanced by the pump, whilst the spin correlations are most strongly suppressed at nearest- and next-nearest-neighbor spacings. After the pump the system is in a nonthermal correlated metallic state, with gapless charge excitations and approximately equal spin and charge correlation lengths, emphasizing the importance of treating these degrees of freedom on an equal footing in nonequilibrium systems.

36 MATERIALS SCIENCE↗

Vortex precession in trapped superfluids from effective field theory

We apply a recently developed effective string theory for vortex lines to the case of two-dimensional trapped superfluids. We do not assume a perturbative microscopic description for the superfluid, but only a gradient expansion for the long-distance hydrodynamical description and for the trapping potential. For any regular trapping potential, we compute the spatial dependence of the superfluid density and the orbital frequency and trajectory of an off-center vortex. Our results are fully relativistic and in the nonrelativistic limit reduce to known results based on the Gross-Pitaevskii model. In our formalism, the leading effect in the nonrelativistic limit arises from two simple Feynman diagrams in which the vortex interacts with the trapping potential through the exchange of hydrodynamical modes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

OpenStudio® HPXML workflow [SWR-25-13]

OpenStudio-HPXML allows running residential EnergyPlus™ simulations using an HPXML file for the building description. It is intended to be used by user interfaces or other automated software workflows that automatically produce the HPXML file. OpenStudio-HPXML can accommodate a wide range of different building technologies and geometries. End-to-end simulations typically run in 3-10 seconds, depending on complexity, computer platform and speed, etc. For more information on running simulations, generating HPXML files with the appropriate inputs to run EnergyPlus, etc., please visit the documentation linked below. https://openstudio-hpxml.readthedocs.io/en/latest

Horowitz, Scott↗

Imaging Bragg Edge Analysis TooLs for Engineering Structures (iBeatles)

The Spallation Neutron Source (SNS) at Oak Ridge National Laboratory (ORNL) provides pulsed neutrons with energies varying from epithermal to cold. In preparation for VENUS, the neutron imaging beamline to be located at beam port 10, we have performed a series of experiments focused on wavelength-dependent radiography and computed tomography for a broad range of applications, from materials science to biological tissues.One of the time-of-flight (TOF) techniques that is of interest to the scientific community is the 2-dimensional mapping of phases and average crystalline plane orientation in samples both ex-situ and during applied stresses such as tensile loading and heating. This technique is known as Bragg edgeimaging and relies on the identification of changes of transmission values, fitting of the edge to measure its displacement, and thus identify the shift in lattice parameter due to stresses. One of the challenges of TOF imaging measurements is the amount of data and the inability to observe Bragg edge shifts in real time during an experiment. Thus, we have been focusing on creating a Python-based interface that allows fast data processing and instantaneous mapping and fitting of the Bragg edges, and their evolution through time. Python libraries and Jupyter notebooks have been implemented to facilitate decision making during an experiment. The advantage of the notebooks is the possibility to guide an experiment as they can quickly process and display Bragg edge data. These notebooks can be used independently, or can be combined in a Python Graphical User Interface (GUI) tool called iBeatles. This interface permits visualization and fitting of the Bragg edges, and ultimately back-projects the fitting results onto the radiographs to display a strain map. Assuming data collection has sufficient statistics, the strain mapping analysis can be performed on a pixel-by-pixel basis. This development is a step forward toward a better user experience at the future VENUS beamline in terms of live feedback and productivity. Analysis that used to take days of switching between different applications can now be done in minutes within the

Bilheux, JeanChristophe [Oak Ridge National Labora↗

ESnet-JLab FPGA Accelerated Transport (data plane) [EJFAT (udplb)] v1.0

The ESnet-JLab FPGA Accelerated Transport system is a solution for streaming high-speed scientific measurement data from Data Acquisition Systems (DAQs) to high-performance computing facilties. It is generally compatible with many science workflows, and makes no assumptions about the specifics of any particular experiment. This program (udplb) implements the data plane portion of the EJFAT system. It is an FPGA design that rewrites and forwards data packets from a UDP-based scientific workflow to high-performance compute nodes. It depends on another program (udplbd, disclosed separately) to implement the control system.

Bengough, Peter [Malleable Networks, Inc.]↗

Experimental energy levels and partition function of the {sup 12}C{sub 2} molecule

The carbon dimer, the {sup 12}C{sub 2} molecule, is ubiquitous in astronomical environments. Experimental-quality rovibronic energy levels are reported for {sup 12}C{sub 2}, based on rovibronic transitions measured for and among its singlet, triplet, and quintet electronic states, reported in 42 publications. The determination utilizes the Measured Active Rotational-Vibrational Energy Levels (MARVEL) technique. The 23,343 transitions measured experimentally and validated within this study determine 5699 rovibronic energy levels, 1325, 4309, and 65 levels for the singlet, triplet, and quintet states investigated, respectively. The MARVEL analysis provides rovibronic energies for six singlet, six triplet, and two quintet electronic states. For example, the lowest measurable energy level of the a{sup 3}Π{sub u} state, corresponding to the J = 2 total angular momentum quantum number and the F {sub 1} spin-multiplet component, is 603.817(5) cm{sup −1}. This well-determined energy difference should facilitate observations of singlet–triplet intercombination lines, which are thought to occur in the interstellar medium and comets. The large number of highly accurate and clearly labeled transitions that can be derived by combining MARVEL energy levels with computed temperature-dependent intensities should help a number of astrophysical observations as well as corresponding laboratory measurements. The experimental rovibronic energy levels, augmented, where needed, with ab initio variational ones based on empirically adjusted and spin–orbit coupled potential energy curves obtained using the Duo code, are used to obtain a highly accurate partition function, and related thermodynamic data, for {sup 12}C{sub 2} up to 4000 K.

79 ASTRONOMY AND ASTROPHYSICS↗

Neural units with time-dependent functionality

We show that the time-resolved dynamics of an underdamped harmonic oscillator can be used to do multifunctional computation, performing distinct computations at distinct times within a single dynamical trajectory. We consider the amplitude of an oscillator whose inputs influence its frequency. The activity of the oscillator at fixed times is a nonmonotonic function of its inputs, so it can solve problems such as XOR that are not linearly separable. The activity of the oscillator at fixed input is a nonmonotonic function of time, so it is multifunctional in a temporal sense, and able to carry out distinct nonlinear computations at distinct times within the same dynamical trajectory. We show that a single oscillator, observed at different times, can act as all of the elementary logic gates and perform binary addition, the latter usually implemented in hardware using five logic gates. We show that a set of n oscillators, observed at different times, can perform an arbitrary number of analog-to-n-bit digital conversions. We also show that oscillators can be trained by gradient descent to perform distinct classification tasks at distinct times. Computing with time-dependent functionality can be done in or out of equilibrium, and suggests a way of reducing the number of parameters or devices required to do nonlinear computations.

97 MATHEMATICS AND COMPUTING↗

Performance Analysis of Scientific Computing Workloads on Trusted Execution Environments

Scientific computing sometimes involves computation on sensitive data. Depending on the data and the execution environment, the HPC (high-performance computing) user or data provider may require confidentiality and/or integrity guarantees. To study the applicability of hardware-based trusted execution environments (TEEs) to enable secure scientific computing, we deeply analyze the performance impact of AMD SEV and Intel SGX for diverse HPC benchmarks including traditional scientific computing, machine learning, graph analytics, and emerging scientific computing workloads. We observe three main findings: 1) SEV requires careful memory placement on large scale NUMA machines (1x -3.4x slowdown without and 1x -1.15x slowdown with NUMA aware placement), 2) virtualization - a prerequisite for SEV - results in performance degradation for workloads with irregular memory accesses and large working sets (1x -4x slowdown compared to native execution for graph applications) and 3) SGX is inappropriate for HPC given its limited secure memory size and inflexible programming model (1.2x -126x slowdown over unsecure execution). Finally, we discuss forthcoming new TEE designs and their potential impact on scientific computing.

97 MATHEMATICS AND COMPUTING↗

Size-dependent error of the density functional theory ionization potential in vacuum and solution

Density functional theory is often the method of choice for modeling the energetics of large molecules and including explicit solvation effects. It is preferable to use a method that treats systems of different sizes and with different amounts of explicit solvent on equal footing. However, recent work suggests that approximate density functional theory has a size-dependent error in the computation of the ionization potential. We here investigate the lack of size-intensivity of the ionization potential computed with approximate density functionals in vacuum and solution. We show that local and semi-local approximations to exchange do not yield a constant ionization potential for an increasing number of identical isolated molecules in vacuum. Instead, as the number of molecules increases, the total energy required to ionize the system decreases. Rather surprisingly, we find that this is still the case in solution, whether using a polarizable continuum model or with explicit solvent that breaks the degeneracy of each solute, and we find that explicit solvent in the calculation can exacerbate the size-dependent delocalization error. We demonstrate that increasing the amount of exact exchange changes the character of the polarization of the solvent molecules; for small amounts of exact exchange the solvent molecules contribute a fraction of their electron density to the ionized electron, but for larger amounts of exact exchange they properly polarize in response to the cationic solute. In vacuum and explicit solvent, the ionization potential can be made size-intensive by optimally tuning a long-range corrected hybrid functional.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Theoretical study of orbital ordering induced structural phase transition in iron pnictides

We attribute the structural phase transition (SPT) in the parent compounds of the iron pnictides to orbital ordering. Due to anisotropy of the d{sub xz} and d{sub yz} orbitals in the xy plane, orbital ordering makes the orthorhombic structure more favorable and thus inducing the SPT. We consider a one band model Hamiltonian consisting of first and second-nearest-neighbor hopping of the electrons. We introduce Jahn-Tellar (JT) distortion in the system arising due to the orbital ordering present in this system. We calculate the electron Green’s function by using Zuvareb’s Green’s function technique and hence calculate an expression for the temperature dependent lattice strain which is computed numerically and self-consistently. The temperature dependent electron specific heat is calculated by minimizing the free energy of the system. The lattice strain is studied by varying the JT coupling and elastic constant of the system. The structural anomaly is studied through the electron occupation number and the specific heat by varying the physical parameters like JT coupling, lattice constant, chemical potential and hopping integrals of the system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Ab Initio Reaction Kinetics of CH 3 O$\dot{C}$(=O) and $\dot{C}$H 2 OC(=O)H Radicals

The dissociation and isomerization kinetics of the methyl ester combustion intermediates methoxycarbonyl radical (CH3O$\dot{C}$(=O)) and (formyloxy)methyl radical ($\dot{C}$H2OC(=O)H) are investigated theoretically using high-level ab initio methods and Rice–Ramsperger–Kassel–Marcus (RRKM)/master equation (ME) theory. Geometries obtained at the hybrid density functional theory (DFT) and coupled cluster singles and doubles with perturbative triples correction (CCSD(T)) levels of theory are found to be similar. We employ high-level ab initio wave function methods to refine the potential energy surface: CCSD(T), multireference singles and doubles configuration interaction (MRSDCI) with the Davidson–Silver (DS) correction, and multireference averaged coupled-pair functional (MRACPF2) theory. MRSDCI+DS and MRACPF2 capture the multiconfigurational character of transition states (TSs) and predict lower barrier heights than CCSD(T). The temperature- and pressure-dependent rate coefficients are computed using RRKM/ME theory in the temperature range 300–2500 K and a pressure range of 0.01 atm to the high-pressure limit, which are then fitted to modified Arrhenius expressions. Dissociation of CH3O$\dot{C}$(=O) to $\dot{C}$H3 and CO2 is predicted to be much faster than dissociating to CH3$\dot{O}$ and CO, consistent with its greater exothermicity. Isomerization between CH3O$\dot{C}$(=O) and $\dot{C}$H2OC(=O)H is predicted to be the slowest among the studied reactions and rarely happens even at high temperature and high pressure, suggesting the decomposition pathways of the two radicals are not strongly coupled. The predicted rate coefficients and branching fractions at finite pressures differ significantly from the corresponding high-pressure-limit results, especially at relatively high temperatures. Finally, because it is one of the most important CH3$\dot{O}$ removal mechanisms under atmospheric conditions, the reaction kinetics of CH3$\dot{O}$ + CO was also studied along the PES of CH3O$\dot{C}$(=O); the resulting kinetics predictions are in remarkable agreement with experiments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Shape‐Dependent Multi‐Weight Magnetic Artificial Synapses for Neuromorphic Computing

Abstract In neuromorphic computing, artificial synapses provide a multi‐weight (MW) conductance state that is set based on inputs from neurons, analogous to the brain. Herein, artificial synapses based on magnetic materials that use a magnetic tunnel junction (MTJ) and a magnetic domain wall (DW) are explored. By fabricating lithographic notches in a DW track underneath a single MTJ, 3–5 stable resistance states that can be repeatably controlled electrically using spin‐orbit torque are achieved. The effect of geometry on the synapse behavior is explored, showing that a trapezoidal device has asymmetric weight updates with high controllability, while a rectangular device has higher stochasticity, but with stable resistance levels. The device data is input into neuromorphic computing simulators to show the usefulness of application‐specific synaptic functions. Implementing an artificial neural network (NN) applied to streamed Fashion‐MNIST data, the trapezoidal magnetic synapse can be used as a metaplastic function for efficient online learning. Implementing a convolutional NN for CIFAR‐100 image recognition, the rectangular magnetic synapse achieves near‐ideal inference accuracy, due to the stability of its resistance levels. This work shows MW magnetic synapses are a feasible technology for neuromorphic computing and provides design guidelines for emerging artificial synapse technologies.

Leonard, Thomas↗

The Role of Nuclear Data Sensitivities in Prompt α-Eigenvalue Predictions of Delayed Critical Benchmarks

Alpha (α) eigenvalues, which describe the logarithmic time derivative of the neutron population in a multiplying system, are integral to time-dependent behavior and diagnostic applications. However, uncertainties in the evaluated nuclear data can significantly impact the accuracy of transport simulations for such quantities. This work explores the use of machine learning models to predict two key outputs, α-eigenvalues and keff bias, using input features derived from α-eigenvalue sensitivities to nuclear data. The criticality safety benchmark models used in this study come from the International Handbook of Evaluated Criticality Safety Benchmark Experiments. Three models, random forest, XGBoost, and NGBoost, are trained on both energy-resolved and energy-summed α sensitivities. For the α-eigenvalue bias prediction, NGBoost achieved the highest R 2 (0.9476) using energy-resolved features, while XGBoost performed best using summed sensitivities. In contrast, when predicting the keff bias, all the models showed moderate predictive capability (best R 2 ≈ 0.72), as the mapping from the static α-sensitivities to the static keff bias was less direct. SHAP (SHapley Additive exPlanations) analysis was used to interpret the model predictions. Across both prediction tasks, the features associated with neutron capture [H-1 (n, γ)], uranium scattering reactions (such as 235 U elastic/inelastic), and actinide capture/fission reactions (such as 239 Pu and 234 U) were consistently identified as the most impactful. This highlights the key role of specific nuclear reactions and energy ranges in shaping both time-dependent and steady-state criticality behavior. These results demonstrated that α-sensitivities, despite being computed for time-dependent metrics, can provide valuable insights for predicting both α-eigenvalues and the keff bias. Moreover, machine learning models offer a promising pathway for uncovering important nuclear data dependencies and guiding future data evaluation efforts.

Nuclear data↗

Efficient loading of reduced data ensembles produced at ORNL SNS/HFIR neutron time-of-flight facilities

We present algorithmic improvements to the loading operations of certain reduced data ensembles produced from neutron scattering experiments at Oak Ridge National Laboratory (ORNL) facilities. Ensembles from multiple measurements are required to cover a wide range of the phase space of a sample material of interest. They are stored using the standard NeXus schema on individual HDF5 files. This makes it a scalability challenge, as the number of experiments stored increases in a single ensemble file. The present work follows up on our previous efforts on data management algorithms, to address identified input output (I/O) bottlenecks in Mantid, an open-source data analysis framework used across several neutron science facilities around the world. We reuse an in-memory binary-tree metadata index that resembles data access patterns, to provide a scalable search and extraction mechanism. In addition, several memory operations are refactored and optimized for the current common use cases, ranging most frequently from 10 to 180, and up to 360 separate measurement configurations. Results from this work show consistent speed ups in wall-clock time on the Mantid LoadMD routine, ranging from 19% to 23% on average, on ORNL production computing systems. The latter depends on the complexity of the targeted instrument-specific data and the system I/O and compute variability for the shared computational resources available to users of ORNL’s Spallation Neutron Source (SNS) and the High Flux Isotope Reactor (HFIR) instruments. Nevertheless, we continue to highlight the need for more research to address reduction challenges as experimental data volumes, user time and processing costs increase.

Godoy, William↗

Autoignition of Premixed Liquefied Petroleum Gas in a Rapid Compression Machine: Experimental Results and Chemical Kinetic Mechanism Reduction

Liquefied petroleum gas (LPG) has many properties that make it an attractive alternative fuel such as lower cost than conventional fuels and an established distribution infrastructure. The development of high efficiency, spark ignited LPG engines is currently limited by engine knock and misfire. The knock and misfire limits are further complicated by the wide range of chemical reactivity in LPG, particularly in international markets. In this study, a rapid compression machine (RCM) was used to characterize the effects of variation in LPG fuel reactivity, equivalence ratio, and exhaust gas recirculation (EGR) on autoignition of LPG/oxidizer/inert/EGR blends. Experiments were conducted with 100% propane (C3H8) and blends of propane with propene, ethane, isobutane, and n-butane. EGR was simulated with mixtures of Ar, CO2, CO, and NO at substitution percentages from 0 to 30 mass percent. Equivalence ratio was varied from 0.75 to 1.5. Ignition delay period under homogeneous autoignition conditions was measured at compressed pressures and temperatures of 23 to 25 bar and 701 to 921 K, respectively. Zero-dimensional simulations of the RCM experiments were performed using Chemkin with several detailed chemical kinetic mechanisms to determine their suitability at predicting ignition delay periods. Multiple reduced chemical kinetic mechanisms were created from the NUIGMech1.1 mechanism to determine the optimal balance between accuracy and computational efficiency for future three-dimensional, time-dependent spark-ignited engine computations.

42 ENGINEERING↗