Search NASA⌕ Search

SEARCH · Search NASA

Results for “Contracts”

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 793 records · Page 44

TAMM: Tensor algebra for many-body methods

Tensor algebra operations such as contractions in computational chemistry consume a significant fraction of the computing time on large-scale computing platforms. The widespread use of tensor contractions between large multi-dimensional tensors in describing electronic structure theory has motivated the development of multiple tensor algebra frameworks targeting heterogeneous computing platforms. In this paper, we present Tensor Algebra for Many-body Methods (TAMM), a framework for productive and performance-portable development of scalable computational chemistry methods. TAMM decouples the specification of the computation from the execution of these operations on available high-performance computing systems. With this design choice, the scientific application developers (domain scientists) can focus on the algorithmic requirements using the tensor algebra interface provided by TAMM, whereas high-performance computing developers can direct their attention to various optimizations on the underlying constructs, such as efficient data distribution, optimized scheduling algorithms, and efficient use of intra-node resources (e.g., graphics processing units). The modular structure of TAMM allows it to support different hardware architectures and incorporate new algorithmic advances. We describe the TAMM framework and our approach to the sustainable development of scalable ground- and excited-state electronic structure methods. We present case studies highlighting the ease of use, including the performance and productivity gains compared to other frameworks.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

High-performance strategies for the recent MRSF-TDDFT in GAMESS

Multiple ERI (Electron Repulsion Integral) tensor contractions (METC) with several matrices are ubiquitous in quantum chemistry. In response theories, the contraction operation, rather than ERI computations, can be the major bottleneck, as its computational demands are proportional to the multiplicatively combined contributions of the number of excited states and the kernel pre-factors. Here, this paper presents several high-performance strategies for METC. Optimal approaches involve either the data layout reformations of interim density and Fock matrices, the introduction of intermediate ERI quartet buffer, and loop-reordering optimization for a higher cache hit rate. The combined strategies remarkably improve the performance of the MRSF (mixed reference spin flip)-TDDFT (time-dependent density functional theory) by nearly 300%. The results of this study are not limited to the MRSF-TDDFT method and can be applied to other METC scenarios.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Models of resistive wall tearing mode disruptions

Disruptions are a serious issue in tokamaks. In a disruption, the thermal energy is lost by means of an instability which could be a resistive wall tearing mode (RWTM). During precursors to a disruption, the plasma edge region cools, causing the current to contract. Model sequences of contracted current equilibria are given, and their stability is calculated. A linear stability study shows that there is a maximum value of edge q a ≈ 3 for RWTMs to occur. Further, this also implies a minimum rational surface radius normalized to plasma radius from RWTMs to be unstable. Nonlinear simulations are performed using a similar model sequence derived from an equilibrium reconstruction. There is a striking difference in the results, depending on whether the wall is ideal or resistive. With an ideal wall, the perturbations saturate at moderate amplitude, causing a minor disruption without a thermal quench. With a resistive wall, there is a major disruption with a thermal quench, if the edge q a ≤ 3. There is a sharp transition in nonlinear behavior at q a = 3. This is consistent with the linear model and with experiments. If disruptions are caused by RWTMs, then devices with highly conducting walls, such as the International Tokamak Experimental Reactor will experience much milder, tolerable, disruptions than presently predicted.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The myosin interacting-heads motif present in live tarantula muscle explains tetanic and posttetanic phosphorylation mechanisms

Striated muscle contraction involves sliding of actin thin filaments along myosin thick filaments, controlled by calcium through thin filament activation. In relaxed muscle, the two heads of myosin interact with each other on the filament surface to form the interacting-heads motif (IHM). A key question is how both heads are released from the surface to approach actin and produce force. Here we used time-resolved synchrotron X-ray diffraction to study tarantula muscle before and after tetani. The patterns showed that the IHM is present in live relaxed muscle. Tetanic contraction produced only a very small backbone elongation, implying that mechanosensing—proposed in vertebrate muscle—is not of primary importance in tarantula. Rather, thick filament activation results from increases in myosin phosphorylation that release a fraction of heads to produce force, with the remainder staying in the ordered IHM configuration. After the tetanus, the released heads slowly recover toward the resting, helically ordered state. During this time the released heads remain close to actin and can quickly rebind, enhancing the force produced by posttetanic twitches, structurally explaining posttetanic potentiation. Taken together, these results suggest that, in addition to stretch activation in insects, two other mechanisms for thick filament activation have evolved to disrupt the interactions that establish the relaxed helices of IHMs: one in invertebrates, by either regulatory light-chain phosphorylation (as in arthropods) or Ca 2+ -binding (in mollusks, lacking phosphorylation), and another in vertebrates, by mechanosensing.

59 BASIC BIOLOGICAL SCIENCES↗

Structural complexity in the f -block: small deviations of the complexation of lanthanides by O,Oʹ -diethylmonothiophosphate

Dithiophosphinic acids undergo radiolytic degradation during the extraction of actinides in used nuclear fuel. These will degrade into monothiophosphinic acids and then to phosphinic acids. To elucidate how the complexes that are formed during these radioactive separations change as the ligand degrades, the mixed donor ligand O,O′- diethylmonothiophosphate is chosen as an analog for the monothiophosphinic intermediate. Herein, the monothiophosphate complexes Ln 2 (OPS(OEt) 2 ) 6 (H 2 O) 8 (Ln = La) (La 2 L 6 H 2 O), Ln 2 (OPS(OEt) 2 ) 6 (EtOH) 4 (Ln = La) (La 2 L 6 EtOH), K 2 [Ln(OPS(OEt) 2 ) 5 (H 2 O) 2 ]·H 2 O·CH 2 Cl 2 , (Ln = Ce) (CeL 5 -α), K 2 [Ln(OPS(OEt) 2 ) 5 (H 2 O) 2 ]·H 2 O·CH 2 Cl 2 , (Ln = Pr) (PrL 5 -β), K[Ln(OPS(OEt) 2 ) 4 (H 2 O) 3 ], (Ln = Pr, Sm-Er) (ML 4 ), and K 3 [Ln(OPS(OEt) 2 ) 6 ], (Ln = Dy) (DyL 6 ) were synthesized and characterized using single-crystal X-ray diffraction and optical spectroscopy. Although the lanthanides contract in a nearly linear fashion, the structural changes observed as the f-block is traversed in these compounds are not necessarily a hard line but more so a blend of different structure types possible for each f-element. Furthermore, comparison of the Ln−O bond lengths shows a nearly linear contraction, but the Ln−S bond lengths do not monotonically decrease because of the hard Lewis acidity of the Ln 3+ cations.

37 - INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL C↗

Plasma dynamics, instabilities and OH generation in a pulsed atmospheric pressure plasma with liquid cathode: a diagnostic study

Abstract While plasma–liquid interactions have been an important focus in the plasma research community, the impact of the strong coupling between plasma and liquid on plasma properties and processes remains not fully understood. In this work, we report on the impact of the applied voltage, pulse width and liquid conductivity on the plasma morphology and the OH generation for a positive pulsed DC atmospheric pressure plasma jet with He–0.1% H 2 O mixture interacting with a liquid cathode. We adopted diagnostic techniques of fast imaging, 2D laser induced fluorescence of OH and Thomson scattering spectroscopy. We show that plasma instabilities and enhanced evaporation occur and have a significant impact on the OH generation. At elevated plasma energies, it is found that the plasma contracts due to a thermal instability through Ohmic heating and the contraction coincides with a depletion in the OH density in the core due to electron impact dissociation. For lower plasma energies, the instability is suppressed/delayed by the equivalent series resistor of the liquid electrode. An estimation of the energy flux from the plasma to the liquid shows that the energy flux of the ions released into the liquid by positive ion hydration is dominant, and significantly larger than the energy needed to evaporate sufficient amount of water to account for the measured H 2 O concentration increase near the plasma–liquid interface.

Physics↗

Cosmic topology. Part I. Limits on orientable Euclidean manifolds from circle searches

The Einstein field equations of general relativity constrain the local curvature at every point in spacetime, but say nothing about the global topology of the Universe. Cosmic microwave background anisotropies have proven to be the most powerful probe of non-trivial topology since, within ΛCDM, these anisotropies have well-characterized statistical properties, the signal is principally from a thin spherical shell centered on the observer (the last scattering surface), and space-based observations nearly cover the full sky. The most generic signature of cosmic topology in the microwave background is pairs of circles with matching temperature and polarization patterns. No such circle pairs have been seen above noise in the WMAP or Planck temperature data, implying that the shortest non-contractible loop around the Universe through our location is longer than 98.5% of the comoving diameter of the last scattering surface. We translate this generic constraint into limits on the parameters that characterize manifolds with each of the nine possible non-trivial orientable Euclidean topologies, and provide a code which computes these constraints. In all but the simplest cases, the shortest non-contractible loop in the space can avoid us, and be shorter than the diameter of the last scattering surface by a factor ranging from 2 to at least 6. Further, this result implies that a broader range of manifolds is observationally allowed than widely appreciated. Probing these manifolds will require more subtle statistical signatures than matched circles, such as off-diagonal correlations of harmonic coefficients.

79 ASTRONOMY AND ASTROPHYSICS↗

Collapsing universe before time

Abstract In this note, we prove analytic bounds on the equation of state of a cosmological fluid composed of an arbitrary number of canonical scalars evolving in a negative multi-exponential potential. Because of the negative energy, the universe is contracting and eventually undergoes a big crunch. A contracting universe is a fundamental feature of models of ekpyrosis and cyclic universes, which are a proposed alternative to cosmic inflation to solve the flatness and horizon problems. Our analytic bounds set quantitative constraints on general effective theories of ekpyrosis. When applied to specific top-down constructions, our bounds can be used to determine whether ekpyrosis could in principle be realized. We point out some possible sources of tension in realizing the ekpyrotic universe in controlled constructions of string theory.

Astronomy & Astrophysics↗

Particle production by time-varying dark energy and the end of cosmic expansion

We consider various possible consequences of time-varying dark energy due to a quintessence scalar field whose energy density is partially converted to particles as the field evolves down its potential. This particle production acts as a source of thermal friction on the field that can make it difficult to distinguish whether dark energy is due to a radiating field rolling down a steep potential, a purely self-interacting field moving down a flatter potential, or a cosmological constant. By reducing the acceleration of the scalar field, thermal friction increases the amount of accelerated expansion and can cause a sizable bump in the quintessence equation of state. We take special interest in the case where a steep potential rapidly changes from positive to negative as the field evolves, resulting in the end of cosmic expansion and the beginning of contraction. Even in this case, we find that thermal friction lengthens the period of accelerated expansion and consequently delays the end of cosmic expansion, making it challenging to detect the impending transition to contraction using conventional cosmological tests. However, particle production can also provide alternative avenues for detection by generating a background of thermal dark radiation, partly comprised of neutrinos or other particles, whose energy density exceeds the remnant photon energy density.

cosmological neutrinos↗

Climate change and the future productivity and distribution of crab in the Bering Sea

Crab populations in the eastern Bering Sea support some of the most valuable fisheries in the United States, but their future productivity and distribution are uncertain. Here, we explore observed changes in the productivity and distribution for snow crab, Tanner crab, and Bristol Bay red king crab. We link historical indices of environmental variation and predator biomass with observed time series of centroids of abundance and extent of crab stock distribution; we also fit stock–recruit curves including environmental indices for each stock. We then project these relationships under forcing from global climate models to forecast potential productivity and distribution scenarios. Our results suggest that the productivity of snow crab is negatively related to the Arctic Oscillation (AO) and positively related to ice cover; Tanner crab’s productivity and distribution are negatively associated with cod biomass and sea surface temperature. Aspects of red king crab distribution and productivity appear to be related to bottom temperature, ice cover, the AO, and/or cod biomass. Projecting these relationships forward with available forecasts suggests that Tanner crab may become more productive and shift further offshore, red king crab distribution may contract and move north, and productivity may decrease for snow crab as the population contracts northward.

54 ENVIRONMENTAL SCIENCES↗

The structural OFF and ON states of myosin can be decoupled from the biochemical super- and disordered-relaxed states

Abstract There is a growing awareness that both thick-filament and classical thin-filament regulations play central roles in modulating muscle contraction. Myosin ATPase assays have demonstrated that under relaxed conditions, myosin may reside either in a high-energy-consuming disordered-relaxed (DRX) state available for binding actin to generate force or in an energy-sparing super-relaxed (SRX) state unavailable for actin binding. X-ray diffraction studies have shown that the majority of myosin heads are in a quasi-helically ordered OFF state in a resting muscle and that this helical ordering is lost when myosin heads are turned ON for contraction. It has been assumed that myosin heads in SRX and DRX states are equivalent to the OFF and ON states, respectively, and the terms have been used interchangeably. In this study, we use X-ray diffraction and ATP turnover assays to track the structural and biochemical transitions of myosin heads, respectively, induced with either omecamtiv mecarbil (OM) or piperine in relaxed porcine myocardium. We find that while OM and piperine induce dramatic shifts of myosin heads from the OFF to the ON state, there are no appreciable changes in the population of myosin heads in the SRX and DRX states in both unloaded and loaded preparations. Our results show that biochemically defined SRX and DRX can be decoupled from structurally defined OFF and ON states. In summary, while SRX/DRX and OFF/ON transitions can be correlated in some cases, these two phenomena are measured using different approaches, reflect different properties of the thick filament, and should be investigated and interpreted separately.

Jani, Vivek P. (ORCID:0000000238114973)↗

Nonlinear elasticity and short-range mechanical coupling govern the rate and symmetry of mouth opening in Hydra

Hydra has a tubular bilayered epithelial body column with a dome-shaped head on one end and a foot on the other. Hydra lacks a permanent mouth: its head epithelium is sealed. Upon neuronal activation, a mouth opens at the apex of the head which can exceed the body column diameter in seconds, allowing Hydra to ingest prey larger than itself. While the kinematics of mouth opening are well characterized, the underlying mechanism is unknown. We show that Hydra mouth opening is generated by independent local contractions that require tissue-level coordination. We model the head epithelium as an active viscoelastic nonlinear spring network. The model reproduces the size, timescale and symmetry of mouth opening. It shows that radial contractions, travelling inwards from the outer boundary of the head, pull the mouth open. Nonlinear elasticity makes mouth opening larger and faster, contrary to expectations. Here, the model correctly predicts changes in mouth shape in response to external forces. By generating innervated : nerve-free chimera in experiments and simulations, we show that nearest-neighbor mechanical signaling suffices to coordinate mouth opening. Hydra mouth opening shows that in the absence of long-range chemical or neuronal signals, short-range mechanical coupling is sufficient to produce long-range order in tissue deformations.

59 BASIC BIOLOGICAL SCIENCES↗

Stoichiometry dependent properties of cerium hydride: An active learning developed interatomic potential study

Cerium hydride has a variety of interesting properties, including a known lattice contraction and densification with increasing hydrogen content. However, precise stoichiometric control is not experimentally straightforward and ab initio approaches are not computationally feasible for many properties such as melting and low temperature diffusion. Therefore, we develop a machine-learned interatomic potential for cerium hydride that is valid for H to Ce ratios from 2.0 to 3.0. A query-by-committee active learning approach is used to develop the training set. Leveraging classical molecular dynamics simulations, we assess a range of properties and provide fundamental mechanisms for the trends with stoichiometry. Finally, a majority of the properties follow the trend of lattice contraction, being governed by the stronger lattice binding induced by adding octahedral atoms.

36 MATERIALS SCIENCE↗

Cosmological Bounces Induced by a Fermion Condensate

Here, we show that it is possible for fermion condensation of the Nambu-Jona-Lasinio type to induce a non-singular bounce that smoothly connects a phase of slow contraction to a phase of expansion. A chiral condensate—a non-zero vacuum expectation value of the spinor bilinear <Ψ¯Ψ> — can form spontaneously after a slow contraction phase smooths and flattens the universe and the Ricci-curvature exceeds a critical value. In this approach, a high density of spin-aligned free fermions is not required, which avoids the problem of generating a large anisotropy and initiating chaotic mixmaster behavior during the bounce phase.

79 ASTRONOMY AND ASTROPHYSICS↗

Tensorized Feature Spaces for Feature Explosion

In this paper 1 1 This manuscript has been authored by UT-Battelle, LLC under Contract No. DE-AC05-000R22725 with the U.S. Department of Energy. The United States Government retains and the publisher, by accepting the article for publication, acknowledges that the United States Government retains a nonexclusive, paid-up, irrevocable, worldwide license to publish or reproduce the published form of this manuscript, or allow others to do so, for United States Government purposes. The Department of Energy will provide public access to these results of federally sponsored research in accordance with the DOE Public Access Plan (http://energy.gov/downloads/doe-public-access-plan) This research used resources of the Oak Ridge Leadership Computing Facility, which is a DOE Office of Science User Facility supported under Contract DE-AC05-000R22725., we present a novel framework that uses tensor factorization to generate richer feature spaces for pixel classification in hyperspectral images. In particular, we assess the performance of different tensor rank decomposition methods as compared to the traditional kernel-based approaches for the hyperspectral image classification problem. We propose Orion, which takes as input a hyperspectral image tensor and a rank and outputs an enhanced feature space from the factor matrices of the decomposed tensor. Our method is a feature explosion technique that inherently maps low dimensional input space in $\mathbb{R}^{K}$ to high dimensional space in $\mathbb{R}^{R}$ , where $R\gg K$ , say in the order of 1000x, like a kernel. We show how the proposed method exploits the multi-linear structure of hyperspectral three dimensional tensor. We demonstrate the effectiveness of our method with experiments on three publicly available hyperspectral datasets with labeled pixels and compare their classification performance against traditional linear and non-linear supervised learning methods such as SVM with Linear, Polynomial, RBF kernels, and the Multi-Layer Perceptron model. Finally, we explore the relationship between the rank of the tensor decomposition and the classification accuracy using several hyperspectral datasets with ground truth.

Pasricha, Ravdeep Singh↗

MICCO: An Enhanced Multi-GPU Scheduling Framework for Many-Body Correlation Functions

Calculation of many-body correlation functions is one of the critical kernels utilized in many scientific computing areas, especially in Lattice Quantum Chromodynamics (Lattice QCD). It is formalized as a sum of a large number of contraction terms each of which can be represented by a graph consisting of vertices describing quarks inside a hadron node and edges designating quark propagations at specific time intervals. Due to its computation- and memory-intensive nature, real-world physics systems (e.g., multi-meson or multi-baryon systems) explored by Lattice QCD prefer to leverage multi-GPUs. Different from general graph processing, many-body correlation function calculations show two specific features: a large number of computation-/data-intensive kernels and frequently repeated appearances of original and intermediate data. The former results in expensive memory operations such as tensor movements and evictions. The latter offers data reuse opportunities to mitigate the data-intensive nature of many-body correlation function calculations. However, existing graph-based multi-GPU schedulers cannot capture these data-centric features, thus resulting in a sub-optimal performance for many-body correlation function calculations. To address this issue, this paper presents a multi-GPU scheduling framework, MICCO, to accelerate contractions for correlation functions particularly by taking the data dimension (e.g., data reuse and data eviction) into account. This work first performs a comprehensive study on the interplay of data reuse and load balance, and designs two new concepts: local reuse pattern and reuse bound to study the opportunity of achieving the optimal trade-off between them. Based on this study, MICCO proposes a heuristic scheduling algorithm and a machine-learning-based regression model to generate the optimal setting of reuse bounds. Specifically, MICCO is integrated into a real-world Lattice QCD system, Redstar, for the first time running on multiple GPUs. The evaluation demonstrates MICCO outperforms other state-of-art works, achieving up to 2.25× speedup in synthesized datasets, and 1.49× speedup in real-world correlation functions.

Wang, Qihan↗

A PEPS Plugin for TNQVM

This work introduces an extension to the Tensor Network Quantum Virtual Machine (TNQVM) tool, enhancing the existing stack of ExaScale Tensor Network (ExaTN), ExaScale Accelerator (XACC), and TNQVM. It features a new plugin that enables efficient simulation of a Projected Entangled Pair State (PEPS), a 2D tensor network. To improve simulation efficiency for PEPS, we have implemented the snake boundary contraction algorithm. By integrating this capability into the existing stack, we enhance the overall functionality and versatility of the framework. We tested this new PEPS topology for a simple GHZ bell-pair generation quantum circuit and saw that its runtime is very close to that of the MPS topology. We estimate that the real potential of the PEPS topology becomes discernible when quantum circuits with multidimensional entanglement are simulated using tensor networks. In such cases, 1D tensor networks fail to represent or contract them efficiently.

Chundury, Srikar↗