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,369 records · Page 76

Cryogenic “Iodide-Tagging” Photoelectron Spectroscopy: A Sensitive Probe for Specific Binding Sites of Amino Acids

This work showcases cryogenic and temperature-dependent “iodide-tagging” photoelectron spectroscopy to probe specific binding sites of amino acids using the glycine-iodide complex (Gly?I-) as a case study. Multiple Gly?I- isomers were generated from ambient electrospray ionization and kinetically isolated in a cryogenic ion trap. These structures were characterized with temperature-dependent “iodide-tagging” negative ion photoelectron spectroscopy (NIPES) where iodide was used as the “messenger” to interpret electronic energetics and structural information of various Gly?I- isomers. Accompanied by theoretical computations and Franck-Condon simulations, a total of five cluster structures have been identified along with their various binding motifs. This work demonstrates that “iodide-tagging” NIPES is a powerful general means for probing specific binding interactions in biological molecules of interest.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Correlated Anion Disorder in Heteroanionic Cubic TiOF 2

Resolving anion configurations in heteroanionic materials is crucial for understanding and controlling their properties. For anion-disordered oxyfluorides, conventional Bragg diffraction cannot fully resolve the anionic structure, necessitating alternative structure determination methods. We have investigated the anionic structure of anion-disordered cubic (ReO 3 -type) TiOF 2 using X-ray pair distribution function (PDF), 19 F MAS NMR analysis, density functional theory (DFT), cluster expansion modeling, and genetic-algorithm structure prediction. Our computational data predict short-range anion ordering in TiOF 2 , characterized by predominant cis-[O 2 F 4 ] titanium coordination, resulting in correlated anion disorder at longer ranges. To validate our predictions, we generated partially disordered supercells using genetic-algorithm structure prediction and computed simulated X-ray PDF data and 19 F MAS NMR spectra, which we compared directly to experimental data. To construct our simulated 19 F NMR spectra, we derived new transformation functions for mapping calculated magnetic shieldings to predicted magnetic chemical shifts in titanium (oxy)fluorides, obtained by fitting DFT-calculated magnetic shieldings to previously published experimental chemical shift data for TiF 4 . We find good agreement between our simulated and experimental data, which supports our computationally predicted structural model and demonstrates the effectiveness of complementary experimental and computational techniques in resolving anionic structure in anion-disordered oxyfluorides. From additional DFT calculations, we predict that increasing anion disorder makes lithium intercalation more favorable by, on average, up to 2 eV, highlighting the significant effect of variations in short-range order on the intercalation properties of anion-disordered materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Electronic and chemical state of aluminum from the single- (K) and double-electron excitation (KL II&III , KL I ) x-ray absorption near-edge spectra of α-alumina, sodium aluminate, aqueous Al³⁺•(H₂O)₆, and aqueous Al(OH)₄⁻

We probe, at high energy resolution, the double electron excitation (KL II&II ) x-ray absorption region that lies approximately 115 eV above the main Al K-edge (1566 eV) of α-alumina and sodium aluminate. The two solid standards, α-alumina (octahedral) and sodium aluminate (tetrahedral) are compared to aqueous species that have the same Al coordination symmetries, Al³⁺•6H₂O (octahedral) and Al(OH)₄⁻ (tetrahedral). For the octahedral species, the edge height of the KL II&III -edge is approximately 10% of the main K-edge however the edge height is much weaker (3% of K-edge height) for Al species with tetrahedral symmetry. For the α-alumina and aqueous Al³⁺•6H₂O the KL II&III spectra contain white line features and extended absorption fine structure (EXAFS) that mimics the K-edge spectra. The KL II&III -edge feature interferes with an important region of the extended-XAFS region of the spectra for the K-edge of the crystalline and aqueous standards. The K-edge spectra and K-edge positions are predicted using time-dependent density functional theory (TDDFT). The TDDFT calculations for the K-edge XANES spectra reproduce the observed transitions in the experimental spectra of the four Al species. The KL II&III and KL I onsets and their corresponding chemical shifts for the four standards are estimated using the delta self-consistent field (ΔSCF) method. Research by JLF, NG, EJB, AV, TDS was supported by U.S. Department of Energy’s (DOE), Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences, and Biosciences. NG thanks Amity Andersen for help with the α-Al₂O₃ and tetrahedral sodium aluminate (NaAlO₂) clusters. All the calculations were performed using the Molecular Science Computing Capability at EMSL, a national scientific user facility sponsored by the U.S. Department of Energy’s Office of Biological and Environmental Research and located at Pacific Northwest National Laboratory (PNNL). PNNL is a multi-program national laboratory operated for DOE by Battelle under Contract # AC05-76RL01830.

aluminum↗

A comparison of neutral and charged species of one- and two-dimensional models of graphene nanoribbons using multireference theory

This study examines the dependence of the polyradical character of charged quasi-linear n-acenes and two-dimensional periacenes used as models for graphene nanoribbons in comparison to the corresponding neutral compounds. For this purpose, high-level ab initio calculations have been performed using the multireference averaged quadratic coupled cluster theory. Vertical ionization energies and electron affinities have been computed. Systematic tests show that the dependence on chain length of these quantities can be obtained from a consideration of the π system only and that remaining contributions coming from the σ orbitals or extended basis sets remain fairly constant. Using best estimate values, the experimental values for the ionization energy of the acene series can be reproduced within 0.1 eV and the experimental electron affinities within 0.4 V. The analysis of the natural orbital occupations and related unpaired electron densities shows that the ionic species exhibit a significant decrease in polyradical character and thus an increased chemical stability as compared to the neutral state.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Application of vibrational correlation formalism to internal conversion rate: Case study of Cu{sub n} (n = 3, 6, and 9) and H{sub 2}/Cu{sub 3}

This work reports non-radiative internal conversion (IC) rate constants obtained for Cu{sub n} with n = 3, 6, and 9 and H{sub 2} on Cu{sub 3}. The Time-Dependent Density Functional Theory (TDDFT) method was employed with three different functionals in order to investigate the electronic structures and the absorption spectra. The performance of the generalized gradient approximation of Perdew, Burke and Ernzerhof (PBE) and the hybrid B3LYP and PBE0 exchange correlation functionals in combination with the SVP and the def2-TZVP basis sets was examined. TDDFT results were used as input data to compute internal conversion rate constants. For this purpose, we have developed a program package. A description of the theoretical background used in our numerical implementation and the program input file is presented. In view of future applications of this program package in photoinduced catalysis, we present the analysis of the IC rate processes for the photodissociation of H{sub 2} on Cu{sub 3}. These results showed the applicability of the method and the computational program to identify the vibrational modes in transition metal clusters giving rise to the largest IC rate constant due to their interactions with the excited electronic states occurring in the hot-electron induced dissociation phenomena.

74 ATOMIC AND MOLECULAR PHYSICS↗

Coupled Cluster Green's function formulations based on the effective Hamiltonians

In this work, we demonstrate that the effective Hamiltonians obtained with the downfolding procedure based on double unitary coupled cluster (DUCC) ansatz can be used in the context of Green’s function coupled cluster (GFCC) formalism to calculate spectral functions of molecular systems. This combined approach (DUCC-GFCC) provides a significant reduction of numerical effort and good agreement with the corresponding all-orbital GFCC methods in energy windows that are consistent with the choice of active space. These features are demonstrated on the example of two benchmark systems: H 2 O and N 2 , where DUCC-GFCC calculations were performed for active spaces of various sizes.

74 ATOMIC AND MOLECULAR PHYSICS↗

Generalized representative structures for atomistic systems

A new method is presented to generate atomic structures that reproduce the essential characteristics of arbitrary material systems, phases, or ensembles. Previous methods allow one to reproduce the essential characteristics (e.g. the chemical disorder) of a large random alloy within a small crystal structure. The ability to generate small representations of random alloys, along with the restriction to crystal systems, results from using the fixed-lattice cluster correlations to describe structural characteristics. A more general description of the structural characteristics of atomic systems is obtained using complete sets of atomic environment descriptors. These are used within for generating representative atomic structures without restriction to fixed lattices. A general data-driven approach is provided here utilizing the atomic cluster expansion (ACE) basis. The N-body ACE descriptors are a complete set of atomic environment descriptors that span both chemical and spatial degrees of freedom and are used within for describing atomic structures. The generalized representative structure (GRS) method presented within generates small atomic structures that reproduce ACE descriptor distributions corresponding to arbitrary structural and chemical complexity. It is shown that systematically improvable representations of crystalline systems on fixed parent lattices, amorphous materials, liquids, and ensembles of atomic structures may be produced efficiently through optimization algorithms. With the GRS method, we highlight reduced representations of atomistic machine-learning training datasets that contain similar amounts of information and small 40–72 atom representations of liquid phases. The ability to use GRS methodology as a driver for informed novel structure generation is also demonstrated. The advantages over other data-driven methods and state-of-the-art methods restricted to high-symmetry systems are highlighted.

atomic cluster expansion↗

Comparison of weak lensing by NFW and Einasto halos and systematic errors

Recent N-body simulations have shown that Einasto radial profiles provide the most accurate description of dark matter halos. Predictions based on the traditional NFW functional form may fail to describe the structural properties of cosmic objects at the percent level required by precision cosmology. We computed the systematic errors expected for weak lensing analyses of clusters of galaxies if one wrongly models the lens density profile. Even though the NFW fits of observed tangential shear profiles can be excellent, viral masses and concentrations of very massive halos (∼> 10{sup 15}M{sub ⊙}/h) can be over- and underestimated by 0∼ 1 per cent, respectively. Misfitting effects also steepen the observed mass-concentration relation, as observed in multi-wavelength observations of galaxy groups and clusters. Based on shear analyses, Einasto and NFW halos can be set apart either with deep observations of exceptionally massive structures (∼> 2×10{sup 15}M{sub ⊙}/h) or by stacking the shear profiles of thousands of group-sized lenses (∼> 10{sup 14}M{sub ⊙}/h)

79 ASTRONOMY AND ASTROPHYSICS↗

Large-scale 16S gene assembly using metagenomics shotgun sequences

Combining a 16S rRNA (16S) gene database with metagenomic shotgun sequences allows unbiased identification of known and novel microbes.To attain this, we herein report reference-based ribosome assembly (RAMBL), a computational pipeline, which integrates taxonomic tree search and Dirichlet process clustering to reconstruct full-length 16S gene sequences from metagenomic sequencing data with high accuracy. By benchmarking against the synthetic and real shotgun sequences, we demonstrated that full-length 16S gene assemblies of RAMBL were a good proxy for known and putative microbes, including Candidate Phyla Radiation. Here, we found that 30-40% of bacteria genera in the terrestrial and intestinal biomes have no closely related genome sequences. We also observed that RAMBL was able to generate a more accurate determination of environmental microbial diversity and yield better disease classification, suggesting that full-length 16S gene assemblies are a powerful alternative to marker gene set and 16S short reads. RAMBL first realizes the access to full-length 16S gene sequences in the near-terabase-scale metagenomic shotgun sequences, which markedly improve metagenomic data analysis and interpretation.

59 BASIC BIOLOGICAL SCIENCES↗

Partitioning of Large-Scale Power Electronics-Based Power Systems for Small-Signal Stability Analysis

The nodal admittance matrix (NAM)-based approach is suitable for analyzing the small-signal stability of large-scale power electronics-based power systems (PEPSs) as it preserves the system structure by utilizing the admittance matrix. Previously, NAM-based area partition has been proposed, which divides the system into various subareas and interconnections for easier analysis of the low-dimension matrix compared to the entire system-based high-dimension matrix. However, no partition algorithm has been presented for the NAM-based area partition method. This paper focuses on implementing the spectral partitioning algorithm for partitioning large-scale PEPSs into a low-dimension matrix to reduce the computation complexity of the analysis. These spectral components facilitate data transformation into a new space, enabling the application of traditional clustering methods like k-means. To evaluate the performance of the partitioning method, the subareas and interconnections obtained from the spectral clustering algorithm are incorporated into the NAM-based area partition method for a large system with 140 buses. The computational times of the original method, where the NAM-based criterion is directly applied to the entire system, are compared with those of the NAM-based partition method in MATLAB. PSCAD simulations of the whole system and the obtained subareas are conducted to validate the effectiveness of the proposed algorithm.

Nupur, Nupur↗

Designing Catalysts for Dehydrogenation of Methane for Reducing Greenhouse Gas during Natural Gas Extraction

Catalytic conversion of methane (CH 4 ) into useful products is critical for maximizing the utility of natural gas output and for reducing green house gas release associated with flaring (burning off CH 4 at natural gas extraction sites). One particular useful technique is methane dry reforming (DRM), which involves the chemical reaction of CH 4 with carbon dioxide (CO 2 ) to generate carbon monoxide (CO), hydrogen gas (H 2 ), and subsequently other useful products. New and improved catalysts are required to facilitate efficient dry methane reforming. In this report, we apply the Density Functional Theory (DFT) computational technique to investigate a catalyst consisting of small nickel clusters (Ni n , n < 10) on ceria (Ce0 2 (111) surfaces) support. One main thrust of this project is to study the initial CH 4 and CO 2 reactions with the catalyst. We find that CH 4 exhibits barrierless reactive adsorption on to the catalyst. In order words, this step is likely not the rate-determining step. A second thrust is to perform detailed studies of the catalyst itself and examine the role of oxygen vacancies. Using a specific DFT method and a hypothesis about the absence of the Ce(III) redox state, we obtain predictions about oxygen vacancies in good agreement with experimental observations.

03 NATURAL GAS↗

Nuclear Theory from First Principles to Forefront Experiments (Final Report for DE-SC0018638)

The Lee research group is a part of the Nuclear Lattice Effective Field Theory Collaboration, which has developed and performed ab initio lattice simulations of nuclear structure, scattering, reactions, and many-body systems. The other senior members of the collaboration include Ulf-G. Meißner at Bonn/Julich, Evgeny Epelbaum and Hermann Krebs at Bochum, Timo Lahde and Thomas Luu at Julich, and Gautam Rupak at Mississippi State. Our letter “Ab initio alpha-alpha scattering” was featured in a Nature News and Views article. Another letter “Nuclear binding near a quantum phase transition” was highlighted in a Viewpoint article in the online APS journal Physics as well as a news article in the IOP publication Physics World (September 21, 2016). Our letter “Eigenvector continuation with subspace learning” was also highlighted a Synopsis article in Physics. There are many promising ab initio approaches being used to calculate the properties of few-and many-nucleon systems. This includes the no-core shell model, symmetry-adapted no-core shell model quantum Monte Carlo, auxiliary-field diffusion Monte Carlo, self-consistent Green’s functions, many-body perturbation theory, in-medium similarity renormalization group, and coupled cluster methods. Lattice effective field theory is another ab initio approach which combines the framework of effective field theory with lattice Monte Carlo methods to allow favorable scaling from few- to many-body systems. Perhaps the most important aspect of lattice effective field theory is that its strengths and weaknesses are orthogonal to that of other ab initio methods. For example, lattice effective field theory has little difficulty in probing cluster structures and collectivity, while such features are much more difficult using other methods. Furthermore it can be used to compute superfluid condensate fractions as well as the phase diagram of strongly matter and the density and temperature dependence of clustering. Lattice effective field theory was first used in simulations of infinite nuclear matter and infinite neutron matter with pions and without pions. In addition to the efforts by our collaboration, there have been recent efforts by other groups as well.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Nuclear Computational Low Energy Initiative (NUCLEI)

The NUCLEI project, as defined by the scope of work, developed, implemented and run codes for large-scale computations of many topics in low-energy nuclear physics. Physics studied include the properties of nuclei and nuclear decays, nuclear structure and reactions, and the properties of nuclear matter. The computational techniques used include Quantum Monte Carlo, Configuration Interaction, Coupled Cluster, and Density Functional methods. The research program emphasized areas of high interest to current and possible future DOE nuclear physics facilities, including ATLAS and FRIB (nuclear structure and reactions, and nuclear astrophysics), TJNAF (neutron distributions in nuclei, few body systems, and electroweak processes), NIF (thermonuclear reactions), MAJORANA and FNPB (neutrinoless double-beta decay and physics beyond the Standard Model), and LANSCE (fission studies).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Interconversion of diborane(4) isomers

Highly correlated electronic structure computations using many-body perturbation theory and coupled-cluster gradient techniques are used to study the reaction pathway that links the two forms (C2u and D2d) of diborane(4). The results obtained indicate that a low-energy pathway exists for interconversion of the two low-lying isomers of diborane(4). The proposed mechanism consists of a single concerted but nonsynchronous rotation of the BH2 groups. The pathway first follows an idealized reaction coordinate which preserves C2 symmetry, but then bifurcates at a branch point, leading to two equivalent transition states which lack nontrivial elements of symmetry.

Stanton, John F.↗

Experimental Program to Stimulate Competitive Research (EPSCoR)

Report includes: (1) CLUSTER: "Studies in Macromolecular Behavior in Microgravity Environment": The Role of Protein Oligomers in Protein Crystallization; Phase Separation Phenomena in Microgravity; Traveling Front Polymerizations; Investigating Mechanisms Affecting Phase Transition Response and Changes in Thermal Transport Properties in ER-Fluids under Normal and Microgravity Conditions. (2) CLUSTER: "Computational/Parallel Processing Studies": Flows in Local Chemical Equilibrium; A Computational Method for Solving Very Large Problems; Modeling of Cavitating Flows.

Dingerson, Michael R.↗

Nitride Fuel Potential for Power Uprates

This work presents the methodology and current results for design analysis of grid supported pressurized water reactor cores using uranium nitride (UN) fuel in square arrays. Uranium nitride has a higher theoretical density and higher thermal conductivity than uranium oxide, allowing for lower fuel temperatures, longer cycle length and higher power densities with similar enrichment. In order to avoid the UN and high temperature water reaction, silicon carbide cladding has been proposed in place of Zircaloy cladding. The analysis aims at finding the maximum thermal power over a range of bundle geometries (rod diameter and pitch). The same methodology is applied to the design of UO{sub 2} oxide fueled cores to provide a fair comparison of the achievable power rating between the two fuel types. Steady-state limits are considered, including fuel bundle pressure drop, minimum departure from nucleate boiling ratio (MDNBR), fuel temperature, fuel rod vibrations and pin fretting and sliding wear. Transients such as loss of flow and overpower accident will be considered. A steady state analysis was undertaken to determine the maximum achievable power that can be sustained by square array fuel assemblies loaded with nitride and oxide fuels with same or better safety performances as the reference design. Steady-state design limits were applied including MDNBR, fuel bundle pressure drop, fuel temperature, cumulative sliding and fretting wear. The major conclusions of this study are as follows: 1) It is possible to achieve significant power increase for UN and UO{sub 2} compared to the reference core power output of 3400 MWth with smaller but more numerous rods and higher pitch over diameter ratios. A power uprate of 25% can be reached according to the VIPRE-MATLAB interface for full core analysis with the same number of spacer grids, and 32% according to the new MATLAB script for single hot channel analysis, with five more spacer grids. For this latter core, a higher enrichment would be required. 2) The nitride fuel temperature constraint from dissociation doesn't influence steady state performances as other design constraints such as MDNBR and pressure loss are more restrictive. The oxide fuel temperature constraint from fission gas release is reached in part of the study domain, but not in the region that allows for possible uprates. 3) As a result, the thermal hydraulic performance of UN and UO{sub 2} are similar. Differences in possible power uprates will arise from fretting and sliding wear. Wear depends of the time spent by the fuel in the core, which depends on burnup and heavy metal inventory. UN being denser, it can allow for power uprates or increased cycle length. There is a trade off between power uprates and cycle length with regards to fuel pin wear. Reducing fuel pin wear by reducing cycle length can allow for higher thermal power. This affects the plants' economics as well. This work will be completed by a transient analysis using CASMO4E, SIMULATE3 and SIMULATE3k to model the neutronics of a nitride fueled core. The transient performances will be used to determine the maximum power uprates while performing as well or better during transients as the reference design. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

The Clustering of High-redshift (2.9 ≤ $z$ ≤ 5.1) Quasars in SDSS Stripe 82

We present a measurement of the two-point autocorrelation function of photometrically selected high-$z$ quasars over ~100 deg 2 on the Sloan Digital Sky Survey Stripe 82 field. Selection is performed using three machine-learning algorithms in a six-dimensional optical/mid-infrared color space. Optical data from the Sloan Digital Sky Survey are combined with overlapping deep mid-infrared data from the $Spitzer$ IRAC Equatorial Survey and the $Spitzer$-HETDEX Exploratory Large-Area survey. Our selection algorithms are trained on the colors of known high-$z$ quasars. The selected quasar sample consists of 1378 objects and contains both spectroscopically confirmed quasars and photometrically selected quasar candidates. These objects span a redshift range of 2.9 ≤ $z$ ≤ 5.1 and are generally fainter than $i$ = 20.2, a regime that has lacked sufficient number density to perform autocorrelation function measurements of photometrically classified quasars. We compute the angular correlation function of these data, marginally detecting quasar clustering. We fit a single power law with an index of $δ$ = 1.39 ± 0.618 and amplitude of $θ_0$ = 0$^{\prime}_{.}$71 ± 0$^{\prime}_{.}$546 . A dark matter model is fit to the angular correlation function to estimate the linear bias. At the average redshift of our survey ($\langle z\rangle =3.38$), the bias is $b$ = 6.78 ± 1.79. Using this bias, we calculate a characteristic dark matter halo mass of 1.70–9.83$\times {10}^{12}{h}^{-1}\,{M}_{\odot }$. Our bias estimate suggests that quasar feedback intermittently shuts down the accretion of gas onto the central supermassive black hole at early times. If confirmed, these results hint at a level of luminosity dependence in the clustering of quasars at high-$z$.

79 ASTRONOMY AND ASTROPHYSICS↗

Efficient calculation of nuclear spin-rotation constants from auxiliary density functional theory

The computation of the spin-rotation tensor within the framework of auxiliary density functional theory (ADFT) in combination with the gauge including atomic orbital (GIAO) scheme, to treat the gauge origin problem, is presented. For the spin-rotation tensor, the calculation of the magnetic shielding tensor represents the most demanding computational task. Employing the ADFT-GIAO methodology, the central processing unit time for the magnetic shielding tensor calculation can be dramatically reduced. In this work, the quality of spin-rotation constants obtained with the ADFT-GIAO methodology is compared with available experimental data as well as with other theoretical results at the Hartree-Fock and coupled-cluster level of theory. It is found that the agreement between the ADFT-GIAO results and the experiment is good and very similar to the ones obtained by the coupled-cluster single-doubles-perturbative triples-GIAO methodology. With the improved computational performance achieved, the computation of the spin-rotation tensors of large systems or along Born-Oppenheimer molecular dynamics trajectories becomes feasible in reasonable times. Three models of carbon fullerenes containing hundreds of atoms and thousands of basis functions are used for benchmarking the performance. Furthermore, a theoretical study of temperature effects on the structure and spin-rotation tensor of the H{sup 12}C–{sup 12}CH–DF complex is presented. Here, the temperature dependency of the spin-rotation tensor of the fluorine nucleus can be used to identify experimentally the so far unknown bent isomer of this complex. To the best of our knowledge this is the first time that temperature effects on the spin-rotation tensor are investigated.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗