Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 685 records · Page 38

Structural transition in sputter-deposited amorphous germanium films by aging at ambient temperature

The structure of amorphous Ge (a-Ge) films prepared by sputter-deposition and the effects of aging at ambient temperature and pressure were studied by pair-distribution-function (PDF) analysis from electron scattering and molecular dynamics simulations. The PDFs of the as-deposited and aged samples for 3–13 months showed that the major peaks for Ge-Ge bonds decrease in intensity and broaden with aging for up to 7 months. In the PDFs of a-Ge of molecular dynamics simulation obtained by quenching liquid at different rates, the major peak intensities of a slowly cooled model are higher than those of a rapidly cooled model. Analyses on short- and medium-range configurations show that the slowly cooled model includes a certain amount of medium-range ordered (MRO) clusters, while the rapidly cooled model includes liquid-like configurations rather than MRO clusters. The similarity between experimental and computational PDFs implies that as-deposited films are similar in structure to the slowly cooled model, whereas the fully aged films are similar to the rapidly cooled model. It is assumed that as they undergo room-temperature aging, the MRO clusters disintegrate and transform into liquid-like regions in the same matrix. This transition in local configurations is discussed in terms of instability and the non-equilibrium of nanoclusters produced by a vapor-deposition process.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Excited-state downfolding using ground-state formalisms

Downfolding coupled cluster (CC) techniques are powerful tools for reducing the dimensionality of many-body quantum problems. This work investigates how ground-state downfolding formalisms can target excited states using non-Aufbau reference determinants, paving the way for applications of quantum computing in excited-state chemistry. This study focuses on doubly excited states for which canonical equation-of-motion CC approaches struggle to describe unless one includes higher-than-double excitations. The downfolding technique results in state-specific effective Hamiltonians that, when diagonalized in their respective active spaces, provide ground- and excited-state total energies (and therefore excitation energies) comparable to high-level CC methods. The performance of this procedure is examined with doubly excited states of H 2 , Methylene, Formaldehyde, and Nitroxyl.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Shadow Molecular Dynamics and Atomic Cluster Expansions for Flexible Charge Models

Here, a shadow molecular dynamics scheme for flexible charge models is presented where the shadow Born–Oppenheimer potential is derived from a coarse-grained approximation of range-separated density functional theory. The interatomic potential, including the atomic electronegativities and the charge-independent short-range part of the potential and force terms, is modeled by the linear atomic cluster expansion (ACE), which provides a computationally efficient alternative to many machine learning methods. The shadow molecular dynamics scheme is based on extended Lagrangian (XL) Born–Oppenheimer molecular dynamics (BOMD). XL-BOMD provides stable dynamics while avoiding the costly computational overhead associated with solving an all-to-all system of equations, which normally is required to determine the relaxed electronic ground state prior to each force evaluation. To demonstrate the proposed shadow molecular dynamics scheme for flexible charge models using atomic cluster expansion, we emulate the dynamics generated from self-consistent charge density functional tight-binding (SCC-DFTB) theory using a second-order charge equilibration (QEq) model. The charge-independent potentials and electronegativities of the QEq model are trained for a supercell of uranium oxide (UO 2 ) and a molecular system of liquid water. The combined ACE+XL-QEq molecular dynamics simulations are stable over a wide range of temperatures both for the oxide and for the molecular systems and provide a precise sampling of the Born–Oppenheimer potential energy surfaces. Accurate ground Coulomb energies are produced by the ACE-based electronegativity model during an NVE simulation of UO 2 , predicted to be within 1 meV of those from SCC-DFTB on average during comparable simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Comparative study of water reactivity with Mo{sub 2}O{sub y}{sup −} and W{sub 2}O{sub y}{sup −} clusters: A combined experimental and theoretical investigation

A computational investigation of the Mo{sub 2}O{sub y}{sup −} + H{sub 2}O (y = 4, 5) reactions as well as a photoelectron spectroscopic probe of the deuterated Mo{sub 2}O{sub 6}D{sub 2}{sup −} product have been carried out to understand a puzzling question from a previous study: Why is the rate constant determined for the Mo{sub 2}O{sub 5}{sup −} + H{sub 2}O/D{sub 2}O reaction, the terminal reaction in the sequential oxidation of Mo{sub 2}O{sub y}{sup −} by water, higher than the W{sub 2}O{sub 5}{sup −} + H{sub 2}O/D{sub 2}O reaction? This disparity was intriguing because W{sub 3}O{sub y}{sup −} clusters were found to be more reactive toward water than their Mo{sub 3}O{sub y}{sup −} analogs. A comparison of molecular structures reveals that the lowest energy structure of Mo{sub 2}O{sub 5}{sup −} provides a less hindered water addition site than the W{sub 2}O{sub 5}{sup −} ground state structure. Several modes of water addition to the most stable molecular and electronic structures of Mo{sub 2}O{sub 4}{sup −} and Mo{sub 2}O{sub 5}{sup −} were explored computationally. The various modes are discussed and compared with previous computational studies on W{sub 2}O{sub y}{sup −} + H{sub 2}O reactions. Calculated free energy reaction profiles show lower barriers for the initial Mo{sub 2}O{sub y}{sup −} + H{sub 2}O addition, consistent with the higher observed rate constant. The terminal Mo{sub 2}O{sub y}{sup −} sequential oxidation product predicted computationally was verified by the anion photoelectron spectrum of Mo{sub 2}O{sub 6}D{sub 2}{sup −}. Based on the computational results, this anion is a trapped dihydroxide intermediate in the Mo{sub 2}O{sub 5}{sup −} + H{sub 2}O/D{sub 2}O → Mo{sub 2}O{sub 6}{sup −} + H{sub 2}/D{sub 2} reaction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Shadow molecular dynamics and atomic cluster expansions for flexible charge models

A shadow molecular dynamics scheme for flexible charge models is presented, where the shadow Born-Oppenheimer potential is derived from a coarse-grained approximation of range-separated density functional theory. The interatomic potential, including the atomic electronegativities and the charge-independent short-range part of the potential and force terms, are modeled by the linear atomic cluster expansion (ACE), which provides a computationally efficient alternative to many machine learning methods. The shadow molecular dynamics scheme is based on extended Lagrangian (XL) Born-Oppenheimer molecular dynamics (BOMD) [Eur. Phys. J. B 94, 164 (2021)]. XL-BOMD provides a stable dynamics, while avoiding the costly computational overhead associated with solving an all-to-all system of equations, which normally is required to determine the relaxed electronic ground state prior to each force evaluation. To demonstrate the proposed shadow molecular dynamics scheme for flexible charge models using the atomic cluster expansion, we emulate the dynamics generated from self-consistent charge density functional tight-binding (SCC-DFTB) theory using a second-order charge equilibration (QEq) model. The charge-independent potentials and electronegativities of the QEq model are trained for a supercell of uranium oxide (UO2) and a molecular system of liquid water. The combined ACE + XL-QEq dynamics are stable over a wide range of temperatures both for the oxide and the molecular systems, and provide a precise sampling of the Born-Oppenheimer potential energy surfaces. Accurate ground Coulomb energies are produced by the ACE-based electronegativity model during an NVE simulation of UO 2 , predicted to be within 1 meV of those from SCC-DFTB on average during comparable simulations.

74 ATOMIC AND MOLECULAR PHYSICS↗

Understanding the Angular Momentum Evolution of T Tauri and Herbig Ae/Be Stars

We investigate a sample of six Herbig Ae/Be stars belonging to the Orion OB1 association, as well as 73 low-mass objects, members of the σ Orionis cluster, in order to explore the angular momentum evolution at early stages of evolution, and its possible connection with main-sequence Ap/Bp magnetic stars. Using FIES and HECTOCHELLE spectra, we obtain projected rotational velocities through two independent methods. Individual masses, radii, and ages are computed using evolutionary models, distance, and cluster extinction. Under the assumption that similar physical processes operate in both T Tauri and Herbig Ae/Be stars, we construct snapshots of the protostar’s rotation against mass during the first 10 Myr with the aid of a rotational model that includes a variable disk lifetime, changes in the stellar moment of inertia, a dipolar magnetic field with variable strength, and angular momentum loss through stellar winds powered by accretion. We use these snapshots, as well as the rotational data, to infer a plausible scenario for the angular momentum evolution. We find that magnetic field strengths of a few kilo-Gauss at 3 Myr are required to match the rotational velocities of both groups of stars. Models with masses between 2 M {sub ⊙} and 3 M {sub ⊙} display larger angular momentum values by a factor of ~3, in comparison to stars of similar spectral types on the main sequence. Even though some quantitative estimates on this dramatic decrease with age for Ap/Bp magnetic main-sequence stars are presented, the results obtained for the angular momentum evolution do not explain their low rotation rates.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Clusters, molecular layers, and 3D crystals of water on Ni(111)

We examined the growth and stability of ice layers on Ni(111) up to ∼7 molecular layers (ML) thick using scanning tunneling microscopy. At low coverage, films were comprised of ∼1 nm wide two-dimensional (2D) clusters. Only above ∼0.5 ML did patches of continuous 2D layers emerge, coexisting with the clusters until the first ML was complete. The structure of the continuous layer is clearly different from that of the 2D clusters. Subsequently, a second molecular layer grew on top of the first. 3D crystallites started to form only after this 2nd ML was complete. 2D clusters re-appeared when thicker films were partially evaporated, implying that these clusters represent the equilibrium configuration at low coverage. Binding energies and image simulations computed with density functional theory suggest that the 2D clusters are partially dissociated and surrounded by H adatoms. The complete 2D layer contains only intact water molecules because of the lack of favorable binding sites for H atoms. We propose molecular structures for the 2D layer that are composed of the same pentagon-heptagon binding motif and water density observed on Pt(111). The similarity of the water structures on Pt and Ni suggests a general prescription for generating low-energy configurations on close-packed metal substrates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A new scheme for perturbative triples correction to (0,1) sector of Fock space multi-reference coupled cluster method: Theory, implementation, and examples

We propose a new elegant strategy to implement third order triples correction in the light of many-body perturbation theory to the Fock space multi-reference coupled cluster method for the ionization problem. The computational scaling as well as the storage requirement is of key concerns in any many-body calculations. Our proposed approach scales as N{sup 6} does not require the storage of triples amplitudes and gives superior agreement over all the previous attempts made. This approach is capable of calculating multiple roots in a single calculation in contrast to the inclusion of perturbative triples in the equation of motion variant of the coupled cluster theory, where each root needs to be computed in a state-specific way and requires both the left and right state vectors together. The performance of the newly implemented scheme is tested by applying to methylene, boron nitride (B{sub 2}N) anion, nitrogen, water, carbon monoxide, acetylene, formaldehyde, and thymine monomer, a DNA base.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

TAMM: Tensor algebra for many-body methods

Tensor algebra operations such as contractions in computational chemistry consume a significant fraction of the computing time on large-scale computing platforms. The widespread use of tensor contractions between large multi-dimensional tensors in describing electronic structure theory has motivated the development of multiple tensor algebra frameworks targeting heterogeneous computing platforms. In this paper, we present Tensor Algebra for Many-body Methods (TAMM), a framework for productive and performance-portable development of scalable computational chemistry methods. TAMM decouples the specification of the computation from the execution of these operations on available high-performance computing systems. With this design choice, the scientific application developers (domain scientists) can focus on the algorithmic requirements using the tensor algebra interface provided by TAMM, whereas high-performance computing developers can direct their attention to various optimizations on the underlying constructs, such as efficient data distribution, optimized scheduling algorithms, and efficient use of intra-node resources (e.g., graphics processing units). The modular structure of TAMM allows it to support different hardware architectures and incorporate new algorithmic advances. We describe the TAMM framework and our approach to the sustainable development of scalable ground- and excited-state electronic structure methods. We present case studies highlighting the ease of use, including the performance and productivity gains compared to other frameworks.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effect of Coulomb Correlation on the Magnetic Properties of Mn Clusters

In spite of decades of research, a fundamental understanding of the unusual magnetic behavior of small Mn clusters remains a challenge. Experiments show that Mn 2 is antiferromagnetic while small clusters containing up to five Mn atoms are ferromagnetic with magnetic moments of 5 μB/atom and become ferrimagnetic as they grow further. Theoretical studies based on density functional theory (DFT), however, find Mn2 to be ferromagnetic, with ferrimagnetic order setting in at different sizes that depend upon the computational methods used. While quantum chemical techniques correctly account for the antiferromagnetic ground state of Mn 2 , they are computationally too demanding to treat larger clusters, making it difficult to understand the evolution of magnetism. These studies clearly point to the importance of correlation and the need to find ways to treat it effectively for larger clusters and nanostructures. In this report we show that the DFT+U method can be used to account for strong correlation. We determine the on-site Coulomb correlation, Hubbard U self-consistently by using the linear response theory and study its effect on the magnetic coupling of Mn clusters containing up to five atoms. With a calculated U value of 4.8 eV, we show that the ground state of Mn 2 is antiferromagnetic with a Mn–Mn distance of 3.34 Å, which agrees well with the electron spin resonance experiment. Equally important, we show that on-site Coulomb correlation also plays an important role in the evolution of magnetic coupling in larger clusters, as the results differ significantly from standard DFT calculations. We conclude that for a proper understanding of magnetism of Mn nanostructures (clusters, chains, and layers) one must take into account the effect of strong correlation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Nuclear Schiff Moments and CP Violation

This article reviews the calculation of nuclear Schiff moments, which one must know in order to interpret experiments that search for time-reversal-violating electric dipole moments in certain atoms and molecules. After briefly reviewing the connection between dipole moments and CP violation in and beyond the Standard Model of particle physics; Schiff's theorem, which concerns the screening of nuclear electric dipole moments by electrons; Schiff moments; and experiments to measure dipole moments in atoms and molecules, this review examines attempts to compute Schiff moments in nuclei such as 199 Hg and octupole-deformed isotopes such as 225 Ra, which are particularly useful in experiments. It then turns to ab initio nuclear-structure theory, describing ways in which both the in-medium similarity renormalization group and coupled-cluster theory can be used to compute important Schiff moments more accurately than the less controlled methods that have been applied so far.

CP violation↗

Monte Carlo neutronics analysis of the SNS STS monolith - 25684

Neutronics analysis for the Oak Ridge National Laboratory's Spallation Neutron Source second target station monolith conceptual design was performed entirely with the Monte Carlo method. It has been demonstrated that even though the monolith provides ∼peta-scale attenuation of the dose rates, Monte Carlo simulations are within reach of computing capabilities of relatively inexpensive Linux clusters. Several monolith configurations achieved or exceeded design requirements and provide basis for the final design. (authors)

43 PARTICLE ACCELERATORS↗

Vibrational energies for HFCO using a neural network sum of exponentials potential energy surface

A six-dimensional potential energy surface (PES) for formyl fluoride (HFCO) is fit in a sum-of-products form using neural network exponential fitting functions. The ab initio data upon which the fit is based were computed at the explicitly correlated coupled cluster with single, double, and perturbative triple excitations [CCSD(T)-F12]/cc-pVTZ-F12 level of theory. The PES fit is accurate (RMSE = 10 cm{sup −1}) up to 10 000 cm{sup −1} above the zero point energy and covers most of the experimentally measured IR data. The PES is validated by computing vibrational energies for both HFCO and deuterated formyl fluoride (DFCO) using block improved relaxation with the multi-configuration time dependent Hartree approach. The frequencies of the fundamental modes, and all other vibrational states up to 5000 cm{sup −1} above the zero-point energy, are more accurate than those obtained from the previous MP2-based PES. The vibrational frequencies obtained on the PES are compared to anharmonic frequencies at the MP2/aug-cc-pVTZ and CCSD(T)/aug-cc-pVTZ levels of theory obtained using second-order vibrational perturbation theory. The new PES will be useful for quantum dynamics simulations for both HFCO and DFCO, e.g., studies of intramolecular vibrational redistribution leading to unimolecular dissociation and its laser control.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Aerodynamic optimization of supersonic compressor cascade using differential evolution on GPU

Differential Evolution (DE) is a powerful stochastic optimization method. Compared to gradient-based algorithms, DE is able to avoid local minima but requires at the same time more function evaluations. In turbomachinery applications, function evaluations are performed with time-consuming CFD simulation, which results in a long, non affordable, design cycle. Modern High Performance Computing systems, especially Graphic Processing Units (GPUs), are able to alleviate this inconvenience by accelerating the design evaluation itself. In this work we present a validated CFD Solver running on GPUs, able to accelerate the design evaluation and thus the entire design process. An achieved speedup of 20x to 30x enabled the DE algorithm to run on a high-end computer instead of a costly large cluster. The GPU-enhanced DE was used to optimize the aerodynamics of a supersonic compressor cascade, achieving an aerodynamic loss minimization of 20%.

97 MATHEMATICS AND COMPUTING↗

Execute BEE workflows on private cloud infrastructure-2.3.6.01 - LANL ATDM ST / STNS01-22 Milestone Completion Documentation (BEE-FY21 P6-2) [Slides]

This work involves the creation of the Cloud Launcher, a new subcomponent of BEE, and the extension of the BEETaskManager to run on Cloud systems. BEE will be able to interact with the Google Compute Engine and OpenStack cloud APIs to set up simple Cloud clusters for launching HPC job scripts. BEE will use existing functionality to launch jobs that previously could only be launched on HPC systems. The BEETaskManager will handle launching tasks on the Cloud cluster.

97 MATHEMATICS AND COMPUTING↗

ASCR Workshop on In Situ Data Management: Enabling Scientific Discovery from Diverse Data Sources

In January 2019, the U.S. Department of Energy, Office of Science program in Advanced Scientific Computing Research, convened a workshop to identify priority research directions for in situ data management (ISDM). The workshop defined ISDM as the practices, capabilities, and procedures to control the organization of data and enable the coordination and communication among heterogeneous tasks, executing simultaneously in a high-performance computing system, cooperating toward a common objective. The workshop revealed two primary, interdependent motivations for processing and managing data in situ. The first motivation is that the in situ methodology enables scientific discovery from a broad range of data sources over a wide scale of computing platforms: leadership-class systems, clusters, clouds, workstations, and embedded devices at the edge. The successful development of ISDM capabilities will benefit real-time decision-making, design optimization, and data-driven scientific discovery. The second motivation is the need to decrease data volumes. ISDM can make critical contributions to managing large data volumes from computations and experiments to minimize data movement, save storage space, and boost resource efficiency, often while simultaneously increasing scientific precision. A fundamental finding of this workshop is that the methodologies used to manage data among a variety of tasks in situ can be used to facilitate scientific discovery from many different data sources—simulation, experiment, and sensors, for example—and that being able to do so at numerous computing scales will benefit real-time decision-making, design optimization, and data-driven scientific discovery across the Office of Science mission space. Applications wanting to use the in situ capabilities include those where data analysis feeds back to the simulation, decisions are made autonomously, big data or machine learning is among the tasks to be coordinated, and computations need to be completed in real time. The workshop identified six priority research directions that highlight the components and capabilities needed for ISDM to be successful for the wide variety of applications discussed: making ISDM capabilities more pervasive, controllable, composable, and transparent, with a focus on greater coordination with the software stack and a diversity of fundamentally new data algorithms.

97 MATHEMATICS AND COMPUTING↗

ASCR Workshop on In Situ Data Management

In January 2019, the U.S. Department of Energy, Office of Science program in Advanced Scientific Computing Research, convened a workshop to identify priority research directions for in situ data management (ISDM). The workshop defined ISDM as the practices, capabilities, and procedures to control the organization of data and enable the coordination and communication among heterogeneous tasks, executing simultaneously in a high-performance computing system, cooperating toward a common objective. The workshop revealed two primary, interdependent motivations for processing and managing data in situ. The first motivation is that the in situ methodology enables scientific discovery from a broad range of data sources over a wide scale of computing platforms: leadership-class systems, clusters, clouds, workstations, and embedded devices at the edge. The successful development of ISDM capabilities will benefit real-time decision-making, design optimization, and data-driven scientific discovery. The second motivation is the need to decrease data volumes. ISDM can make critical contributions to managing large data volumes from computations and experiments to minimize data movement, save storage space, and boost resource efficiency, often while simultaneously increasing scientific precision.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗

Hydration dynamics in water clusters via quantum molecular dynamics simulations

We have investigated the hydration dynamics in size selected water clusters with n = 66, 104, 200, 500, and 1000 water molecules using molecular dynamics simulations. To study the most fundamental aspects of relaxation phenomena in clusters, we choose one of the simplest, still realistic, quantum mechanically treated test solute, an excess electron. The project focuses on the time evolution of the clusters following two processes, electron attachment to neutral equilibrated water clusters and electron detachment from an equilibrated water cluster anion. The relaxation dynamics is significantly different in the two processes, most notably restoring the equilibrium final state is less effective after electron attachment. Nevertheless, in both scenarios only minor cluster size dependence is observed. Significantly different relaxation patterns characterize electron detachment for interior and surface state clusters, interior state clusters relaxing significantly faster. This observation may indicate a potential way to distinguish surface state and interior state water cluster anion isomers experimentally. A comparison of equilibrium and non-equilibrium trajectories suggests that linear response theory breaks down for electron attachment at 200 K, but the results converge to reasonable agreement at higher temperatures. Relaxation following electron detachment clearly belongs to the linear regime. Cluster relaxation was also investigated using two different computational models, one preferring cavity type interior states for the excess electron in bulk water, while the other simulating non-cavity structure. While the cavity model predicts appearance of several different hydrated electron isomers in agreement with experiment, the non-cavity model locates only cluster anions with interior excess electron distribution. The present simulations show that surface isomers computed with the cavity predicting potential show similar dynamical behavior to the interior clusters of the non-cavity type model. Relaxation associated with cavity collapse presents, however, unique dynamical signatures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗