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 1,351 records · Page 75

Magnetic field induced phase transitions of the triangular spin model

In this work, motivated by a recent inelastic neutron scattering experiment on YbMgGaO 4 , I study the homogeneous spin model on the triangular lattice using the cluster mean-field theory and exact diagonalization. Phase diagrams with and without magnetic field are computed. The results show that the magnetic field can induce the phase transition from the spin liquid state (or the stripe state) to the 120 ° antiferromagnetic state. These phase transitions are suppressed by the next-nearest-neighbor exchange interaction J 2 / J 1 and vanish as J 2 / J 1 > 0.11 . I analyze a parameter space at J 2 / J 1 = 0.1 , where a field induced phase transition can occur, and find that the deviation of the theoretical spin excitation energy from the experimental data is only about 0.054. These results imply that an effective homogeneous spin model still applies to YbMgGaO 4 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

LACEwING: A New Moving Group Analysis Code

We present a new nearby young moving group (NYMG) kinematic membership analysis code, LocAting Constituent mEmbers In Nearby Groups (LACEwING), a new Catalog of Suspected Nearby Young Stars, a new list of bona fide members of moving groups, and a kinematic traceback code. LACEwING is a convergence-style algorithm with carefully vetted membership statistics based on a large numerical simulation of the Solar Neighborhood. Given spatial and kinematic information on stars, LACEwING calculates membership probabilities in 13 NYMGs and three open clusters within 100 pc. In addition to describing the inputs, methods, and products of the code, we provide comparisons of LACEwING to other popular kinematic moving group membership identification codes. As a proof of concept, we use LACEwING to reconsider the membership of 930 stellar systems in the Solar Neighborhood (within 100 pc) that have reported measurable lithium equivalent widths. We quantify the evidence in support of a population of young stars not attached to any NYMGs, which is a possible sign of new as-yet-undiscovered groups or of a field population of young stars.

79 ASTRONOMY AND ASTROPHYSICS↗

On the geometric dependence of the molecular dipole polarizability in water: A benchmark study of higher-order electron correlation, basis set incompleteness error, core electron effects, and zero-point vibrational contributions

In this work, we investigate how geometric changes influence the static dipole polarizability ($α$) of a water molecule by explicitly computing the corresponding dipole polarizability surface (DPS) across 3125 total (1625 symmetry-unique) geometries using linear response coupled cluster theory including single, double, and triple excitations (LR-CCSDT) and the doubly augmented triple-$ζ$ basis set (d-aug-cc-pVTZ). Analytical formulae based on power series expansions of this $ab initio$ surface are generated using linear least-squares analysis and provide highly accurate estimates of this quantity as a function of molecular geometry (i.e., bond and angle variations) in a computationally tractable manner. An additional database, which consists of 25 representative molecular geometries and incorporates a more thorough treatment of both basis sets and core electron effects, is provided as a current benchmark for this quantity and the corresponding leading-order $C_6$ dispersion coefficient. This database has been utilized to assess the importance of these effects as well as the relative accuracy that can be obtained using several quantum chemical methods and a library of density functional approximations. In addition to high-level electron correlation methods (like CCSD) and our analytical least-squares formulae, we find that the SCAN0, PBE0, MN15, and B97-2 hybrid functionals yield the most accurate descriptions of the molecular polarizability tensor in H 2 O. Using first-order perturbation theory, we compute the zero-point vibrational correction to $α$ at the CCSDT/d-aug-cc-pVTZ level and find that this correction contributes approximately 3% to the isotropic ($α$ iso ) and nearly 50% to the anisotropic ($α$ aniso ) polarizability values. In doing so, we find that $α$ iso = 9.8307 bohr 3 , which is in excellent agreement with the experimental value of 9.83 ± 0.02 bohr 3 provided by Russell and Spackman. The DPS reported herein provides a benchmark-quality quantum mechanical estimate of this fundamental quantity of interest and should find extensive use in the development (and assessment) of next-generation force fields and machine-learning based approaches for modeling water in complex condensed-phase environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Pairing correlations in the cuprates: A numerical study of the three-band Hubbard model

We study the three-band Hubbard model for the copper oxide plane of the high-temperature superconducting cuprates using determinant quantum Monte Carlo and the dynamical cluster approximation (DCA) and provide a comprehensive view of the pairing correlations in this model using these methods. Specifically, we compute the pair-field susceptibility and study its dependence on temperature, doping, interaction strength, and charge-transfer energy. Using the DCA, we also solve the Bethe-Salpeter equation for the two-particle Green's function in the particle-particle channel to determine the transition temperature to the superconducting phase on smaller clusters. Our calculations reproduce many aspects of the cuprate phase diagram and indicate that there is an “optimal” value of the charge-transfer energy for the model where T c is largest. Finally, these results have implications for our understanding of superconductivity in both the cuprates and other doped charge-transfer insulators.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

The spin and orbital contributions to the total magnetic moments of free Fe, Co, and Ni clusters

We present size dependent spin and orbital magnetic moments of cobalt (Co{sub n}{sup +}, 8 ≤ n ≤ 22), iron (Fe{sub n}{sup +}, 7 ≤ n ≤ 17), and nickel cluster (Ni{sub n}{sup +}, 7 ≤ n ≤ 17) cations as obtained by X-ray magnetic circular dichroism (XMCD) spectroscopy of isolated clusters in the gas phase. The spin and orbital magnetic moments range between the corresponding atomic and bulk values in all three cases. We compare our findings to previous XMCD data, Stern-Gerlach data, and computational results. We discuss the application of scaling laws to the size dependent evolution of the spin and orbital magnetic moments per atom in the clusters. We find a spin scaling law “per cluster diameter,” ∼n{sup −1/3}, that interpolates between known atomic and bulk values. In remarkable contrast, the orbital moments do likewise only if the atomic asymptote is exempt. A concept of “primary” and “secondary” (induced) orbital moments is invoked for interpretation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Relativistic Few-Body Hadronic Physics Calculations

The goal of this research proposal was to use ``few-body'' methods to understand the structure and reactions of systems of interacting hadrons (neutrons, protons, mesons, quarks) over a broad range of energy scales. Realistic mathematical models of few-hadron systems have the advantage that they are sufficiently simple that they can be solved with mathematically controlled errors. These systems are also simple enough that it is possible to perform complete accurate experimental measurements on these systems. Comparison between theory and experiment puts strong constraints on the structure of the models. Even though these systems are ``simple'', both the experiments and computations push the limits of technology. The important property of ``few-body'' systems is that the ``cluster property'' implies that the interactions that appear in few-body systems are identical to the interactions that appear in complicated many-body systems. Of particular interest are models that correctly describe physics at distance scales that are sensitive to the internal structure of the individual nucleons. The Heisenberg uncertainty principle implies that in order to be sensitive to physics on distance scales that are a fraction of the proton or neutron radius, a relativistic treatment of quantum mechanics is necessary. The research supported by this grant involved 30 years of effort devoted to studying all aspects of interacting two and three-body systems. Realistic interactions were used to compute bound states of two- and three-nucleon, and two- and three-quark systems. Scattering observables for these systems were computed for a broad range of energies - from zero energy scattering to few GeV scattering, where experimental evidence of sub-nucleon degrees of freedom is beginning to appear. Benchmark calculations were produced, which when compared with calculations of other groups provided an essential check on these complicated calculations. In addition to computing bound state properties and scattering cross section, we also computed electron scattering cross sections in few-nucleon and few-quark systems, which are sensitive to the electric currents in these systems. We produced the definitive review on article on relativistic quantum mechanics, which and been used by many groups. In addition we developed and tested many computational techniques are used by other groups. Many of these techniques have applications in other areas of physics. The research benefited by collaborations with physicists from many different institutions and countries. It also involved working with seventeen undergraduate and graduate students.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Particle clustering effects on damage mechanisms in elastomeric syntactic foams

To elucidate the damage mechanisms in syntactic foams with hollow glass microballoon (GMB) reinforcement and elastomer matrices, in situ X-ray computed tomography mechanical testing was performed on syntactic foams with increasing GMB volume fraction. Image processing and digital volume correlation techniques identified very different damage mechanisms compared to syntactic foams with brittle matrices. Specifically, the prevailing mechanism transitioned from dispersed GMB collapse at low volume fraction to clustered GMB collapse at high volume fraction. Moreover, damage initiated and propagated earlier in closely-packed GMBs for all specimens. Both of these trends were attributed to increased interaction between closely-packed GMBs. This was validated by statistical analysis of GMB damage, which identified a consistent, inverse relationship among the probability of survival and the local coordination number (N neighbor ) across all specimens.

42 ENGINEERING↗

First-principles study of hydrogen diffusion and self-clustering below tungsten surfaces

The diffusion and self-clustering nucleation behavior of hydrogen (H) without vacancies below tungsten (W) surfaces are important for understanding the retention of deuterium (D) in W crystals exposed to high-flux and low-energy D plasma. The H migration energy and binding energy of H to self-clusters near the W(100), W(110), and W(111) surfaces have been investigated by first-principles computer simulations using density functional theory. H diffusion from adsorption sites on the W(100), W(110), and W(111) surfaces into solute sites in the bulk requires energies of at least 1.21, 1.78, and 1.80 eV, respectively, while 0.27, 0.31, and 0.24 eV for the reverse process. In addition, the lateral diffusion of H between two subsurface layers below the W surfaces has been investigated. Two H atoms at a depth of 0.08 nm below the W(110) surface have the highest binding energy, followed by H atom pairs below the W(111) and W(110) surfaces. The nucleation and stability of H clusters depend on the surface orientation. A planar configuration between the first nearest neighbor {100} planes is energetically favorable for H self-clustering below the W(100), W(110), and W(111) surfaces. In conclusion, the thermal stability of a platelet containing 16 H atoms below the W surfaces at 300 and 600 K was also studied by ab initio molecular dynamics simulations, which indicate that the H platelet below a W(111) surface is more stable than that below either the W(100) or the W(110) surface.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Benchmarks and Tests of a Multidimensional Cluster Dynamics Model of Helium Implantation in Tungsten

Here, we present a hierarchical multiscale modeling study of implanted helium (He) segregation near grain boundaries (GBs) of tungsten. We extend our spatially dependent cluster dynamics model to two spatial dimensions in order to take into account the biased drift of mobile He clusters toward the GBs observed in atomic-scale simulations. We are also able to reproduce the results from large-scale molecular dynamics simulations near and away from the GBs at low fluence with the extended cluster dynamics model. We suggest and verify that the sink (surface and GB) strengths are attenuated by the increasing concentration of He clusters at high fluence. This cluster dynamics model continues to set the stage for development of fully atomistically informed, coarse-grained models for computationally efficient predictions of He retention and surface morphological evolution, advancing progress toward the goal of efficient and optimal design of plasma-facing components.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Single-Walled Carbon Nanotubes with Confined Chalcogens as the Catalysts and Electrodes for Oxygen Reduction Reaction in Fuel Cells

The goal of this project is to synthesize and characterize a new non-metal electro-catalyst for oxygen reduction reaction (ORR) for fuel cell applications. The intended catalyst is a composite material composed of sulfur chains encapsulated in narrow diameter single-walled carbon nanotubes (S@SWNTs). S@SWNTs were successfully synthesized through sulfur vapor infusion method, and validated with Raman spectroscopy. However, the electrochemical analysis on the ORR catalytic activity of S@SWNTs indicated that it had low ORR catalytic activity. Our theoretical study based on density functional theory (DFT) revealed that the poor oxygen adsorption (low binding energy) on the surface of S@SWNTs was the bottleneck of the entire catalytic reaction. The focus of the project was subsequently pivoted to the development of a new non-noble metal ORR catalyst that could provide durability in acidic electrolyte. The hypothesis was to encapsulate small iron (Fe) clusters in SWNTs (Fe@SWNTs) can provide ORR electro-catalytic activity and long durability in acidic environment. DFT-based computational studies were carried out to explore the feasibility of the Fe@SWNTs catalyst. The theoretical study indicated that Fe@SWNTs indeed could catalyze the ORR with lower theoretical overpotential than platinum (Pt). However, its weaker bonding energy with oxygen was the bottleneck of the overall reaction. On the other hand, Fe clusters (composed of 7 Fe atoms) encapsulated in nitrogen-doped SWNTs (Fe 7 @N 4 WSNTs) showed proper oxygen adsorption by the Fe cluster and low theoretical overpotential of ORR. Comparing to the Fe single atom catalyst on N-doped SWNTs (Fe-N 4 SWNTs), which is one the best non-noble metal ORR catalysts reported in the literature, the Fe 7 @N 4 WSNTs showed lower overpotential and better resistance to acidic environment. Fe encapsulated N-doped SWNTs were synthesized with ferrocene as the Fe precursors through vapor infusion method, and experimental validate is underway. This study theoretically demonstrated the feasibility of a new type of non-noble metal electro-catalyst for ORR that could have high catalytic activity and long durability.

08 HYDROGEN↗

Understanding the Existence of a Na 2 Dimer in a High-Spin State

The recent observation of a high-spin Na 2 dimer formed on the surface of liquid helium nanodroplets raises some fundamental questions, as the ground state of Na 2 is known to have zero spin. Is it protected against spontaneous dissociation? What is its binding energy and interatomic distance? Is it stable at a higher temperature? Using calculations based on density functional theory (with and without long-range interaction) and coupled cluster methods, CCSD(T), we show that the bonding in the high-spin Na 2 dimer is governed by van der Waals interaction with binding energy (bond length) varying between −0.030 eV (5.108 Å) and −0.192 eV (4.231 Å), depending on the computational method used. Thus, the experimental method used by Kresin and coworkers can be very useful to study larger metastable high-spin clusters such as Li 4 which was predicted in 1985 to have a tetrahedral structure carrying a magnetic moment of 2 μ B , while its ground state is planar and nonmagnetic.

Basis sets↗

Cluster formation in precompound nuclei in the time-dependent framework

Background: Modern applications of nuclear time-dependent density functional theory (TDDFT) are often capable of providing quantitative description of heavy ion reactions. However, the structures of precompound (preequilibrium, prefission) states produced in heavy ion reactions are difficult to assess theoretically in TDDFT as the single-particle density alone is a weak indicator of shell structure and cluster states. Purpose: We employ the time-dependent nucleon localization function (NLF) to reveal the structure of precompound states in nuclear reactions involving light and medium-mass ions. We primarily focus on spin saturated systems with N = Z . Furthermore, we study reactions with oxygen and carbon ions, for which some experimental evidence for α clustering in precompound states exists. Method: We utilize the symmetry-free TDDFT approach with the Skyrme energy density functional UNEDF1 and compute the time-dependent NLFs to describe 16 O + 16 O, 40 Ca + 16 O, 40 Ca + 40 Ca , and 16,18 O + 12 C collisions at energies above the Coulomb barrier. Results: We show that NLFs reveal a variety of time-dependent modes involving cluster structures. For instance, the 16 O + 16 O collision results in a vibrational mode of a quasimolecular α - 12 C - 12 C- α state. For heavier ions, a variety of cluster configurations are predicted. For the collision of 16,18 O + 12 C, we showed that the precompound system has a tendency to form α clusters. This result supports the experimental findings that the presence of cluster structures in the projectile and target nuclei gives rise to strong entrance channel effects and enhanced α emission. Conclusion: The time-dependent nucleon localization measure is a very good indicator of cluster structures in complex precompound states formed in heavy-ion fusion reactions. Finally, the localization reveals the presence of collective vibrations involving cluster structures, which dominate the initial dynamics of the fusing system.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Global Transformations for Legacy Parallel Applications via Structural Analysis and Rewriting

Performance and scalability optimization of large HPC applications is currently a labor-intensive, manual process with very low productivity. Major difficulties come from the disaggregated environment for HPC application development: the compiler is only involved in local decisions (core or multithreaded domain), while a library-based, communication-oriented programming model realizes whole-machine parallelism. Realizing any major global change in such a disaggregated environment is very difficult and involves changing large portions of the source code. We present semi-automated techniques, based on structural analysis and rewriting, for performing global transformations on an HPC application source code. We present two case studies using the Self-Consistent Field (SCF) standalone benchmark as well as the Coupled Cluster (CCSD) module (2.9 million lines of Fortran code), a key module of the NWChem computational chemistry application. We demonstrate how structural rewriting techniques can be used to automate transformations that affect multiple sections of the application's source code. We show that the transformations can be applied in a systematic fashion across the source code bases with minimal manual effort. These transformations improve the scalability of the SCF benchmark by more than two orders of magnitude and the performance of the full CCSD module by a factor of 4.0.

legacy parallel applications↗

Green’s Function Coupled-Cluster Approach: Simulating Photoelectron Spectra for Realistic Molecular Systems

In this paper we present an efficient implementation for the analytical energy- dependent coupled cluster Green’s function with singles and doubles (GFCCSD) approach with our first practice being computing spectral functions of real molecular systems. Due to its algebraic structure, the presented method is highly scalable, and is capable of computing spectral function for a given molecular system in any energy region. Several typical examples have been given to demonstrate its capability of computing spectral functions not only in the valence band, but also in the core-level energy region. Satellite peaks have been observed in the inner valence band and core-level energy region where many-body effect becomes significant and single particle picture of ionization often breaks down. The accuracy test has been carried out by extensively comparing our GFCCSD results with experimental results as well as the theoretical results from other methods. It turns out GFCCSD method is able to provide qualitative or semi-quantitative level of description of ionization processes in both the core and valence regimes. To significantly improve the GFCCSD results for the main states, larger basis set can usually be employed, whereas the improvement of the GFCCSD results for the satellite states needs higher order many-body terms to be included in the GFCC implementation.

electronic structure theory↗

A map of the inorganic ternary metal nitrides

Exploratory synthesis in new chemical spaces is the essence of solid-state chemistry. However, uncharted chemical spaces can be difficult to navigate, especially when materials synthesis is challenging. Nitrides represent one such space, where stringent synthesis constraints have limited the exploration of this important class of functional materials. Here, we employ a suite of computational materials discovery and informatics tools to construct a large stability map of the inorganic ternary metal nitrides. Our map clusters the ternary nitrides into chemical families with distinct stability and metastability, and highlights hundreds of promising new ternary nitride spaces for experimental investigation - from which we experimentally realized seven new Zn- and Mg-based ternary nitrides. As a result, by extracting the mixed metallicity, ionicity and covalency of solid-state bonding from the density functional theory (DFT)-computed electron density, we reveal the complex interplay between chemistry, composition and electronic structure in governing large-scale stability trends in ternary nitride materials.

36 MATERIALS SCIENCE↗

Membrane lipids drive formation of KRAS4b-RAF1 RBDCRD nanoclusters on the membrane

The oncogene RAS, extensively studied for decades, presents persistent gaps in understanding, hindering the development of effective therapeutic strategies due to a lack of precise details on how RAS initiates MAPK signaling with RAF effector proteins at the plasma membrane. Recent advances in X-ray crystallography, cryo-EM, and super-resolution fluorescence microscopy offer structural and spatial insights, yet the molecular mechanisms involving protein-protein and protein-lipid interactions in RAS-mediated signaling require further characterization. This study utilizes single-molecule experimental techniques, nuclear magnetic resonance spectroscopy, and the computational Machine-Learned Modeling Infrastructure (MuMMI) to examine KRAS4b and RAF1 on a biologically relevant lipid bilayer. MuMMI captures long-timescale events while preserving detailed atomic descriptions, providing testable models for experimental validation. Both in vitro and computational studies reveal that RBDCRD binding alters KRAS lateral diffusion on the lipid bilayer, increasing cluster size and decreasing diffusion. RAS and membrane binding cause hydrophobic residues in the CRD region to penetrate the bilayer, stabilizing complexes through β-strand elongation. These cooperative interactions among lipids, KRAS4b, and RAF1 are proposed as essential for forming nanoclusters, potentially a critical step in MAP kinase signal activation.

59 BASIC BIOLOGICAL SCIENCES↗

Calculation of the Gibbs Free Energy of Solvation and Dissociation of HCl in Water via Monte Carlo Simulations and Continuum Solvation Models

The free energy of solvation and dissociation of hydrogen chloride in water is calculated through a combined molecular simulation quantum chemical approach at four temperatures between T = 300 and 450 K. The free energy is first decomposed into the sum of two components: the Gibbs free energy of transfer of molecular HCl from the vapor to the aqueous liquid phase and the standard-state free energy of acid dissociation of HCl in aqueous solution. The former quantity is calculated using Gibbs ensemble Monte Carlo simulations using either Kohn-Sham density functional theory or a molecular mechanics force field to determine the system’s potential energy. The latter free energy contribution is computed using a continuum solvation model utilizing either experimental reference data or micro-solvated clusters. The predicted combined solvation and dissociation free energies agree very well with available experimental data. CJM was supported by the US Department of Energy,Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences & Biosciences. Pacific Northwest National Laboratory is operated by Battelle for the US Department of Energy.

McGrath, Matthew↗

Sensitivity of CoSi{sub 2} precipitation in silicon to extra-low dopant concentrations. II. First-principles calculations

The paper is the second part of the study on the influence of very low dopant content in silicon on CoSi{sub 2} precipitation during high-temperature cobalt ion implantation into transmission electron microscope samples. It deals with the computational justification of various assumptions used in Paper I when rationalizing the kinetics of cobalt clustering in ion-implanted intrinsic silicon (both undoped and containing low concentrations of phosphorus atoms). In particular, it is proven that divacancies are efficient nucleation centers for the new Co-Si phase. It is shown that the capture of vacancies and divacancies on phosphorus atoms increases their lifetime in silicon matrix, but practically does not affect the mechanism of their interaction with interstitial cobalt atoms. Finally, it is demonstrated that the mobility of phosphorus interstitials at temperatures of our experiment is orders of magnitude higher than might be expected from the published literature data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗