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 811 records · Page 45

Brassica rapa Domestication: Untangling Wild and Feral Forms and Convergence of Crop Morphotypes

The study of domestication contributes to our knowledge of evolution and crop genetic resources. Human selection has shaped wild Brassica rapa into diverse turnip, leafy, and oilseed crops. Despite its worldwide economic importance and potential as a model for understanding diversification under domestication, insights into the number of domestication events and initial crop(s) domesticated in B. rapa have been limited due to a lack of clarity about the wild or feral status of conspecific noncrop relatives. To address this gap and reconstruct the domestication history of B. rapa, we analyzed 68,468 genotyping-by-sequencing-derived single nucleotide polymorphisms for 416 samples in the largest diversity panel of domesticated and weedy B. rapa to date. To further understand the center of origin, we modeled the potential range of wild B. rapa during the mid-Holocene. Our analyses of genetic diversity across B. rapa morphotypes suggest that noncrop samples from the Caucasus, Siberia, and Italy may be truly wild, whereas those occurring in the Americas and much of Europe are feral. Clustering, tree-based analyses, and parameterized demographic inference further indicate that turnips were likely the first crop type domesticated, from which leafy types in East Asia and Europe were selected from distinct lineages. These findings clarify the domestication history and nature of wild crop genetic resources for B. rapa, which provides the first step toward investigating cases of possible parallel selection, the domestication and feralization syndrome, and novel germplasm for Brassica crop improvement.

59 BASIC BIOLOGICAL SCIENCES↗

Software for the frontiers of quantum chemistry: An overview of developments in the Q-Chem 5 package

This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange-correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear-electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an "open teamware" model and an increasingly modular design.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modeling the Performance of an Algebraic Multigrid Cycle Using Hybrid MPI/OpenMP

The rise of multicore cluster architectures has led to intense interest in using a combination of MPI and OpenMP to more effectively program these machines. We present a performance model for hybrid implementation of the solve cycle of algebraic multigrid (AMG), a popular iterative solver for large sparse linear systems and a key component of many scientific simulations. We validate the model on two leading parallel platforms, and discuss implications for applications programmed in a hybrid model on future machines.

Gahvari, Hormozd↗

Summary of VERA Core Simulator Performance Improvements

This report describes the performance improvements made to the VERA Core Simulator (VERA CS) during FY2016. At the beginning of the year, a timing study was performed to identify areas of the code where the run-times could be improved. This study identified six broad areas for improvement: • Cross section lookup and CMFD stabilization, • CTF parallel partitioning, • MOC methodology, • Subgroup methodology, • MOC parallel partitioning, and • CMFD Krylov solution techniques. The goal of this work was to reduce the VERA-CS run-times so a complete cycle depletion could be run with 1,000 cores overnight (i.e., the “1,000 core depletion” criterion). This goal was set by our industry partners so they could achieve reasonable turnaround on medium-sized computer clusters. The results presented in this report show a code speed-up by a factor of 5, which successfully satisfies the goal set by industry.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Validation of astronaut psychological select-in criteria

An optional astronaut selection strategy would select-in individuals on the basis of personality attributes associated with superior performance. Method: A test battery, the Astronaut Personal Characteristics Inventory (ASTROPCI) was developed which assesses positive and negative components of achievement, motivation, and interpersonal orientations and skills. The battery was administered to one hundred three astronaut candidates and sixty-six current U.S. Shuttle astronauts. To determine performance, a series of conceptual areas related to space flight performance were defined. Astronauts rated their peers on each of these dimensions. Ratings were obtained on all eighty-four current astronauts (excluding those selected in 1990). In addition to peer ratings, supervisor assessments of the same dimensions were obtained for each astronaut. Results: Cluster and factor analysis techniques were employed to isolate subgroups of astronauts. Those astronauts with both high achievement needs and interpersonal skills were most often rated among the top five by their peers and least often rated among the lowest five. A number of scales discriminated between astronauts rated high and low on one or more performance dimensions. Conclusions: The results parallel findings from the personality assessment of individuals in other demanding professions, including aircraft pilots and research scientists, suggesting that personality factors are significant determinants to performance in the space environment.

Rose, R. M.↗

Scalable Multiprocessor for High-Speed Computing in Space

A report discusses the continuing development of a scalable multiprocessor computing system for hard real-time applications aboard a spacecraft. "Hard realtime applications" signifies applications, like real-time radar signal processing, in which the data to be processed are generated at "hundreds" of pulses per second, each pulse "requiring" millions of arithmetic operations. In these applications, the digital processors must be tightly integrated with analog instrumentation (e.g., radar equipment), and data input/output must be synchronized with analog instrumentation, controlled to within fractions of a microsecond. The scalable multiprocessor is a cluster of identical commercial-off-the-shelf generic DSP (digital-signal-processing) computers plus generic interface circuits, including analog-to-digital converters, all controlled by software. The processors are computers interconnected by high-speed serial links. Performance can be increased by adding hardware modules and correspondingly modifying the software. Work is distributed among the processors in a parallel or pipeline fashion by means of a flexible master/slave control and timing scheme. Each processor operates under its own local clock; synchronization is achieved by broadcasting master time signals to all the processors, which compute offsets between the master clock and their local clocks.

Lux, James↗

Computational Inference of Vibratory System with Incomplete Modal Information Using Parallel, Interactive and Adaptive Markov Chains

Inverse analysis of vibratory system is an important subject in fault identification, model updating, and robust design and control. It is challenging subject because 1) the problem is oftentimes underdetermined while the measurements are limited and/or incomplete; 2) many combinations of parameters may yield results that are similar with respect to actual response measurements; and 3) uncertainties inevitably exist. The aim of this research is to leverage upon computational intelligence through statistical inference to facilitate an enhanced, probabilistic framework using incomplete modal response measurement. This new framework is built upon efficient inverse identification through optimization, whereas Bayesian inference is employed to account for the effect of uncertainties. To overcome the computational cost barrier, we adopt Markov chain Monte Carlo (MCMC) to characterize the target function/distribution. Instead of using single Markov chain in conventional Bayesian approach, we develop a new sampling theory with multiple parallel, interactive and adaptive Markov chains and incorporate it into Bayesian inference. This can harness the collective power of these Markov chains to realize the concurrent search of multiple local optima. The number of required Markov chains and their respective initial model parameters are automatically determined via Monte Carlo simulation-based sample pre-screening followed by K-means clustering analysis. These enhancements can effectively address the aforementioned challenges in finite element inverse analysis. The validity of this framework is systematically demonstrated through case studies.

K Zhou↗

Massively parallel X-ray scattering simulations

Although present X-ray scattering techniques can provide tremendous information on the nano-structural properties of materials that are valuable in the design and fabrication of energy-relevant nano-devices, a primary challenge remains in the analyses of such data. In this paper we describe a high-performance, flexible, and scalable Grazing Incidence Small Angle X-ray Scattering simulation algorithm and codes that we have developed on multi-core/CPU and many-core/GPU clusters. We discuss in detail our implementation, optimization and performance on these platforms. Our results show speedups of ~125x on a Fermi-GPU and ~20x on a Cray-XE6 24-core node, compared to a sequential CPU code, with near linear scaling on multi-node clusters. To our knowledge, this is the first GISAXS simulation code that is flexible to compute scattered light intensities in all spatial directions allowing full reconstruction of GISAXS patterns for any complex structures and with highresolutions while reducing simulation times from months to minutes.

Sarje, Abhinav↗

Linear-Scaling Coupled Cluster with Perturbative Triple Excitations: The Divide–Expand–Consolidate CCSD(T) Model

Here, we propose a reformulation of the traditional (T) triples correction to the coupled cluster singles and doubles (CCSD) energy in terms of local Hartree–Fock (HF) orbitals such that its structural form aligns with our recently developed linear-scaling divide–expand–consolidate (DEC) coupled cluster family of local correlation methods. In a DEC-CCSD(T) calculation, a basis of local occupied and virtual HF orbitals is used to partition the correlated calculation on the full system into a number of independent atomic fragment and pair fragment calculations, each performed within a truncated set of the complete orbital space. In return, this leads to a massively parallel algorithm for the evaluation of the DEC-CCSD(T) correlation energy, which formally scales linearly with the size of the full system and has a tunable precision with respect to a conventional CCSD(T) calculation via a single energy-based input threshold. The theoretical developments are supported by proof of concept DEC-CCSD(T) calculations on a series of medium-sized molecular systems.

Algorithms↗

A quantum processor based on coherent transport of entangled atom arrays

The ability to engineer parallel, programmable operations between desired qubits within a quantum processor is key for building scalable quantum information systems. In most state-of-the-art approaches, qubits interact locally, constrained by the connectivity associated with their fixed spatial layout. Here we demonstrate a quantum processor with dynamic, non-local connectivity, in which entangled qubits are coherently transported in a highly parallel manner across two spatial dimensions, between layers of single- and two-qubit operations. Our approach makes use of neutral atom arrays trapped and transported by optical tweezers; hyperfine states are used for robust quantum information storage, and excitation into Rydberg states is used for entanglement generation. We use this architecture to realize programmable generation of entangled graph states, such as cluster states and a seven-qubit Steane code state. Furthermore, we shuttle entangled ancilla arrays to realize a surface code state with thirteen data and six ancillary qubits and a toric code state on a torus with sixteen data and eight ancillary qubits. Finally, we use this architecture to realize a hybrid analogue–digital evolution and use it for measuring entanglement entropy in quantum simulations, experimentally observing non-monotonic entanglement dynamics associated with quantum many-body scars. Realizing a long-standing goal, these results provide a route towards scalable quantum processing and enable applications ranging from simulation to metrology.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

An Investigation of Parallel Programming Techniques Applied to Monte Carlo Simulations for Post-Flight Reconstruction of Spacecraft Trajectory

Parallelizing software to execute on multi-core central processing units (CPUs) and graphics processing units (GPUs) can be challenging. For some fields outside of Computer Science, this transition comes with new issues. For example, memory limitations can require modifications to code not initially developed to run on GPUs. This work applies the Open Multi-Processing (OpenMP) and Open Accelerators (OpenACC) directive-based parallelization strategies on a Monte Carlo simulation approach for trajectory reconstruction enabling it to run on multi-core CPUs and GPUs. Large matrix operations are the most common use of GPUs, which are not present in this algorithm; however, the natural parallelism of independent trajectories in Monte Carlo simulations is exploited. Benchmarking data are presented comparing execution times of the software for single-thread CPUs, multi-thread CPUs with OpenMP, and multi-thread GPUs using OpenACC. These data were collected using nodes with Intel® Xeon® E5-2670 (Sandy Bridge) CPUs enhanced with NVIDIA® Tesla® K40 GPUs on the Pleiades Supercomputer cluster at the National Aeronautics and Space Administration (NASA) Ames Research Center (ARC) and a local Intel® Xeon Phi™ node at NASA Langley Research Center (LaRC).

Williams, R. Anthony↗

Water–Hydrocarbon Interactions in Anionic Pyrene Monohydrate

Interactions between water and polycyclic aromatic hydrocarbons are essential in many aspects of chemistry, from interstellar and atmospheric processes to interfacial hydrophobicity and wetting phenomena. Despite their growing importance, the intermolecular potentials of the water-hydrocarbon interactions are underdeveloped compared to water-water potentials, and there are similarly few experimental probes that are sensitive to the details of the water-hydrocarbon potential. We present a combined experimental and computational study of anionic pyrene monohydrate, one of the simplest water/hydrocarbon clusters. The action spectrum in the OH region of the mass-selected cluster ion provides a rigorous benchmark for intermolecular potentials and computational methodologies. We identify missing intermolecular interactions and shortcomings in conventional dynamics calculations by comparing experimental data to density functional theory and classical molecular dynamics calculations. Kinetic trapping is prevalent, even for one water molecule and one pyrene molecule, leading to slow equilibration in conventional molecular dynamics calculations, even on nanosecond timescales and at low temperatures (50 K). At constant energy, temperature fluctuations for the pair of molecules are substantial. Immersing the system in a bath of soft spheres and employing parallel tempering alleviates kinetic trapping and dampens temperature fluctuations, bringing the system closer to the thermodynamic limit. With such augmented sampling, a simple, flexible water model reproduces the linewidth and the asymmetric broadening of the symmetric OH stretching mode, which we assign to spectral diffusion. In the OH stretching region, dynamics calculations predict a more intense antisymmetric peak than experiments observe but do not predict the bimodal split symmetric peak that the experiments show. Furthermore, our work suggests that electronic polarization, missing in the empirical force field, is responsible for the first discrepancy and that quantum nuclear effects, captured neither in density functional theory nor in classical dynamics, may be responsible for the second.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A k-mer based approach for classifying viruses without taxonomy identifies viral associations in human autism and plant microbiomes

Viruses are an underrepresented taxa in the study and identification of microbiome constituents; however, they play an essential role in health, microbiome regulation, and transfer of genetic material. Only a few thousand viruses have been isolated, sequenced, and assigned a taxonomy, which limits the ability to identify and quantify viruses in the microbiome. Additionally, the vast diversity of viruses represents a challenge for classification, not only in constructing a viral taxonomy, but also in identifying similarities between a virus’ genotype and its phenotype. However, the diversity of viral sequences can be leveraged to classify their sequences in metagenomic and metatranscriptomic samples, even if they do not have a taxonomy. To identify and quantify viruses in transcriptomic and genomic samples, we developed a dynamic programming algorithm for creating a classification tree out of 715,672 metagenome viruses. To create the classification tree, we clustered proportional similarity scores generated from the k-mer profiles of each of the metagenome viruses to create a database of metagenomic viruses. The resulting Kraken2 database of the metagenomic viruses can be found here: https://www.osti.gov/biblio/1615774 and is compatible with Kraken2. We then integrated the viral classification database with databases created with genomes from NCBI for use with ParaKraken (a parallelized version of Kraken provided in Supplemental Zip 1), a metagenomic/transcriptomic classifier. To illustrate the breadth of our utility for classifying metagenome viruses, we analyzed data from a plant metagenome study identifying genotypic and compartment specific differences between two Populus genotypes in three different compartments. We also identified a significant increase in abundance of eight viral sequences in post mortem brains in a human metatranscriptome study comparing Autism Spectrum Disorder patients and controls. We also show the potential accuracy for classifying viruses by utilizing both the JGI and NCBI viral databases to identify the uniqueness of viral sequences. Finally, we validate the accuracy of viral classification with NCBI databases containing viruses with taxonomy to identify pathogenic viruses in known COVID-19 and cassava brown streak virus infection samples. Our method represents the compulsory first step in better understanding the role of viruses in the microbiome by allowing for a more complete identification of sequences without taxonomy. Better classification of viruses will improve identifying associations between viruses and their hosts as well as viruses and other microbiome members. Despite the lack of taxonomy, this database of metagenomic viruses can be used with any tool that utilizes a taxonomy, such as Kraken, for accurate classification of viruses.

59 BASIC BIOLOGICAL SCIENCES↗

Parallel Evolution of Chemokine Binding by Structurally Related Herpesvirus Decoy Receptors

A wide variety of pathogens targets chemokine signaling networks in order to disrupt host immune surveillance and defense. Here, we report a structural and mutational analysis of rodent herpesvirus Peru encoded R17, a potent chemokine inhibitor that sequesters CC and C chemokines with high affinity. R17 consists of a pair of β-sandwich domains linked together by a bridging sheet, which form an acidic binding cleft for the chemokine CCL3 on the opposite face of a basic surface cluster that binds glycosaminoglycans. R17 promiscuously engages chemokines primarily through the same N-loop determinants used for host receptor recognition while residues located in the chemokine 40s loop drive kinetically stable complex formation. Here, the core fold adopted by R17 is unexpectedly similar to that of the M3 chemokine decoy receptor encoded by MHV-68, although, strikingly, neither the location of ligand engagement nor the stoichiometry of binding is conserved, suggesting that their functions evolved independently.

59 BASIC BIOLOGICAL SCIENCES↗

Necessary and Sufficient Conditions for Attitude Estimation in Fractionated Spacecraft Systems

This paper addresses the problem of attitude estimation in fractionated spacecraft clusters. Each module in the cluster may have either a star-tracker, a relative attitude sensor, or both. Using results in nonlinear ob- servability theory, we provide graph-theoretic sufficient conditions for the attitude of every module to be observable. In particular we show that the attitude of every module in the cluster can be observed if every module has either a star tracker with non-collinear stars, or there is a path through the sensing network from a module with a star tracker to the module without a star tracker, and each of the relative measurements along the path has either multiple non-collinear beacons or a single beacon that is not parallel to the rotation vector of the target module.

Blackmore, Lars↗

Exact quantum scattering calculations of transport properties for the H{sub 2}O–H system

Transport properties for collisions of water with hydrogen atoms are computed by means of exact quantum scattering calculations. For this purpose, a potential energy surface (PES) was computed for the interaction of rigid H{sub 2}O, frozen at its equilibrium geometry, with a hydrogen atom, using a coupled-cluster method that includes all singles and doubles excitations, as well as perturbative contributions of connected triple excitations. To investigate the importance of the anisotropy of the PES on transport properties, calculations were performed with the full potential and with the spherical average of the PES. We also explored the determination of the spherical average of the PES from radial cuts in six directions parallel and perpendicular to the C{sub 2} axis of the molecule. Finally, the computed transport properties were compared with those computed with a Lennard-Jones 12-6 potential.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Interactions of model biomolecules. Benchmark CC calculations within MOLCAS

We present results using the OVOS approach (Optimized Virtual Orbitals Space) aimed at enhancing the effectiveness of the Coupled Cluster calculations. This approach allows to reduce the total computer time required for large-scale CCSD(T) calculations about ten times when the original full virtual space is reduced to about 50% of its original size without affecting the accuracy. The method is implemented in the MOLCAS computer program. When combined with the Cholesky decomposition of the two-electron integrals and suitable parallelization it allows calculations which were formerly prohibitively too demanding. We focused ourselves to accurate calculations of the hydrogen bonded and the stacking interactions of the model biomolecules. Interaction energies of the formaldehyde, formamide, benzene, and uracil dimers and the three-body contributions in the cytosine – guanine tetramer are presented. Other applications, as the electron affinity of the uracil affected by solvation are also shortly mentioned.

74 ATOMIC AND MOLECULAR PHYSICS↗

Catalog of narrow Mg II absorption lines in the Baryon Oscillation Spectroscopic Survey

Using the Data Release 9 Quasar spectra from the Baryonic Oscillation Spectroscopic Survey, which does not include quasar spectra from the Sloan Digital Sky Survey Data Release 7, we detect narrow Mg ii λλ2796, 2803 absorption doublets in the spectral data redward of 1250 Å (quasar rest frame) until the red wing of the Mg ii λ2800 emission line. Our survey is limited to quasar spectra with a median signal-to-noise ratio 〈S/N〉 ⩾4 pixel{sup −1} in the surveyed spectral region, resulting in a sample that contains 43,260 quasars. We have detected a total of 18,598 Mg ii absorption doublets with 0.2933 ≤ z{sub abs} ≤ 2.6529. About 75% of absorbers have an equivalent width at rest frame of W{sub r}{sup λ2796}⩾1 A-ring . About 75% of absorbers have doublet ratios (DR=W{sub r}{sup λ2796}/W{sub r}{sup λ2803}) in the range of 1 ≤ DR ≤ 2, and about 3.2% lie outside the range of 1 − σ{sub DR} ≤ DR ≤ 2 + σ{sub DR}. We characterize the detection false positives/negatives by the frequency of detected Mg ii absorption doublets in the limits of the S/N of the spectral data. The S/N = 4.5 limit is assigned a completeness fraction of 53% and tends to be complete when the S/N is greater than 4.5. The redshift number densities of all of the detected Mg ii absorbers moderately increase from z ≈ 0.4 to z ≈ 1.5, which parallels the evolution of the cosmic star formation rate density. Limiting our investigation to those quasars whose emission redshift can be determined from narrow emission lines, the relative velocities (β) of Mg ii absorbers have a complex distribution which probably consists of three classes of Mg ii absorbers: (1) cosmologically intervening absorbers; (2) environmental absorbers that reside within the quasar host galaxies or galaxy clusters; (3) quasar outflow absorbers. After subtracting contributions from cosmologically intervening absorbers and environmental absorbers, the β distribution of the Mg iiabsorbers might mainly be contributed by the quasar outflow absorbers and peaks at υ ≈ 1500 km s{sup −1}. This peak velocity is lower than the value of 2000 km s{sup −1} found in statistical analysis of C iv absorbers.

79 ASTRONOMY AND ASTROPHYSICS↗