Search NASA⌕ Search

SEARCH · Search NASA

Results for “parallel clustering”

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 757 records · Page 42

A simple method for obtaining a three-dimensional proton distribution function from Voyager plasma data

The main sensor of the Vogager plasma experiment consists of a cluster of three, modulated-grid Faraday cups whose normals are arranged symmetrically about the symmetry axis of the cluster at an angle of 20 degrees to that axis. In interplanetary space, each cup explores the positive ion distribution by accepting particles from contiguous slices in velocity space. The slices are narrow in the direction of the normal to the modulating grid but are broad in planes parallel to that grid. The resulting three sets of measurements can be combined to yield the three-dimensional distribution function in the following way: the distribution function is assumed to be gyrotropic. For each value of speed in a frame of reference moving with the bulk velocity of the solar wind, the variation of the distribution function with angle from the field direction is represented by a series of Legendre polynomials. Effects such as double-streaming and heat flow can be well represented by using only the first three terms of the series which are fully specified by the measurements. Examples of the use of this method in the analysis of Voyager data are shown.

Olbert, S.↗

Investigations on the usefulness of the Massively Parallel Processor for study of electronic properties of atomic and condensed matter systems

The usefulness of the Massively Parallel Processor (MPP) for investigation of electronic structures and hyperfine properties of atomic and condensed matter systems was explored. The major effort was directed towards the preparation of algorithms for parallelization of the computational procedure being used on serial computers for electronic structure calculations in condensed matter systems. Detailed descriptions of investigations and results are reported, including MPP adaptation of self-consistent charge extended Hueckel (SCCEH) procedure, MPP adaptation of the first-principles Hartree-Fock cluster procedure for electronic structures of large molecules and solid state systems, and MPP adaptation of the many-body procedure for atomic systems.

Das, T. P.↗

Low-Dimensional Water on Ru(0001)Model System for X-ray Absorption Spectroscopy Studies of Liquid Water

We present an x-ray absorption spectroscopy results for fully broken to a complete H-bond network of water molecules on Ru(0001) by varying the morphology from isolated water molecules via two-dimensional clusters to a fully covered monolayer as probed by scanning tunneling microscopy. The sensitivity of x-ray absorption to the symmetry of H-bonding is further elucidated for the amino (-NH{sub 2}) group in glycine adsorbed on Cu(110) where the E-vector is parallel either to the NH donating an H-bond or to the non-H-bonded NH. The results give further evidence for the interpretation of the various spectral features of liquid water and for the general applicability of x-ray absorption spectroscopy to analyze H-bonded systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Supercomputer-Based Ensemble Docking Drug Discovery Pipeline with Application to Covid-19

In this work, we present a supercomputer-driven pipeline for in silico drug discovery using enhanced sampling molecular dynamics (MD) and ensemble docking. Ensemble docking makes use of MD results by docking compound databases into representative protein binding-site conformations, thus taking into account the dynamic properties of the binding sites. We also describe preliminary results obtained for 24 systems involving eight proteins of the proteome of SARS-CoV-2. The MD involves temperature replica exchange enhanced sampling, making use of massively parallel supercomputing to quickly sample the configurational space of protein drug targets. Using the Summit supercomputer at the Oak Ridge National Laboratory, more than 1 ms of enhanced sampling MD can be generated per day. We have ensemble docked repurposing databases to 10 configurations of each of the 24 SARS-CoV-2 systems using AutoDock Vina. Comparison to experiment demonstrates remarkably high hit rates for the top scoring tranches of compounds identified by our ensemble approach. We also demonstrate that, using Autodock-GPU on Summit, it is possible to perform exhaustive docking of one billion compounds in under 24 h. Finally, we discuss preliminary results and planned improvements to the pipeline, including the use of quantum mechanical (QM), machine learning, and artificial intelligence (AI) methods to cluster MD trajectories and rescore docking poses.

60 APPLIED LIFE SCIENCES↗

Performance analysis of distributed symmetric sparse matrix vector multiplication algorithm for multi-core architectures

In this article, sparse matrix vector multiply (SpMVM) is an important kernel that frequently arises in high performance computing applications. Due to its low arithmetic intensity, several approaches have been proposed in literature to improve its scalability and efficiency in large scale computations. In this paper, our target systems are high end multi-core architectures and we use messaging passing interface + open multiprocessing hybrid programming model for parallelism. We analyze the performance of recently proposed implementation of the distributed symmetric SpMVM, originally developed for large sparse symmetric matrices arising in ab initio nuclear structure calculations. We also study important features of this implementation and compare with previously reported implementations that do not exploit underlying symmetry. Our SpMVM implementations leverage the hybrid paradigm to efficiently overlap expensive communications with computations. Our main comparison criterion is the "CPU core hours" metric, which is the main measure of resource usage on supercomputers. We analyze the effects of topology-aware mapping heuristic using simplified network load model. Furthermore, we have tested the different SpMVM implementations on two large clusters with 3D Torus and Dragonfly topology. Our results show that the distributed SpMVM implementation that exploits matrix symmetry and hides communication yields the best value for the "CPU core hours" metric and significantly reduces data movement overheads.

97 MATHEMATICS AND COMPUTING↗

Coupling LAMMPS and the vl3 Framework for Co-Visualization of Atomistic Simulations

In this work we perform in-situ visualization of molecular dynamics simulations, which can help scientists to visualize simulation output on-the-fly, without incurring storage overheads. We present a case study to couple LAMMPS, the large-scale molecular dynamics simulation code with vl3, our parallel framework for large-scale visualization and analysis. Our motivation is to identify effective approaches for covisualization and exploration of large-scale atomistic simulations at interactive frame rates.We propose a system of coupled libraries and describe its architecture, with an implementation that runs on GPU-based clusters. We present the results of strong and weak scalability experiments, as well as future research avenues based on our results.

co-visualization↗

Observations of How Magnetofluid Turbulence Dissipates at Small Scales

The solar wind is a turbulent magneto fluid that can be studied intensively at multiple scales. Investigations using single spacecraft have revealed much about the properties of the solar wind throughout the heliosphere (from 0.3 AU to 100 AU). More recently, data from multiple spacecraft have provided further details of both the statistical properties of the turbulence and its small-scale structure. In particular, high time resolution magnetic field measurements from the four Cluster spacecrafl have led to the conclusion that at spatial scales of order the proton inertial length and smaller, the turbulence becomes strongly anisotropic and the power in fluctuations that are perpendicular to the (local) magnetic field is measured to be much larger than that in fluctuations that are parallel to the magnetic field. As the spatial scales approach the electron inertial length, the power is almost completely dissipated. Various analysis techniques and theoretical ideas have been put forward to account for the properties of those measurements. The talk will describe the current state of observations, theory and simulations.

Goldstein, M. L.↗

Self-optimized construction of transition rate matrices from accelerated atomistic simulations with Bayesian uncertainty quantification

Here, a massively parallel method to build large transition rate matrices from temperature-accelerated molecular dynamics trajectories is presented. Bayesian Markov model analysis is used to estimate the expected residence time in the known state space, providing crucial uncertainty quantification for higher-scale simulation schemes such as kinetic Monte Carlo or cluster dynamics. The estimators are additionally used to optimize where exploration is performed and the degree of temperature acceleration on the fly, giving an autonomous, optimal procedure to explore the state space of complex systems. The method is tested against exactly solvable models and used to explore the dynamics of C15 interstitial defects in iron. Our uncertainty quantification scheme allows for accurate modeling of the evolution of these defects over timescales of several seconds.

36 MATERIALS SCIENCE↗

A Study of Parallel Scalability and Dynamic Workload Balancing in GlennICE

The Glenn Icing Computational Environment (GlennICE) is a computational tool designed to calculate ice growth on complex three-dimensional geometries. It utilizes user-supplied computational fluid dynamics solutions for the geometry of interest. Key developments include advancements in convergence of collection efficiency, trajectory optimization, and refinement methodology. These improvements have significantly enhanced GlennICE’s efficiency for practical engineering applications. A recent study focused on benchmarking GlennICE’s scalability in a parallel environment using static scheduling. Findings indicated a potential twofold increase in efficiency through workload balance enhancements. This paper presents an analysis of the solver’s new workload balancing improvements, incorporating shared memory and dynamic scheduling routines. Results demonstrate a highly efficient and consistent algorithm across high-performance computing clusters.

Computational Icing↗

A Study of Parallel Scalability and Dynamic Workload Balancing in GlennICE

The Glenn Icing Computational Environment (GlennICE) is a computational tool designed to calculate ice growth on complex three-dimensional geometries. It utilizes user-supplied computational fluid dynamics solutions for the geometry of interest. Key developments include advancements in convergence of collection efficiency, trajectory optimization, and refinement methodology. These improvements have significantly enhanced GlennICE’s efficiency for practical engineering applications. A recent study focused on benchmarking GlennICE’s scalability in a parallel environment using static scheduling. Findings indicated a potential twofold increase in efficiency through workload balance enhancements. This paper presents an analysis of the solver’s new workload balancing improvements, incorporating shared memory and dynamic scheduling routines. Results demonstrate a highly efficient and consistent algorithm across high-performance computing clusters.

Computational Icing↗

Structural evidence for intermediates during O 2 formation in photosystem II

In natural photosynthesis, the light-driven splitting of water into electrons, protons and molecular oxygen forms the first step of the solar-to-chemical energy conversion process. The reaction takes place in photosystem II, where the Mn 4 CaO 5 cluster first stores four oxidizing equivalents, the S 0 to S 4 intermediate states in the Kok cycle, sequentially generated by photochemical charge separations in the reaction center and then catalyzes the O–O bond formation chemistry. Here, we report room temperature snapshots by serial femtosecond X-ray crystallography to provide structural insights into the final reaction step of Kok’s photosynthetic water oxidation cycle, the S 3 →[S 4 ]→S 0 transition where O 2 is formed and Kok’s water oxidation clock is reset. Our data reveal a complex sequence of events, which occur over micro- to milliseconds, comprising changes at the Mn 4 CaO 5 cluster, its ligands and water pathways as well as controlled proton release through the hydrogen-bonding network of the Cl1 channel. Importantly, the extra O atom O x , which was introduced as a bridging ligand between Ca and Mn1 during the S 2 →S 3 transition disappears or relocates in parallel with Y z reduction starting at approximately 700 μs after the third flash. The onset of O 2 evolution, as indicated by the shortening of the Mn1–Mn4 distance, occurs at around 1,200 μs, signifying the presence of a reduced intermediate, possibly a bound peroxide.

59 BASIC BIOLOGICAL SCIENCES↗

Evidence of the ion cyclotron resonance at proton kinetic scales in the solar wind

The k-filtering and wave polarization techniques are applied to Cluster magnetic field data at four intervals of fast free (not connected to Earth's foreshock) solar wind. In addition to the commonly observed population of magnetic field fluctuations propagating at quasi-perpendicular angles to the global mean field B{sub 0}, a population of fluctuations propagating at quasi-parallel angles are also observed when no local plasma instability can be identified. At low wavenumbers (kv{sub A}/Ω{sub p} ≤ 0.6) both components are present, and have powers of similar strength, while at higher wavenumbers (kv{sub A}Ω{sub p} > 0.6) only the perpendicular component can be identified. Here v{sub A} and Ω{sub p} are the Alfvén speed and the proton (angular) gyration frequency, k denotes the wavenumber. The dispersion curve obtained (and the ratio of the magnitudes of left and right hand polarized fluctuations) shows a sharp transition (and decrease) at (kv{sub A}/Ω{sub p}) ∼ 0.6, which corresponds to a spacecraft frequency f {sub sc} = 0.32 Hz. At higher wavenumbers the measurements agree with previous studies indicating the presence of kinetic Alfvén wave (KAW) turbulence or a mixture of KAW turbulence and convected structures. The parallel component displays a scaling of k{sub ∥}{sup −δ} where the spectral index δ ranges between about 5/3 and 2 although the number of data points does not warrant an accurate determination. The origins of the observed parallel component are unclear. The presence of a parallel component is consistent with ion cyclotron waves, and we interpret the sharp drop of this component at higher wavenumbers as evidence that the ion cyclotron resonance is occurring in the solar wind intervals studied.

79 ASTRONOMY AND ASTROPHYSICS↗

In Vitro Experimental Model to Investigate the Biological Effects across the Bragg Curve of High-LET Radiation

The space environment consists of a varying field of radiation particles including high energy ions, with a spacecrafts shielding material providing the only major protection to astronauts from harmful exposure. Unlike lowLET gamma or Xrays, the presence of shielding does not always reduce the radiation risks for energetic charged particle exposure since the dose delivered by the charged particle increases sharply as the particle approaches the end of its range, a position known as the Bragg peak and the correlating spatial dose distribution identified as the Bragg curve. The Bragg curve does not necessarily represent the biological damage along the particle traversal since biological effects are influenced by the track structure of both primary and secondary particles. Therefore, the biological Bragg curve is dependent of the energy and the type of the primary particle, and may vary for different biological endpoints. Here we describe a unique irradiation geometry and experimental system to measure the biological response across the Bragg curve in one consistent biological sample. Polyethylene shielding was used to achieve a Bragg curve distribution with the beam geometry parallel to a monolayer of fibroblast cells. We present data that highlights the differential formation of DNA double strand breaks (DSBs) and chromosomal deletions across the Bragg curve in human fibroblasts irradiated with 600 MeV/nucleon iron ion beams. Qualitative analyses of gammaH2AX fluorescence, a known marker of DSBs, indicated potentially increased clustering of DNA damage before the Bragg peak, enhanced homogenous distribution at the peak, and provided visual evidence of high linear energy transfer (LET) particle traversal of cells beyond the Bragg peak in agreement with one-dimensional transport approximations. A biological response curve generated for micronuclei induction across the Bragg curve for 600 MeV/n Fe ions did not reveal an increase in the yield of micronuclei at the Bragg peak location. Assessment of such biological parameters employing the described in vitro experimental system may provide improved platforms to measure a number of biological consequences of shielding materials across the Bragg curve for high charge and energy (HZE) ions.

Desai, N.↗

A Case Study of LLVM-Based Analysis for Optimizing SIMD Code Generation

This paper presents a methodology for using LLVM-based tools to tune the DCA++ (dynamical cluster approximation) application that targets the new ARM A64FX processor. The goal is to describe the changes required for the new architecture and generate efficient single instruction/multiple data (SIMD) instructions that target the new Scalable Vector Extension instruction set. During manual tuning, the authors used the LLVM tools to improve code parallelization by using OpenMP SIMD, refactored the code and applied transformation that enabled SIMD optimizations, and ensured that the correct libraries were used to achieve optimal performance. By applying these code changes, code speed was increased by 1.98× and 78 GFlops were achieved on the A64FX processor. The authors aim to automatize parts of the efforts in the OpenMP Advisor tool, which is built on top of existing and newly introduced LLVM tooling.

Huber, Joseph↗

Structural modulation and luminescent properties of four Cd{sup II} coordination architectures based on 3-(pyridin-4-yl)-5-(pyrazin-2-yl)-1H-1,2,4-triazole and flexible/rigid dicarboxylate ligands

To systematically investigate the influence of the flexible or rigid auxiliary ligands on the structures and properties of transition metal compounds, we synthesized four new d{sup 10} coordination polymers (CPs) from 3-(pyridin-4-yl)-5-(pyrazin-2-yl)-1H-1,2,4-triazole (4-Hpzpt) and flexible/rigid dicarboxylate ligands, [Cd(4-pzpt){sub 2}]{sub n} (1), [Cd{sub 3}(4-pzpt){sub 2}(suc){sub 2}]{sub n} (2), [Cd{sub 2}(4-Hpzpt)(nbc){sub 2}(H{sub 2}O)]{sub n} (3) and ([Cd{sub 2}(4-pzpt){sub 2}(tfbdc)(H{sub 2}O){sub 4}]·H{sub 2}O){sub n} (4) (H{sub 2}suc=1,2-ethanedicarboxylic acid, H{sub 2}nbc=hthalene-1,4-dicarboxylic acid, H{sub 2}tfbdc =2,3,5,6-tetrafluoroterephthalic acid). Single crystal X-ray analysis indicates that compound 1 shows a 4{sup 4}-sql layer, which is extended to a 3D network via nonclassical C–H{sup …}N hydrogen bonds. Compound 2 possesses a 6-connected pcu-4{sup 12}0.6{sup 3} net composed of trinuclear Cd{sup II}-clusters. Compound 3 represents a rare 3D (3,4,4,5)-connected topology with a Schläfli symbol of (4·6·7)(4·5{sup 3}·7{sup 2})(5{sup 3}·6·7·9)(4{sup 2}·5{sup 5}·6·7{sup 2}). Compound 4 exhibits a 2D+2D→2D parallel interpenetrated 6{sup 3}-hcb network. The adjacent 2D networks are interdigitated with each other to form the resulting 3D supramolecular architecture through classical O–H{sup …}N and O–H{sup …}O hydrogen bonds. Structural diversities indicate that the nature of flexible/rigid-dicarboxlates plays crucial roles in modulating structures of these compounds. Moreover, the luminescent properties of them have been briefly investigated. - Graphical abstract: Four new Cd{sup II} coordination architectures constructed from the primary ligand 4-Hpzpt and flexible/rigid dicarboxylate coligands. Structural diversities indicate that the nature of flexible/rigid-dicarboxlates plays crucial roles in modulating structures of these compounds. And more, the thermal stability and luminescence are discussed. - Highlights: • Four new Cd{sup II} compounds with 4-Hpzpt and flexible/rigid dicarboxylate coligands. • Structural analysis of all compounds. • luminescent property of all compounds.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Direct conversion of methane to aromatics and hydrogen via a heterogeneous trimetallic synergistic catalyst

Abstract Non-oxidative methane dehydro-aromatization reaction can co-produce hydrogen and benzene effectively on a molybdenum-zeolite based thermochemical catalyst, which is a very promising approach for natural-gas upgrading. However, the low methane conversion and aromatics selectivity and weak durability restrain the realistic application for industry. Here, a mechanism for enhancing catalysis activity on methane activation and carbon-carbon bond coupling has been found to promote conversion and selectivity simultaneously by adding platinum–bismuth alloy cluster to form a trimetallic catalyst on zeolite (Pt-Bi/Mo/ZSM-5). This bimetallic alloy cluster has synergistic interaction with molybdenum: the formed CH 3 * from Mo 2 C on the external surface of zeolite can efficiently move on for C-C coupling on the surface of Pt-Bi particle to produce C 2 compounds, which are the key intermediates of oligomerization. This pathway is parallel with the catalysis on Mo inside the cage. This catalyst demonstrated 18.7% methane conversion and 69.4% benzene selectivity at 710 °C. With 95% methane/5% nitrogen feedstock, it exhibited robust stability with slow deactivation rate of 9.3% after 2 h and instant recovery of 98.6% activity after regeneration in hydrogen. The enhanced catalytic activity is strongly associated with synergistic interaction with Mo and ligand effects of alloys by extensive mechanism studies and DFT calculation.

03 NATURAL GAS↗

CRAGE enables rapid activation of biosynthetic gene clusters in undomesticated bacteria

Advances in the prediction of biosynthetic gene clusters (BGCs) from genomic and metagenomic data have made the large-scale identification of candidate secondary metabolite biosynthesis pathways routine. However, most predicted BGCs cannot be successfully expressed in the limited range of expression hosts currently available, which prevents the identification and characterization of their respective products. Here we show that parallelized use of a phylogenetically and physiologically diverse panel of purpose-engineered chassis strains offers a versatile strategy for the rapid activation of BGCs. We domesticated an ensemble of 25 diverse ?-proteobacteria from 11 different genera by integrating a unified landing pad that enables efficient introduction of large BGCs into their genome. We deployed this chassis-independent recombinase-assisted genome engineering (CRAGE) strategy to express several previously posited but experimentally elusive non-ribosomal peptide synthase (NRPS) and NRPSpolyketide synthase (PKS) hybrid BGCs from Photorhabdus luminescens. Successful activation of 6 BGCs identified 21 products, including several with likely roles in host-microbe interactions. Taken together, we illustrate how the regulatory and physiological diversity naturally present across bacteria can be harnessed for the exploration of novel bioproducts and provide a sizable panel of domesticated chassis strains to enable this strategy.

applied microbiology, genetic engineering, industr↗