Search NASA⌕ Search

SEARCH · Search NASA

Results for “Cluster Computing”

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 379 records · Page 21

Post-processing interstitialcy diffusion from molecular dynamics simulations

An algorithm to rigorously trace the interstitialcy diffusion trajectory in crystals is developed. The algorithm incorporates unsupervised learning and graph optimization which obviate the need to input extra domain specific information depending on crystal or temperature of the simulation. The algorithm is implemented in a flexible framework as a post-processor to molecular dynamics (MD) simulations. We describe in detail the reduction of interstitialcy diffusion into known computational problems of unsupervised clustering and graph optimization. We also discuss the steps, computational efficiency and key components of the algorithm. Using the algorithm, thermal interstitialcy diffusion from low to near-melting point temperatures is studied. We encapsulate the algorithms in a modular framework with functionality to calculate diffusion coefficients, migration energies and other trajectory properties. The study validates the algorithm by establishing the conformity of output parameters with experimental values and provides detailed insights for the interstitialcy diffusion mechanism. The algorithm along with the help of supporting visualizations and analysis gives convincing details and a new approach to quantifying diffusion jumps, jump-lengths, time between jumps and to identify interstitials from lattice atoms. -- Graphical abstract:.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Technologies for node-degree based clustering of data sets

Technologies for node-degree based clustering include a computing device to construct a graph that includes multiple vertices corresponding to the data points of a data set. The computing device inserts an edge between each pair of vertices that has a corresponding similarity metric that meets a predetermined threshold similarity metric. The computing device determines a node degree for each vertex in the graph and initializes a cutoff node degree as the lowest node degree of the vertices. The computing device selects a test subset of the graph that includes vertices having a node degree less than or equal to the cutoff node degree. The computing device determines whether the test subset covers the graph and if not increases the cutoff node degree. If the test subset covers the graph, the data points corresponding to the vertices of the test subset are the representative cluster. Other embodiments are described and claimed.

Sitik, Ahmet Can↗

Automated Classification of Vehicle Movements at Signalized Intersections Using Vehicle Trajectories

Accurate vehicle movement classification through signalized intersections is of paramount importance to the analysis of intersection performance and the optimization of traffic control strategies. Conventional techniques for tracking vehicle turning movements depend on infrastructure-based strategies like human counts, loop detectors, and video analytics, all of which are costly, prone to errors, and spatially constrained. High-frequency trajectory data can be utilized to determine vehicle movement patterns in a scalable and infrastructure-independent method due to the adoption of connected vehicles (CVs). In recent years, several studies have utilized connected vehicle data to generate performance measures. Most of the trajectory-based performance measures approaches, however, require map matching-i.e., extracting geospatial references from maps to identify the movements that individual vehicles make at a signalized intersection. These approaches are often time-consuming and hinder scalability since geographic features need to be provided for an analysis to be conducted. Map matching methods are prone to errors as different map versions change these geographic features. This research presents a novel automatic classification pipeline that uses CV trajectory data to classify vehicle movements at signalized crossings, specifically pass-through left-turn and right-turn maneuvers. The process starts by filtering trips that cross a spatial bounding box that has been defined at the target intersection. Approach and departure headings for each trajectory crossing the boundary are computed and are clustered together to identify dominant movements. The proposed algorithm is used to classify the movement of vehicles at 10 intersections in the state of California, and the results indicate that the algorithm can classify movements at these intersections with varying traffic volumes and road network configurations, all in a map-less framework with no need for conflation of vehicle trajectories to a digital base map.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Reduced‐cost sparsity‐exploiting algorithm for solving coupled‐cluster equations

We present an algorithm for reducing the computational work involved in coupled‐cluster (CC) calculations by sparsifying the amplitude correction within a CC amplitude update procedure. We provide a theoretical justification for this approach, which is based on the convergence theory of inexact Newton iterations. We demonstrate by numerical examples that, in the simplest case of the CCD equations, we can sparsify the amplitude correction by setting, on average, roughly 90% nonzero elements to zeros without a major effect on the convergence of the inexact Newton iterations.

Brabec, Jiri↗

An ab initio exciton model for singlet fission

Here we present an ab initio exciton model that extends the Frenkel exciton model and includes valence, charge-transfer, and multiexcitonic excited states. It serves as a general, parameter-free, yet computationally efficient and scalable approach for simulation of singlet fission processes in multichromophoric systems. A comparison with multiconfigurational methods confirms that our exciton model predicts consistent energies and couplings for the pentacene dimer and captures the correct physics. Calculations of larger pentacene clusters demonstrate the computational scalability of the exciton model and suggest that the mixing between local and charge-transfer excitations narrows the gap between singlet and multiexcitonic states. Local vibrations of pentacene molecules are found to facilitate singlet–multiexcitonic state-crossing and hence are important for understanding singlet fission. The exciton model developed in this work also sets the stage for further implementation of the nuclear gradients and nonadiabatic couplings needed for first principles nonadiabatic quantum molecular dynamics simulations of singlet fission.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Superfluid response of two-dimensional parahydrogen clusters in confinement

We study by computer simulations the effect of confinement on the superfluid properties of small two-dimensional (2D) parahydrogen clusters. For clusters of fewer than twenty molecules, the superfluid response in the low temperature limit is found to remain comparable in magnitude to that of free clusters, within a rather wide range of depth and size of the confining well. The resilience of the superfluid response is attributable to the “supersolid” character of these clusters. We investigate the possibility of establishing a bulk 2D superfluid “cluster crystal” phase of p-H{sub 2}, in which a global superfluid response would arise from tunnelling of molecules across adjacent unit cells. The computed energetics suggests that for clusters of about ten molecules, such a phase may be thermodynamically stable against the formation of the equilibrium insulating crystal, for values of the cluster crystal lattice constant possibly allowing tunnelling across adjacent unit cells.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Spin–spin interactions in defects in solids from mixed all-electron and pseudopotential first-principles calculations

Abstract Understanding the quantum dynamics of spin defects and their coherence properties requires an accurate modeling of spin-spin interaction in solids and molecules, for example by using spin Hamiltonians with parameters obtained from first principles calculations. We present a real-space approach based on density functional theory for the calculation of spin-Hamiltonian parameters, where only selected atoms are treated at the all-electron level, while the rest of the system is described with the pseudopotential approximation. Our approach permits calculations for systems containing more than 1000 atoms, as demonstrated for defects in diamond and silicon carbide. We show that only a small number of atoms surrounding the defect needs to be treated at the all-electron level, in order to obtain an overall all-electron accuracy for hyperfine and zero-field splitting tensors. We also present results for coherence times, computed with the cluster correlation expansion method, highlighting the importance of accurate spin-Hamiltonian parameters for quantitative predictions of spin dynamics.

36 MATERIALS SCIENCE↗

Recent Advances on Computational Modeling of Supported Single-Atom and Cluster Catalysts: Characterization, Catalyst–Support Interaction, and Active Site Heterogeneity

To satisfy the need for catalyst materials with high activity, selectivity, and stability for energy conversion, material design and discovery guided by theoretical insights are a necessity. In the past decades, the rise in theoretical investigations into the properties of catalyst materials, reaction mechanisms, and catalyst design principles has shed light on the catalysis field. Quantitative structure–activity relationships have been developed through incorporating spectroscopic simulations, electronic structure calculations, and reaction mechanistic studies. Here, in this review, we report the state-of-the-art computational approaches to catalyst materials characterization for supported single-atom and cluster catalysts utilizing spectroscopic simulations, i.e., XANES simulation, and material properties investigation via electronic-structure calculations. Furthermore, approaches regarding reaction mechanisms, focusing on active site heterogeneity, are also discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Electronic Structure Theory Calculations Using Modern Architectures: KNL vs Haswell

The time to solution and parallel efficiency of several commonly used electronic structure methods (Hartree–Fock, density functional theory, second order perturbation theory, resolution of the identity second order perturbation theory, coupled cluster) are evaluated on both the Intel Xeon Haswell and the Intel Xeon Phi Knights Landing (KNL) architectures. The Haswell completes the benchmark calculations with a faster time to solution than the KNL for all molecules and methods tested. While the Haswell exhibits an average speedup of at least 3.5 relative to the KNL for all nonthreaded computations, the KNL has a better parallel efficiency than the Haswell with increasing core counts. The architectures are further tested using a more computationally costly coupled cluster method on a transition state reaction. The Haswell appears to be the best choice to minimize the time to solution, though for very large systems and high levels of theory that require memory intensive processes the superior memory hierarchy and larger on node memory of the KNL can make it a better choice. These results are used to showcase aspects of novel architectures that will increase efficiency for quantum chemistry applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Computational Study of the Resonance Enhancement of Raman Signals of Ligands Adsorbed to CdSe Clusters through Photoexcitation of the Cluster

This paper describes density-functional-theory-based computations of resonance Raman (RR) spectra of ligand molecules adsorbed to the surface of a Cd 16 Se 13 cluster. Signals from asymmetric vibrational modes of ligand binding groups, such as the asymmetric O–C–O stretching modes of carboxylates, are enhanced relative to the symmetric vibrational modes when the excitation energy is on-resonance with the excitonic energy of the cluster. Certain ligand molecules have frontier orbitals with the correct energies and symmetries to mix with the orbitals of the CdSe cluster, and as a result, the wave functions of the electron and the hole delocalize from the cluster onto the ligand molecules; experimentally, this delocalization results in a bathochromic shift of the band edge excitonic absorption. Increased excitonic delocalization results in greater vibronic coupling between the exciton and the ligand vibrations and, on average, preferential enhancements in the RR signals of those vibrations. This work suggests that the use of exciton-delocalizing ligands to optimize electronic coupling between neighboring CdSe nanoparticles may, at the same time, enhance the rates of nonradiative exciton decay by coupling the exciton and ligand vibrational modes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Poplar: a phylogenomics pipeline

Motivation Generating phylogenomic trees from the genomic data is essential in understanding biological systems. Each step of this complex process has received extensive attention and has been significantly streamlined over the years. Given the public availability of data, obtaining genomes for a wide selection of species is straightforward. However, analyzing that data to generate a phylogenomic tree is a multistep process with legitimate scientific and technical challenges, often requiring a significant input from a domain-area scientist. Results We present Poplar, a new, streamlined computational pipeline, to address the computational logistical issues that arise when constructing the phylogenomic trees. It provides a framework that runs state-of-the-art software for essential steps in the phylogenomic pipeline, beginning from a genome with or without an annotation, and resulting in a species tree. Running Poplar requires no external databases. In the execution, it enables parallelism for execution for clusters and cloud computing. The trees generated by Poplar match closely with state-of-the-art published trees. The usage and performance of Poplar is far simpler and quicker than manually running a phylogenomic pipeline. Availability and implementation Freely available on GitHub at https://github.com/sandialabs/poplar. Implemented using Python and supported on Linux.

Koning, Elizabeth [Sandia National Laboratories (S↗

Pressure Effects on the Relaxation of an Excited Ethane Molecule in High-Pressure Bath Gases

Here, we use molecular dynamics to calculate the rotational and vibrational energy relaxation of C 2 H 6 in Ar, Kr, and Xe bath gases over a pressure range of 10 to 400 atm and at temperatures of 300 K and 800 K. The C 2 H 6 is instantaneously excited by 80 kcal/mol randomly distributed into both vibrational and rotational modes. The computed relaxation rates show little sensitivity to the identity of the noble gas in the bath. Vibrational relaxation rates show a non-linear pressure dependence at 300 K. At 800 K the reduced range of bath gas densities covered by the range of pressures do not yet show any non-linearity in the pressure dependence. Rotational relaxation is characterized with two relaxation rates. The slower rate is comparable to the vibrational relaxation rate. The faster rate has a linear pressure dependence at 300 K but an irregular, nonlinear pressure dependence at 800 K. To understand this, a model was developed based on approximating the periodic box used in the molecular dynamics simulations by an equal-volume collection of cubes where each cube is sized to allow only single occupancy by the noble gas or the molecule. Combinatorial statistics then leads to a pressure and temperature dependent analytic distribution of the bath gas species the molecule encounters in a collision. This distribution, the dissociation energy of molecule/bath gas complexes and bath gas clusters, and the computed energy release per collision combine to show that only at 300 K is the energy release sufficient to dissociate likely complexes and clusters. This suggests that persistent and pressure-dependent clusters and complexes at 800 K may be responsible for the non-linear pressure dependence of rotational relaxation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Detecting CAN Masquerade Attacks with Signal Clustering Similarity

Vehicular Controller Area Networks (CANs) are susceptible to cyber attacks of different levels of sophistication. Fabrication attacks are the easiest to administer—an adversary simply sends (extra) frames on a CAN—but also the easiest to detect because they disrupt frame frequency. To overcome time-based detection methods, adversaries must administer masquerade attacks by sending frames in lieu of (and therefore at the expected time of) benign frames but with malicious payloads. Research efforts have proven that CAN attacks, and masquerade attacks in particular, can affect vehicle functionality. Examples include causing unintended acceleration, deactivation of vehicle’s brakes, as well as steering the vehicle. We hypothesize that masquerade attacks modify the nuanced correlations of CAN signal time series and how they cluster together. Therefore, changes in cluster assignments should indicate anomalous behavior. We confirm this hypothesis by leveraging our previously developed capability for reverse engineering CAN signals (i.e., CAN-D [Controller Area Network Decoder]) and focus on advancing the state of the art for detecting masquerade attacks by analyzing time series extracted from raw CAN frames. Specifically, we demonstrate that masquerade attacks can be detected by computing time series clustering similarity using hierarchical clustering on the vehicle’s CAN signals (time series) and comparing the clustering similarity across CAN captures with and without attacks. We test our approach in a previously collected CAN dataset with masquerade attacks (i.e., the ROAD dataset) and develop a forensic tool as a proof of concept to demonstrate the potential of the proposed approach for detecting CAN masquerade attacks.

Moriano Salazar, Pablo↗

Limits on clustering and smooth quintessence from the EFTofLSS

We apply the Effective Field Theory of Large-Scale Structure (EFTofLSS) to analyze cosmological models with clustering quintessence, which allows us to consistently describe the parameter region in which the quintessence equation of state w < -1. First, we extend the description of biased tracers in redshift space to the presence of clustering quintessence, and compute the one-loop power spectrum. We solve the EFTofLSS equations using the exact time dependence, which is relevant to obtain unbiased constraints. Then, fitting the full shape of BOSS pre-reconstructed power spectrum measurements, the BOSS post-reconstruction BAO measurements, BAO measurements from 6DF/MGS and eBOSS, the Supernovae from Pantheon, and a prior from BBN, we bound the clustering quintessence equation of state parameter w = -1.011 -0.048 +0.053 at 68% C.L. Further combining with Planck, we obtain w = -1.028 -0.030 +0.037 at 68% C.L. We also obtain constraints on smooth quintessence, in the physical regime w ≥ -1: combining all datasets, we get -1 ≤ w < -0.979 at 68% C.L. These results strongly support a cosmological constant.

79 ASTRONOMY AND ASTROPHYSICS↗

The effect of hydrogen-bonding on flavin’s infrared absorption spectrum

Cluster and continuum solvation computational models are employed to model the effect of hydrogen bonding interactions on the vibrational modes of lumiflavin. Calculated spectra were compared to experimental Fourier-transform infrared (FTIR) spectra in the diagnostic 1450–1800 cm –1 range, where intense V C=C , V C=N , V C 2 =O , and V C 4 =O stretching modes of flavin’s isoalloxazine ring are found. Local mode analysis is used to describe the strength of hydrogen-bonding in cluster models. The computations indicate that V C=C and V C=N mode frequencies are relatively insensitive to intermolecular interactions while the V C 2 =O and V C 4 =O modes are sensitive to direct (and also indirect for V C 2 =O ) hydrogen-bonding interactions. Although flavin is neutral, basis sets without the diffuse functions provide incorrect relative frequencies and intensities. The 6-31+G* basis set is found to be adequate for this system, and there is limited benefit to considering larger basis sets. Calculated vibrational mode frequencies agree with experimentally determined frequencies in solution when cluster models with multiple water molecules are used. Accurate simulation of relative FTIR band intensities, on the other hand, requires a continuum (or possibly quantum mechanical/molecular mechanical) model that accounts for long-range electrostatic effects. Lastly, an experimental peak at ca. 1624 cm –1 that is typically assigned to the V C 2 =O vibrational stretching mode has a complicated shape that suggests multiple underlying contributions. Our calculations show that this band has contributions from both the C 6 –C 7 and C 2 = O stretching vibrations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The halo model in a massive neutrino cosmology

We provide a quantitative analysis of the halo model in the context of massive neutrino cosmologies. We discuss all the ingredients necessary to model the non-linear matter and cold dark matter power spectra and compare with the results of N-body simulations that incorporate massive neutrinos. Our neutrino halo model is able to capture the non-linear behavior of matter clustering with a ∼20% accuracy up to very non-linear scales of k = 10 h/Mpc (which would be affected by baryon physics). The largest discrepancies arise in the range k = 0.5 – 1 h/Mpc where the 1-halo and 2-halo terms are comparable and are present also in a massless neutrino cosmology. However, at scales k < 0.2 h/Mpc our neutrino halo model agrees with the results of N-body simulations at the level of 8% for total neutrino masses of < 0.3 eV. We also model the neutrino non-linear density field as a sum of a linear and clustered component and predict the neutrino power spectrum and the cold dark matter-neutrino cross-power spectrum up to k = 1 h/Mpc with ∼30% accuracy. For masses below 0.15 eV the neutrino halo model captures the neutrino induced suppression, casted in terms of matter power ratios between massive and massless scenarios, with a 2% agreement with the results of N-body/neutrino simulations. Finally, we provide a simple application of the halo model: the computation of the clustering of galaxies, in massless and massive neutrinos cosmologies, using a simple Halo Occupation Distribution scheme and our halo model extension.

79 ASTRONOMY AND ASTROPHYSICS↗

Spatiotemporal modeling of node temperatures in supercomputers

Los Alamos National Laboratory (LANL) is home to many large supercomputing clusters. These clusters require an enormous amount of power (~500-2000 kW each), and most of this energy is converted into heat. Thus, cooling the components of the supercomputer becomes a critical and expensive endeavor. Recently a project was initiated to investigate the effect that changes to the cooling system in a machine room had on three large machines that were housed there. Coupled with this goal was the aim to develop a general good-practice for characterizing the effect of cooling changes and monitoring machine node temperatures in this and other machine rooms. This paper focuses on the statistical approach used to quantify the effect that several cooling changes to the room had on the temperatures of the individual nodes of the computers. The largest cluster in the room has 1,600 nodes that run a variety of jobs during general use. Since extremes temperatures are important, a Normal distribution plus generalized Pareto distribution for the upper tail is used to model the marginal distribution, along with a Gaussian process copula to account for spatio-temporal dependence. A Gaussian Markov random field (GMRF) model is used to model the spatial effects on the node temperatures as the cooling changes take place. This model is then used to assess the condition of the node temperatures after each change to the room. The analysis approach was used to uncover the cause of a problematic episode of overheating nodes on one of the supercomputing clusters. Lastly, this same approach can easily be applied to monitor and investigate cooling systems at other data centers, as well.

97 MATHEMATICS AND COMPUTING↗

Properties of the water to boron nitride interaction: From zero to two dimensions with benchmark accuracy

Molecular adsorption on surfaces plays an important part in catalysis, corrosion, desalination, and various other processes that are relevant to industry and in nature. As a complement to experiments, accurate adsorption energies can be obtained using various sophisticated electronic structure methods that can now be applied to periodic systems. The adsorption energy of water on boron nitride substrates, going from zero to 2-dimensional periodicity, is particularly interesting as it calls for an accurate treatment of polarizable electrostatics and dispersion interactions, as well as posing a practical challenge to experiments and electronic structure methods. Here, we present reference adsorption energies, static polarizabilities, and dynamic polarizabilities, for water on BN substrates of varying size and dimension. Adsorption energies are computed with coupled cluster theory, fixed-node quantum Monte Carlo (FNQMC), the random phase approximation, and second order Møller-Plesset theory. These wavefunction based correlated methods are found to agree in molecular as well as periodic systems. The best estimate of the water/h-BN adsorption energy is -107 ± 7 meV from FNQMC. In addition, the water adsorption energy on the BN substrates could be expected to grow monotonically with the size of the substrate due to increased dispersion interactions, but interestingly, this is not the case here. This peculiar finding is explained using the static polarizabilities and molecular dispersion coefficients of the systems, as computed from time-dependent density functional theory (DFT). Dynamic as well as static polarizabilities are found to be highly anisotropic in these systems. Finally, the many-body dispersion method in DFT emerges as a particularly useful estimation of finite size effects for other expensive, many-body wavefunction based methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗