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,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Hilbert space multireference coupled cluster tailored by matrix product states

In the past decade, the quantum chemical version of the density matrix renormalization group method has established itself as the method of choice for strongly correlated molecular systems. However, despite its favorable scaling, in practice, it is not suitable for computations of dynamic correlation. Several approaches to include that in post-DMRG methods exist; in our group, we focused on the tailored coupled cluster (TCC) approach. This method works well in many situations; however, in exactly degenerate cases (with two or more determinants of equal weight), it exhibits a bias toward the reference determinant representing the Fermi vacuum. Although sometimes it is possible to use a compensation scheme to avoid this bias for energy differences, it is certainly a drawback. In order to overcome this bias of the TCC method, we have developed a Hilbert-space multireference version of tailored CC, which can treat several determinants on an equal footing. We have implemented and compared the performance of three Hilbert-space multireference coupled cluster (MRCC) variants—the state universal one and the Brillouin–Wigner and Mukherjee’s state specific ones. We have assessed these approaches on the cyclobutadiene and tetramethyleneethane molecules, which are both diradicals with exactly degenerate determinants at a certain geometry. We have also investigated the sensitivity of the results on the orbital rotation of the highest occupied and lowest unoccupied molecular orbital (HOMO–LUMO) pair, as it is well known that Hilbert-space MRCC methods are not invariant to such transformations.

Chemistry↗
DOE OSTI2024

Jacobian-scaled K-means clustering for physics-informed segmentation of reacting flows

This work introduces Jacobian-scaled K-means (JSK-means) clustering, which is a physicsinformed clustering strategy centered on the K-means framework. The method allows for the injection of underlying physical knowledge into the clustering procedure through a distance function modification: instead of leveraging conventional Euclidean distance vectors, the JSKmeans procedure operates on distance vectors scaled by matrices obtained from dynamical system Jacobians evaluated at the cluster centroids. The goal of this work is to show how the JSKmeans algorithm - without modifying the input dataset - produces clusters that capture regions of dynamical similarity, in that the clusters are redistributed towards high-sensitivity regions in phase space and are described by similarity in the source terms of samples instead of the samples themselves. The algorithm is demonstrated on a complex reacting flow simulation dataset (a channel detonation configuration), where the dynamics in the thermochemical composition space are known through the highly nonlinear and stiff Arrhenius-based chemical source terms. Interpretations of cluster partitions in both physical space and composition space reveal how JSK-means shifts clusters produced by standard K-means towards regions of high chemical sensitivity (e.g., towards regions of peak heat release rate near the detonation reaction zone). Furthermore, the findings presented here illustrate the benefits of utilizing Jacobian-scaled distances in clustering techniques, and the JSK-means method in particular displays promising potential for improving former partition-based modeling strategies in reacting flow (and other multi-physics) applications.

Clustering↗
DOE OSTI2015

Clustering high dimensional data using RIA

Clustering may simply represent a convenient method for organizing a large data set so that it can easily be understood and information can efficiently be retrieved. However, identifying cluster in high dimensionality data sets is a difficult task because of the curse of dimensionality. Another challenge in clustering is some traditional functions cannot capture the pattern dissimilarity among objects. In this article, we used an alternative dissimilarity measurement called Robust Influence Angle (RIA) in the partitioning method. RIA is developed using eigenstructure of the covariance matrix and robust principal component score. We notice that, it can obtain cluster easily and hence avoid the curse of dimensionality. It is also manage to cluster large data sets with mixed numeric and categorical value.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Global Optimization of Adsorbate Covered Supported Cluster Catalysts: The Case of Pt 7 H 10 CH 3 on α-Al 2 O 3

Determining structures of supported sub-nanometer catalytic clusters is important for finding binding sites, active sites, and reaction mechanisms. In realistic conditions, the cluster may be covered with adsorbed reactants and intermediates, which strongly affect the structure of the cluster compared to the low coverage regime. The global optimization of supported adsorbate-covered clusters is not possible with existing methods. In this paper, we present a customized Basin Hopping (BH) approach, with a core (metal cluster) – shell (adsorbates) separation scheme, and perform first applications to the search of low-energy isomers of Pt 7 /α-Al 2 O 3 at high hydrogen coverage (10H). The presence of adsorbates significantly alters the low-energy isomers compared to bare Pt 7 /α-Al 2 O 3 . Molecular dynamics simulations additionally show the different time-scales associated with the mobility of H and Pt atoms. The revised BH method with core-shell separation presents key advantages for the structural exploration of adsorbate-covered supported catalysts, relevant to realistic catalytic conditions.

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

Water interactions with hydrophobic groups: Assessment and recalibration of semiempirical molecular orbital methods

In this work, we present a study of the ability of different semiempirical methods to describe intermolecular interactions in water solution. In particular, we focus on methods based on the Neglect of Diatomic Differential Overlap approximation. Significant improvements of these methods have been reported in the literature in the past years regarding the description of non-covalent interactions. In particular, a broad range of methodologies has been developed to deal with the properties of hydrogen-bonded systems, with varying degrees of success. In contrast, the interactions between water and a molecule containing hydrophobic groups have been little analyzed. Indeed, by considering the potential energy surfaces obtained using different semiempirical Hamiltonians for the intermolecular interactions of model systems, we found that none of the available methods provides an entirely satisfactory description of both hydrophobic and hydrophilic interactions in water. In addition, a vibrational analysis carried out in a model system for these interactions, a methane clathrate cluster, showed that some recent methods cannot be used to carry out studies of vibrational properties. Following a procedure established in our group [M. I. Bernal-Uruchurtu, M. T. C. Martins-Costa, C. Millot, and M. F. Ruiz-López, J. Comput. Chem. 21, 572 (2000); W. Harb, M. I. Bernal-Uruchurtu, and M. F. Ruiz-López, Theor. Chem. Acc. 112, 204 (2004)], we developed new parameters for the core-core interaction terms based on fitting potential energy curves obtained at the MP2 level for our model system. We investigated the transferability of the new parameters to describe a system, having both hydrophilic and hydrophobic groups, interacting with water. We found that only by introducing two different sets of parameters for hydrophilic and hydrophobic hydrogen atom types we are able to match the features of the ab initio calculated properties. Once this assumption is made, a good agreement with the MP2 reference is achieved. The results reported in this work provide therefore a direction for future developments of semiempirical approaches that are still required to investigate chemical processes in biomolecules and in large disordered systems.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS2021

Atacama Cosmology Telescope measurements of a large sample of candidates from the Massive and Distant Clusters of WISE Survey: Sunyaev-Zeldovich effect confirmation of MaDCoWS candidates using ACT

Context. Galaxy clusters are an important tool for cosmology, and their detection and characterization are key goals for current and future surveys. Using data from the Wide-field Infrared Survey Explorer (WISE), the Massive and Distant Clusters of WISE Survey (MaDCoWS) located 2839 significant galaxy overdensities at redshifts 0.7 . z . 1.5, which included extensive follow-up imaging from the Spitzer Space Telescope to determine cluster richnesses. Concurrently, the Atacama Cosmology Telescope (ACT) has produced large area millimeter-wave maps in three frequency bands along with a large catalog of Sunyaev-Zeldovich (SZ)-selected clusters as part of its Data Release 5 (DR5). Aims. We aim to verify and characterize MaDCoWS clusters using measurements of, or limits on, their thermal SZ effect signatures. We also use these detections to establish the scaling relation between SZ mass and the MaDCoWS-defined richness. Methods. Using the maps and cluster catalog from DR5, we explore the scaling between SZ mass and cluster richness. We do this by comparing cataloged detections and extracting individual and stacked SZ signals from the MaDCoWS cluster locations. We use complementary radio survey data from the Very Large Array, submillimeter data from Herschel, and ACT 224 GHz data to assess the impact of contaminating sources on the SZ signals from both ACT and MaDCoWS clusters. We use a hierarchical Bayesian model to fit the mass-richness scaling relation, allowing for clusters to be drawn from two populations: one, a Gaussian centered on the mass-richness relation, and the other, a Gaussian centered on zero SZ signal. Results. We find that MaDCoWS clusters have submillimeter contamination that is consistent with a gray-body spectrum, while the ACT clusters are consistent with no submillimeter emission on average. Additionally, the intrinsic radio intensities of ACT clusters are lower than those of MaDCoWS clusters, even when the ACT clusters are restricted to the same redshift range as the MaDCoWS clusters. We find the best-fit ACT SZ mass versus MaDCoWS richness scaling relation has a slope of p1 = 1.84+0.15 −0.14, where the slope is defined as M ∝ λ p1 15 and λ15 is the richness. We also find that the ACT SZ signals for a significant fraction (∼57%) of the MaDCoWS sample can statistically be described as being drawn from a noise-like distribution, indicating that the candidates are possibly dominated by low-mass and unvirialized systems that are below the mass limit of the ACT sample. Further, we note that a large portion of the optically confirmed ACT clusters located in the same volume of the sky as MaDCoWS are not selected by MaDCoWS, indicating that the MaDCoWS sample is not complete with respect to SZ selection. Finally, we find that the radio loud fraction of MaDCoWS clusters increases with richness, while we find no evidence that the submillimeter emission of the MaDCoWS clusters evolves with richness. Conclusions. We conclude that the original MaDCoWS selection function is not well defined and, as such, reiterate the MaDCoWS collaboration’s recommendation that the sample is suited for probing cluster and galaxy evolution, but not cosmological analyses. We find a best-fit mass-richness relation slope that agrees with the published MaDCoWS preliminary results. Additionally, we find that while the approximate level of infill of the ACT and MaDCoWS cluster SZ signals (1–2%) is subdominant to other sources of uncertainty for current generation experiments, characterizing and removing this bias will be critical for next-generation experiments hoping to constrain cluster masses at the sub-percent level.

large↗
DOE OSTI2015

Structural evolution in the crystallization of rapid cooling silver melt

The structural evolution in a rapid cooling process of silver melt has been investigated at different scales by adopting several analysis methods. The results testify Ostwald’s rule of stages and Frank conjecture upon icosahedron with many specific details. In particular, the cluster-scale analysis by a recent developed method called LSCA (the Largest Standard Cluster Analysis) clarified the complex structural evolution occurred in crystallization: different kinds of local clusters (such as ico-like (ico is the abbreviation of icosahedron), ico-bcc like (bcc, body-centred cubic), bcc, bcc-like structures) in turn have their maximal numbers as temperature decreases. And in a rather wide temperature range the icosahedral short-range order (ISRO) demonstrates a saturated stage (where the amount of ico-like structures keeps stable) that breeds metastable bcc clusters. As the precursor of crystallization, after reaching the maximal number bcc clusters finally decrease, resulting in the final solid being a mixture mainly composed of fcc/hcp (face-centred cubic and hexagonal-closed packed) clusters and to a less degree, bcc clusters. This detailed geometric picture for crystallization of liquid metal is believed to be useful to improve the fundamental understanding of liquid–solid phase transition. - Highlights: • A comprehensive structural analysis is conducted focusing on crystallization. • The involved atoms in our analysis are more than 90% for all samples concerned. • A series of distinct intermediate states are found in crystallization of silver melt. • A novelty icosahedron-saturated state breeds the metastable bcc state.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Measuring metallicities with Hubble space telescope/wide-field camera 3 photometry

We quantified and calibrated the metallicity and temperature sensitivities of colors derived from nine Wide-Field Camera 3 filters on board the Hubble Space Telescope using Dartmouth isochrones and Kurucz atmosphere models. The theoretical isochrone colors were tested and calibrated against observations of five well studied galactic clusters, M92, NGC 6752, NGC 104, NGC 5927, and NGC 6791, all of which have spectroscopically determined metallicities spanning –2.30 < [Fe/H] <+0.4. We found empirical corrections to the Dartmouth isochrone grid for each of the following color-magnitude diagrams (CMDs): (F555W-F814W, F814W), (F336W-F555W, F814W), (F390M-F555W, F814W), and (F390W-F555W, F814W). Using empirical corrections, we tested the accuracy and spread of the photometric metallicities assigned from CMDs and color-color diagrams (which are necessary to break the age-metallicity degeneracy). Testing three color-color diagrams [(F336W-F555W),(F390M-F555W),(F390W-F555W), versus (F555W-F814W)], we found the colors (F390M-F555W) and (F390W-F555W) to be the best suited to measure photometric metallicities. The color (F390W-F555W) requires much less integration time, but generally produces wider metallicity distributions and, at very low metallicity, the metallicity distribution function (MDF) from (F390W-F555W) is ∼60% wider than that from (F390M-F555W). Using the calibrated isochrones, we recovered the overall cluster metallicity to within ∼0.1 dex in [Fe/H] when using CMDs (i.e., when the distance, reddening, and ages are approximately known). The measured MDF from color-color diagrams shows that this method measures metallicities of stellar clusters of unknown age and metallicity with an accuracy of ∼0.2-0.5 dex using F336W-F555W, ∼0.15-0.25 dex using F390M-F555W, and ∼0.2-0.4 dex with F390W-F555W, with the larger uncertainty pertaining to the lowest metallicity range.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Inferring Viral Dynamics in Chronically HCV Infected Patients from the Spatial Distribution of Infected Hepatocytes

Chronic liver infection by hepatitis C virus (HCV) is a major public health concern. Despite partly successful treatment options, several aspects of intrahepatic HCV infection dynamics are still poorly understood, including the preferred mode of viral propagation, as well as the proportion of infected hepatocytes. Answers to these questions have important implications for the development of therapeutic interventions. In this study, we present methods to analyze the spatial distribution of infected hepatocytes obtained by single cell laser capture microdissection from liver biopsy samples of patients chronically infected with HCV. By characterizing the internal structure of clusters of infected cells, we are able to evaluate hypotheses about intrahepatic infection dynamics. We found that individual clusters on biopsy samples range in size from 4-50 infected cells. In addition, the HCV RNA content in a cluster declines from the cell that presumably founded the cluster to cells at the maximal cluster extension. These observations support the idea that HCV infection in the liver is seeded randomly (e.g. from the blood) and then spreads locally. Assuming that the amount of intracellular HCV RNA is a proxy for how long a cell has been infected, we estimate based on models of intracellular HCV RNA replication and accumulation that cells in clusters have been infected on average for less than a week. Further, we do not find a relationship between the cluster size and the estimated cluster expansion time. Lastly, our method represents a novel approach to make inferences about infection dynamics in solid tissues from static spatial data.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1993

Designers' models of the human-computer interface

Understanding design models of the human-computer interface (HCI) may produce two types of benefits. First, interface development often requires input from two different types of experts: human factors specialists and software developers. Given the differences in their backgrounds and roles, human factors specialists and software developers may have different cognitive models of the HCI. Yet, they have to communicate about the interface as part of the design process. If they have different models, their interactions are likely to involve a certain amount of miscommunication. Second, the design process in general is likely to be guided by designers' cognitive models of the HCI, as well as by their knowledge of the user, tasks, and system. Designers do not start with a blank slate; rather they begin with a general model of the object they are designing. The author's approach to a design model of the HCI was to have three groups make judgments of categorical similarity about the components of an interface: human factors specialists with HCI design experience, software developers with HCI design experience, and a baseline group of computer users with no experience in HCI design. The components of the user interface included both display components such as windows, text, and graphics, and user interaction concepts, such as command language, editing, and help. The judgments of the three groups were analyzed using hierarchical cluster analysis and Pathfinder. These methods indicated, respectively, how the groups categorized the concepts, and network representations of the concepts for each group. The Pathfinder analysis provides greater information about local, pairwise relations among concepts, whereas the cluster analysis shows global, categorical relations to a greater extent.

Gillan, Douglas J.↗
DOE OSTI2025

Slimmer Geminals For Accurate F12 Electronic Structure Models

The Slater-type F12 geminal length scales originally tuned for the second-order Mo̷ller-Plesset F12 method are too large for higher-order F12 methods formulated using the SP (diagonal fixed-coefficient spin-adapted) F12 ansatz. The new geminal parameters reported herein reduce the basis set incompleteness errors (BSIEs) of absolute coupled-cluster singles and doubles F12 correlation energies by a significant─and increase with the cardinal number of the basis─margin. The effect of geminal reoptimization is especially pronounced for the cc-pVXZ-F12 basis sets (specifically designed for use with F12 methods) relative to their conventional aug-cc-pVXZ counterparts. The BSIEs of relative energies are less affected, but substantial reductions can be obtained, especially for atomization energies and ionization potentials with the cc-pVXZ-F12 basis sets. The new geminal parameters are therefore recommended for all applications of high-order F12 methods, such as coupled-cluster F12 methods and transcorrelated F12 methods.

Powell, Samuel R. [Virginia Polytechnic Inst. and ↗
DOE OSTI2019

Identification of stable adsorption sites and diffusion paths on nanocluster surfaces: an automated scanning algorithm

The diverse coordination environments on the surfaces of discrete, three-dimensional (3D) nanoclusters contribute significantly to their unique catalytic properties. Identifying the numerous adsorption sites and diffusion paths on these clusters is however tedious and time-consuming, especially for large, asymmetric nanoclusters. Here, we present a simple, automated method for constructing approximate 2D potential energy surfaces for the adsorption of atomic species on the surfaces of 3D nanoclusters with minimal human intervention. These potential energy surfaces fully characterize the important adsorption sites and diffusion paths on the nanocluster surfaces with accuracies similar to current approaches and at comparable computational cost. Our method can treat complex nanoclusters, such as alloy nanoclusters, and accounts for cluster relaxation and adsorbate-induced reconstruction, important for obtaining accurate energetics. Moreover, its highly parallelizable nature is ideal for modern supercomputer architectures. We showcase our method using two clusters: Au 18 and Pt 55 . For Au 18 , diffusion of atomic hydrogen between the most stable sites occurs via non-intuitive paths, underlining the necessity of exploring the complete potential energy surface. By enabling the rapid and unbiased assessment of adsorption and diffusion on large, complex nanoclusters, which are particularly difficult to handle manually, our method will help advance materials discovery and the rational design of catalysts.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Accurate ab initio potential energy surface, thermochemistry, and dynamics of the F{sup −} + CH{sub 3}F S{sub N}2 and proton-abstraction reactions

We develop a full-dimensional global analytical potential energy surface (PES) for the F{sup −} + CH{sub 3}F reaction by fitting about 50 000 energy points obtained by an explicitly correlated composite method based on the second-order Møller–Plesset perturbation-F12 and coupled-cluster singles, doubles, and perturbative triples-F12a methods and the cc-pVnZ-F12 [n = D, T] basis sets. The PES accurately describes the (a) back-side attack Walden inversion mechanism involving the pre- and post-reaction (b) ion-dipole and (c) hydrogen-bonded complexes, the configuration-retaining (d) front-side attack and (e) double-inversion substitution pathways, as well as (f) the proton-abstraction channel. The benchmark quality relative energies of all the important stationary points are computed using the focal-point analysis (FPA) approach considering electron correlation up to coupled-cluster singles, doubles, triples, and perturbative quadruples method, extrapolation to the complete basis set limit, core-valence correlation, and scalar relativistic effects. The FPA classical(adiabatic) barrier heights of (a), (d), and (e) are −0.45(−0.61), 46.07(45.16), and 29.18(26.07) kcal mol{sup −1}, respectively, the dissociation energies of (b) and (c) are 13.81(13.56) and 13.73(13.52) kcal mol{sup −1}, respectively, and the endothermicity of (f) is 42.54(38.11) kcal mol{sup −1}. Quasiclassical trajectory computations of cross sections, scattering (θ) and initial attack (α) angle distributions, as well as translational and internal energy distributions are performed for the F{sup −} + CH{sub 3}F(v = 0) reaction using the new PES. Apart from low collision energies (E{sub coll}), the S{sub N}2 excitation function is nearly constant, the abstraction cross sections rapidly increase with E{sub coll} from a threshold of ∼40 kcal mol{sup −1}, and retention trajectories via double inversion are found above E{sub coll} = ∼ 30 kcal mol{sup −1}, and at E{sub coll} = ∼ 50 kcal mol{sup −1}, the front-side attack cross sections start to increase very rapidly. At low E{sub coll}, the indirect mechanism dominates (mainly isotropic backward-forward symmetric θ distribution and translationally cold products) and significant long-range orientation effects (isotropic α distribution) and barrier recrossings are found. At higher E{sub coll}, the S{sub N}2 reaction mainly proceeds with direct rebound mechanism (backward scattering and hot product translation)

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

Observational constraints of an anisotropic boost due to the projection effects using redMaPPer clusters

Optical clusters identified from red-sequence galaxies suffer from projection effects, where interloper galaxies along the line of sight to a cluster are mistaken as genuine members of the cluster. In the previous study, we found that the projection effects cause the boost on the amplitudes of clustering and lensing on large scale compared to the expected amplitudes in the absence of any projection effects. These boosts are caused by preferential selections of filamentary structure aligned to the line of sight due to distance uncertainties in photometric surveys. We model the projection effects with two simple assumptions and develop a novel method to quantify the size of the boost using cluster-galaxy cross-correlation functions. We validate our method using mock cluster catalogues built from cosmological N-body simulations and find that we can obtain unbiased constraints on the boost parameter with our model. We then apply our analysis on the SDSS redMaPPer clusters and find that the size of the boost is roughly 20 per cent for all the richness bins except the cluster sample with the richness bin λ ∈ [30, 40]. This is the first study to constrain the boost parameter independent from cluster cosmology studies and provides a self-consistency test for the projection effects.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Stochastic Formulation of the Resolution of Identity: Application to Second Order Møller–Plesset Perturbation Theory

A stochastic orbital approach to the resolution of identity (RI) approximation for 4-index electron repulsion integrals (ERIs) is presented. The stochastic RI-ERIs are then applied to second order Møller–Plesset perturbation theory (MP2) utilizing a multiple stochastic orbital approach. The introduction of multiple stochastic orbitals results in an O(N AO 3 ) scaling for both the stochastic RI-ERIs and stochastic RI-MP2, NAO being the number of basis functions. For a range of water clusters we demonstrate that this method exhibits a small prefactor and observed scalings of O(N e 2.4 ) for total energies and O(N e 3.1 ) for forces (N e being the number of correlated electrons), outperforming MP2 for clusters with as few as 21 water molecules.

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

Anion photoelectron spectroscopy and chemical bonding of ThS 2 – and ThSO –

Anion photoelectron spectra of ThSO – and ThS 2 – were recorded using the third (355 nm) harmonic of an Nd-YAG laser; these provided the measured vertical detachment energies of each anion. The experiments are supported by extensive coupled cluster calculations on ThSO, ThSO – , ThS 2 , and ThS 2 – , as well as the oxygen congeners ThO 2 and ThO 2 – . The ab initio calculations, which included complete basis set extrapolations, spin–orbit effects using four-component coupled cluster, and higher-order correlation contributions through CCSDT(Q), yielded an adiabatic electron affinity for ThO 2 – that was within 0.02 eV of the previously determined experimental value. The singly occupied molecular orbital (SOMO) in all three anions corresponds primarily to the 7s orbital on Th. Successive substitution of S for each O in ThO 2 leads to larger electron affinities and smaller bond angles in the neutral molecules, but larger angles in the anions. As demonstrated by Franck–Condon simulations of the spectra using the CCSD(T) spectroscopic constants, substitution of O by S significantly complicates the resulting detachment spectra due to the lower vibrational frequencies in the sulfur species. Altogether the calculated vertical detachment energies are in very good agreement with the experiment. In addition to the adiabatic electron affinities of each species, atomization energies and heats of formation have also been determined via the FPD approach with expected uncertainties of 1–2 kcal/mol.

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

Beyond chemical accuracy in the heavy p-block: The first ionization potentials and electron affinities of Ga-Kr, In-Xe, and Tl-Rn

A relativistic coupled-cluster version of the Feller-Peterson-Dixon (FPD) composite method was previously used to accurately calculate the first ionization potentials (IP) and electron affinities (EAs) of the post-d, p-block elements Ga–Rn. Complete basis set (CBS) extrapolations including outer-core correlation at the CCSD(T) level of theory were combined with contributions from higher order electron correlation up to CCSDTQ, quantum electrodynamic effects (Lamb shift), and spin-orbit (SO) coupling including the Gaunt contribution. Several methods for including SO were investigated, in which all involved the four-component (4c) Dirac-Coulomb (DC) Hamiltonian. Here, the treatment of SO coupling was the contribution that limited the final accuracy of the present results. In the cases where 4c-DC-CCSD(T) could be reliably used for the SO contributions, the final composite IPs and EAs agreed with the available experimental values to within an unsigned average error of just 0.16 and 0.20 kcal/mol, respectively. In all cases the final IPs and EAs were within 1 kcal/mol of the available experimental values, except for the EAs of the group 13 elements (Ga, In, Tl), where the currently accepted experimental values appear to be too large by as much as 4 kcal/mol. The values predicted in this work, which have estimated uncertainties of ±0.5 kcal/mol, are 5.25 (Ga), 7.69 (In), and 7.39 (Tl) kcal/mol. For the EAs of Po and At, which do not have experimental values, the current calculations predict values of 34.2 and 55.8 kcal/mol with estimated uncertainties of ±0.6 and ±0.3 kcal/mol, respectively.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗