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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Computer simulation structure and vibrations of small metal cluster on the Cu (111) surface

Vibrational properties of the small tetrahedral cluster of Co on the Cu (111) surface are studied by using tight-binding second moment approximation interatomic interaction potentials. It was shown that interaction of the clusters with substrate leads to arising of frustrated translation and frustrated rotation in-plane polarized vibrational modes localized on the cluster atoms. The Co{sub 4} cluster on the surface the high frequency modes remain strongly localized and mixed with the nearest neighbor atoms vibrations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

An experimental and computational view of the photoionization of diol–water clusters

In the interstellar medium, diols and other prebiotic molecules adsorb onto icy mantles surrounding dust grains. Water in the ice may affect the reactivity and photoionization of these diols. Ethylene glycol (EG), 1,2-propylene glycol, and 1,3-propylene glycol clusters with water clusters were used as a proxy to study these interactions. The diol–water clusters were generated in a continuous supersonic molecular beam, photoionized by synchrotron-based vacuum ultraviolet light from the Advanced Light Source, and subsequently detected by reflectron time-of-flight mass spectrometry. The appearance energies for the detected clusters were determined from the mass spectra, collected at increasing photon energy. Clusters of both diol fragments and unfragmented diols with water were detected. The lowest energy geometry optimized conformers for the observed EG–water clusters and EG fragment–water clusters have been visualized using density functional theory (DFT), providing insight into hydrogen bonding networks and how these affect fragmentation and appearance energy. As the number of water molecules clustered around EG fragments (m/z 31 and 32) increased, the appearance energy for the cluster decreased, indicating a stabilization by water. This trend was supported by DFT calculations. Fragment clusters from 1,2-propylene glycol exhibited a similar trend, but with a smaller energy decrease, and no trend was observed from 1,3-propylene glycol. In conclusion, we discuss and suggest that the reactivity and photoionization of diols in the presence of water depend on the size of the diol, the location of the hydroxyl group, and the number of waters clustered around the diol.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Spatially Resolved Mapping of Three-Dimensional Molecular Orientations with ~2 nm Spatial Resolution through Tip-Enhanced Raman Scattering

We record local optical field images of silver nanocubes (75 nm) using tip-enhanced Raman (TER) spectral imaging. The images that we observe are consistent with several recent reports from our group, but here, we demonstrate sub-2 nm spatial resolution in local optical field nanoimaging under ambient laboratory conditions. This is achieved by scanning the substrate (nanocube on Si) relative to a 4-thiobenzonitrile (TBN)-functionalized Ag-coated TER probe. The spatial resolution that we obtain necessitates that only a few molecules govern the recorded optical response; molecular orientation becomes an important consideration in such measurements. We model the orientation through geometry optimization of a TBN molecule chemisorbed onto an Ag79 cluster (sphere with a ~1 nm diameter). Using the computed orientation of the cluster-bound molecule, we then model the optical response using formalism that accounts for the orientation of the molecule relative to vector components of the local optical fields. We find optimal agreement between experiment and theory. In effect, this work reveals the parallels between single-molecule Raman scattering and high-spatial-resolution TER spectroscopy, even when the images themselves cannot be used to visualize a single molecule in real space.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

Sodium Superionic Conductors Based on Clusters

Owing to the high abundance and low cost of sodium (Na), Na-based rechargeable batteries hold great potential for large-scale applications in the future energy industry. However, as key component of the battery electrolyte, only a few Na-based superionic conductors can reach the ionic conductivity comparable to that of liquid or gel electrolytes. Here, we provide a guideline for the development of cluster-based Na-rich antiperovskite superionic conductors using computational studies. With a selected cluster ion BCl 4 – , we are able to achieve high-room-temperature Na-ionic conductivity over 10 –3 S/cm and low activation energies below 0.2 eV in the antiperovskite crystals Na 3 S(BCl 4 ) and Na 3 S(BCl 4 ) 0.5 I 0.5 . In addition, these materials have large bandgaps and favorable mechanical properties. Here, a comprehensive study of the stability and formation energy of these materials further illustrates possible routes for their synthesis. New insights into the conduction mechanism of these cluster-based superionic conductors are provided, including the cooperative motion of Na ions and the significant reduction of the migration barrier due to the changing orientation of the cluster ion.

25 ENERGY STORAGE↗
DOE OSTI2023

Proper orthogonal descriptors for efficient and accurate interatomic potentials

Here, we present the proper orthogonal descriptors for efficient and accuracy representation of the potential energy surface. The potential energy surface is represented as a many-body expansion of parametrized potentials in which the potentials are functions of atom positions and parameters. The proper orthogonal decomposition is employed to decompose the parametrized potentials into a set of proper orthogonal descriptors (PODs). Because of the rapid convergence of the proper orthogonal decomposition, relevant snapshots can be sampled exhaustively to represent the atomic neighborhood environment accurately with a small number of descriptors. The proper orthogonal descriptors are used to develop interatomic potentials by using a linear expansion of the descriptors and determining the expansion coefficients from a weighted least-squares regression against a density functional theory (DFT) training set. We present a comprehensive evaluation of the POD potentials on previously published DFT data sets comprising Li, Mo, Cu, Ni, Si, Ge, and Ta elements. The data sets represent a diverse pool of metals, transition metals, and semiconductors. The accuracy of the POD potentials are comparable to that of state-of-the-art machine learning potentials such as the spectral neighbor analysis potential (SNAP) and the atomic cluster expansion (ACE).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

A hybrid computational strategy to address WGS variant analysis in >5000 samples

Abstract Background The decreasing costs of sequencing are driving the need for cost effective and real time variant calling of whole genome sequencing data. The scale of these projects are far beyond the capacity of typical computing resources available with most research labs. Other infrastructures like the cloud AWS environment and supercomputers also have limitations due to which large scale joint variant calling becomes infeasible, and infrastructure specific variant calling strategies either fail to scale up to large datasets or abandon joint calling strategies. Results We present a high throughput framework including multiple variant callers for single nucleotide variant (SNV) calling, which leverages hybrid computing infrastructure consisting of cloud AWS, supercomputers and local high performance computing infrastructures. We present a novel binning approach for large scale joint variant calling and imputation which can scale up to over 10,000 samples while producing SNV callsets with high sensitivity and specificity. As a proof of principle, we present results of analysis on Cohorts for Heart And Aging Research in Genomic Epidemiology (CHARGE) WGS freeze 3 dataset in which joint calling, imputation and phasing of over 5300 whole genome samples was produced in under 6 weeks using four state-of-the-art callers. The callers used were SNPTools, GATK-HaplotypeCaller, GATK-UnifiedGenotyper and GotCloud. We used Amazon AWS, a 4000-core in-house cluster at Baylor College of Medicine, IBM power PC Blue BioU at Rice and Rhea at Oak Ridge National Laboratory (ORNL) for the computation. AWS was used for joint calling of 180 TB of BAM files, and ORNL and Rice supercomputers were used for the imputation and phasing step. All other steps were carried out on the local compute cluster. The entire operation used 5.2 million core hours and only transferred a total of 6 TB of data across the platforms. Conclusions Even with increasing sizes of whole genome datasets, ensemble joint calling of SNVs for low coverage data can be accomplished in a scalable, cost effective and fast manner by using heterogeneous computing platforms without compromising on the quality of variants.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Connected three-body terms in single-reference unitary many-body theories: Iterative and perturbative approximations

This work introduces various approaches to include connected three-body terms in unitary many-body theories, focusing on the driven similarity renormalization group (DSRG). Starting from the least approximate method—the linearized DSRG truncated to one-, two-, and three-body operators [LDSRG(3)]—we develop several approximate LDSRG(3) models with reduced computational cost. Through a perturbative analysis, we motivate a family of iterative LDSRG(3)-n and -n' (n = 1, 2, 3, 4) methods that contain a subset of the LDSRG(3) diagrams. Among these variants, the LDSRG(3)-2 scheme has the same computational complexity of coupled cluster theory with singles, doubles, and triples (CCSDT), but it outperforms CCSDT in the accuracy of the predicted correlation energies. We also propose and implement two perturbative triples corrections based on the linearized DSRG truncated to one- and two-body operators augmented with recursive semi-quadratic commutators [qDSRG(2)]. Overall, the resulting qDSRG(2)+(T) approach matches the accuracy of the “gold-standard” coupled cluster theory with singles, doubles, and perturbative triples model on the energetics of twenty-eight closed-shell atoms and small molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

Simulation of circularly polarized luminescence spectra using coupled cluster theory

We report the first computations of circularly polarized luminescence (CPL) rotatory strengths at the equation-of-motion coupled cluster singles and doubles (EOM-CCSD) level of theory. Using a test set of eight chiral ketones, we compare both dipole and rotatory strengths for absorption (electronic circular dichroism) and emission to the results from time-dependent density-functional theory (TD-DFT) and available experimental data for both valence and Rydberg transitions. For two of the compounds, we obtained optimized geometries of the lowest several excited states using both EOM-CCSD and TD-DFT and determined that structures and EOM-CCSD transition properties obtained with each structure were sufficiently similar that TD-DFT optimizations were acceptable for the remaining test cases. Agreement between EOM-CCSD and the Becke three-parameter exchange function and Lee-Yang-Parr correlation functional (B3LYP) corrected using the Coulomb attenuating method (CAM-B3LYP) is typically good for most of the transitions, though agreement with the uncorrected B3LYP functional is significantly worse for all reported properties. The choice of length vs. velocity representation of the electric dipole operator has little impact on the EOM-CCSD transition strengths for nearly all of the states we examined. For a pair of closely related β, γ-enones, (1R)-7-methylenebicyclo[2.2.1]heptan-2-one and (1S)-2-methylenebicyclo[2.2.1]heptan-7-one, we find that EOM-CCSD and CAM-B3LYP agree with the energetic ordering of the two possible excited-state conformations, resulting in good agreement with experimental rotatory strengths in both absorption and emission, whereas B3LYP yields a qualitatively incorrect result for the CPL signal of (1S)-2-methylenebicyclo[2.2.1]heptan-7-one. Finally, we predict that one of the compounds considered here, trans-bicyclo[3.3.0]octane-3,7-dione, is unique in that it exhibits an achiral ground state and a chiral first excited state, leading to a strong CPL signal but a weak circular dichroism signal.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

Machine learning based job status prediction in scientific clusters

Large high-performance computing systems are built with increasing number of components with more CPU cores, more memory, and more storage space. At the same time, scientific applications have been growing in complexity. Together, they are leading to more frequent unsuccessful job statuses on HPC systems. From measured job statuses, 23.4% of CPU time was spent to the unsuccessful jobs. Here, we set out to study whether these unsuccessful job statuses could be anticipated from known job characteristics. To explore this possibility, we have developed a job status prediction method for the execution of jobs on scientific clusters. The Random Forests algorithm was applied to extract and characterize the patterns of unsuccessful job statuses. Experimental results show that our method can predict the unsuccessful job statuses from the monitored ongoing job executions in 99.8% the cases with 83.6% recall and 94.8% precision. Lastly, this prediction accuracy can be sufficiently high that it can be used to mitigation procedures of predicted failures.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Cluster Morphology-Polymer Dynamics Correlations in Sulfonated Polystyrene Melts: Computational Study

Reaching exceptionally long times up to 500 ns in equilibrium and nonequilibrium molecular dynamics simulations studies, we have attained a fundamental molecular understanding of the correlation of ionomer clusters structure and multiscale dynamics, providing new insight into one critical, long-standing challenge in ionic polymer physics. The cluster structure in melts of sulfonated polystyrene with Na + and Mg 2+ counterions are resolved and correlated with the dynamics on multiple length and time scales extracted from measurements of the dynamic structure factor and shear rheology. We find that as the morphology of the ionic clusters changes from ladderlike for Na + to disordered structures for Mg 2+ , the dynamic structure factor is affected on the length scale corresponding to the ionic clusters. Lastly, rheology studies show that the viscosity for Mg 2+ melts is higher than for Na + ones for all shear rates, which is well correlated with the larger ionic clusters’ size for the Mg 2+ melts.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2016

Emergent properties of nuclei from ab initio coupled-cluster calculations

Emergent properties such as nuclear saturation and deformation, and the effects on shell structure due to the proximity of the scattering continuum and particle decay channels are fascinating phenomena in atomic nuclei. In recent years, ab initio approaches to nuclei have taken the first steps towards tackling the computational challenge of describing these phenomena from Hamiltonians with microscopic degrees of freedom. Our endeavor is now possible due to ideas from effective field theories, novel optimization strategies for nuclear interactions, ab initio methods exhibiting a soft scaling with mass number, and ever-increasing computational power. We review some of the recent accomplishments. We also present new results. The recently optimized chiral interaction NNLO${}_{{\rm{sat}}}$ is shown to provide an accurate description of both charge radii and binding energies in selected light- and medium-mass nuclei up to 56Ni. We derive an efficient scheme for including continuum effects in coupled-cluster computations of nuclei based on chiral nucleon–nucleon and three-nucleon forces, and present new results for unbound states in the neutron-rich isotopes of oxygen and calcium. Finally, the coupling to the continuum impacts the energies of the ${J}^{\pi }=1/{2}^{-},3/{2}^{-},7/{2}^{-},3/{2}^{+}$ states in ${}^{\mathrm{17,23,25}}$O, and—contrary to naive shell-model expectations—the level ordering of the ${J}^{\pi }=3/{2}^{+},5/{2}^{+},9/{2}^{+}$ states in ${}^{\mathrm{53,55,61}}$Ca.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2015

THE POWER OF IMAGING: CONSTRAINING THE PLASMA PROPERTIES OF GRMHD SIMULATIONS USING EHT OBSERVATIONS OF Sgr A*

Recent advances in general relativistic magnetohydrodynamic simulations have expanded and improved our understanding of the dynamics of black-hole accretion disks. However, current simulations do not capture the thermodynamics of electrons in the low density accreting plasma. This poses a significant challenge in predicting accretion flow images and spectra from first principles. Because of this, simplified emission models have often been used, with widely different configurations (e.g., disk- versus jet-dominated emission), and were able to account for the observed spectral properties of accreting black holes. Exploring the large parameter space introduced by such models, however, requires significant computational power that exceeds conventional computational facilities. In this paper, we use GRay, a fast graphics processing unit (GPU) based ray-tracing algorithm, on the GPU cluster El Gato, to compute images and spectra for a set of six general relativistic magnetohydrodynamic simulations with different magnetic field configurations and black-hole spins. We also employ two different parametric models for the plasma thermodynamics in each of the simulations. We show that, if only the spectral properties of Sgr A* are used, all 12 models tested here can fit the spectra equally well. However, when combined with the measurement of the image size of the emission using the Event Horizon Telescope, current observations rule out all models with strong funnel emission, because the funnels are typically very extended. Our study shows that images of accretion flows with horizon-scale resolution offer a powerful tool in understanding accretion flows around black holes and their thermodynamic properties.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

FORMATION AND PROPERTIES OF ASTROPHYSICAL CARBONACEOUS DUST. I. AB-INITIO CALCULATIONS OF THE CONFIGURATION AND BINDING ENERGIES OF SMALL CARBON CLUSTERS

The binding energies of n < 100 carbon clusters are calculated using the ab initio density functional theory code Quantum Espresso. Carbon cluster geometries are determined using several levels of classical techniques and further refined using density functional theory. The resulting energies are used to compute the work of cluster formation and the nucleation rate in a saturated, hydrogen-poor carbon gas. Compared to classical calculations that adopt the capillary approximation, we find that nucleation of carbon clusters is enhanced at low temperatures and depressed at high temperatures. This difference is ascribed to the different behavior of the critical cluster size. We find that the critical cluster size is at n = 27 or n = 8 for a broad range of temperatures and saturations, instead of being a smooth function of such parameters. The results of our calculations can be used to follow carbonaceous cluster/grain formation, stability, and growth in hydrogen-poor environments, such as the inner layers of core-collapse supernovae and supernova remnants.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

The HEPCloud Facility: elastic computing for High Energy Physics – The NOvA Use Case

The need for computing in the HEP community follows cycles of peaks and valleys mainly driven by conference dates, accelerator shutdown, holiday schedules, and other factors. Because of this, the classical method of provisioning these resources at providing facilities has drawbacks such as potential overprovisioning. As the appetite for computing increases, however, so does the need to maximize cost efficiency by developing a model for dynamically provisioning resources only when needed. To address this issue, the HEPCloud project was launched by the Fermilab Scientific Computing Division in June 2015. Its goal is to develop a facility that provides a common interface to a variety of resources, including local clusters, grids, high performance computers, and community and commercial Clouds. Initially targeted experiments include CMS and NOvA, as well as other Fermilab stakeholders. In its first phase, the project has demonstrated the use of the “elastic” provisioning model offered by commercial clouds, such as Amazon Web Services. In this model, resources are rented and provisioned automatically over the Internet upon request. In January 2016, the project demonstrated the ability to increase the total amount of global CMS resources by 58,000 cores from 150,000 cores - a 25 percent increase - in preparation for the Recontres de Moriond. In March 2016, the NOvA experiment has also demonstrated resource burst capabilities with an additional 7,300 cores, achieving a scale almost four times as large as the local allocated resources and utilizing the local AWS s3 storage to optimize data handling operations and costs. NOvA was using the same familiar services used for local computations, such as data handling and job submission, in preparation for the Neutrino 2016 conference. In both cases, the cost was contained by the use of the Amazon Spot Instance Market and the Decision Engine, a HEPCloud component that aims at minimizing cost and job interruption. This paper describes the Fermilab HEPCloud Facility and the challenges overcome for the CMS and NOvA communities.

Fuess, S.↗
DOE OSTI2022

Critical points of the random cluster model with Newman–Ziff sampling

Here, we present a method for computing transition points of the random cluster model using a generalization of the Newman–Ziff algorithm, a celebrated technique in numerical percolation, to the random cluster model. The new method is straightforward to implement and works for real cluster weight q > 0. Furthermore, results for an arbitrary number of values of q can be found at once within a single simulation. Because the algorithm used to sweep through bond configurations is identical to that of Newman and Ziff, which was conceived for percolation, the method loses accuracy for large lattices when q > 1. However, by sampling the critical polynomial, accurate estimates of critical points in two dimensions can be found using relatively small lattice sizes, which we demonstrate here by computing critical points for non-integer values of q on the square lattice, to compare with the exact solution, and on the unsolved non-planar square matching lattice. The latter results would be much more difficult to obtain using other techniques.

97 MATHEMATICS AND COMPUTING↗