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.
Origin of mean-field behavior in an elastic Ising model
Simple elastic models of spin-crossover compounds are known empirically to exhibit classical critical behavior. Here, we demonstrate how the long-range interactions responsible for this behavior arise naturally upon integrating out mechanical fluctuations of such a model. A mean-field theory applied to the resulting effective Hamiltonian quantitatively accounts for both thermodynamics and kinetics observed in computer simulations, including a barrier to magnetization reversal that grows extensively with system size. For nanocrystals, which break translational symmetry, a straightforward extension of mean-field theory yields similarly accurate results.
pnnl/arena
CFA ARENA is a novel programming model with the support of a runtime targeting asynchronous data-centric execution paradigm in a distributed system. All the machine nodes in ARENA are connected by a ring network to bring the specialized computation to the data rather than the reverse to minimize data movement. The programming interfaces are implemented using C++
Towards reverse mode automatic differentiation of Kokkos-based codes
Derivative computation is a key component of optimization, sensitivity analysis, uncertainty quantification, and the solving of nonlinear problems. Automatic differentiation (AD) is a powerful technique for evaluating such derivatives, and in recent years, has been integrated into programming environments such as Jax, PyTorch, and TensorFlow to support derivative computations needed for training of machine learning models, facilitating wide-spread use of these technologies. The C++ language has become the de facto standard for scientific computing due to numerous factors, yet language complexity has made the wide-spread adoption of AD technologies for C++ difficult, hampering the incorporation of powerful differentiable programming approaches into C++ scientific simulations. This is exacerbated by the increasing emergence of architectures, such as GPUs, with limited memory capabilities and requiring massive thread-level concurrency. C++ AD tools must effectively use these environments to bring novel scientific simulations to next-generation DOE experimental and observational facilities. In this project, we investigated source transformation-based automatic differentiation using LLVM compiler infrastructure to automatically generate portable and efficient gradient computations of Kokkos-based code. We have demonstrated that our proposed strategy is feasible by investigating the usage of a prototype LLVM-based source transformation tool to generate gradients of simple functions made of sequences of simple Kokkos parallel regions. Speedups of up to 500x compared to Sacado were observed on NVIDIA V100 GPU.
RanCompute: Computational Security in Embedded Devices via Random Input and Output Encodings
An embedded device in an insecure environment is subject to additional security risk through capture and reverse-engineering by a capable adversary. If this device contains a microchip performing sensitive computations, capture of the chip may leak functionality to an adversary. In this paper we propose a novel method in which we randomly encode the input operands and the outputs of a computation, thus not revealing the arithmetic operations being performed. The operations are sequenced in a graph representing the overall application. Once the initialization values are overwritten and lost, the results of these computations are indecipherable by the device performing the calculations as well as by any adversary. The result is transmitted back to a secure server which has stored the initialization values and so can decode the results which appear random to the adversary.
Numerical investigation of initial condition effects on Rayleigh-Taylor instability with acceleration reversals
The influence of initial conditions on miscible incompressible baroclinically driven Rayleigh-Taylor instability undergoing nonuniform acceleration is explored computationally using an implicit large eddy simulation (ILES) technique. We consider here the particular case of evolution during multiple reversals of acceleration direction, where the flow is alternately statically stable or unstable. In the unstable phase, the flow is driven by the baroclinic release of potential energy, whereas in the stable phase, work is done against the density stratification with the energy exchange taking place by wavelike mechanisms. These dynamics are fundamentally different; here, we track the evolution of volume-averaged turbulent statistics that are most sensitive to changes in the distribution of spectral power and bandwidth of the initial conditions as the flow alternates between dynamical regimes due to acceleration reversal.
Engineering and characterization of carbohydrate‐binding modules for imaging cellulose fibrils biosynthesis in plant protoplasts
Abstract Carbohydrate binding modules (CBMs) are noncatalytic domains that assist tethered catalytic domains in substrate targeting. CBMs have therefore been used to visualize distinct polysaccharides present in the cell wall of plant cells and tissues. However, most previous studies provide a qualitative analysis of CBM‐polysaccharide interactions, with limited characterization of engineered tandem CBM designs for recognizing polysaccharides like cellulose and limited application of CBM‐based probes to visualize cellulose fibrils synthesis in model plant protoplasts with regenerating cell walls. Here, we examine the dynamic interactions of engineered type‐A CBMs from families 3a and 64 with crystalline cellulose‐I and phosphoric acid swollen cellulose. We generated tandem CBM designs to determine various characteristic properties including binding reversibility toward cellulose‐I using equilibrium binding assays. To compute the adsorption ( nk on ) and desorption ( k off ) rate constants of single versus tandem CBM designs toward nanocrystalline cellulose, we employed dynamic kinetic binding assays using quartz crystal microbalance with dissipation. Our results indicate that tandem CBM3a exhibited the highest adsorption rate to cellulose and displayed reversible binding to both crystalline/amorphous cellulose, unlike other CBM designs, making tandem CBM3a better suited for live plant cell wall biosynthesis imaging applications. We used several engineered CBMs to visualize Arabidopsis thaliana protoplasts with regenerated cell walls using confocal laser scanning microscopy and wide‐field fluorescence microscopy. Lastly, we also demonstrated how CBMs as probe reagents can enable in situ visualization of cellulose fibrils during cell wall regeneration in Arabidopsis protoplasts.
Reverse-Mode Algorithmic Differentiation of an OpenMP-Parallel Compressible Flow Solver
y Reverse-mode algorithmic differentiation (AD) is an established method for obtaining adjoint derivatives of computer simulation applications. In computational fluid dynamics (CFD), adjoint derivatives of a cost function output such as drag or lift with respect to design parameters such as surface coordinates or geometry control points are a key ingredient for shape optimization, uncertainty quantification and flow control. The computational cost of CFD applications and their derivatives makes it essential to use high-performance computing hardware efficiently, including multi- and many-core architectures. Nevertheless, OpenMP is not supported in most AD tools, and previously shown methods achieve poor scalability of the derivative code. We present the AD of an OpenMP-parallelized finite volume compressible flow solver for unstructured meshes. Our approach enables us to reuse the parallelization of the original code in the computation of adjoint derivatives. The method works by identifying code segments that can be differentiated in reverse-mode without changing their memory access pattern. The OpenMP parallelization is integrated into the derivative code during the build process in a way that is robust to modifications of the original code and independent of the OpenMP support of the differentiation tool. We show the scalability of our adjoint CFD solver on test cases ranging from thousands to millions of finite volume mesh cells on CPUs with up to 16 threads as well as on an Intel XeonPhi card with 236 threads. We demonstrate that our approach is more practical to implement for production-sized CFD codes and produces more efficient adjoint derivative code than previously shown AD methods.
Issues in measure-preserving three dimensional flow integrators: Self-adjointness, reversibility, and non-uniform time stepping
Properties of integration schemes for solenoidal fields in three dimensions are studied, with a focus on integrating magnetic field lines in a plasma using adaptive time stepping. It is shown that implicit midpoint (IM) and a scheme we call three-dimensional leapfrog (LF) can do a good job (in the sense of preserving KAM tori) of integrating fields that are reversible, or (for LF) have a “special divergence-free” (SDF) property. We review the notion of a self-adjoint scheme, showing that such schemes are at least second order accurate and can always be formed by composing an arbitrary scheme with its adjoint. We also review the concept of reversibility, showing that a reversible but not exactly volume-preserving scheme can lead to a fractal invariant measure in a chaotic region, although this property may not often be observable. We also show numerical results indicating that the IM and LF schemes can fail to preserve KAM tori when the reversibility property (and the SDF property for LF) of the field is broken. We discuss extensions to measure preserving flows, the integration of magnetic field lines in a plasma and the integration of rays for several plasma waves. The main new result of this paper relates to non-uniform time stepping for volume-preserving flows. We investigate two potential schemes, both based on the general method of Feng and Shang [Numer. Math. 71, 451 (1995)], in which the flow is integrated in split time steps, each Hamiltonian in two dimensions. The first scheme is an extension of the method of extended phase space, a well-proven method of symplectic integration with non-uniform time steps. This method is found not to work, and an explanation is given. The second method investigated is a method based on transformation to canonical variables for the two split-step Hamiltonian systems. This method, which is related to the method of non-canonical generating functions of Richardson and Finn [Plasma Phys. Controlled Fusion 54, 014004 (2012)], appears to work very well.
Radiative Heat Transfer Calculation on 16384 GPUs Using a Reverse Monte Carlo Ray Tracing Approach with Adaptive Mesh Refinement
Modeling thermal radiation is computationally challenging in parallel due to its all-to-all physical and resulting computational connectivity, and is also the dominant mode of heat transfer in practical applications such as next-generation clean coal boilers, being modeled by the Uintah framework. However, a direct all-to-all treatment of radiation is prohibitively expensive on large computers systems whether homogeneous or heterogeneous. DOE Titan and the planned DOE Summit and Sierra machines are examples of current and emerging GPU-based heterogeneous systems where the increased processing capability of GPUs over CPUs exacerbates this problem. These systems require that computational frameworks like Uintah leverage an arbitrary number of on-node GPUs, while simultaneously utilizing thousands of GPUs within a single simulation. We show that radiative heat transfer problems can be made to scale within Uintah on heterogeneous systems through a combination of reverse Monte Carlo ray tracing (RMCRT) techniques combined with AMR, to reduce the amount of global communication. In particular, significant Uintah infrastructure changes, including a novel lock and contention-free, thread-scalable data structure for managing MPI communication requests and improved memory allocation strategies were necessary to achieve excellent strong scaling results to 16384 GPUs on Titan.
Utilizing Essential Symmetry Breaking in Auxiliary-Field Quantum Monte Carlo: Application to the Spin Gaps of the C 36 Fullerene and an Iron Porphyrin Model Complex
We present three distinct examples where phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) can be reliably performed with a single-determinant trial wave function with essential symmetry breaking. Essential symmetry breaking was first introduced by Lee and Head-Gordon [ Phys. Chem. Chem. Phys. 2019, 21, 4763-4778, 10.1039/C8CP07613H]. We utilized essential complex and time-reversal symmetry breaking with ph-AFQMC to compute the triplet–singlet energy gap in the TS12 set. We found statistically better performance of ph-AFQMC with complex-restricted orbitals than with spin-unrestricted orbitals. We then showed the utilization of essential spin symmetry breaking when computing the singlet–triplet gap of a known biradicaloid, C 36 . ph-AFQMC with spin-unrestricted Hartree–Fock (ph-AFQMC+UHF) fails catastrophically even with spin-projection and predicts no biradicaloid character. With approximate Brueckner orbitals obtained from regularized orbital-optimized second-order Møller–Plesset perturbation theory (κ-OOMP2), ph-AFQMC quantitatively captures strong biradicaloid character of C 36 . Lastly, we applied ph-AFQMC to the computation of the quintet–triplet gap in a model iron porphyrin complex where brute-force methods with a small active space fail to capture the triplet ground state. We show unambiguously that neither triplet nor quintet is strongly correlated using UHF, κ-OOMP2, and coupled-cluster with singles and doubles (CCSD) performed on UHF and κ-OOMP2 orbitals. There is no essential symmetry breaking in this problem. By virtue of this, we were able to perform UHF+ph-AFQMC reliably with a cc-pVTZ basis set and predicted a triplet ground state for this model geometry. The largest ph-AFQMC in this work correlated 186 electrons in 956 orbitals. Our work highlights the utility, scalability, and accuracy of ph-AFQMC with a single-determinant trial wave function with essential symmetry breaking for systems mainly dominated by dynamical correlation with little static correlation.
Statistics of time delay and scattering correlation functions in chaotic systems. II. Semiclassical approximation
We consider S-matrix correlation functions for a chaotic cavity having M open channels, in the absence of time-reversal invariance. Relying on a semiclassical approximation, we compute the average over E of the quantities Tr[S{sup †}(E − ϵ) S(E + ϵ)]{sup n}, for general positive integer n. Our result is an infinite series in ϵ, whose coefficients are rational functions of M. From this, we extract moments of the time delay matrix Q = − iħS{sup †}dS/dE and check that the first 8 of them agree with the random matrix theory prediction from our previous paper [M. Novaes, J. Math. Phys. 56, 062110 (2015)].
PtychoNet: Fast and High Quality Phase Retrieval for Ptychography
Ptychography is a coherent diffractive imaging method that captures multiple diffraction patterns of a sample with a set of shifted localized illuminations (“probes”). The reconstruction problem, known as “phase retrieval”, is typically solved by iterative algorithms. In this paper, we propose PtychoNet, a deep learning based method to perform phase retrieval for ptychography in a non-iterative manner. We devise a generative network to encode a full ptychography scan, reverse the diffractions at each scanning point and compute the amplitude and phase of the object. We demonstrate successful reconstructions using PtychoNet as well as recovering fine features in the case of extreme sparse scanning where conventional iterative methods fail to give recognizable features.
De novo design of modular protein hydrogels with programmable intra- and extracellular viscoelasticity
Relating the macroscopic properties of protein-based materials to their underlying component microstructure is an outstanding challenge. Here, we exploit computational design to specify the size, flexibility, and valency of de novo protein building blocks, as well as the interaction dynamics between them, to investigate how molecular parameters govern the macroscopic viscoelasticity of the resultant protein hydrogels. We construct gel systems from pairs of symmetric protein homo-oligomers, each comprising 2, 5, 24, or 120 individual protein components, that are crosslinked either physically or covalently into idealized step-growth biopolymer networks. Through rheological assessment, we find that the covalent linkage of multifunctional precursors yields hydrogels whose viscoelasticity depends on the crosslink length between the constituent building blocks. In contrast, reversibly crosslinking the homo-oligomeric components with a computationally designed heterodimer results in viscoelastic biomaterials exhibiting fluid-like properties under rest and low shear, but solid-like behavior at higher frequencies. Exploiting the unique genetic encodability of these materials, we demonstrate the assembly of protein networks within living mammalian cells and show via fluorescence recovery after photobleaching (FRAP) that mechanical properties can be tuned intracellularly in a manner similar to formulations formed extracellularly. We anticipate that the ability to modularly construct and systematically program the viscoelastic properties of designer protein-based materials could have broad utility in biomedicine, with applications in tissue engineering, therapeutic delivery, and synthetic biology.
A new method to compute periodic orbits in general symplectic maps
The search of high-order periodic orbits has been typically restricted to problems with symmetries that help to reduce the dimension of the search space. Well-known examples include reversible maps with symmetry lines. The present work proposes a new method to compute high-order periodic orbits in twist maps without the use of symmetries. The method is a combination of the parameterization method in Fourier space and a Newton–Gauss multiple shooting scheme. The parameterization method has been successfully used in the past to compute quasi-periodic invariant circles. However, this is the first time that this method is used in the context of periodic orbits. Numerical examples are presented showing the accuracy and efficiency of the proposed method. Furthermore, the method is also applied to verify the renormalization prediction of the residues’ convergence at criticality (extensively studied in reversible maps) in the relatively unexplored case of maps without symmetries.
Interpreting the Operando X-ray Absorption Near-Edge Structure of Supported Cu and CuPd Clusters in Conditions of Oxidative Dehydrogenation of Propane: Dynamic Changes in Composition and Size
Supported subnano-cluster catalysts are highly dynamic, developing true active sites only under the pressures and temperatures of reaction conditions. Operando X-ray absorption near-edge structure (XANES) spectroscopy can track changes in the oxidation state and the local environment of cluster atoms, providing insight into the development of these active sites. While bulk metal, oxide, and hydroxide standards are often used for fitting experimental XANES spectra to obtain average oxidation states, we recently showed that computed cluster standards of relevant compositions are a more suitable basis, producing more accurate fits. Here, we theoretically interpret the operando XANES of supported Cu 3 Pd and Cu 4 clusters during temperature-programmed reaction (TPRx) of oxidative dehydrogenation of propane. We use an expanded basis set including both globally optimized computed clusters and bulk standards. Not only can we track reversible composition/oxidation state change with temperature, but also the irreversible growth of the bulk fraction upon heating, which we attribute to cluster sintering. This has important implications for the mechanism of the catalyzed reaction and the nature of the available active sites. Here, we propose that operando XANES provides most significant insight into the nature of supported cluster catalysts in reaction conditions when interpreted using mixed computed cluster and bulk standards.
Reversible Control of Interfacial Magnetism through Ionic-Liquid-Assisted Polarization Switching
The ability to control magnetism of materials via electric field enables a myriad of technological innovations in information storage, sensing, and computing. In this paper, we use ionic-liquid-assisted ferroelectric switching to demonstrate reversible modulation of interfacial magnetism in a multiferroic heterostructure composed of ferromagnetic (FM) La 0.8 Sr 0.2 MnO 3 and ferroelectric (FE) PbZr 0.2 Ti 0.8 O 3 . It is shown that ionic liquids can be used to persistently and reversibly switch a large area of a FE film. Finally, this is a prerequisite for polarized neutron reflectometry (PNR) studies that are conducted to directly probe magnetoelectric coupling of the FE polarization to the interfacial magnetization.
Up‐And‐Coming Advances in Optical and Microwave Nonreciprocity: From Classical to Quantum Realm
Reciprocity is a fundamental physical principle that roots in the time‐reversal symmetry of physical laws. It allows making predictions on any arbitrary complex system's response and operation and hence simplifies the analysis. However, there are many practical situations in which it is advantageous to break reciprocity, e.g., isolators preventing wave scattering back to lasers and generators, full‐duplex systems for multiplexing transmission and receiving in the same channel, nonreciprocal cavity excitation, and protection of fragile states of superconductor quantum computers from thermal noise. The most widespread approach to time‐reversal symmetry breaking and nonreciprocity based on magnetic field biasing suffers from bulkiness, cost ineffectiveness, and loss, motivating researchers and engineers to search for more practical approaches. Herein, the up‐and‐coming advances in optical nonreciprocity, including new materials (Weyl semimetals, topological insulators, metasurfaces), active structures, time‐modulation, parity‐time (PT)‐symmetry breaking, nonlinearity combined with a structural asymmetry, quantum nonlinearity, unidirectional gain and loss, chiral quantum states and valley polarization are overviewed. A general description of nonreciprocal systems is provided and the pros and cons of the mentioned approaches to nonreciprocity are discussed.