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 505 records · Page 28

Simplifying computational workflows with the Multiscale Atomic Zeolite Simulation Environment (MAZE)

Zeolites, an important class of 3-dimensional nanoporous materials, have been widely explored for a variety of applications including gas storage, separations, and catalysis. As the properties of these aluminosilicate materials depend on a number of factors (e.g., framework topology, Si/Al ratio, extra-framework cations etc.), detailed experiments (e.g., catalytic properties, adsorption capacities etc.) are often limited to only a handful of materials. Computational methods have played an important role in (1) providing molecular level insights to rationalize experimental observations, and (2) screening large libraries of zeolites to identify promising candidates for experimental synthesis and validation. Different levels of theory and computational chemistry codes are necessary to describe the range of relevant phenomena such as adsorption (e.g., grand canonical Monte Carlo), diffusion (e.g., molecular dynamics), and chemical reactions (e.g., density functional theory). Manipulation of atomic structures, handling of input files, and developing robust workflows becomes quite cumbersome. To mitigate these challenges, we describe the development of the Multiscale Atomic Zeolite Simulation Environment (MAZE) – a Python package that simplifies zeolite-specific calculation workflows by providing a user-friendly interface for systematically manipulating zeolite structures

97 MATHEMATICS AND COMPUTING↗

DONKEY: A Flexible and Accurate Algorithm for Clustering

We propose an accurate clustering algorithm suitable for the varied and multidimensional data sets that correspond to temporal snapshots from on-the-fly nonadiabatic trajectory-based simulations of photoexcited dynamics. The algorithm approximates the underlying probability density function using variable kernel density estimation, with local maxima corresponding to cluster centers. Each data point is then assigned to one of the maxima by employing a maximization procedure. Finally, clusters artificially separated by minor fluctuations in the probability density are merged. The algorithm does not require parameter tuning, which ensures flexibility and reduces the risk of bias. It is tested on several synthetic data sets, where it consistently outperforms conventional clustering algorithms. As a final example, the algorithm is applied to the excited dynamics of the norbornadiene ⇌ quadricyclane (C 7 H 8 ) molecular photoswitch, demonstrating how distinct reaction pathways can be identified.

algorithms↗

LeaPP: Learning Pathways to Polymorphs through Machine Learning Analysis of Atomic Trajectories

Understanding the mechanisms underlying crystal nucleation and growth is crucial for many technological applications. Due to the short length and time scales involved, crystal nucleation is often studied using molecular simulations. Most existing approaches to extract the nucleation mechanism from simulations focus on the analysis of static snapshots of the configurations, potentially overlooking subtle local fluctuations and the history of the particles involved in the formation of solid nuclei. Here, in this work, we propose a novel methodology called LeaPP that categorizes nucleation trajectories based on the temporal information of their constituent particles. We leverage the time evolution of the local environment of the crystallizing particles to encapsulate the relationship between the structure and dynamics and distinguish between different evolving particle paths. Identification of the distinct particle paths further enables characterizing the nucleation trajectories into different pathways. Collectively, LeaPP provides a more nuanced understanding of nucleation through an unsupervised approach with lesser dependence on traditional order parameters. Furthermore, the pathways identified by LeaPP are predictive of the resulting polymorph. We demonstrate LeaPP on three different systems─Lennard-Jones-like particles, Ni 3 Al, and water on surfaces. The general methodology underlying LeaPP─considering the time evolution of the building blocks─applies to a wide range of self-assembly problems.

36 MATERIALS SCIENCE↗

Determining Atomistic SAXS Models of Tri-Ubiquitin Chains from Bayesian Analysis of Accelerated Molecular Dynamics Simulations

Small-angle X-ray scattering (SAXS) has become an increasingly popular technique for characterizing the solution ensemble of flexible biomolecules. However, data resulting from SAXS is typically low-dimensional and is therefore difficult to interpret without additional structural knowledge. In theory, molecular dynamics (MD) trajectories can provide this information, but conventional simulations rarely sample the complete ensemble. Here, we demonstrate that accelerated MD simulations can be used to produce higher quality models in shorter time scales than standard simulations, and we present an iterative Bayesian Monte Carlo method that is able to identify multistate ensembles without overfitting. This methodology is applied to several ubiquitin trimers to demonstrate the effect of linkage type on the solution states of the signaling protein. Here, we observe that the linkage site directly affects the solution flexibility of the trimer and theorize that this difference in plasticity contributes to their disparate roles in vivo.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

First-Principles Simulations of CuCl in High-Temperature Water Vapor

Experimental data suggest that the solubility of copper in high-temperature water vapor is controlled by the formation of hydrated clusters of the form CuCl(H 2 O) n , where the average number of water molecules in the cluster generally increases with increasing density [Migdisov, A. A.; et al. Geochim. Cosmochim. Acta2014, 129, 33–53]. However, the precise nature of these clusters is difficult to probe experimentally. Moreover, there are some discrepancies between experimental estimates of average cluster size and prior simulation work [Mei, Y. Geofluids2018, 2018, 4279124]. Here, we have performed first-principles Monte Carlo (MC) and molecular dynamics (MD) simulations to explore these clusters in finer detail. We find that molecular dynamics is not the most appropriate technique for studying aggregation in vapor phases, even at relatively high temperatures. Specifically, our MD simulations exhibit substantial problems in adequately sampling the equilibrium cluster size distribution. In contrast, MC simulations with specialized cluster moves are able to accurately sample the phase space of hydrogen-bonding vapors. At all densities, we find a stable, slightly distorted linear H 2 O–Cu–Cl structure, which is in agreement with the earlier simulations, surrounded by a variable number of water molecules. The surrounding water molecules do not form a well-defined second solvation shell but rather a loose network of hydrogen-bonded water with molecular CuCl on the outside edge of the water cluster. We also find a broad distribution of hydration numbers, especially at higher densities. In contrast to previous simulation work but in agreement with experimental data, we find that the average hydration number substantially increases with increasing density. Moreover, the value of the hydration number depends on the choice of cluster definition.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Crystal Nucleation Kinetics and Mechanism: Influence of Interaction Potential

Modulating liquid-to-solid transitions and the resulting crystalline structure for tailored properties is much desired. Colloidal systems are exemplary to this end, but the fundamental knowledge gaps in relating the influence of intermolecular interactions to crystallization behavior continue to hinder progress. In this study, we address this knowledge gap by studying nucleation and growth in systems with modified Lennard-Jones potential. Specifically, we study the commonly used 12-6 potential and a softer 7-6 potential. The thermodynamic state point for the study is chosen such that both systems are investigated at the same level of supercooling and pressure. Under these conditions, we find that the nucleation rate for both systems is comparable. Interestingly, the nucleation pathways and resulting crystal structures are different. In the 12-6 system, nucleation and growth occur predominantly through the FCC structure. Softening the potential alters the critical nucleus composition and introduces two distinct nucleation pathways. One pathway predominantly leads to the nucleus with a body-centered cubic (BCC) structure, while the other favors the face-centered cubic (FCC) arrangement. Our study illustrates that polymorph selection can be achieved through modifications to intermolecular interactions without impacting nucleation kinetics. The results have significant implications in designing approaches for polymorph selection and modulating self-assembly mechanisms.

36 MATERIALS SCIENCE↗

Reagent-Triggered Isomerization of Fluxional Cluster Catalyst via Dynamic Coupling

Metallic cluster catalysts have many thermodynamically accessible isomers with diverse active sites and low reaction barriers, and lately a strong hypothesis emerged that the many catalyst states collectively drive the catalysis. Furthermore, it remained a hypothesis that catalyst isomerization is actually kinetically feasible under the current reaction conditions. Using high-temperature dynamics simulations and sampling, a range of orientations, and vibrational energy distributions, we probe how thermal effects and molecular events affect cluster catalyst dynamics. We show that even such a delicate affair as the dissociation or scattering of a methane molecule on the heavy and thus slow Pt 13 cluster triggers substantial isomerization of the catalyst, far beyond thermal at 700 K. A kinetic coupling between the methane activity and cluster catalyst dynamics is observed. In return, the thermal dynamics of the cluster affects the methane reaction and scattering probabilities. Hence, molecular events at the surfaces of fluxional cluster catalysts should facilitate the population of an ensemble of catalyst states under the current reaction conditions, with implications for available active sites, reaction mechanisms, and apparent rates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Formation of Exotic Networks of Water Clusters in Helium Droplets Facilitated by the Presence of Neon Atoms

Water clusters are formed in helium droplets via the sequential capture of monomers. One or two neon atoms are added to each droplet prior to the addition of water. The infrared spectrum of the droplet ensemble reveals several signatures of polar, water tetramer clusters having dipole moments between 2D and 3D. Comparison with ab initio computations supports the assignment of the cluster networks to noncyclic “3+1” clusters, which are ~5.3 kcal/mol less stable than the global minimum nonpolar cyclic tetramer. The (H2O)3Ne + H2O ring insertion barrier is sufficiently large, such that evaporative helium cooling is capable of kinetically quenching the nonequilibrium tetramer system prior to its rearrangement to the lower energy cyclic species. To this end, the reported process results in the formation of exotic water cluster networks that are either higher in energy than the most stable gas-phase analogs or not even stable in the gas phase.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Deducing subnanometer cluster size and shape distributions of heterogeneous supported catalysts

Abstract Infrared (IR) spectra of adsorbate vibrational modes are sensitive to adsorbate/metal interactions, accurate, and easily obtainable in-situ or operando. While they are the gold standards for characterizing single-crystals and large nanoparticles, analogous spectra for highly dispersed heterogeneous catalysts consisting of single-atoms and ultra-small clusters are lacking. Here, we combine data-based approaches with physics-driven surrogate models to generate synthetic IR spectra from first-principles. We bypass the vast combinatorial space of clusters by determining viable, low-energy structures using machine-learned Hamiltonians, genetic algorithm optimization, and grand canonical Monte Carlo calculations. We obtain first-principles vibrations on this tractable ensemble and generate single-cluster primary spectra analogous to pure component gas-phase IR spectra. With such spectra as standards, we predict cluster size distributions from computational and experimental data, demonstrated in the case of CO adsorption on Pd/CeO 2 (111) catalysts, and quantify uncertainty using Bayesian Inference. We discuss extensions for characterizing complex materials towards closing the materials gap.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Water confinement in small polycyclic aromatic hydrocarbons

The confinement of water molecules is vital in fields from biology to nanotechnology. The conditions allowing confinement in small finite polycyclic aromatic hydrocarbons (PAHs) are unclear, yet are crucial for understanding confinement in larger systems. Here, we report a computational study of water cluster confinement within PAHs dimers. Our results serve as a model for larger carbon allotropes and for understanding molecular interactions in confined systems. We identified size and structural motifs allowing confinement and demonstrated the motifs in various PAHs systems. Here we show that optimal OH∙∙∙π interactions between water clusters and the PAH dimer permit optimal confinement to occur. However, the lack of such interactions leads to the formation of CH∙∙∙O interactions, resulting in less ideal confinement. Confinement of layered clusters is also possible, provided that the optimal OH∙∙∙π interactions are conserved.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Laboratory study of isocyanic acid ions: Rotational spectroscopy of NCO{sup −}, H{sub 2}NCO{sup +}, and HNCOH{sup +}

We report detection of protonated isocyanic acid in two isomeric forms, H{sub 2}NCO{sup +} and HNCOH{sup +}, by high-resolution spectroscopy. The two ions were first observed at centimeter wavelengths by Fourier Transform (FT) microwave spectroscopy, in a discharge through HNCO heavily diluted in hydrogen in the throat of a supersonic nozzle. Spectroscopic constants derived from the two lowest rotational transitions of both isomers agree very well with those derived from theoretical structures computed at the coupled cluster level of theory. In the same molecular beam, the fundamental rotational transition of NCO{sup −} was observed with well-resolved nitrogen quadrupole hyperfine structure. Detection of NCO{sup −} and H{sub 2}NCO{sup +} in our beam was subsequently confirmed by observation of several millimeter-wave transitions in a low pressure discharge through cyanogen and water. The spectroscopic constants of NCO{sup −} obtained earlier by infrared laser spectroscopy are in good agreement with the highly accurate constants derived here. Owing to the high abundance of HNCO in many galactic molecular sources, both ions are excellent candidates for astronomical detection in the radio band.

74 ATOMIC AND MOLECULAR PHYSICS↗

Merging ab initio theory and few-body approach for (d, p) reactions

A framework for A(d, p)B reactions is introduced by merging the microscopic approach to computing the properties of the nucleon-target systems and the three-body n + p + A reaction formalism, thus providing a consistent link between the reaction cross sections and the underlying microscopic structure. In this first step toward a full microscopic description, we focus on the inclusion of the neutron-target microscopic properties. The properties of the neutron-target subsystem are encapsulated in the Green's function which is computed with the coupled Cluster theory using a chiral nucleon-nucleon and three-nucleon interactions. Subsequently, this many-body information is introduced in the few-body Green's Function Transfer approach to (d, p) reactions. Our benchmarks on stable targets 40,48 Ca show an excellent agreement with the data. We then proceed to make specific predictions for (d, p) on neutron rich 52,54 Ca isotopes. These predictions are directly relevant to testing the new magic numbers N = 32, 34 and are expected to be feasible in the first campaign of the projected FRIB facility.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Collisional excitation of CH 2 rotational/fine-structure levels by helium

Accurate determination of the abundance of CH 2 in interstellar media relies on both radiative and collisional rate coefficients. Here, we investigate here the rotational/fine-structure excitation of CH 2 induced by collisions with He. We employ a recoupling technique to generate fine-structure-resolved cross-sections and rate coefficients from close coupling spin-free scattering calculations. The calculations are based on a recent, high-accuracy CH 2 –He potential energy surface computed at the coupled clusters level of theory. The collisional cross-section calculations are performed for all fine-structure transitions among the first 22 and 24 energy levels of ortho- and para-CH 2 , respectively, and for temperatures up to 300 K. As a first application, we simulate the excitation of CH 2 in typical molecular clouds. The excitation temperatures of the CH 2 lines are found to be small at typical densities of molecular clouds, showing that the non-local thermodynamic equilibrium approach has to be used to analyse interstellar spectra. We also found that the fine-structure lines connected with the 404 – 313 and 505 – 414 rotational transitions show possible maser emissions so that they can be easily seen in emission. In conclusion, these calculations show that CH 2 may have to be detected mainly through absorption spectra.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Coherent elastic neutrino-nucleus scattering on Ar 40 from first principles

Coherent elastic neutrino scattering on the 40 Ar nucleus is computed with coupled-cluster theory based on nuclear Hamiltonians inspired by effective field theories of quantum chromodynamics. Our approach is validated by calculating the charge form factor and comparing it to data from electron scattering. We make predictions for the weak form factor, the neutron radius, and the neutron skin and estimate systematic uncertainties. The neutron-skin thickness of 40 Ar is consistent with results from density-functional theory. Precision measurements from coherent elastic neutrino-nucleus scattering could potentially be used to extract these observables and help to constrain nuclear models.

39 ≤ A ≤ 58↗

Exactness of the normal-ordered two-body truncation of three-nucleon forces

Reference-state-based many-body methods start from Hamiltonians that are normal ordered with respect to the reference state. In low-energy nuclear physics applications, normal-ordered Hamiltonians consisting of two- and three-nucleon forces are usually truncated at the two-body rank with residual three-nucleon operators being discarded. Benchmark computations have shown that this truncation is accurate, but we lack an understanding about why it works. Here, we show that the normal-ordered two-body truncation is exact for zero-range three-body forces when nuclei are computed using the coupled cluster with singles and doubles method. As the nuclear three-nucleon force is short ranged and a three-body contact is a leading term in effective field theories of quantum chromodynamics, our result provides an analytical basis for the popular normal-ordered two-body approximation.

Rothman, Maxwell [Univ. of Tennessee, Knoxville, T↗

A Visual Analytics Paradigm Enabling Trillion-Edge Graph Exploration

We present a visual analytics paradigm and a system prototype for exploring web-scale graphs. A web-scale graph is described as a graph with ~one trillion edges and ~50 billion vertices. While there is an aggressive R&D effort in processing and exploring web-scale graphs among internet vendors such as Facebook and Google, visualizing a graph of that scale still remains an underexplored R&D area. The paper describes a nontraditional peek-and-filter strategy that facilitates the exploration of a graph database of unprecedented size for visualization and analytics. We demonstrate that our system prototype can 1) preprocess a graph with ~25 billion edges in less than two hours and 2) support database query and visualization on the processed graph database afterward. Based on our computational performance results, we argue that we most likely will achieve the one trillion edge mark (a computational performance improvement of 40 times) for graph visual analytics in the near future.

graph visualization↗