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 217 records · Page 12

A conformational equilibrium in the nitrogenase MoFe protein with an α-V70I amino acid substitution illuminates the mechanism of H 2 formation

Study of α-V70I-substituted nitrogenase MoFe protein identified Fe6 of FeMo-cofactor (Fe 7 S 9 MoC-homocitrate) as a critical N 2 binding/reduction site. Freeze-trapping this enzyme during Ar turnover captured the key catalytic intermediate in high occupancy, denoted E 4 (4H), which has accumulated 4[e — /H + ] as two bridging hydrides, Fe2–H–Fe6 and Fe3–H–Fe7, and protons bound to two sulfurs. E 4 (4H) is poised to bind/reduce N 2 as driven by mechanistically-coupled H 2 reductive-elimination of the hydrides. This process must compete with ongoing hydride protonation (HP), which releases H 2 as the enzyme relaxes to state E 2 (2H), containing 2[e — /H + ] as a hydride and sulfur-bound proton; accumulation of E 4 (4H) in α-V70I is enhanced by HP suppression. EPR and 95 Mo ENDOR spectroscopies now show that resting-state α-V70I enzyme exists in two conformational states, both in solution and as crystallized, one with wild type (WT)-like FeMo-co and one with perturbed FeMo-co. These reflect two conformations of the Ile residue, as visualized in a reanalysis of the X-ray diffraction data of α-V70I and confirmed by computations. EPR measurements show delivery of 2[e — /H + ] to the E 0 state of the WT MoFe protein and to both α-V70I conformations generating E 2 (2H) that contains the Fe3–H–Fe7 bridging hydride; accumulation of another 2[e — /H + ] generates E 4 (4H) with Fe2–H–Fe6 as the second hydride. E 4 (4H) in WT enzyme and a minority α-V70I E 4 (4H) conformation as visualized by QM/MM computations relax to resting-state through two HP steps that reverse the formation process: HP of Fe2–H–Fe6 followed by slower HP of Fe3–H–Fe7, which leads to transient accumulation of E 2 (2H) containing Fe3–H–Fe7. In the dominant α-V70I E 4 (4H) conformation, HP of Fe2–H–Fe6 is passively suppressed by the positioning of the Ile sidechain; slow HP of Fe3–H–Fe7 occurs first and the resulting E 2 (2H) contains Fe2–H–Fe6. Importantly, it is this HP suppression in E 4 (4H) that enables α-V70I MoFe to accumulate E 4 (4H) in high occupancy. In addition, HP suppression in α-V70I E 4 (4H) kinetically unmasks hydride reductive-elimination without N 2 -binding, a process that is precluded in WT enzyme.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Phase diagram of the TIP4P/Ice water model by enhanced sampling simulations

In this work, we studied the phase diagram for the TIP4P/Ice water model using enhanced sampling molecular dynamics simulations. Our approach is based on the calculation of ice–liquid free energy differences from biased coexistence simulations that reversibly sample the melting and growth of layers of ice. We computed a total of 19 melting points for five different ice polymorphs, which are in excellent agreement with the melting lines obtained from the integration of the Clausius–Clapeyron equation. For proton-ordered and fully proton-disordered ice phases, the results are in very good agreement with previous calculations based on thermodynamic integration. For the partially proton-disordered ice III, we find a large increase in stability that is in line with previous observations using direct coexistence simulations for the TIP4P/2005 model. This issue highlights the robustness of the approach employed here for ice polymorphs with diverse degrees of proton disorder. Our approach is general and can be applied to the calculation of other complex phase diagrams.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DYNAMO ACTIVITIES DRIVEN BY MAGNETOROTATIONAL INSTABILITY AND THE PARKER INSTABILITY IN GALACTIC GASEOUS DISKS

We carried out global three-dimensional magnetohydrodynamic simulations of dynamo activities in galactic gaseous disks without assuming equatorial symmetry. Numerical results indicate the growth of azimuthal magnetic fields non-symmetric to the equatorial plane. As the magnetorotational instability (MRI) grows, the mean strength of magnetic fields is amplified until the magnetic pressure becomes as large as 10% of the gas pressure. When the local plasma {beta} (=p {sub gas}/p {sub mag}) becomes less than 5 near the disk surface, magnetic flux escapes from the disk by the Parker instability within one rotation period of the disk. The buoyant escape of coherent magnetic fields drives dynamo activities by generating disk magnetic fields with opposite polarity to satisfy the magnetic flux conservation. The flotation of the azimuthal magnetic flux from the disk and the subsequent amplification of disk magnetic field by the MRI drive quasi-periodic reversal of azimuthal magnetic fields on a timescale of 10 rotation periods. Since the rotation speed decreases with radius, the interval between the reversal of azimuthal magnetic fields increases with radius. The rotation measure computed from the numerical results shows symmetry corresponding to a dipole field.

79 ASTRONOMY AND ASTROPHYSICS↗

Rupture Dynamics Simulations of Shallow Crustal Earthquakes on Reverse Slip Faults

The main goal of this study was to perform dynamic rupture simulations in order constrain the characterization of kinematic rupture models implemented in techniques for simulating strong ground motion for crustal earthquakes. First, we investigated the rupture process of the 2008 Iwate-Miyagi Nairiku earthquake by developing a dynamic rupture model on a reverse fault using a trial-and-error technique that produced a slip distribution and near-fault ground motion that matched the recorded ones. The simulations were performed in the frequency range 0-2 Hz, using a 3D staggered grid finite-difference method and a linear slip weakening friction law. Constrained by the observed slip distribution and consistent with dynamic rupture models, the derived kinematic rupture model of the Iwate-Miyagi Nairiku earthquake contains areas with relatively large slip rate, representing strong motion generation areas, set against lower amplitude heterogeneous background slip, and a relatively low slip rate in the weak zone of the top-most crust (upper 3km). Second, we performed rupture dynamics modeling to constrain shallow slip characterization in rupture models for strike-slip crustal earthquakes. The objective was to establish general rules about the characterization of slip rate function and slip at shallow depths and in the strong motion generation areas (SMGAs). The simulations of spontaneous rupture were performed in the frequency range 0-2.0Hz, using a 3D staggered-grid finite-difference method and a layer over half-space 1D crustal velocity model with a minimum shear-wave velocity of 2.8 km/s. In order to account for changes in material ductility and reduction of stress drop, observed in the shallow crust (upper 3-5 km), and the transition from ductile state to brittle state in the upper seismogenic zone, in our stress models we included a shallow weak zone (<4km). In this zone the stress drop was set to zero at the free surface and gradually increased with depth, while the slip weakening distance was set to 75 cm at the free surface and decreased to 50 cm at the base of the weak zone. As in the case of reverse faulting during the Iwate-Miyagi Nairiku earthquake, from these computations we found a systematic change in the shape of the slip-rate function, from Kostrov-type in the deeper part of the fault to a more symmetric cosine-type in the upper few kilometers, near the free-surface. Moreover, the average slip duration in the weak zone, with respect to slip duration in the deeper parts of the fault, increases by at most a factor of 1.5. We found a systematic and gradual change in the shape of the slip-rate function from Kostrov-type in the asperity areas (SMGAs), to more symmetric cosine-type in the upper few km near the free-surface, on the long period motion generation area (LMGA). The effective rise time in the LMGA, with respect to that in the SMGAs, increases by at most a factor of 2. Effective rise time is the time difference between the time at which the slip rate drops to a level that is equal to 25% of its peak and the onset time of the signal. In addition, the slip in the LMGA located above the SMGAs is almost the same as the one in the SMGAs, and about 1.5 times larger than the average slip.

58 GEOSCIENCES↗

Computational and experimental search for potential polyanionic K-ion cathode materials

Discovering high-energy cathode materials is critical to construct K-ion batteries for practical applications. Owing to the great success of layered oxides in Li- and Na-ion systems, K layered cathodes have also been investigated in recent years. However, the much larger size of K + compared to Li or Na introduces strong K + –K + interaction within the layer, which results in a sloped voltage profile thereby limiting the specific capacity and operating voltage. In contrast, polyanionic materials with a three-dimensional K + arrangement can effectively mitigate K + –K + interaction. In this work, ten K polyanionic compounds with theoretical capacity >100 mA h g –1 are screened from the Inorganic Crystal Structure Database as potential cathode materials for K-ion batteries. Among the ten proposed compounds, K 2 MnP 2 O 7 , K 2 Mn 2 P 2 O 7 F 2 , K 2 Fe 2 P 2 O 7 F 2 , and K 6 V 2 (PO 4 ) 4 with average voltage <4.5 V are synthesized and evaluated electrochemically. While the re-insertion of K into these compounds is not fully reversible, it may be related to the very high migration barrier that we compute for K ions. In addition, we show the successful synthesis of a series of K 3 V 3–x Cr x (PO 4 ) 4 (x = 0, 1, 2, 3) compounds. Among these, K 3 V 2 Cr(PO 4 ) 4 exhibits the largest reversible capacity, as revealed by the in situ investigation. Lastly, we find that the redox couples in many of these compounds sit at remarkably high potential, even higher than in equivalent Li compounds, which brings both opportunities and challenges in the future research of K polyanion cathodes.

25 ENERGY STORAGE↗

Design, synthesis, and biological testing of biphenylmethyloxazole inhibitors targeting HIV-1 reverse transcriptase

In this article, we report non-nucleoside inhibitors of HIV-1 reverse transcriptase (NNRTIs) using a biphenylmethyloxazole pharmacophore. A crystal structure of benzyloxazole 1 was obtained and suggested the potential viability of biphenyl analogues. In particular, 6a, 6b, and 7 turned out to be potent NNRTIs with low-nanomolar activity in enzyme inhibition and infected T-cell assays, and with low cytotoxicity. Though modeling further suggested that analogues with fluorosulfate and epoxide warheads might provide covalent modification of Tyr188, synthesis and testing did not find evidence for this outcome.

60 APPLIED LIFE SCIENCES↗

Reversible Olefin Addition to Extended Lattices of a Nickel–Selenium Framework

We performed density functional theory computations to study the structural and electronic properties as the basis of ethylene addition activity for [Ni(XC)4]n (X = Se, S)-extended lattices. We demonstrated that the mechanism of ethylene cycloaddition to a periodic [Ni(SeC)4]n two-dimensional (2D) network is analogous to that previously described for [Ni(SC)4]n 2D sheets and similar to the metal bis(dithiolene) molecular complexes [M(S2C2R2)2] (M = Ni, Pd, Pt, Co, Cu). These nanosheet materials avoid decomposition upon olefin addition, which is one of the main limitations of the molecular metal bis(dithiolene) complexes, as we find the decomposition processes to be thermodynamically unfavorable. Our calculations also suggest that the preferred conformation of the [Ni(SeC)4]n bilayer lattice is parallel displaced, with the Se atoms positioned above the Ni atoms, which is different from the eclipsed conformation found for [Ni(SC)4]n. We also managed to optimize an adduct of [Ni(SC)4]n in the bilayer form, which exceed the ethylene coverage of molecular complexes. We calculate that the preferred three-dimensional geometry of the stacked sheets is eclipsed because of strong van der Waals interactions. Such an arrangement of the sheets indicates that these materials should be highly porous, pointing to the high capacity for olefin bindings. Indeed, a few moderately stable ethylene adducts have been located. Owing to their unique structures and chemical reactivity, these newly predicted materials can be potentially developed as electrocatalysts for olefin purification. © 2018 American Chemical Society.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Thermodynamics and kinetics of precipitation and austenite reversion during aging of Ti-free grade 300 maraging steel manufactured by laser powder bed fusion (LPBF)

Austenite reversion and intermetallic precipitation are two key phase transformations in maraging steels which determine the strength and ductility. Usually, an optimum aging temperature/time is determined using trial and error approach with inputs from thermodynamic databases and a-posteriori microscopy/hardness evaluations. In this study, a thermo-kinetic model is developed to predict the precipitation and austenite reversion kinetics during various post-fabrication heat treatments in a Ti-free variant of grade 300 maraging steel manufactured using laser powder bed fusion. The thermo-kinetic model was based on the principles of Scheil solidification and the classical nucleation theory. The predictions from the thermo-kinetic model were verified using suitable microscopy techniques, and the results from the thermo-kinetic model agree well with the experimental measurements. Overall, results highlight the fact that an Integrated Computational Materials Engineering approach can be used to predict the precipitation and austenite reversion kinetics, and be used as an input for an efficient design of heat treatments for achieving optimum mechanical properties.

36 MATERIALS SCIENCE↗

Sequence optimizations in a high-performance computing environment

Embodiments are directed to techniques to determine dataflow graph instructions comprising one or more pick/switch instruction pairs and generate a reverse static single assignment graph based on the dataflow graph instructions, the reverse static single assignment graph comprising strongly connected components, each of the strongly connected components associated with at least one of the one or more pick/switch instruction pairs. Embodiments also include traversing the reverse static single assignment graph depth-first, and replace pick/switch instructions associated with strongly connected components having configuration values with compound instructions.

97 MATHEMATICS AND COMPUTING↗

Statistical effects related to low numbers of reacting molecules analyzed for a reversible association reaction A + B = C in ideally dispersed systems: An apparent violation of the law of mass action

Theoretical analysis and computer simulations (Monte Carlo and numerical integration of differential equations) show that the statistical effect of a small number of reacting molecules depends on a way the molecules are distributed among the small volume nano-reactors (droplets in this study). A simple reversible association A + B = C was chosen as a model reaction, enabling to observe both thermodynamic (apparent equilibrium constant) and kinetic effects of a small number of reactant molecules. When substrates are distributed uniformly among droplets, all containing the same equal number of substrate molecules, the apparent equilibrium constant of the association is higher than the chemical one (observed in a macroscopic—large volume system). The average rate of the association, being initially independent of the numbers of molecules, becomes (at higher conversions) higher than that in a macroscopic system: the lower the number of substrate molecules in a droplet, the higher is the rate. This results in the correspondingly higher apparent equilibrium constant. A quite opposite behavior is observed when reactant molecules are distributed randomly among droplets: the apparent association rate and equilibrium constants are lower than those observed in large volume systems, being the lower, the lower is the average number of reacting molecules in a droplet. The random distribution of reactant molecules corresponds to ideal (equal sizes of droplets) dispersing of a reaction mixture. Our simulations have shown that when the equilibrated large volume system is dispersed, the resulting droplet system is already at equilibrium and no changes of proportions of droplets differing in reactant compositions can be observed upon prolongation of the reaction time.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Broken time-reversal symmetry in the topological superconductor UPt 3

Topological properties of materials are of fundamental as well as practical importance. Of particular interest are unconventional superconductors that break time-reversal symmetry, for which the superconducting state is protected topologically and vortices can host Majorana fermions with potential use in quantum computing. However, in striking contrast to the unconventional A phase of superfluid 3 He where chiral symmetry was directly observed, identification of broken time-reversal symmetry of the superconducting order parameter, a key component of chiral symmetry, has presented a challenge in bulk materials. The two leading candidates for bulk chiral superconductors are UPt 3 and Sr 2 RuO 4 , although evidence for broken time-reversal symmetry comes largely from surface-sensitive measurements. A long-sought demonstration of broken time-reversal symmetry in bulk Sr 2 RuO 4 is the observation of edge currents, which has so far not been successful. The situation for UPt 3 is not much better. In this work, we use vortices to probe the superconducting state in ultraclean crystals of UPt 3 . Using small-angle neutron scattering, a strictly bulk probe, we demonstrate that the vortices possess an internal degree of freedom in one of its three superconducting phases, providing direct evidence for bulk broken time-reversal symmetry in this material.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Exact Optimal Power Dispatch in Unbalanced Distribution Systems With High PV Penetration

Smart inverters provide additional control capability to help optimize the operation of distribution systems. This paper proposes a framework for exact optimal active and reactive power dispatch of distributed photovoltaic (PV) generation, switched capacitors, and tap changers in large multi-phase unbalanced distribution systems. The objectives of the optimal dispatch are minimization of the energy loss, PV active power curtailment, and capacitor and tap changer switching operations, in addition to elimination of voltage violations and reverse power flow. The optimization problem is formulated in rectangular coordinates as a nonlinear, nonconvex problem. Effective computational strategies are proposed to allow the application of predictorcorrector primal-dual interior point method to solve optimization problems in real-time with a large number of constraints and variables, including discrete variables corresponding to switched capacitors and tap changers. The accuracy of the numerical solution and the ability to implement the proposed framework are validated using the unbalanced multi-phase IEEE 34-bus and EPRI 2,998-bus distribution systems with 15-minute load and PV data. The results show a significant loss reduction and complete removal of both voltage violations and reverse power flow.

Nguyen, Quan (ORCID:0000000347769557)↗

Local Explosion Detection and Infrasound Localization by Reverse Time Migration Using 3-D Finite-Difference Wave Propagation

Infrasound data are routinely used to detect and locate volcanic and other explosions, using both arrays and single sensor networks. However, at local distances (< km) topography often complicates acoustic propagation, resulting in inaccurate acoustic travel times leading to biased source locations when assuming straight-line propagation. Here we present a new method, termed Reverse Time Migration-Finite-Difference Time Domain (RTM-FDTD), that integrates numerical modeling into the standard RTM back-projection process. Travel time information is computed across the entire potential source grid via FDTD modeling to incorporate the effects of topography. The waveforms are then back-projected and stacked at each grid point, with the stack maximum corresponding to the likely source. We apply our method to three volcanoes with different network configurations, source-receiver distances, and topography. At Yasur Volcano, Vanuatu, RTM-FDTD locates explosions within ~20 m of the source and differentiates between multiple vents. RTM-FDTD produces a more accurate location for the two Yasur subcraters than standard RTM and doubles the number of detected events. At Sakurajima Volcano, Japan, RTM-FDTD locates the source within 50 m of the active vent despite notable topographic blocking. The RTM-FDTD location is similar to that from the Time Reversal Mirror method, but is more computationally efficient. Lastly, at Shishaldin Volcano, Alaska, RTM and RTM-FDTD both produce realistic source locations (<50 m) for ground-coupled airwaves recorded on a four-station seismic network. We show that RTM is an effective method to detect and locate infrasonic sources across a variety of scenarios, and by integrating numerical modeling, RTM-FDTD produces more accurate source locations and increases the detection capability.

58 GEOSCIENCES↗

Recent Advances for Improving the Accuracy, Transferability, and Efficiency of Reactive Force Fields

Reactive force fields provide an affordable model for simulating chemical reactions at a fraction of the cost of quantum mechanical approaches. However, classically accounting for chemical reactivity often comes at the expense of accuracy and transferability, while computational cost is still large relative to nonreactive force fields. Here, we summarize recent efforts for improving the performance of reactive force fields in these three areas with a focus on the ReaxFF theoretical model. To improve accuracy, we describe recent reformulations of charge equilibration schemes to overcome unphysical long-range charge transfer, new ReaxFF models that account for explicit electrons, and corrections for energy conservation issues of the ReaxFF model. To enhance transferability we also highlight new advances to include explicit treatment of electrons in the ReaxFF and hybrid nonreactive/reactive simulations that make it possible to model charge transfer, redox chemistry, and large systems such as reverse micelles within the framework of a reactive force field. To address the computational cost, we review recent work in extended Lagrangian schemes and matrix preconditioners for accelerating the charge equilibration method component of ReaxFF and improvements in its software performance in LAMMPS.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

String Method for Protein–Protein Binding Free-Energy Calculations

A powerful computational strategy to determine the equilibrium association constant of two macromolecules with explicit-solvent molecular dynamics (MD) simulations is the "geometric route", which considers the reversible physical separation of the bound complex in solution. Nonetheless, multiple challenges remain to render this type of methodology reliable and computationally efficient in practice. In particular, in one, formulation of the geometric route relies on the potential of mean force (PMF) for physically separating the two binding partners restrained along a straight axis, which must be selected prior to the calculation. However, practical applications indicate that the calculation of the separation PMF along the predefined rectilinear pathway may be suboptimal and slowly convergent. Recognizing that a rectilinear straight separation pathway is generally not representative of how the protein complex physically separates in solution, we put forth a novel theoretical framework for binding free-energy calculations, leaning on the optimal curvilinear minimum free-energy path (MFEP) determined from the string method. The proposed formalism is validated by comparing the results obtained using both rectilinear and curvilinear pathways for a prototypical host-guest complex formed by cucurbit[7]uril (CB[7]) binding benzene, and for the barnase-barstar protein complex. On the basis of multi-microsecond MD calculations, we find that the calculations following the traditional rectilinear pathway and the string-based curvilinear pathway agree quantitatively, but convergence is faster with the latter.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Efficient Bayesian inference with latent Hamiltonian neural networks in No-U-Turn Sampling

When sampling for Bayesian inference, one popular approach in the computational field is to use Hamiltonian Monte Carlo (HMC) and specifically the No-U-Turn Sampler (NUTS), which automatically decides the end time of the Hamiltonian trajectory. However, HMC and NUTS can require numerous numerical gradients of the target density and can prove slow in practice when relying on computationally expensive forward models. We propose Latent Hamiltonian neural networks (L-HNNs) with HMC and NUTS for solving Bayesian inference problems. Once trained, L-HNNs do not require numerical gradients of the target density during sampling, and hence numerous evaluations of the forward computational model. Moreover, L-HNNs satisfy important properties such as perfect time reversibility and Hamiltonian conservation, making them well-suited for use within HMC and NUTS because stationarity can be shown. We also propose the integration of L-HNNs in an online error monitoring scheme, in which numerical gradients of the target density are used for a few samples whenever the L-HNNs prediction errors are large. This online error monitor scheme prevents sample degeneracy in regions of low probability density and ensures robust uncertainty quantification. We demonstrate L-HNNs in NUTS with online error monitoring on several analytical examples involving complex, heavy-tailed, and high-local-curvature probability densities. We then demonstrate the applicability of L-HNNs in NUTS to two computational case studies, namely the Allen-Cahn stochastic partial differential equation and an elliptic partial differential equation with 25 and 50 inference parameters, respectively. Overall, the L-HNNs in NUTS with online error monitoring satisfactorily inferred these probability densities. In conclusion, compared to traditional NUTS, L-HNNs in NUTS with online error monitoring required 1–2 orders of magnitude fewer numerical gradients of the target density and improved the effective sample size (ESS) per gradient (which is a measure of both the sampling quality and the computational expense) by an order of magnitude.

97 MATHEMATICS AND COMPUTING↗

Ultrahigh pressure compaction-resistant thin film crosslinked composite reverse osmosis membranes

In this study, we present a class of thin-film crosslinked (TFX) composite reverse osmosis (RO) membranes that resist physical compaction at ultrahigh pressures (up to 200 bar). Since RO membranes experience compaction at virtually all pressure ranges, the ability to resist compaction has widespread implications for RO membrane technology. The process described herein involves crosslinking a phase inverted porous polyimide (PI) support membrane followed by interfacial polymerization of a polyamide layer, thereby forming a fully thermoset composite membrane structure. We explore a range of phase inversion membrane formation parameters such as PI concentration, solvent-cosolvent ratios, coagulation bath composition, and crosslinking methods in addition to interfacial polymerization reaction chemistry and conditions. Overall, TFX membranes exhibit significantly less compaction compared to hand-cast and commercial high-pressure RO membranes, experiencing less than 10% decline in water permeance and maintaining salt rejection over 99% for NaCl solutions up to 180,000 mg/L with 200 bar applied pressure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗