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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
NASA NTRS1998

Singlet Excited States of Cl and Br Molecules: New Theories Applied to the -XO and -XO2 (X=C1 and Br) Chromophores

Electronic excitation energies are determined using single-reference based theories derived from response equations involving perturbation theory and coupled-cluster theory. These methods are applied to the singlet manifold of excited electronic states of the HClO, HBrO, HOClO, HOBrO, HClO2, and HBrO2 molecules. The reliability of the various perturbation theory approaches is assessed by comparison to the linear-response singles and doubles coupled-cluster (LRCCSD) method. The excitation energies for the Y-XO compounds are compared and contrasted for Y=H and HO, and X=Cl and Br. A similar comparison is performed for the H-XO2 compounds.

Lee, Timothy J.↗
DOE OSTI2019

Excited states via coupled cluster theory without equation-of-motion methods: Seeking higher roots with application to doubly excited states and double core hole states

Here, in this work, we revisited the idea of using the coupled-cluster (CC) ground state formalism to target excited states. Our main focus was targeting doubly excited states and double core hole states. Typical equation-of-motion (EOM) approaches for obtaining these states struggle without higher-order excitations than doubles. We showed that by using a non-Aufbau determinant optimized via the maximum overlap method, the CC ground state solver can target higher energy states. Furthermore, just with singles and doubles (i.e., CCSD), we demonstrated that the accuracy of ΔCCSD and ΔCCSD(T) (triples) far surpasses that of EOM-CCSD for doubly excited states. The accuracy of ΔCCSD(T) is nearly exact for doubly excited states considered in this work. For double core hole states, we used an improved ansatz for greater numerical stability by freezing core hole orbitals. The improved methods, core valence separation (CVS)-ΔCCSD and CVS-ΔCCSD(T), were applied to the calculation of the double ionization potential of small molecules. Even without relativistic corrections, we observed qualitatively accurate results with CVS-ΔCCSD and CVS-ΔCCSD(T). Remaining challenges in ΔCC include the description of open-shell singlet excited states with the single-reference CC ground state formalism as well as excited states with genuine multireference character. The tools and intuition developed in this work may serve as a stepping stone toward directly targeting arbitrary excited states using ground state CC methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Non-local corrections to the typical medium theory of Anderson localization

We use the recently developed finite cluster typical medium approach to study the Anderson localization transition in three dimensions. Applying our method to the box and binary alloy disorder distributions, we find a fast convergence with the cluster size. We demonstrate the importance of the typical medium environment and the non-local spatial correlations for the proper characterization of the localization transition. As the cluster size increases, our typical medium cluster method recovers the correct critical disorder strength for the transition. Our findings highlight the importance of the non-local cluster corrections for capturing the localization behavior of the mobility edge trajectories. Our results demonstrate that the typical medium cluster approach developed here provides a consistent and systematic description of the Anderson localization transition in the framework of the effective medium embedding schemes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Benchmarking third-order cluster perturbation theory for electronically excited states

In this study, we investigate the reliability of cluster perturbation (CP) theory applied to the calculation of electronically excited states through a comprehensive benchmark. In CP theory, perturbative corrections are added to the properties of a parent excitation space, which converge toward the properties of a target excitation space. For the CPS(D-n) model, perturbative corrections through order n are added to the coupled cluster singles (CCS) excitation energies to target the coupled cluster singles and doubles (CCSD) excitation energies. Through a comparative analysis of excitation energy calculations across a diverse set of molecules and wavefunction methods, we present a comprehensive evaluation of the accuracy of the third-order CPS(D) model, CPS(D-3), in calculating excitation energies. Further, our findings demonstrate that CPS(D-3) is a reliable alternative to established methods, particularly CCSD, while systematically overestimating the excitation energies compared to high-level coupled cluster methods such as CC3. These results highlight the strengths and limitations of CPS(D-3), as well as the promising directions for its future development.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

CASTELO: clustered atom subtypes aided lead optimization—a combined machine learning and molecular modeling method

Background: Drug discovery is a multi-stage process that comprises two costly major steps: pre-clinical research and clinical trials. Among its stages, lead optimization easily consumes more than half of the pre-clinical budget. We propose a combined machine learning and molecular modeling approach that partially automates lead optimization workflow in silico, providing suggestions for modification hot spots. Results: The initial data collection is achieved with physics-based molecular dynamics simulation. Contact matrices are calculated as the preliminary features extracted from the simulations. To take advantage of the temporal information from the simulations, we enhanced contact matrices data with temporal dynamism representation, which are then modeled with unsupervised convolutional variational autoencoder (CVAE). Finally, conventional and CVAE-based clustering methods are compared with metrics to rank the submolecular structures and propose potential candidates for lead optimization. Conclusion: With no need for extensive structure-activity data, our method provides new hints for drug modification hotspots which can be used to improve drug potency and reduce the lead optimization time. It can potentially become a valuable tool for medicinal chemists.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2007

Monitoring by Use of Clusters of Sensor-Data Vectors

The inductive monitoring system (IMS) is a system of computer hardware and software for automated monitoring of the performance, operational condition, physical integrity, and other aspects of the health of a complex engineering system (e.g., an industrial process line or a spacecraft). The input to the IMS consists of streams of digitized readings from sensors in the monitored system. The IMS determines the type and amount of any deviation of the monitored system from a nominal or normal ( healthy ) condition on the basis of a comparison between (1) vectors constructed from the incoming sensor data and (2) corresponding vectors in a database of nominal or normal behavior. The term inductive reflects the use of a process reminiscent of traditional mathematical induction to learn about normal operation and build the nominal-condition database. The IMS offers two major advantages over prior computational monitoring systems: The computational burden of the IMS is significantly smaller, and there is no need for abnormal-condition sensor data for training the IMS to recognize abnormal conditions. The figure schematically depicts the relationships among the computational processes effected by the IMS. Training sensor data are gathered during normal operation of the monitored system, detailed computational simulation of operation of the monitored system, or both. The training data are formed into vectors that are used to generate the database. The vectors in the database are clustered into regions that represent normal or nominal operation. Once the database has been generated, the IMS compares the vectors of incoming sensor data with vectors representative of the clusters. The monitored system is deemed to be operating normally or abnormally, depending on whether the vector of incoming sensor data is or is not, respectively, sufficiently close to one of the clusters. For this purpose, a distance between two vectors is calculated by a suitable metric (e.g., Euclidean distance) and "sufficiently close" signifies lying at a distance less than a specified threshold value. It must be emphasized that although the IMS is intended to detect off-nominal or abnormal performance or health, it is not necessarily capable of performing a thorough or detailed diagnosis. Limited diagnostic information may be available under some circumstances. For example, the distance of a vector of incoming sensor data from the nearest cluster could serve as an indication of the severity of a malfunction. The identity of the nearest cluster may be a clue as to the identity of the malfunctioning component or subsystem. It is possible to decrease the IMS computation time by use of a combination of cluster-indexing and -retrieval methods. For example, in one method, the distances between each cluster and two or more reference vectors can be used for the purpose of indexing and retrieval. The clusters are sorted into a list according to these distance values, typically in ascending order of distance. When a set of input data arrives and is to be tested, the data are first arranged as an ordered set (that is, a vector). The distances from the input vector to the reference points are computed. The search of clusters from the list can then be limited to those clusters lying within a certain distance range from the input vector; the computation time is reduced by not searching the clusters at a greater distance.

Iverson, David L.↗
DOE OSTI2022

Methods to Calculate Electronic Excited-State Dynamics for Molecules on Large Metal Clusters with Many States: Ensuring Fast Overlap Calculations and a Robust Choice of Phase

Here, we present an efficient set of methods for propagating excited-state dynamics involving a large number of configuration interaction singles (CIS) or Tamm-Dancoff approximation (TDA) single-reference excited states. Specifically, (i) following Head-Gordon et al., we implement an exact evaluation of the overlap of singly-excited CIS/TDA electronic states at different nuclear geometries using a biorthogonal basis and (ii) we employ a unified protocol for choosing the correct phase for each adiabat at each geometry. For many-electron systems, the combination of these techniques significantly reduces the computational cost of integrating the electronic Schrodinger equation and imposes minimal overhead on top of the underlying electronic structure calculation. As a demonstration, we calculate the electronic excited-state dynamics for a hydrogen molecule scattering off a silver metal cluster, focusing on high-lying excited states, where many electrons can be excited collectively and crossings are plentiful. Interestingly, we find that the high-lying, plasmon-like collective excitation spectrum changes with nuclear dynamics, highlighting the need to simulate non-adiabatic nuclear dynamics and plasmonic excitations simultaneously. In the future, the combination of methods presented here should help theorists build a mechanistic understanding of plasmon-assisted charge transfer and excitation energy relaxation processes near a nanoparticle or metal surface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Gas-Phase Stability of Large Lanthanide:Ligand Clusters Evaluated Using Collision-Induced Dissociation

Introduction In the reprocessing of f-elements present in used nuclear fuels, a variety of diglycolamides (DGA’s) are used as extractants for actinide partitioning. In particular, the Actinide-Lanthanide Separation (ALSEP) process typically utilizes either the N,N,N’,N’-tetraoctyl diglycolamide (TODGA) or N,N,N',N'-tetra-2-ethylhexyl diglycolamide (T2EHDGA) extractant ligands following the partitioning of uranium and plutonium from used nuclear fuel. To better understand fundamental interactions in these processes, covalent bonding of several f-elements with diglycolamides, primarily TODGA, is investigated in the gas phase using nanospray ionization and a quadrupole time-of-flight mass spectrometer. Further, analysis of the identity and relative strength of the cluster is enabled by MS2 isolation and collision induced dissociation. Methods Metal ion cluster analysis was completed using a Bruker (Billerica, MA, USA) mircOTOF-Q II quadrupole time-of-flight mass spectrometer with a CaptiveSpray nanospray ion source. Detection was accomplished using positive ionization mode. Metal: ligand solutions were assembled as 30 µM europium nitrate, samarium nitrate, cerium nitrate, or holmium nitrate and 3 µM DGA in acetonitrile or a 50:50 mixture of acetonitrile: isopropanol. Preliminary data The samarium cluster experiments yielded clusters with a samarium:TODGA ratio of up to 1:7 able to be isolated and evidence of greater ratios present in the mass spectrum. This is surprising, as metal clusters are not expected to have a coordination space able to accommodate this many TODGA ligands, due to its size and tridenticity. Collisional activation of [Sm(TODGA)3]3+ suggested loss of a TODGA radical cation, in addition to ligand fragmentation. In contrast, activation of clusters with higher Sm:TODGA ratios resulted in loss of entire ligands, with no evidence of fragmentation. A lower collision energy was required to remove ligands as the number of bound TODGAs increased, suggesting that in larger clusters, ligands are more delicately complexed to the metal. In addition, several clusters were observed with the composition [Sm(NO3)x(TODGA)n x]+3 x. With a single nitrate ion, clusters with up to six TODGAs were able to be isolated. In a similar pattern to the samarium clusters containing only TODGA, less collision energy was required to eliminate one or more TODGAs with increasing size. Clusters with composition [Sm(NO3)(TODGA)n-1]2+ appeared in lower abundance and were more collisionally stable than [Sm(TODGA)n]3+ clusters. With two nitrate ions, only clusters with a single TODGA were able to be isolated. Analogous europium experiments resulted in similar clusters. Ratios of up to 1:7 Eu:TODGA and clusters with one nitrate and up to five TODGAs were isolated. In clusters with two nitrate ions, only one or two TODGAs were observed to be bound. Similar to samarium, MS2 experiments with the Eu clusters suggested that larger clusters required less collision energy to eliminate TODGA. Europium clusters with the composition [Eu(NO3)(TODGA)n-1]2+ were observed in greater abundance and with greater stability than the equivalent cluster with the composition [Eu(TODGA)n]3+. Novel Aspect These are the first reported Ln:TODGA clusters, allowing us to begin to investigate intrinsic complexation of lanthanides with process-relevant ligands.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Investigating Covalent Bonding in f-elements using Gas-phase Ion Chemistry

Introduction In the reprocessing of f-elements present in used nuclear fuels, a variety of diglycolamides (DGA’s) are used as extractants for actinide partitioning. In particular, the Actinide-Lanthanide Separation (ALSEP) process typically utilizes either the N,N,N’,N’-tetraoctyl diglycolamide (TODGA) or N,N,N',N'-tetra-2-ethylhexyl diglycolamide (T2EHDGA) extractant ligands following the partitioning of uranium and plutonium from used nuclear fuel. To better understand fundamental interactions in these processes, covalent bonding of several f-elements with diglycolamides, primarily TODGA, is investigated in the gas phase using nanospray ionization and a quadrupole time-of-flight mass spectrometer. Further, analysis of the identity and relative strength of the cluster is enabled by MS2 isolation and collision induced dissociation. Methods Metal ion cluster analysis was completed using a Bruker mircOTOF-Q II quadrupole time-of-flight mass spectrometer equipped with a CaptiveSpray nanospray ion source. Metal:ligand solutions were prepared as 30 µM europium nitrate, samarium nitrate, cerium nitrate, or holmium nitrate and 3 µM DGA in acetonitrile or a 50:50 mixture of acetonitrile: isopropanol. Cluster mass spectra and collision-induced dissociation experiments were conducted in positive mode. Preliminary data To examine the patterns and relative strength of lanthanide: DGA interactions, MS2 experiments were completed with each lanthanide species listed above. Preliminary analyses of samarium and europium TODGA clusters suggest several combinations of TODGA and nitrate forming. The samarium cluster experiments yielded Sm(TODGA)x clusters with a samarium:TODGA ratio of up to 1:7 able to be isolated and evidence of greater ratios present in the mass spectrum. This is surprising, as metal clusters are not expected to have a coordination space able to accommodate this many ligands as large as TODGA. MS2 experiments show that, at higher ratios and with sufficient collision energy, entire TODGA ligands are removed instead of being fragmented. These experiments show that a lower collision energy is required to remove ligands as the number of bound TODGA’s increases, suggesting that in larger clusters, ligands are more delicately complexed to the metal. In addition to Sm(TODGA)x, several clusters were observed with nitrate ions bound to the metal in addition to TODGA. With a single nitrate ion, clusters with up to six TODGA’s were able to be isolated. In a similar pattern to the samarium clusters with only TODGA, less collision energy is required to eliminate one or more TODGA’s with increasing size. MS2 experiments suggest clusters with one nitrate appear to be of an equivalent or greater stability to clusters which replace the nitrate with a TODGA, as more collision energy is required to remove a TODGA ligand. These species with one nitrate are also in a higher abundance than the equivalent TODGA only cluster. With two nitrate ions, only clusters with a single TODGA were able to be isolated. Analogous europium experiments resulted in very similar clusters. Ratios of up to 1:7 Eu:TODGA were able to be isolated, and clusters with one nitrate and up to five TODGAs were isolated. In clusters with two nitrate ions, one or two TODGA’s could also be bound to the metal. MS2 experiments suggested, similarly to samarium, that larger clusters required less collision energy to eliminate TODGA. Europium clusters with one nitrate are in greater abundance and are stronger than the equivalent cluster which replaces the nitrate with TODGA. Similar analysis with cerium and holmium is ongoing, as well as analysis with other DGA ligands to compare relative strengths of the lanthanide metals with various extractant ligands.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Experiments to validate Thermodynamics and Transport models of Strongly Coupled Dusty Plasma Matter (Final Technical Report for DE-SC0023416)

The goal of this two-year grant is to provide access to the PI to dusty plasma experimental facilities at the DOE-funded Collaborative Research Facility Magnetized Plasma Research Laboratory, Auburn University to become a user of that facility, to obtain experimental data to support another ongoing grant DE-SC0021146 (an Early Career Award to the PI that is focused on modeling of dusty plasma thermodynamics and transport processes), to generate experimental data for funding proposals, and to provide exposure to University of Memphis students to advanced experimental techniques. The following technical accomplishments were made: 1. Development of a novel Bidirectional Electrode Control Arms Assembly (BECAA) for producing perfect 2D grain layers for complex plasma experimentation. BECAA uses movable electrode arms to tilt or move the electrode in a RF discharge from outside the chamber, allowing for the manipulation of grain clouds without needing to change the plasma parameters or gas pressure. This work addresses a longstanding gap in the literature for a method to produce clusters of selectable number of grains and that are perfectly two dimensional as opposed to being only quasi-2D. 2. Experimental investigation of the structural properties of finite-N clusters with N=2 to 50. Individual particle behavior in clusters could vary from grain to grain and this study measured systematically produced clusters for two different grain sizes. Analysis (funded by another grant DE-SC0021146) is currently underway to quantify the differences between grains that are found on the surface vs. the interior of clusters, the shell structure, and the decay of correlations in position, velocity, and kinetic energy. 3. An experimental method to measure the structural entropy of clusters was developed by observing the self-induced structural transitions between various possible arrangements. In a series of heating and cooling cycles, the number of times each possible arrangement was attained was experimentally observed and used to compute the probability of existence of that arrangement, and subsequently the configurational entropy of the cluster. Analysis (funded by another grant DE-SC0021146) is currently underway to produce the entropy of clusters as a function of the number of grains and use the same to compute thermodynamic state variables for 2D complex plasma/grain clusters. 4. A preliminary experimental study of multibody collisions between N grains (N=2 – 10) was conducted. The clusters were produced using the BECAA technique and velocities were imparted to the grains using manipulation laser pulses. Analysis (funded by another grant DE-SC0021146) is currently underway to develop a theoretical framework to describe multibody collisions analogous to classical two-body interactions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Local Embedding and Effective Downfolding in the Auxiliary-Field Quantum Monte Carlo Method

A local embedding and effective downfolding scheme has been developed and implemented in the auxiliary-field quantum Monte Carlo (AFQMC) method. A local cluster in which electrons are fully correlated is defined, and the frozen orbital method is used on the remainder of the system to construct an effective Hamiltonian, which operates within the local cluster. Local embedding, which involves only the occupied sector, has previously been employed in the context of Co/graphene. Here, the methodology is extended in order to allow for effective downfolding of the virtual sector, thus allowing for significant reduction in the computational effort required for AFQMC calculations. The system size, which can be feasibly treated with AFQMC, is therefore greatly extended as only a single local cluster is explicitly correlated at the AFQMC level of theory. The approximation is controlled by the separate choice of the spatial size of the active occupied region (Ro) and of the active virtual region (Rv). In conclusion, the systematic dependence of the AFQMC energy on Ro and Rv is investigated, and it is found that relative AFQMC energies of physical and chemical interest converge rapidly to the full AFQMC treatment (i.e., using no embedding or downfolding).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Assessment of DFT functionals for a minimal nitrogenase [Fe(SH) 4 H]- model employing state-of-the-art ab initio methods

We have designed a [Fe(SH)4H]- model with the fifth proton binding either to Fe or S. We show that the energy difference between these two isomers (ΔE) is hard to estimate with quantum-mechanical (QM) methods. For example, different density functional theory (DFT) methods give ΔE estimates that vary by almost 140 kJ/mol, mainly depending on the amount of exact Hartree–Fock included (0%–54%). The model is so small that it can be treated by many high-level QM methods, including coupled-cluster (CC) and multiconfigurational perturbation theory approaches. With extrapolated CC series (up to fully connected coupled-cluster calculations with singles, doubles, and triples) and semistochastic heat-bath configuration interaction methods, we obtain results that seem to be converged to full configuration interaction results within 5 kJ/mol. Our best result for ΔE is 101 kJ/mol. With this reference, we show that M06 and B3LYP-D3 give the best results among 35 DFT methods tested for this system. Brueckner doubles coupled cluster with perturbaitve triples seems to be the most accurate coupled-cluster approach with approximate triples. CCSD(T) with Kohn–Sham orbitals gives results within 4–11 kJ/mol of the extrapolated CC results, depending on the DFT method. Single-reference CC calculations seem to be reasonably accurate (giving an error of ~5 kJ/mol compared to multireference methods), even if the D1 diagnostic is quite high (0.25) for one of the two isomers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Vertical detachment energies of ammonia cluster anions using self-interaction-corrected methods

Systems with weakly bound extra electrons impose great challenges to semilocal density functional approximations (DFAs), which suffer from self-interaction errors. Small ammonia clusters are one such example of weakly bound anions where the extra electron is weakly bound. We applied two self-interaction correction (SIC) schemes, viz., the well-known Perdew–Zunger and the recently developed locally scaled SIC (LSIC) with the local spin density approximation (LSDA), Perdew–Burke–Ernzerhof (PBE) generalized gradient approximation (GGA), and the SCAN meta-GGA functionals to calculate the vertical detachment energies (VDEs) of small ammonia cluster anions (NH3)n–. Our results show that the LSIC significantly reduces the errors in calculations of VDE with LSDA and PBE-GGA functionals leading to better agreement with the reference values calculated with coupled cluster singles and doubles with perturbative triples [CCSD(T)]. Accurate prediction of VDE as an absolute of the highest occupied molecular orbital (HOMO) is challenging for DFAs. Our results show that VDEs estimated from the negative of HOMO eigenvalues with the LSIC-LSDA and Perdew–Zunger SIC-PBE are within 11 meV of the reference CCSD(T) results. Here, the LSIC method performs consistently well for the VDE estimates, from both the total energy differences and the absolute HOMO eigenvalues.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Interpreting the Operando XANES of Surface-Supported Subnanometer Clusters: When Fluxionality, Oxidation State, and Size Effect Fight.

X-ray absorption near edge structure (XANES) spectroscopy is widely used for operando catalyst characterization. We show that, for highly fluxional supported nanoclusters, the customary extraction of the oxidation state of the metal from the XANES data by fitting to the bulk standards is highly questionable. The XANES signatures as well as the apparent oxidation state for such clusters arise from a complex combination of many factors, and not only from the chemical composition in reaction conditions (e.g., oxygen content in oxidizing atmosphere). The thermally accessible isomerization and population of several structurally distinct cluster forms, cluster-support interaction, and intrinsic size effects all impact the metal oxidation state and XANES signal. We demonstrate this on copper oxide clusters with different compositions, Cu4Ox(x = 2-5) and Cu5Oy (y = 3, 5), deposited on amorphous alumina and ultrananocrystalline diamond, for which we computed the XANES spectra and compare the results to the experiment. We show in addition that fitting the experimental spectrum to calculated spectra of supported clusters can, in contrast, provide good agreement and insight into the spectrum-composition-structure relation. Experimental XANES interpreted using the proposed fitting scheme shows the partial reduction of Cu oxide clusters at rising temperatures, and pinpoints the specific stoichiometries that dominate in the ensemble of cluster states as the temperature changes.

Zandkarimi, Borna↗
DOE OSTI2019

Revisiting surface core-level shifts for ionic compounds

The shifts of core-level binding energies can provide powerful information about the electronic structure of a material. Understanding the physical origin of these shifts for catalytically relevant oxides may provide important insight into their properties. This requires reliable theoretical methods which are able to relate the binding energy shifts to the electronic structure. In order to establish such a methodology, the CaO(100) surface to bulk core-level binding energy shifts have been studied with Hartree-Fock and density-functional theory methods using both cluster and periodic slab models. The shifts obtained from the different theoretical methods are compared with each other and with data from synchrotron-based x-ray photoelectron spectroscopy (XPS) measurements. With a common approximation for the slab model treatment of XPS, the predicted binding energy shifts are seriously in error. In this work, the origin of the error is identified as arising from a flawed treatment of the surface atom binding energies, and a method for correcting the failure is presented.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1997

Sulfur in Cometary Dust

The computer-intensive project consisted of the analysis and synthesis of existing data on composition of comet Halley dust particles. The main objective was to obtain a complete inventory of sulfur containing compounds in the comet Halley dust by building upon the existing classification of organic and inorganic compounds and applying a variety of statistical techniques for cluster and cross-correlational analyses. A student hired for this project wrote and tested the software to perform cluster analysis. The following tasks were carried out: (1) selecting the data from existing database for the proposed project; (2) finding access to a standard library of statistical routines for cluster analysis; (3) reformatting the data as necessary for input into the library routines; (4) performing cluster analysis and constructing hierarchical cluster trees using three methods to define the proximity of clusters; (5) presenting the output results in different formats to facilitate the interpretation of the obtained cluster trees; (6) selecting groups of data points common for all three trees as stable clusters. We have also considered the chemistry of sulfur in inorganic compounds.

Fomenkova, M. N.↗