Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reversible computation”

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 253 records · Page 14

Carboxyl intermediate formation via an in situ-generated metastable active site during water-gas shift reaction

Definitive experimental proof for catalytic pathways and active sites during the low-temperature water-gas shift reaction remains elusive. Herein, we combine spectroscopic, kinetic, and computational analyses to address the decades-long mechanistic controversy by studying the reverse water-gas shift over Pd/Al2O3. Isotopic transient kinetic analysis clearly established the minor role of the formate (HCOO) intermediate. Moreover, hydrogen titration experiments unambiguously confirmed the intermediacy of carboxyl (COOH) for the first time. The ability to decouple the parallel formate and carboxyl pathways led to the identification of a distinct active site exhibiting regio- and chemoselective H-addition to CO2 to yield carboxyl. The metastable active site is formed in situ, resulting in hydroxylation of the metal-support interface and electronic restructuring of Pd via reduced coordination. Parallel atomistic simulations of the active site electronic structure and mechanistic landscape provided a framework consistent with experimental observations. Our study highlights the dynamic creation of a coordinatively unsaturated metal site caused by substrate adsorption on an adjacent support site.

Nelson, Nicholas C.↗

Near-equilibrium measurements of nonequilibrium free energy

A central endeavor of thermodynamics is the measurement of free energy changes. Regrettably, although we can measure the free energy of a system in thermodynamic equilibrium, typically all we can say about the free energy of a nonequilibrium ensemble is that it is larger than that of the same system at equilibrium. Herein, we derive a formally exact expression for the probability distribution of a driven system, which involves path ensemble averages of the work over trajectories of the time-reversed system. From this we find a simple near-equilibrium approximation for the free energy in terms of an excess mean time-reversed work, which can be experimentally measured on real systems. With analysis and computer simulation, we demonstrate the accuracy of our approximations for several simple models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Efficient 3-D velocity model building using joint inline and crossline plane-wave wave-equation migration velocity analyses

SUMMARY Wave-equation migration velocity analysis (WEMVA) is an image-domain inversion method for velocity model building. Automatic plane-wave WEMVA (PWEMVA) calculates the moveouts of plane-wave common-image gathers (CIGs) by searching a best-fitting parabola with semblance analysis and backprojects residual CIG moveouts into wavefield wave paths with a reflection tomographic kernel. However, 3-D PWEMVA is very computationally expensive because 3-D reflection tomographic inversion requires at least five 3-D reverse-time migrations per iteration and stores two types of source wavefields at model boundaries. We develop a joint inline and crossline PWEMVA method for efficient 3-D velocity model building. We alternatively implement the inline and crossline PWEMVAs with a constraint for each other, in which we iteratively construct the 3-D velocity model update through 1-D spline interpolation of 2-D gradients. The inline and crossline joint inversion is practical since PWEMVA only inverts for low-wavenumber velocity perturbations along wave paths, and the method can take less than 1 per cent of the computational cost of full 3-D PWEMVA. To construct unaliased plane waves for our joint inline and crossline PWEMVA, we develop a 3-D data interpolation method in the frequency–wavenumber (FK) domain to recover regularly and randomly missing traces. The method minimizes the misfit on sufficiently localized data subsets with iterative optimal step lengths and a gradient preconditioner that iteratively selects dominant dips along different azimuths. In numerical experiments, we use a 3-D synthetic seismic data set and a land 3-D field seismic data set acquired at the Farnsworth CO2-EOR (enhanced oil recovery) field to demonstrate the efficacy of our velocity model building and data interpolation methods.

Liu, Xuejian↗

Voltage control of unidirectional anisotropy in ferromagnet-multiferroic system

Demonstration of ultralow energy switching mechanisms is imperative for continued improvements in computing devices. Ferroelectric (FE) and multiferroic (MF) order and their manipulation promise an ideal combination of state variables to reach attojoule range for logic and memory (i.e., ~30× lower switching energy than nanoelectronics). In BiFeO 3 (BFO), the coupling between the antiferromagnetic (AFM) and FE order is robust at room temperature, scalable in voltage, stabilized by the FE order, and can be integrated into a fabrication process for a beyond-CMOS (complementary metal-oxide semiconductor) era. The presence of the AFM order and a canted magnetic moment in this system causes exchange interaction with a ferromagnet such as Co 0.9 Fe 0.1 or La 0.7 Sr 0.3 MnO 3 . Previous research has shown that exchange coupling (uniaxial anisotropy) can be controlled with an electric field. However, voltage modulation of unidirectional anisotropy, which is preferred for logic and memory technologies, has not yet been demonstrated. Here, we present evidence for electric field control of exchange bias of laterally scaled spin valves that is exchange coupled to BFO at room temperature. We show that the exchange bias in this bilayer is robust, electrically controlled, and reversible. We anticipate that magnetoelectricity at these scaled dimensions provides a powerful pathway for computing beyond modern nanoelectronics by enabling a new class of nonvolatile, ultralow energy computing elements.

36 MATERIALS SCIENCE↗

ELECTROMAGNETIC EMISSION FROM LONG-LIVED BINARY NEUTRON STAR MERGER REMNANTS. II. LIGHT CURVES AND SPECTRA

Recent observations indicate that in a large fraction of binary neutron star (BNS) mergers a long-lived neutron star (NS) may be formed rather than a black hole. Unambiguous electromagnetic (EM) signatures of such a scenario would strongly impact our knowledge on how short gamma-ray bursts (SGRBs) and their afterglow radiation are generated. Furthermore, such EM signals would have profound implications for multimessenger astronomy with joint EM and gravitational-wave (GW) observations of BNS mergers, which will soon become reality thanks to the ground-based advanced LIGO/Virgo GW detector network. Here we explore such EM signatures based on the model presented in a companion paper, which provides a self-consistent evolution of the post-merger system and its EM emission up to ∼10{sup 7} s. Light curves and spectra are computed for a wide range of post-merger physical properties. We present X-ray afterglow light curves corresponding to the “standard” and the “time-reversal” scenario for SGRBs (prompt emission associated with the merger or with the collapse of the long-lived NS). The light curve morphologies include single and two-plateau features with timescales and luminosities that are in good agreement with Swift observations. Furthermore, we compute the X-ray signal that should precede the SGRB in the time-reversal scenario, the detection of which would represent smoking-gun evidence for this scenario. Finally, we find a bright, highly isotropic EM transient peaking in the X-ray band at ∼10{sup 2}–10{sup 4} s after the BNS merger with luminosities of L{sub X} ∼ 10{sup 46}–10{sup 48} erg s{sup −1}. This signal represents a very promising EM counterpart to the GW emission from BNS mergers.

79 ASTRONOMY AND ASTROPHYSICS↗

The Memory Scaling of Reverse-Mode Differentiation in Particle Accelerator Simulations with Space Charge

The recent development of differentiable simulation codes for particle accelerators has enabled gradient-based workflows that promise finer control and more realistic modeling of accelerator facilities. However, when using reverse-mode automatic differentiation, the memory usage continuously increases during the simulation, and can potentially exceed the available hardware memory - especially when costly space charge computation is included. To study the memory requirements for differentiable simulations, we have implemented space charge in Cheetah, a PyTorch-based beam tracking code that supports reverse-mode differentiation. We find that the memory usage for reverse-mode differentiation grows linearly with the number of macroparticles and cells, and that it is proportional to the number of space charge kicks involved in the simulation. This general scaling can be used to evaluate whether a given differentiable simulation is feasible given hardware memory constraints.

Dhamrait, Arjun↗

Buried Dirac Points in Quantum Spin Hall Insulators: Implications for Majorana Kramers Pair-Based Quantum Computing

Quantum spin Hall insulators (QSHIs) host helical electronic edge states that are protected from backscattering due to time-reversal symmetry (TRS). Despite considerable work investigating QSHI edge states, there is still an open question about their unexpected resilience to large magnetic fields where TRS is undoubtedly broken. In this work, we investigate the transport properties of helical edge states in a QSHI-superconductor (QSHI-SC) junction formed by a In⁢As(15 nm)/Ga⁢Sb(5 nm) double quantum well and a superconducting tantalum (Ta) constriction. We observe a robust conductance plateau up to 2 T, signaling resilient edge-state transport. Using a modified Landauer-Büttiker analysis, we find that the zero-field conductance is consistent with 98% Andreev-reflection probability owing to the high transparency of the (In⁢As/Ga⁢Sb)-Ta interface. Such resilience is consistent with the Dirac point for the edge states being buried in the bulk valence band. We further theoretically show that a buried Dirac point does not affect the robustness of the quasi-one-dimensional topological superconducting phase. We find that a buried Dirac point favors the hybridization of Majorana Kramers pairs (MKPs)—predicted to exist in a QSHI-SC constriction—and fermionic modes in the QSHI vacuum edge resulting in extended MKP states, highlighting the subtle role of buried Dirac points in probing MKPs.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Computational thermodynamics and its applications

Thermodynamics is a science concerning the state of a system, whether it is stable, metastable or un- stable, when interacting with the surroundings. In this overview, fundamentals of thermodynamics are briefly reviewed through the combination of first and second laws of thermodynamics for open and nonequilibrium systems to demonstrate that the reversible equilibrium and irreversible nonequilibrium thermodynamics can be integrated to enhance the power and utilities of thermodynamics. The recent progresses in computational thermodynamics, the remaining challenges, and potential impacts in broad scientific fields are discussed here. It is shown that computational thermodynamics enables the modeling of thermodynamics of a state as a function of both external and internal variables and quantitative calculations of a broad range of properties of a multicomponent system in terms of first and second derivatives of energy, including not only equilibrium states when there are no driving forces for any internal processes and but also non-equilibrium states with driving forces for internal processes. Consequently, external constraints such as fixed strain and internal degree of freedoms such as ordering and defects can be described in a coherent framework and applied to materials design. Furthermore, two important but largely overlooked aspects in thermodynamics will be discussed, i.e. the rigorous application of statistical thermodynamics with the probability of configurations and their contributions to system properties, and the applications of second derivatives of energy with respect to either two extensive variables or two potentials or a mixture of them in terms of understanding and predicting emergent behaviors, critical phenomena, kinetic coefficients, and mechanical properties.

36 MATERIALS SCIENCE↗

De novo design of transmembrane β barrels

Transmembrane β-barrel proteins (TMBs) are of great interest for single-molecule analytical technologies because they can spontaneously fold and insert into membranes and form stable pores, but the range of pore properties that can be achieved by repurposing natural TMBs is limited. We leverage the power of de novo computational design coupled with a “hypothesis, design, and test” approach to determine TMB design principles, notably, the importance of negative design to slow β-sheet assembly. We design new eight-stranded TMBs, with no homology to known TMBs, that insert and fold reversibly into synthetic lipid membranes and have nuclear magnetic resonance and x-ray crystal structures very similar to the computational models. These advances should enable the custom design of pores for a wide range of applications.

Vorobieva, Anastassia A.↗

On the role of topological complexity in spontaneous development of current sheets

The computations presented in this work aim to asses the importance of field line interlacing on spontaneous development of current sheets. From Parker's magnetostatic theorem, such development of current sheets is inevitable in a topologically complex magnetofluid, with infinite electrical conductivity, at equilibrium. Relevant initial value problems are constructed by superposition of two untwisted component fields, each component field being represented by a pair of global magnetic flux surface. The intensity of field line interlacing is then specified by the relative amplitude of the two superposed fields. The computations are performed by varying this relative amplitude. Also to have a direct visualization of current sheet formation, we follow the evolution of flux surfaces instead of the vector magnetic field. An important finding of this paper is in the demonstration that initial field lines having intense interlacing tend to develop current sheets which are distributed throughout the computational domain with no preference for topologically favorable sites like magnetic nulls or field reversal layers. The onsets of these current sheets are attributed to favorable contortions of magnetic flux surfaces where two oppositely directed parts of the same field line or different field lines come to close proximity. However, for less intensely interlaced field lines, the simulations indicate development of current sheets at sites only where the magnetic topology is favorable. These current sheets originate as two sets of anti-parallel complimentary field lines press onto each other.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Computational and Experimental Simulation of Static Memory Cells of Submicron Microcircuits under the Effect of Neutron Fluxes

The simulation of reversible single events in test samples of static memory microcircuits with design norms of 0.5, 0.35, 0.25, and 0.1 μm under the effect of neutron fluxes with various energies is performed. It is shown theoretically and experimentally that reversible single events can occur in modern microelectronics and nanoelectronics products under the effect of a fission-spectrum neutron flux caused by the passage of primary recoil atoms and nuclear reaction products along the microcircuit surface perpendicularly to the electric current lines in the near-drain transistor area. A series of irradiation experiments of static memory circuits with design norms of 0.35 μm is interpreted based on the proposed model.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Optimal Multistage Algorithm for Adjoint Computation

We reexamine the work of Stumm and Walther on multistage algorithms for adjoint computation. We provide an optimal algorithm for this problem when there are two levels of checkpoints, in memory and on disk. Previously, optimal algorithms for adjoint computations were known only for a single level of checkpoints with no writing and reading costs; a well-known example is the binomial checkpointing algorithm of Griewank and Walther. Stumm and Walther extended that binomial checkpointing algorithm to the case of two levels of checkpoints, but they did not provide any optimality results. We bridge the gap by designing the first optimal algorithm in this context. We experimentally compare our optimal algorithm with that of Stumm and Walther to assess the difference in performance.

97 MATHEMATICS AND COMPUTING↗

Schiff Moment of 129 Xe and 199 Hg [Slides]

Using Chicoma in 2025, we advanced DFT-based calculations needed to connect fundamental T violating interactions to measurable Schiff moments in heavy nuclei. We implemented DFT for polarized odd nuclear systems, a first step toward treating the time-reversal-breaking mean fields required for Schiff moment calculations. As a validation test, we computed polarizabilities using a perturbation with the dipole operator (see figure to the right) for the odd systems shown in Table 1; while no odd-system effect is expected for polarizabilities, this capability is essential for Schiff moments. Initial mean-field solutions for 129 Xe Schiff moment calculations were also generated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Numerical Analysis of Regular Material Point Method and its Application to Multiphase Flows

The material point method (MPM) is gaining wide popularity in engineering research to model and simulate complex multiphase flow dynamics. The method relies on solving the governing equations of motion and transport in a Lagrangian framework using particles also known as material points. The fluid and kinematic properties are stored on the material points while the spatial gradient calculation and temporal integration are performed on a background grid. This Lagrangian framework allows for large deformations, easy integration of constitutive models, and direct import of complex geometries as particles. However, despite their increasing popularity, very few studies have addressed the issues of numerical resolution and stability of MPM techniques. The presence of additional factors such as the number of material points-per-cell, the location of the material points, the CFL-like condition used in time update, and the grid shape functions also increase the complexity of the error analysis when compared to other finite element methods. In this presentation, we analyze the various forms of error incurred in the application of MPM to continuum mechanics and multiphase flows. The effect of the previously mentioned factors on the error dynamics is studied. The application of these principles to canonical and industrial problems is also presented.

high pressure reverse osmosis↗

Computationally Predicted High-Throughput Free-Energy Phase Diagrams for the Discovery of Solid-State Hydrogen Storage Reactions

The design of multinary solid-state material systems that undergo reversible phase changes via changes in temperature and pressure provides a potential means of safely storing hydrogen. However, fully mapping the stabilities of known or newly targeted compounds relative to competing phases at reaction conditions has previously required many stringent experiments or computationally demanding calculations of each compound’s change in Gibbs energy with respect to temperature, G(T). Here, we have extended the approach of constructing chemical potential phase diagrams based on ΔG f (T) to enable the analysis of phase stability at non-zero temperatures. We first performed density functional theory calculations to compute the formation enthalpies of binary, ternary, and quaternary compounds within several compositional spaces of current interest for solid-state hydrogen storage. Temperature effects on solid compound stability were then accounted for using our recently introduced machine learned descriptor for the temperature-dependent contribution G δ (T) to the Gibbs energy G(T). From these Gibbs energies, we evaluated each compound’s stability relative to competing compounds over a wide range of conditions and show using chemical potential and composition phase diagrams that the predicted stable phases and H2 release reactions are consistent with experimental observations. This demonstrates that our approach rapidly computes the thermochemistry of hydrogen release reactions for compounds at sufficiently high accuracy relative to experiment to provide a powerful framework for analyzing hydrogen storage materials. This framework based on G(T) enables the accelerated discovery of active materials for a variety of technologies that rely on solid-state reactions involving these materials.

08 HYDROGEN↗

Efficient Discovery of Novel Multicomponent Mixtures for Hydrogen Storage: A Combined Computational/Experimental Approach

The objective of the proposed program is to discover novel mixed hydrides for hydrogen storage, which enable the DOE 2010 system-level goals. Our goal is to find a material that desorbs 8.5 wt.% H 2 or more at temperatures below 85°C. The research program will combine first-principles calculations of reaction thermodynamics and kinetics with material and catalyst synthesis, testing, and characterization. We will combine materials from distinct categories (e.g., chemical and complex hydrides) to form novel multicomponent reactions. Systems to be studied include mixtures of complex hydrides and chemical hydrides [e.g. LiNH 2+ NH 3 BH 3 ] and nitrogen-hydrogen based borohydrides [e.g. Al(BH 4 ) 3 (NH 3 ) 3 ]. The 2010 and 2015 FreedomCAR/DOE targets for hydrogen storage systems are very challenging, and cannot be met with existing materials. The vast majority of the work to date has delineated materials into various classes, e.g., complex and metal hydrides, chemical hydrides, and sorbents. However, very recent studies indicate that mixtures of storage materials, particularly mixtures between various classes, hold promise to achieve technological attributes that materials within an individual class cannot reach. Our project involves a systematic, rational approach to designing novel multicomponent mixtures of materials with fast hydrogenation/dehydrogenation kinetics and favorable thermodynamics using a combination of state-of-the-art scientific computing and experimentation. We will use the accurate predictive power of first-principles modeling to understand the thermodynamic and microscopic kinetic processes involved in hydrogen release and uptake and to design new material/catalyst systems with improved properties. Detailed characterization and atomic-scale catalysis experiments will elucidate the effect of dopants and nanoscale catalysts in achieving fast kinetics and reversibility. And, state-of-the-art storage experiments will give key storage attributes of the investigated reactions, validate computational predictions, and help guide and improve computational methods. In sum, our approach involves a powerful blend of: 1) H2 Storage measurements and characterization, 2) State-of-the-art computational modeling, 3) Detailed catalysis experiments, 4) In-depth automotive perspective.

08 HYDROGEN↗

Improved Variable-Transformation Algorithm for the Hexagonal Variational Nodal Method

Hexagonal Variational Nodal Method (VNM) has been widely employed in fast reactor core calculation codes. It provides the weak form solution of the second order transport equation. It employs a direct three-dimensional expansion instead of the transverse integration technique, with the singular term caused by transverse integration over the hexagonal nodal avoided. However, as the increasing Pn expansion order in angle, the number of moments increases dramatically, which significantly increases the requirements in both memory and computational time. As a consequence, Yang proposed the Variable-Transformation Algorithm. It firstly used the reverse Cuthill-McKee (RCMK) reordering method to transform the response matrices into block-diagonal. Then, for the unreducible matrices, orthogonal transformation was applied, after which the matrices were sparse enough to be reduced by using the RCMK reordering once more. In this summary, based on a newly developed hexagonal VNM code NECP-VioletHex, the Variable Transformation Algorithm was investigated and improved by adopting more constraints to the basis functions to further take advantage of the symmetric properties of the hexagonal-z nodes. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A Reversed Impedance-Based Stability Criterion for IBR Grids: Preprint

The existing impedance-based stability criterion is effective for analyzing local control interactions; however, it is difficult to scale the existing criterion to analyze wide-area control interactions among numerous IBRs through a complex power system network. The scaled version of the existing criterion requires the impedance response of each IBR in the system as well as of the network looking from all the IBRs. It is quite challenging to obtain all these impedance responses because of the computational effort and the requirement of separately scanning the impedance of the network and the IBRs. We propose a reversed criterion for the impedance-based stability analysis to address these problems. In contrast to the existing criterion, the reversed criterion analyzes the stability of a power system when an IBR is disconnected from the system. The reversed criterion estimates the impact of an IBR on the frequency and damping of power system oscillation modes using the impedance scans of only the IBR and the grid at its terminal. It can be sequentially applied at different IBRs to evaluate their impact on the power system stability. In addition to scalability, the reversed criterion gives flexibility to focus only on a few selected IBRs, depending on their rating, the magnitude of oscillations observed at their terminals, and the vendor support available for implementing stabilizing control system updates. The reversed criterion is demonstrated on a 14-bus power system with 100% IBRs.

control interactions↗