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 1,585 records · Page 88

The void spectrum in two-dimensional numerical simulations of gravitational clustering

An algorithm for deriving a spectrum of void sizes from two-dimensional high-resolution numerical simulations of gravitational clustering is tested, and it is verified that it produces the correct results where those results can be anticipated. The method is used to study the growth of voids as clustering proceeds. It is found that the most stable indicator of the characteristic void 'size' in the simulations is the mean fractional area covered by voids of diameter d, in a density field smoothed at its correlation length. Very accurate scaling behavior is found in power-law numerical models as they evolve. Eventually, this scaling breaks down as the nonlinearity reaches larger scales. It is shown that this breakdown is a manifestation of the undesirable effect of boundary conditions on simulations, even with the very large dynamic range possible here. A simple criterion is suggested for deciding when simulations with modest large-scale power may systematically underestimate the frequency of larger voids.

Kauffmann, Guinevere↗

Three-dimensional paganica fault morphology obtained from hypocenter clustering (L'Aquila 2009 seismic sequence, Central Italy)

In seismic modelling, fault planes are normally assumed to be flat due to the lack of data which can constrain fault morphology. However, incorporating 3D fault morphology is important for modelling several phenomena, for example calculating mainshock induced stress changes. Here we utilize a data-analytical method to unveil the 3D rupture morphology of faults using unsupervised clustering techniques applied to earthquake hypocenters in seismic sequences. We apply this method to the 2009 L'Aquila seismic sequence which involved a M W 6.1 mainshock on April 6th. We use a dataset of about 50,000 relocated events, mostly microearthquakes, reaching magnitude of completeness equal to 0.7. Clustering distinguishes the earthquakes as occurring in three main clusters along with other minor fault segments. We then represent the morphology of the main Paganica fault system (responsible for the largest mainshock) using splines. This method shows promise as a step toward robustly and quickly obtaining 3D rupture morphologies where earthquake sequences have been monitored. The 3D model is presented interactively online, and the processing is presented in an interactive Jupyter Notebook (https://bit.ly/2MnCFdj).

58 GEOSCIENCES↗

Unleashed from constrained optimization: quantum computing for quantum chemistry employing generator coordinate inspired method

Hybrid quantum-classical approaches offer potential solutions to quantum chemistry problems, yet they often manifest as constrained optimization problems. Here, we explore the interconnection between constrained optimization and generalized eigenvalue problems through the Unitary Coupled Cluster (UCC) excitation generators. Inspired by the generator coordinate method, we employ these UCC excitation generators to construct non-orthogonal, overcomplete many-body bases, projecting the system Hamiltonian into an effective Hamiltonian, which bypasses issues such as barren plateaus that heuristic numerical minimizers often encountered in standard variational quantum eigensolver (VQE). Diverging from conventional quantum subspace expansion methods, we introduce an adaptive scheme that robustly constructs the many-body basis sets from a pool of the UCC excitation generators. This scheme supports the development of a hierarchical ADAPT quantum-classical strategy, enabling a balanced interplay between subspace expansion and ansatz optimization to address complex, strongly correlated quantum chemical systems cost-effectively, setting the stage for more advanced quantum simulations in chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Selective oxidation of propane to propylene oxide

The invention provides a one pot method for generating propylene oxide, the method having the steps of contacting propane with catalyst clusters no greater than 30 atoms in the presence of oxygen for a time sufficient to directly convert the propane to the propylene oxide. The invented method eliminates the generation of intermediate compounds or intermediate reaction steps.

Vajda, Stefan↗

First Principles Prediction of Structure, Structure Selectivity, and Thermodynamic Stability under Realistic Conditions

Novel materials are often the enabler for new energy technologies. In ab-initio computational materials science, method are developed to predict the behavior of materials starting from the laws of physics, so that properties can be predicted before compounds have to be synthesized and tested. As such, a virtual materials laboratory can be constructed, saving time and money. The objectives of this program were to develop first-principles theory to predict the structure and thermodynamic stability of materials. Since its inception the program focused on the development of the cluster expansion to deal with the increased complexity of complex oxides. This research led to the incorporation of vibrational degrees of freedom in ab-initio thermodynamics, developed methods for multi-component cluster expansions, included the explicit configurational degrees of freedom of localized electrons, developed the formalism for stability in aqueous environments, and culminated in the first ever approach to produce exact ground state predictions of the cluster expansion. Many of these methods have been disseminated to the larger theory community through the Materials Project, pymatgen software, or individual codes. We summarize three of the main accomplishments.

36 MATERIALS SCIENCE↗

Background based Gaussian mixture model lesion segmentation in PET

Purpose: Quantitative {sup 18}F-fluorodeoxyglucose positron emission tomography is limited by the uncertainty in lesion delineation due to poor SNR, low resolution, and partial volume effects, subsequently impacting oncological assessment, treatment planning, and follow-up. The present work develops and validates a segmentation algorithm based on statistical clustering. The introduction of constraints based on background features and contiguity priors is expected to improve robustness vs clinical image characteristics such as lesion dimension, noise, and contrast level. Methods: An eight-class Gaussian mixture model (GMM) clustering algorithm was modified by constraining the mean and variance parameters of four background classes according to the previous analysis of a lesion-free background volume of interest (background modeling). Hence, expectation maximization operated only on the four classes dedicated to lesion detection. To favor the segmentation of connected objects, a further variant was introduced by inserting priors relevant to the classification of neighbors. The algorithm was applied to simulated datasets and acquired phantom data. Feasibility and robustness toward initialization were assessed on a clinical dataset manually contoured by two expert clinicians. Comparisons were performed with respect to a standard eight-class GMM algorithm and to four different state-of-the-art methods in terms of volume error (VE), Dice index, classification error (CE), and Hausdorff distance (HD). Results: The proposed GMM segmentation with background modeling outperformed standard GMM and all the other tested methods. Medians of accuracy indexes were VE <3%, Dice >0.88, CE <0.25, and HD <1.2 in simulations; VE <23%, Dice >0.74, CE <0.43, and HD <1.77 in phantom data. Robustness toward image statistic changes (±15%) was shown by the low index changes: <26% for VE, <17% for Dice, and <15% for CE. Finally, robustness toward the user-dependent volume initialization was demonstrated. The inclusion of the spatial prior improved segmentation accuracy only for lesions surrounded by heterogeneous background: in the relevant simulation subset, the median VE significantly decreased from 13% to 7%. Results on clinical data were found in accordance with simulations, with absolute VE <7%, Dice >0.85, CE <0.30, and HD <0.81. Conclusions: The sole introduction of constraints based on background modeling outperformed standard GMM and the other tested algorithms. Insertion of a spatial prior improved the accuracy for realistic cases of objects in heterogeneous backgrounds. Moreover, robustness against initialization supports the applicability in a clinical setting. In conclusion, application-driven constraints can generally improve the capabilities of GMM and statistical clustering algorithms.

60 APPLIED LIFE SCIENCES↗

Diverse Power Iteration Embeddings and Its Applications

Abstract—Spectral Embedding is one of the most effective dimension reduction algorithms in data mining. However, its computation complexity has to be mitigated in order to apply it for real-world large scale data analysis. Many researches have been focusing on developing approximate spectral embeddings which are more efficient, but meanwhile far less effective. This paper proposes Diverse Power Iteration Embeddings (DPIE), which not only retains the similar efficiency of power iteration methods but also produces a series of diverse and more effective embedding vectors. We test this novel method by applying it to various data mining applications (e.g. clustering, anomaly detection and feature selection) and evaluating their performance improvements. The experimental results show our proposed DPIE is more effective than popular spectral approximation methods, and obtains the similar quality of classic spectral embedding derived from eigen-decompositions. Moreover it is extremely fast on big data applications. For example in terms of clustering result, DPIE achieves as good as 95% of classic spectral clustering on the complex datasets but 4000+ times faster in limited memory environment.

97 MATHEMATICS AND COMPUTING↗

Implementation of density functional embedding theory within the projector-augmented-wave method and applications to semiconductor defect states

We report a new implementation of the density functional embedding theory (DFET) in the VASP code, using the projector-augmented-wave (PAW) formalism. Newly developed algorithms allow us to efficiently perform optimized effective potential optimizations within PAW. The new algorithm generates robust and physically correct embedding potentials, as we verified using several test systems including a covalently bound molecule, a metal surface, and bulk semiconductors. We show that with the resulting embedding potential, embedded cluster models can reproduce the electronic structure of point defects in bulk semiconductors, thereby demonstrating the validity of DFET in semiconductors for the first time. Compared to our previous version, the new implementation of DFET within VASP affords use of all features of VASP (e.g., a systematic PAW library, a wide selection of functionals, a more flexible choice of U correction formalisms, and faster computational speed) with DFET. Furthermore, our results are fairly robust with respect to both plane-wave and Gaussian type orbital basis sets in the embedded cluster calculations. This suggests that the density functional embedding method is potentially an accurate and efficient way to study properties of isolated defects in semiconductors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Aspects of defects in 3d-3d correspondence

In this paper we study supersymmetric co-dimension 2 and 4 defects in the compactification of the 6d (2, 0) theory of type A N-1 on a 3-manifold M . The so-called 3d-3d correspondence is a relation between complexified Chern-Simons theory (with gauge group SL(N,C) ) on M and a 3d N=2 theory T N [M ]. We study this correspondence in the presence of supersymmetric defects, which are knots/links inside the 3-manifold. Our study employs a number of different methods: state-integral models for complex Chern-Simons theory, cluster algebra techniques, domain wall theory T [SU(N )], 5d N=2 SYM, and also supergravity analysis through holography. These methods are complementary and we find agreement between them. In some cases the results lead to highly non-trivial predictions on the partition function. Our discussion includes a general expression for the cluster partition function, which can be used to compute in the presence of maximal and certain class of non-maximal punctures when N > 2. We also highlight the non-Abelian description of the 3d N=2 T N [M ] theory with defect included, when such a description is available. This paper is a companion to our shorter paper, which summarizes our main results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Evaluating gyrochronology on the zero-age-main-sequence: rotation periods in the southern open cluster Blanco 1 from the Kelt-South survey

We report periods for 33 members of Blanco 1 as measured from Kilodegree Extremely Little Telescope-South light curves, the first reported rotation periods for this benchmark zero-age-main-sequence open cluster. The distribution of these stars spans from late-A or early-F dwarfs to mid-K with periods ranging from less than a day to ∼8 days. The rotation period distribution has a morphology similar to the coeval Pleiades cluster, suggesting the universal nature of stellar rotation distributions. Employing two different gyrochronology methods, we find an age of 146{sub −14}{sup +13} Myr for the cluster. Using the same techniques, we infer an age of 134{sub −10}{sup +9} Myr for the Pleiades measured from existing literature rotation periods. These rotation-derived ages agree with independently determined cluster ages based on the lithium depletion boundary technique. Additionally, we evaluate different gyrochronology models and quantify levels of agreement between the models and the Blanco 1/Pleiades rotation period distributions, including incorporating the rotation distributions of clusters at ages up to 1.1 Gyr. We find the Skumanich-like spin-down rate sufficiently describes the rotation evolution of stars hotter than the Sun; however, we find cooler stars rotating faster than predicted by a Skumanich law, suggesting a mass dependence in the efficiency of stellar angular momentum loss rate. Finally, we compare the Blanco 1 and Pleiades rotation period distributions to available nonlinear angular momentum evolution models. We find they require a significant mass dependence on the initial rotation rate of solar-type stars to reproduce the observed range of rotation periods at a given stellar mass and are furthermore unable to predict the observed over-density of stars along the upper envelope of the clusters' rotation distributions.

79 ASTRONOMY AND ASTROPHYSICS↗

A new self-adaptive reconstruction method to identify defects through Wigner–Seitz approach

A new self-adaptive reconstruction method based on local atomic structure at any given molecular dynamics (MD) step has been developed in this article. The method can be used in Wigner–Seitz defect analysis approach to correctly and efficiently explore the information of both point defects and complex defect clusters (e.g. dislocation loops and voids) formed after a displacement cascade where the cascade interacts with grain boundaries and/or dislocations. The algorithm and validation are provided in detail. Results for identification of radiation defects during and after cascades interacting with a dislocation network show that the new method can well recognize all simple and complex defects and defect clusters. Thus, this new method provides a totally new way to explore the density and size of radiation defects at atomic scale after complex MD evolution processes, providing correct information to understand and predict radiation damage in materials through atomic simulations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Zwartia hollandica gen. nov., sp. nov., Jezberella montanilacus gen. nov., sp. nov. and Sheuella amnicola gen. nov., comb. nov., representing the environmental GKS98 (betIII) cluster

We present two strains affiliated with the GKS98 cluster. This phylogenetically defined cluster is representing abundant, mainly uncultured freshwater bacteria, which were observed by many cultivation-independent studies on the diversity of bacteria in various freshwater lakes and streams. Bacteria affiliated with the GKS98 cluster were detected by cultivation-independent methods in freshwater systems located in Europe, Asia, Africa and the Americas. The two strains, LF4-65 T (=CCUG 56422 T =DSM 107630 T ) and MWH-P2sevCIIIb T (=CCUG 56420 T =DSM 107629 T ), are aerobic chemoorganotrophs, both with genome sizes of 3.2 Mbp and G+C values of 52.4 and 51.0 mol%, respectively. Phylogenomic analyses based on concatenated amino acid sequences of 120 proteins suggest an affiliation of the two strains with the family Alcaligenaceae and revealed Orrella amnicola and Orrella marina (= Algicoccus marinus ) as being the closest related, previously described species. However, the calculated phylogenomic trees clearly suggest that the current genus Orrella represents a polyphyletic taxon. Based on the branching order in the phylogenomic trees, as well as the revealed phylogenetic distances and chemotaxonomic traits, we propose to establish the new genus Zwartia gen. nov. and the new species Z. hollandica sp. nov. to harbour strain LF4-65 T and the new genus Jezberella gen. nov. and the new species J. montanilacus sp. nov. to harbour strain MWH-P2sevCIIIb T . Furthermore, we propose the reclassification of the species Orrella amnicola in the new genus Sheuella gen. nov. The new genera Zwartia, Jezberella and Sheuella together represent taxonomically the GKS98 cluster.

Microbiology↗

A new “gold standard”: Perturbative triples corrections in unitary coupled cluster theory and prospects for quantum computing

A major difficulty in quantum simulation is the adequate treatment of a large collection of entangled particles, synonymous with electron correlation in electronic structure theory, with coupled cluster (CC) theory being the leading framework for dealing with this problem. Augmenting computationally affordable low-rank approximations in CC theory with a perturbative account of higher-rank excitations is a tractable and effective way of accounting for the missing electron correlation in those approximations. This is perhaps best exemplified by the “gold standard” CCSD(T) method, which bolsters the baseline CCSD with the effects of triple excitations using considerations from many-body perturbation theory (MBPT). Despite this established success, such a synergy between MBPT and the unitary analog of CC theory (UCC) has not been explored. In this work, we propose a similar approach wherein converged UCCSD amplitudes are leveraged to evaluate energy corrections associated with triple excitations, leading to the UCCSD[T] method. In terms of quantum computing, this correction represents an entirely classical post-processing step that improves the energy estimate by accounting for triple excitation effects without necessitating new quantum algorithm developments or increasing demand for quantum resources. The rationale behind this choice is shown to be rigorous by studying the properties of finite-order UCC energy functionals, and our efforts do not support the addition of the fifth-order contributions as in the (T) correction. We assess the performance of these approaches on a collection of small molecules and demonstrate the benefits of harnessing the inherent synergy between MBPT and UCC theories.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A new calibration of the extragalactic distance scale using Cepheids and RR Lyrae stars

Visual absolute magnitudes of classical Cepheids, metal-poor RR Lyrae stars and short-period type II Cepheids have been determined with very high precision by combining a large number of old and new astrophysical data. Five independent methods (four observational and one theoretical) have been successfully used: (1) secular and statistical parallaxes; (2) moving-group parallaxes; (3) cluster main-sequence fitting; (4) the Baade-Wesselink method and its modifications; and (5) light-curve and velocity-curve fitting (the theoretical method). With these five adopted methods, the zero point of the galactic and extragalactic distance scale can be set on a relatively firm and self-consistent basis. Classical Cepheids and RR Lyrae stars now provide essentially identical distances to nearby galaxies: the distance modulus to the Large Magellanic Cloud is 18.5 and to the Small Magellanic Cloud is 18.8, both with an uncertainty of 0.1. The distance to the center of our Galaxy is 8.6 plus or minus 0.5 kiloparsecs. The major uncertainty in these values lies in the correction for interstellar extinction.

Stothers, R. B.↗

THE DYNAMICS OF MERGING CLUSTERS: A MONTE CARLO SOLUTION APPLIED TO THE BULLET AND MUSKET BALL CLUSTERS

Merging galaxy clusters have become one of the most important probes of dark matter, providing evidence for dark matter over modified gravity and even constraints on the dark matter self-interaction cross-section. To properly constrain the dark matter cross-section it is necessary to understand the dynamics of the merger, as the inferred cross-section is a function of both the velocity of the collision and the observed time since collision. While the best understanding of merging system dynamics comes from N-body simulations, these are computationally intensive and often explore only a limited volume of the merger phase space allowed by observed parameter uncertainty. Simple analytic models exist but the assumptions of these methods invalidate their results near the collision time, plus error propagation of the highly correlated merger parameters is unfeasible. To address these weaknesses I develop a Monte Carlo method to discern the properties of dissociative mergers and propagate the uncertainty of the measured cluster parameters in an accurate and Bayesian manner. I introduce this method, verify it against an existing hydrodynamic N-body simulation, and apply it to two known dissociative mergers: 1ES 0657-558 (Bullet Cluster) and DLSCL J0916.2+2951 (Musket Ball Cluster). I find that this method surpasses existing analytic models-providing accurate (10% level) dynamic parameter and uncertainty estimates throughout the merger history. This, coupled with minimal required a priori information (subcluster mass, redshift, and projected separation) and relatively fast computation ({approx}6 CPU hours), makes this method ideal for large samples of dissociative merging clusters.

79 ASTRONOMY AND ASTROPHYSICS↗

Atomic Cluster Aggregates in Nucleation of Solid Alumina Inclusion in the Aluminum Deoxidation for Liquid Iron

Nucleation plays a decisive role in determining the structure and size distribution of a solid inclusion particle in liquid iron. The literature does not yet provide a clear picture about the pathway leading to solid alumina inclusion from liquid iron. The main confusion surrounds the little-known structures and thermodynamic properties of metastable alumina phase (MAP) that appear in the nucleation of alumina inclusion. In this work, it is suggested that MAP could appear in the Al-deoxidation reaction, based on the summaries and analysis in terms of various equilibrium experimental results and thermodynamic modeling of the Al deoxidation in liquid iron since the 1950s. The structures and thermodynamic properties of a species of MAP, alumina atomic cluster aggregates (AACAs), are calculated by density functional theory (DFT) methods. The thermodynamics on formation and transformation of AACAs shows that AACAs are in equilibrium with Al and O with various contents measured in the Al-deoxidation equilibrium experiments for liquid iron. It is suggested that the residual Al and O, which cannot nucleate, appear in the form of AACAs in Al-deoxidized liquid iron. It can be concluded that AACAs are the metastable intermediates and structural units for alumina inclusion nucleating from liquid iron. Thus, the behaviors of AACAs are important to the control of the nucleation rate of alumina inclusions in the Al deoxidation for liquid iron.

36 MATERIALS SCIENCE↗

A BAYESIAN METHOD FOR THE INTERCALIBRATION OF SPECTRA IN REVERBERATION MAPPING

Flux calibration of spectra in reverberation mapping (RM) is most often performed by assuming the flux constancy of some specified narrow emission lines, which stem from an extended region that is sometimes partially spatially resolved, in contrast to the point-like broad-line region and the central continuum source. The inhomogeneous aperture geometries used among different observation sets in a joint monitoring campaign introduce systematic deviations to the fluxes of broad lines and central continuum, and intercalibration over these data sets is required. As an improvement to the previous empirical correction performed by comparing the (nearly) contemporaneous observation points, we describe a feasible Bayesian method that obviates the need for (nearly) contemporaneous observations, naturally incorporates physical models of flux variations, and fully takes into account the measurement errors. In particular, it fits all the data sets simultaneously regardless of samplings and makes use of all of the information in the data sets. A Markov Chain Monte Carlo implementation is employed to recover the parameters and uncertainties for intercalibration. Application to the RM data sets of NGC 5548 with joint monitoring shows the high fidelity of our method.

79 ASTRONOMY AND ASTROPHYSICS↗

Angular-momentum projection in coupled-cluster theory: Structure of 34 Mg

Single- reference coupled-cluster theory is an accurate and affordable computational method for the nuclear many-body problem. For open-shell nuclei, the reference state typically breaks rotational invariance and angular momentum must be restored as a good quantum number. We perform angular-momentum projection after variation and employ the disentangled coupled-cluster formalism and a Hermitian approach. We compare our results with benchmarks for 8 Be and 20 Ne using a two-nucleon interaction from chiral effective field theory and for pf-shell nuclei within the traditional shell model. We compute the rotational band in the exotic nucleus 34 Mg and find agreement with data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗