Search NASA⌕ Search

SEARCH · Search NASA

Results for “cluster methods”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Theoretical framework for new magnetic materials for quantum computing and information storage. Final report for the Award No. DE-SC0018910

The focus of this grant was on molecular magnetic materials for information storage and quantum computing. We have been developing robust, first-principle methods for computing relevant electronic and magnetic properties of molecular building blocks (SMMs) of novel magnetic materials and quantum computers. These tools enable theoretical modeling of SMMs’ behavior, facilitating the interpretation of experimental studies and aiding the design of novel magnetic materials. Our strategy is based on the spin-flip (SF) approach, which extends the hierarchy of black-box single-reference methods to strongly correlated systems. Specifically, we developed general scalable algorithms and computer codes for calculating molecular properties, with an emphasis on spin-related properties, such as zero-field splittings, hyperfine couplings, and g-tensors. While our primary focus was on SF wave functions and SF-TDDFT, the underlying theory and computer codes were formulated using reduced density matrices, such that these tools are applicable to a broader class of methods. To extend the scope of applicability of wave-function-based SF methods to larger systems, we developed reduced-scaling approaches for the equation-of-motion coupled-cluster (EOM-CC) methods and continue developing libtensor (our open-source general tensor contraction library for many-body methods). We carried out extensive benchmarks and also carried out several applications.

36 MATERIALS SCIENCE↗
NASA NTRS2012

A New Method to Constrain Supernova Fractions Using X-ray Observations of Clusters of Galaxies

Supernova (SN) explosions enrich the intracluster medium (ICM) both by creating and dispersing metals. We introduce a method to measure the number of SNe and relative contribution of Type Ia supernovae (SNe Ia) and core-collapse supernovae (SNe cc) by directly fitting X-ray spectral observations. The method has been implemented as an XSPEC model called snapec. snapec utilizes a single-temperature thermal plasma code (apec) to model the spectral emission based on metal abundances calculated using the latest SN yields from SN Ia and SN cc explosion models. This approach provides a self-consistent single set of uncertainties on the total number of SN explosions and relative fraction of SN types in the ICM over the cluster lifetime by directly allowing these parameters to be determined by SN yields provided by simulations. We apply our approach to XMM-Newton European Photon Imaging Camera (EPIC), Reflection Grating Spectrometer (RGS), and 200 ks simulated Astro-H observations of a cooling flow cluster, A3112.We find that various sets of SN yields present in the literature produce an acceptable fit to the EPIC and RGS spectra of A3112. We infer that 30.3% plus or minus 5.4% to 37.1% plus or minus 7.1% of the total SN explosions are SNe Ia, and the total number of SN explosions required to create the observed metals is in the range of (1.06 plus or minus 0.34) x 10(exp 9), to (1.28 plus or minus 0.43) x 10(exp 9), fromsnapec fits to RGS spectra. These values may be compared to the enrichment expected based on well-established empirically measured SN rates per star formed. The proportions of SNe Ia and SNe cc inferred to have enriched the ICM in the inner 52 kiloparsecs of A3112 is consistent with these specific rates, if one applies a correction for the metals locked up in stars. At the same time, the inferred level of SN enrichment corresponds to a star-to-gas mass ratio that is several times greater than the 10% estimated globally for clusters in the A3112 mass range.

Observations↗
DOE OSTI2021

A Massively Parallel Implementation of the CCSD(T) Method Using the Resolution-of-the-Identity Approximation and a Hybrid Distributed/Shared Memory Parallelization Model

In this work, a parallel algorithm is described for the coupled-cluster singles and doubles method augmented with a perturbative correction for triple excitations [CCSD(T)] using the resolution-of-the-identity (RI) approximation for two-electron repulsion integrals (ERIs). The algorithm bypasses the storage of four-center ERIs by adopting an integral-direct strategy. The CCSD amplitude equations are given in a compact quasi-linear form by factorizing them in terms of amplitude-dressed three-center intermediates. A hybrid MPI/OpenMP parallelization scheme is employed, which uses the OpenMP-based shared memory model for intranode parallelization and the MPI-based distributed memory model for internode parallelization. Parallel efficiency has been optimized for all terms in the CCSD amplitude equations. Two different algorithms have been implemented for the rate-limiting terms in the CCSD amplitude equations that entail and -scaling computational costs, where N O and N V denote the number of correlated occupied and virtual orbitals, respectively. One of the algorithms assembles the four-center ERIs requiring N V 4 and N O 2 N V 2 -scaling memory costs in a distributed manner on a number of MPI ranks, while the other algorithm completely bypasses the assembling of quartic memory-scaling ERIs and thus largely reduces the memory demand. It is demonstrated that the former memory-expensive algorithm is faster on a few hundred cores, while the latter memory-economic algorithm shows a better strong scaling in the limit of a few thousand cores. The program is shown to exhibit a near-linear scaling, in particular for the compute-intensive triples correction step, on up to 8000 cores. The performance of the program is demonstrated via calculations involving molecules with 24–51 atoms and up to 1624 atomic basis functions. As the first application, the complete basis set (CBS) limit for the interaction energy of the π-stacked uracil dimer from the S66 data set has been investigated. This work reports the first calculation of the interaction energy at the CCSD(T)/aug-cc-pVQZ level without local orbital approximation. The CBS limit for the CCSD correlation contribution to the interaction energy was found to be -8.01 kcal/mol, which agrees very well with the value -7.99 kcal/mol reported by Schmitz, Hättig, and Tew [ Phys. Chem. Chem. Phys. 2014 , 16 , 22167-22178]. The CBS limit for the total interaction energy was estimated to be -9.64 kcal/mol.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2011

Improve Data Mining and Knowledge Discovery Through the Use of MatLab

Data mining is widely used to mine business, engineering, and scientific data. Data mining uses pattern based queries, searches, or other analyses of one or more electronic databases/datasets in order to discover or locate a predictive pattern or anomaly indicative of system failure, criminal or terrorist activity, etc. There are various algorithms, techniques and methods used to mine data; including neural networks, genetic algorithms, decision trees, nearest neighbor method, rule induction association analysis, slice and dice, segmentation, and clustering. These algorithms, techniques and methods used to detect patterns in a dataset, have been used in the development of numerous open source and commercially available products and technology for data mining. Data mining is best realized when latent information in a large quantity of data stored is discovered. No one technique solves all data mining problems; challenges are to select algorithms or methods appropriate to strengthen data/text mining and trending within given datasets. In recent years, throughout industry, academia and government agencies, thousands of data systems have been designed and tailored to serve specific engineering and business needs. Many of these systems use databases with relational algebra and structured query language to categorize and retrieve data. In these systems, data analyses are limited and require prior explicit knowledge of metadata and database relations; lacking exploratory data mining and discoveries of latent information. This presentation introduces MatLab(R) (MATrix LABoratory), an engineering and scientific data analyses tool to perform data mining. MatLab was originally intended to perform purely numerical calculations (a glorified calculator). Now, in addition to having hundreds of mathematical functions, it is a programming language with hundreds built in standard functions and numerous available toolboxes. MatLab's ease of data processing, visualization and its enormous availability of built in functionalities and toolboxes make it suitable to perform numerical computations and simulations as well as a data mining tool. Engineers and scientists can take advantage of the readily available functions/toolboxes to gain wider insight in their perspective data mining experiments.

Shaykhian, Gholam Ali↗
NASA NTRS2011

Improve Data Mining and Knowledge Discovery through the use of MatLab

Data mining is widely used to mine business, engineering, and scientific data. Data mining uses pattern based queries, searches, or other analyses of one or more electronic databases/datasets in order to discover or locate a predictive pattern or anomaly indicative of system failure, criminal or terrorist activity, etc. There are various algorithms, techniques and methods used to mine data; including neural networks, genetic algorithms, decision trees, nearest neighbor method, rule induction association analysis, slice and dice, segmentation, and clustering. These algorithms, techniques and methods used to detect patterns in a dataset, have been used in the development of numerous open source and commercially available products and technology for data mining. Data mining is best realized when latent information in a large quantity of data stored is discovered. No one technique solves all data mining problems; challenges are to select algorithms or methods appropriate to strengthen data/text mining and trending within given datasets. In recent years, throughout industry, academia and government agencies, thousands of data systems have been designed and tailored to serve specific engineering and business needs. Many of these systems use databases with relational algebra and structured query language to categorize and retrieve data. In these systems, data analyses are limited and require prior explicit knowledge of metadata and database relations; lacking exploratory data mining and discoveries of latent information. This presentation introduces MatLab(TradeMark)(MATrix LABoratory), an engineering and scientific data analyses tool to perform data mining. MatLab was originally intended to perform purely numerical calculations (a glorified calculator). Now, in addition to having hundreds of mathematical functions, it is a programming language with hundreds built in standard functions and numerous available toolboxes. MatLab's ease of data processing, visualization and its enormous availability of built in functionalities and toolboxes make it suitable to perform numerical computations and simulations as well as a data mining tool. Engineers and scientists can take advantage of the readily available functions/toolboxes to gain wider insight in their perspective data mining experiments.

Shaykahian, Gholan Ali↗
NASA NTRS1993

The distance to the Large Magellanic Cloud - Constraints from Cepheids in Large Magellanic Cloud star clusters

We have used recent observational data for the Cepheids in the rich, young LMC clusters NGC 1866 and NGC 2031 to constrain the cluster distances with the mass-equivalency (ME) method. The basis of this approach is to fix the cluster distances by requiring the Cepheid evolutionary and pulsational masses to be equal. Using evolutionary models incorporating a mild amount of core and envelope overshooting along with recent pulsational models, we derive distance moduli of 18.51 +/- 0.21 and 18.32 +/- 0.20 for NGC 1866 and NGC 2031, respectively. The quoted errors are dominated by the uncertainties in the heavy element abundances of the clusters (assumed to be 0.3 dex for both clusters), with a smaller contribution due to the apparently intrinsic spread in the masses of the Cepheids in each cluster. For the ME method, we find that Delta(m-M)(0)/Delta sub Z(1) = 0.69, where Z(1) = log (Z/0.016). This result implies that the cluster distances can be determined to better than +/- 5 percent if the cluster abundances can be measured to better than about +/- 0.15 dex.

Bertelli, Gianpaolo↗
DOE OSTI2018

Diagnostics for Undersampling and Clustering in Monte Carlo Criticality Calculations

For the past 40 years, it has been known from the theoretical work of Gelbard [1] and Brissenden [2] that the results from Monte Carlo (MC) criticality calculations are biased if the number of neutrons per generation is not sufficiently large. The computed k-effective eigenvalue will be lower than the actual eigenvalue, and the computed eigenfunction will be too low in high-importance regions and too high in low-importance regions. This difficulty is commonly called the undersampling problem for criticality calculations. The biases for undersampling were confirmed for practical, realistic problems with modern MC codes by Brown [3]. The solution to undersampling is simple - use a larger number of neutrons per generation to reduce the biases to levels that are negligible. Over the past few years, the OECD-NEA-WPNCS Expert Group on Advanced Monte Carlo Techniques (EGAMCT) [4,5,6] investigated an extreme form of undersampling called clustering. When the number of neutrons per generation is small, neutrons in successive generations tend to group together (cluster) due to correlations, and the clusters tend to migrate through the problem space. While the problems of undersampling and clustering are well known and well understood theoretically, there were no existing diagnostics for MC codes to determine whether a given calculation is using sufficient neutrons per generation to eliminate these problems. That is, no tools were available to diagnose biases due to undersampling or clustering. The tools described herein resolve this longstanding problem by providing 2 robust methods for detecting undersampling and clustering in MC criticality calculations. The diagnosis of undersampling and clustering is part of a larger effort to thoroughly overhaul the calculational methods for criticality problems, providing for automated sampling of the initial fission distribution, acceleration of the iteration convergence process, automated detection of convergence in the iterations and starting of the active tally cycles, and the diagnosis of undersampling. All of these new methods have been prototyped in a local modified version of mcnp6.2 [7] and tested on a variety of critical systems.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS1991

An efficient formulation and implementation of the analytic energy gradient method to the single and double excitation coupled-cluster wave function - Application to Cl2O2

The analytic energy gradient for the single and double excitation coupled-cluster (CCSD) wave function has been reformulated and implemented in a new set of programs. The reformulated set of gradient equations have a smaller computational cost than any previously published. The iterative solution of the linear equations and the construction of the effective density matrices are fully vectorized, being based on matrix multiplications. The new method has been used to investigate the Cl2O2 molecule, which has recently been postulated as an important intermediate in the destruction of ozone in the stratosphere. In addition to reporting computational timings, the CCSD equilibrium geometries, harmonic vibrational frequencies, infrared intensities, and relative energetics of three isomers of Cl2O2 are presented.

Rendell, Alistair P.↗
NASA NTRS1993

The first detection of a collapsed core globular cluster in M31

We report on the observations of a globular cluster (designated G105 = Bo 343) in the nearby spiral galaxy M31, using the Hubble Space Telescope (HST). Image deconvolutions using three different methods indicate that this cluster has a density profile with the morphology characteristic of stellar systems which have undergone core collapse, i.e., a power-law density cusp near the center. This is the first such detection in a galaxy as distant as M31. This discovery may lead to extensive future statistical studies of the globular cluster system of M31, and thus a better understanding of its evolution.

Bendinelli, O.↗
DOE OSTI2020

Extension of frozen natural orbital approximation to open-shell references: Theory, implementation, and application to single-molecule magnets

Natural orbitals are often used to achieve a more compact representation of correlated wave-functions. Using natural orbitals computed as eigenstates of the virtual–virtual block of the state density matrix instead of the canonical Hartree–Fock orbitals results in smaller errors when the same fraction of virtual space is frozen. This strategy, termed frozen natural orbital (FNO) approach, is effective in reducing the cost of regular coupled-cluster (CC) calculations and some multistate methods, such as EOM-IP-CC (equation-of-motion CC for ionization potentials). This contribution extends the FNO approach to the EOM-SF-CC ansatz (EOM-CC with spin-flip). In contrast to EOM-IP-CCSD, EOM-SF-CCSD relies on high-spin open-shell references. Using FNOs computed for an open-shell reference leads to an erratic behavior of the EOM-SF-CC energies and properties due to an inconsistent truncation of the α and β orbital spaces. A general solution to problems arising in the EOM-CC calculations utilizing open-shell references, termed OSFNO (open-shell FNO), is proposed. Furthermore, by means of singular value decomposition (SVD) of the overlap matrix of the α and β orbitals, the OSFNO algorithm identifies the corresponding orbitals and determines virtual orbitals corresponding to the singly occupied space. This is followed by SVD of the singlet part of the state density matrix in the remaining virtual orbital subspace. The so-computed FNOs preserve the spin purity of the open-shell orbital subspace to the extent allowed by the original reference, thus facilitating a safe truncation of the virtual space. The performance of OSFNO is benchmarked for selected diradicals and triradicals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1999

The Sunyaev-Zel'dovich Effect in Abell 370

We present interferometric measurements of the Sunyaev-Zel'dovich (SZ) effect towards the galaxy cluster Abell 370. These measurements, which directly probe the pressure of the cluster's gas, show the gas is strongly aspherical, on agreement with the morphology revealed by x-ray and gravitational lensing observations. We calculate the cluster's gas mass fraction by comparing the gas mass derived from the SZ measurements to the lensing-derived gravitational mass near the critical lensing radius. We also calculate the gas mass fraction from the SZ data by deriving the total mass under the assumption that the gas is in hydrostatic equilibrium (HSE). We test the assumptions in the HSE method by comparing the total cluster mass implied by the two methods. The Hubble constant derived for this cluster, when the known systematic uncertainties are included, has a very wide range of values and therefore does not provide additional constraints on the validity of the assumptions. We examine carefully the possible systematic errors in the gas fraction measurement. The gas fraction is a lower limit to the cluster's baryon fraction and so we compare the gas mass fraction, calibrated by numerical simulations to approximately the virial radius, to measurements of the global mass fraction of baryonic matter, OMEGA(sub B)/OMEGA(sub matter). Our lower limit to the cluster baryon fraction is f(sub B) = (0.043 +/- 0.014)/h (sub 100). From this, we derive an upper limit to the universal matter density, OMEGA(sub matter) <= 0.72/h(sub 100), and a likely value of OMEGA(sub matter) <= (0.44(sup 0.15, sub -0.12)/h(sub 100).

Grego, Laura↗
DOE OSTI2022

Abundance matching analysis of the emission-line galaxy sample in the extended Baryon Oscillation Spectroscopic Survey

We present the measurements of the small-scale clustering for the emission-line galaxy (ELG) sample from the extended Baryon Oscillation Spectroscopic Survey (eBOSS) in the Sloan Digital Sky Survey IV (SDSS-IV). We use conditional abundance matching method to interpret the clustering measurements from 0.34 to 70 h -1 Mpc. In order to account for the correlation between properties of ELGs and their environment, we add a secondary connection between star formation rate of ELGs and halo accretion rate. Three parameters are introduced to model the ELG [O II ] luminosity and to mimic the target selection of eBOSS ELGs. The parameters in our models are optimized using Markov Chain Monte Carlo (MCMC) method. We find that by conditionally matching star formation rate of galaxies and the halo accretion rate, we are able to reproduce the eBOSS ELG small-scale clustering within 1σ error level. Our best-fitting model shows that the eBOSS ELG sample only consists of ~ 12 per cent of all star-forming galaxies, and the satellite fraction of eBOSS ELG sample is 19.3 per cent. We show that the effect of assembly bias is ~ 20 per cent on the two-point correlation function and ~ 5 per cent on the void probability function at scale of r ~ 20 h -1 Mpc.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2000

Clustered DNA damages induced in isolated DNA and in human cells by low doses of ionizing radiation

Clustered DNA damages-two or more closely spaced damages (strand breaks, abasic sites, or oxidized bases) on opposing strands-are suspects as critical lesions producing lethal and mutagenic effects of ionizing radiation. However, as a result of the lack of methods for measuring damage clusters induced by ionizing radiation in genomic DNA, neither the frequencies of their production by physiological doses of radiation, nor their repairability, nor their biological effects are known. On the basis of methods that we developed for quantitating damages in large DNAs, we have devised and validated a way of measuring ionizing radiation-induced clustered lesions in genomic DNA, including DNA from human cells. DNA is treated with an endonuclease that induces a single-strand cleavage at an oxidized base or abasic site. If there are two closely spaced damages on opposing strands, such cleavage will reduce the size of the DNA on a nondenaturing gel. We show that ionizing radiation does induce clustered DNA damages containing abasic sites, oxidized purines, or oxidized pyrimidines. Further, the frequency of each of these cluster classes is comparable to that of frank double-strand breaks; among all complex damages induced by ionizing radiation, double-strand breaks are only about 20%, with other clustered damage constituting some 80%. We also show that even low doses (0.1-1 Gy) of high linear energy transfer ionizing radiation induce clustered damages in human cells.

NASA Discipline Radiation Health↗
DOE OSTI2024

Performant automatic differentiation of local coupled cluster theories: Response properties and ab initio molecular dynamics

In this work, we introduce a differentiable implementation of the local natural orbital coupled cluster (LNO-CC) method within the automatic differentiation framework of the PySCFAD package. The implementation is comprehensively tuned for enhanced performance, which enables the calculation of first-order static response properties on medium-sized molecular systems using coupled cluster theory with single, double, and perturbative triple excitations [CCSD(T)]. We evaluate the accuracy of our method by benchmarking it against the canonical CCSD(T) reference for nuclear gradients, dipole moments, and geometry optimizations. In addition, we demonstrate the possibility of property calculations for chemically interesting systems through the computation of bond orders and Mössbauer spectroscopy parameters for a [NiFe]-hydrogenase active site model, along with the simulation of infrared spectra via ab initio LNO-CC molecular dynamics for a protonated water hexamer.

Chemistry↗
DOE OSTI2020

Affinity propagation clustering of full-field, high-spatial-dimensional measurements for robust output-only modal identification: A proof-of-concept study

Determination of the model order is a challenging problem in system identification, especially in output-only or operational modal identification where some modes are weakly excited. While existing methods such as the stabilization diagram method (spectral information) are effective, they do not scale to high-dimensional data, which is usually needed for high-fidelity characterization of structural dynamics and has been made available in the emerging full-field measurement techniques using optical methods such as photogrammetry and laser vibrometers. In this proof-of-concept study we present a new non-parametric, data-driven approach for robust output-only identification of high-spatial-dimensional modal parameters of basic structures by efficiently processing and interactively exploiting the full-field measurement (i.e., very dense spatial measurement points). Specifically, we first over-estimate the system model once, producing a pool of candidate modes associated with their modal frequencies and full-field, high-spatial-dimensional mode shapes. This is accomplished by a data-driven method termed affinity propagation clustering (APC), where the active clusters, which are the active modes in our formulations, emerge from the "message-passing" procedure and does not require a pre-determination of the cluster number (mode or model order). Next, rather than using the spectral information to distinguish the physical and spurious modes in the stabilization diagram method, we exploit and visualize the spatial, full-field mode shape associated with each candidate mode to do so. We conduct extensive experiments on basic structural models with comparisons to a few existing methods. The results indicate that the new method is computationally efficient for identifying high-spatial-dimensional modal parameters, and robust to identify weak modes by exploiting the full-field measurement. We also discuss its applicability and limitations for structures with complex geometry (shapes).

42 ENGINEERING↗
DOE OSTI2020

Finite-temperature coupled cluster: Efficient implementation and application to prototypical systems

In this work, we discuss the theory and implementation of the finite temperature coupled cluster singles and doubles (FT-CCSD) method including the equations necessary for an efficient implementation of response properties. Numerical aspects of the method including the truncation of the orbital space and integration of the amplitude equations are tested on some simple systems, and we provide some guidelines for applying the method in practice. The method is then applied to the 1D Hubbard model, the uniform electron gas (UEG) at warm, dense conditions, and some simple materials. The performance of model systems at high temperatures is encouraging: for the one-dimensional Hubbard model, FT-CCSD provides a qualitatively accurate description of finite-temperature correlation effects even at U = 8, and it allows for the computation of systematically improvable exchange–correlation energies of the warm, dense UEG over a wide range of conditions. We highlight the obstacles that remain in using the method for realistic ab initio calculations on materials

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗