Search NASA⌕ Search

SEARCH · Search NASA

Results for “coupled cluster method”

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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Two-fluid simulations of galaxy formation

We investigate the formation of galaxies and larger structure with a simulation modeling two gravitationally coupled fluids representing dark matter and baryons. The baryon gas dynamics are calculated with a smoothed particle hydrodynamics (SPH) method, and the physics modeled includes thermal pressure, shock heating, and radiative cooling. We simulate a 16 Mpc periodic cube with 64(exp 3) particles in each fluid and 10% baryon mass fraction. We confirm, for the first time experimentally, disk formation as a natural consequence of hierarchical clustering in a large-scale cosmological environment. The majority of isolated galaxies exhibit centrifugally supported disks. A power-law relation between cold baryonic mass and maximum rotation velocity is found, M varies as nu(sub rot)(exp alpha) with alpha = 2.5 after correcting for differential numerical resolution. Both the spatial and velocity distributions of the simulated galaxies are biased with respect to the dark matter. A counts-in-cells analysis indicates that an unphysical degree of merging in the central cluster is likely responsible for the antibias signal in the correlation function. A robust, scale-dependent velocity bias is measured. The ratio of galaxy to dark matter pairwise velocity dispersions on a scale of 1 Mpc is 0.7. The amplitude is only mildly dependent on redshift or mass cutoff and scales with separation as r(exp 0.2). The degree to which these results depend on numerical parameters is discussed. Mass resolution plays a key role in controlling the resulting fraction of cold, dense baryons. The mass fraction associated with galaxies decreases by a factor of approximately greater than 3 when the mass per particle is increased by a factor 8. Photoionization and energy input from supernova will have to be included to determine more carefully the fraction of highly dissipated material and the characteristics of the stellar component of galaxies.

Evrard, August E.↗
DOE OSTI2014

Optical and electronic properties of mixed Ag-Au tetramer cations

We present experimental and theoretical studies of the optical response of mixed Ag {sub n} Au {sub 4−n}{sup +} (n=1–3) clusters in the photon energy range ℏω = 1.9–3.5 eV. Absorption spectra are recorded by a newly built longitudinal molecular beam depletion spectroscopy apparatus providing lower limits to absolute photodissociation cross sections. The experimental data are compared to optical response calculations in the framework of long-range corrected time-dependent density functional theory with initial cluster geometries obtained by the unbiased Birmingham Cluster Genetic Algorithm coupled with density functional theory. Experiments and excited state calculations shed light on the structural and electronic properties of the mixed Ag-Au tetramer cations.

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

SU-E-T-523: On the Radiobiological Impact of Lateral Scatter in Proton Beams

Introduction: In proton therapy, justified concern has been voiced with respect to an increased efficiency in cell kill at the distal end of the Bragg peak. This coupled with range uncertainty is a counter indication to use the Bragg peak to define the border of a treated volume with a critical organ. An alternative is to use the lateral edge of the proton beam, obtaining more robust plans. We investigate the spectral and biological effects of the lateral scatter . Methods: A general purpose Monte Carlo simulation engine (MCNPX 2.7c) installed on a Scientific Linux cluster, calculated the dose deposition spectrum of protons, knock on electrons and generated neutrons for a proton beam with maximal kinetic energy of 200MeV. Around the beam at different positions in the beam direction the spectrum is calculated in concentric rings of thickness 1cm. The deposited dose is converted to a double strand break map using an analytical expression.based on micro dosimetric calculations using a phenomenological Monte Carlo code (MCDS). A strict version of RBE is defined as the ratio of generation of double strand breaks in the different modalities. To generate the reference a Varian linac was modelled in MCNPX and the generated electron dose deposition spectrum was used . Results: On a pristine point source 200MeV beam the RBE before the Bragg peak was of the order of 1.1, increasing to 1.7 right behind the Bragg peak. When using a physically more realistic beam of 10cm diameter the effect was smaller. Both the lateral dose and RBE increased with increasing beam depth, generating a dose deposition with mixed biological effect. Conclusions: The dose deposition in proton beams need to be carefully examined because the biological effect will be different depending on the treatment geometry. Deeply penetrating proton beams generate more biologically effective lateral scatter.

07 ISOTOPE AND RADIATION SOURCES↗
NASA NTRS1999

An Intercomparison Between Radar Reflectivity and the IR Cloud Classification Technique for the TOGA-COARE Area

Satellite infrared (IR) and visible (VIS) images from the Tropical Ocean Global Atmosphere - Coupled Ocean Atmosphere Response Experiment (TOGA-COARE) experiment are investigated through the use of Clustering Analysis. The clusters are obtained from the values of IR and VIS counts and the local variance for both channels. The clustering procedure is based on the standardized histogram of each variable obtained from 179 pairs of images. A new approach to classify high clouds using only IR and the clustering technique is proposed. This method allows the separation of the enhanced convection in two main classes: convective tops, more closely related to the most active core of the storm, and convective systems, which produce regions of merged, thick anvil clouds. The resulting classification of different portions of cloudiness is compared to the radar reflectivity field for intensive events. Convective Systems and Convective Tops are followed during their life cycle using the IR clustering method. The areal coverage of precipitation and features related to convective and stratiform rain is obtained from the radar for each stage of the evolving Mesoscale Convective Systems (MCS). In order to compare the IR clustering method with a simple threshold technique, two IR thresholds (Tir) were used to identify different portions of cloudiness, Tir=240K which roughly defines the extent of all cloudiness associated with the MCS, and Tir=220K which indicates the presence of deep convection. It is shown that the IR clustering technique can be used as a simple alternative to identify the actual portion of convective and stratiform rainfall.

Carvalho, L. M. V.↗
DOE OSTI2020

Scalable Heterogeneous Execution of a Coupled-Cluster Model with Perturbative Triples

The CCSD(T) coupled-cluster model with perturbative triples is considered a gold standard for computational modeling of the correlated behavior of electrons in molecular systems. A fundamental constraint is the relatively small global-memory capacity in GPUs compared to the main-memory capacity on host nodes, necessitating relatively smaller tile sizes for high-dimensional tensor contractions in NWChem's GPU-accelerated implementation of the CCSD(T) method. A coordinated redesign is described to address this limitation and associated data movement overheads, including a novel fused GPU kernel for a set of tensor contractions, along with inter-node communication optimization and data caching. The new implementation of GPU-accelerated CCSD(T) improves overall performance by 3.4x. Finally, we discuss the trade-offs in using this fused algorithm on current and future supercomputing platforms.

Kim, Jinsung↗
DOE OSTI2014

Nuclear three-body problem in the complex energy plane: Complex-scaling Slater method

The physics of open quantum systems is an interdisciplinary area of research. The nuclear “openness” manifests itself through the presence of the many-body continuum representing various decay, scattering, and reaction channels. As the radioactive nuclear beam experimentation extends the known nuclear landscape toward the particle drip lines, the coupling to the continuum space becomes exceedingly more important. Of particular interest are weakly bound and unbound nuclear states appearing around particle thresholds. Theories of such nuclei must take into account their open quantum nature. To describe open quantum systems, we introduce a complex-scaling (CS) approach in the Slater basis. We benchmark it with the complex-energy Gamow shell model (GSM) by studying energies and wave functions of the bound and unbound states of the two-neutron halo nucleus 6 He viewed as an $α + n + n$ cluster system. Both CS and GSM approaches are applied to a translationally invariant Hamiltonian with the two-body interaction approximated by the finite-range central Minnesota force. In the CS approach, we use the Slater basis, which exhibits the correct asymptotic behavior at large distances. To extract particle densities from the back-rotated CS solutions, we apply the Tikhonov regularization procedure, which minimizes the ultraviolet numerical noise. In this paper, we show that the CS-Slater method is both accurate and efficient. Its equivalence to the GSM approach has been demonstrated numerically for both energies and wave functions of 6 He. One important technical aspect of our calculation was to fully retrieve the correct asymptotic behavior of a resonance state from the complex-scaled (square-integrable) wave function. While standard applications of the inverse complex transformation to the complex-rotated solution provide unstable results, the stabilization method fully reproduces the GSM benchmark. We also propose a method to determine the smoothing parameter of the Tikhonov regularization. The combined suite of CS-Slater and GSM techniques has many attractive features when applied to nuclear problems involving weakly bound and unbound states. While both methods can describe energies, total widths, and wave functions of nuclear states, the CS-Slater method—if it can be applied—can provide additional information about partial energy widths associated with individual thresholds.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Distributed non-negative matrix factorization with determination of the number of latent features

The holistic analysis and understanding of the latent (that is, not directly observable) variables and patterns buried in large datasets is crucial for data-driven science, decision making and emergency response. Such exploratory analyses require devising unsupervised learning methods for data mining and extraction of the latent features, and non-negative matrix factorization (NMF) is one of the prominent such methods. NMF is based on compute-intense non-convex constrained minimization, which, for large datasets requires fast and distributed algorithms. However, current parallel implementations of NMF fail to estimate the number of latent features. In practice, identifying these features is both difficult and significant for pattern recognition and latent feature analysis, especially for large dense matrices. Here, we introduce a distributed NMF algorithm coupled with distributed custom clustering followed by a stability analysis on dense data, which we call DnMFk, to determine the number of latent variables. The results on synthetic data and the classical Swimmer data set demonstrate the accuracy of model determination while scaling nearly linearly across multiple processors for large data. Further, we employ DnMFk to determine the number of hidden features from a terabyte matrix.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Assessing Electronic Structure Methods for Long-Range Three-Body Dispersion Interactions: Analysis and Calculations on Well-Separated Metal Atom Trimers

Three-body dispersion interactions are much weaker than their two-body counterpart. However, their importance grows quickly as the number of interacting monomers rises. To explore the numerical performance of correlation methods for long-range three-body dispersion, we performed calculations on eight very simple dispersion-dominated model metal trimers: Na 3 , Mg 3 , Zn 3 , Cd 3 , Hg 3 , Cu 3 , Ag 3 , and Au 3 . One encouraging aspect is that relatively small basis sets of augmented triple-ζ size appear to be adequate for three-body dispersion in the long-range. Coupled cluster calculations were performed at high levels to assess MP3, CCSD, CCSD(T), empirical density functional theory dispersion (D3), and the many-body dispersion (MBD) approach. We found that the accuracy of CCSD(T) was generally significantly lower than for two-body interactions, with errors sometimes reaching 20% in the investigated systems, while CCSD and particularly MP3 were generally more erratic. MBD is found to perform better than D3 at large distances, whereas the opposite is true at shorter distances. When computing reference numbers for three-body dispersion, care should be taken to appropriately represent the effect of the connected triple excitations, which are significant in most cases and incompletely approximated by CCSD(T).

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

Gaunt and Breit two-electron contributions to mean-field transformations and fine structure splitting

Materials utilized by novel energy systems are often studied using weakly correlated mean-field theories. However, if these systems incorporate heavy elements, relativistic effects must be included. Therefore, a Kramers unrestricted coupled cluster with singles and doubles excitation formalism within a molecular mean-field exact two-component framework (X2C mmf ) using a four-component Dirac–Hartree–Fock (DHF) reference state is presented. The exact X2C mmf transformed normal-order Hamiltonian incorporates all one-electron and two-electron (2e) contributions from the Coulomb, Gaunt, and Breit operators and is used with the equation of motion method to calculate the excitation energies of the alkali group of elements. Using this framework, the effects of 2e Gaunt and Breit integrals are studied. Results demonstrate growing contributions from these integrals to the generated X2C mmf mean-fields and electronic fine structure calculations with increasing atomic number. Overall, this paper outlines the method, its effect within the X2C mmf approach, and lays the foundation for future theoretical development of relativistic calculations within this framework.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Site occupancy trend of Co in Ni{sub 2}MnIn: Ab initio approach

The trend of site occupation of Co at Ni sites of Ni{sub 2}MnIn system is studied in austenitic phase having L2{sub 1} structure by ab initio density functional theory (DFT) calculation. The Co atoms prefer to be at Ni sites rather than Mn site and are ferromagetically coupled with Ni and Mn. The ground state has tetragonal structure for Ni{sub 1.5}Co{sub 0.5}MnIn and Ni{sub 1.25}Co{sub 0.75}MnIn. The Co tends to form cluster.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

NMR shielding and spin–rotation constants of {sup 175}LuX (X = {sup 19}F, {sup 35}Cl, {sup 79}Br, {sup 127}I) molecules

This presentation demonstrates the relativistic effects on the spin-rotation constants, absolute nuclear magnetic resonance (NMR) shielding constants and shielding spans of {sup 175}LuX (X = {sup 19}F, {sup 35}Cl, {sup 79}Br, {sup 127}I) molecules. The results are obtained from calculations performed using density functional theory (non-relativistic and four-component relativistic) and coupled-cluster calculations. The spin-rotation constants are compared with available experimental values. In most of the molecules studied, relativistic effects make an order of magnitude difference on the NMR absolute shielding constants.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Interpretation of the photoelectron, ultraviolet, and vacuum ultraviolet photoabsorption spectra of bromobenzene by ab initio configuration interaction and DFT computations

New photoelectron, ultraviolet (UV), and vacuum UV (VUV) spectra have been obtained for bromobenzene by synchrotron study with higher sensitivity and resolution than previous work. This, together with use of ab initio calculations with both configuration interaction and time dependent density functional theoretical methods, has led to major advances in interpretation. The VUV spectrum has led to identification of a considerable number of Rydberg states for the first time. The Franck-Condon (FC) analyses including both hot and cold bands lead to identification of the vibrational structure of both ionic and electronically excited states including two Rydberg states. The UV onset has been interpreted in some detail, and an interpretation based on the superposition of FC and Herzberg-Teller contributions has been performed. In a similar way, the 6 eV absorption band which is poorly resolved is analysed in terms of the presence of two ππ* states of {sup 1}A{sub 1} (higher oscillator strength) and {sup 1}B{sub 2} (lower oscillator strength) symmetries, respectively. The detailed analysis of the vibrational structure of the 2{sup 2}B{sub 1} ionic state is particularly challenging, and the best interpretation is based on equation-of-motion-coupled cluster with singles and doubles computations. A number of equilibrium structures of the ionic and singlet excited states show that the molecular structure is less subject to variation than corresponding studies for iodobenzene. The equilibrium structures of the 3b{sub 1}3s and 6b{sub 2}3s (valence shell numbering) Rydberg states have been obtained and compared with the corresponding ionic limit structures.

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

27 Al NMR diffusometry of Al 13 Keggin nanoclusters

Although nanometer-sized aluminum hydroxide clusters (i.e., ϵ-Al 13 , [Al 13 O 4 (OH) 24 (H 2 O) 12 ] 7+ ) command a central role in aluminum ion speciation and transformations between minerals, measurement of their translational diffusion is often limited to indirect methods. Here, in this work, 27 Al pulsed field gradient stimulated echo nuclear magnetic resonance (PFGSTE NMR) spectroscopy has been applied to the AlO 4 core of the ϵ-Al 13 cluster with complementary theoretical simulations of the diffusion coefficient and corresponding hydrodynamic radii from a boundary element-based calculation. The tetrahedral AlO 4 center of the ϵ-Al 13 cluster is symmetric and exhibits only weak quadrupolar coupling, which results in favorable T 1 and T 2 27 Al NMR relaxation coefficients for 27 Al PFGSTE NMR studies. Stokes–Einstein relationship was used to relate the 27 Al diffusion coefficient of the ϵ-Al 13 cluster to the hydrodynamic radius for comparison with theoretical simulations, dynamic light scattering from literature, and previously published 1 H PFGSTE NMR studies of chelated Keggin clusters. This first-of-its-kind observation proves that 27 Al PFGSTE NMR diffusometry can probe symmetric Al environments in polynuclear clusters of greater molecular weight than previously considered.

27Al↗
DOE OSTI2014

Multi-level coupled cluster theory

We present a general formalism where different levels of coupled cluster theory can be applied to different parts of the molecular system. The system is partitioned into subsystems by Cholesky decomposition of the one-electron Hartree-Fock density matrix. In this way the system can be divided across chemical bonds without discontinuities arising. The coupled cluster wave function is defined in terms of cluster operators for each part and these are determined from a set of coupled equations. The total wave function fulfills the Pauli-principle across all borders and levels of electron correlation. We develop the associated response theory for this multi-level coupled cluster theory and present proof of principle applications. The formalism is an essential tool in order to obtain size-intensive complexity in the calculation of local molecular properties.

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

Actinide chemistry using singlet-paired coupled cluster and its combinations with density functionals

Singlet-paired coupled cluster doubles (CCD0) is a simplification of CCD that relinquishes a fraction of dynamic correlation in order to be able to describe static correlation. Combinations of CCD0 with density functionals that recover specifically the dynamic correlation missing in the former have also been developed recently. Here, we assess the accuracy of CCD0 and CCD0+DFT (and variants of these using Brueckner orbitals) as compared to well-established quantum chemical methods for describing ground-state properties of singlet actinide molecules. Finally, the f 0 actinyl series (UO$2+\atop{2}$, NpO$3+\atop{2}$, PuO$4+\atop{2}$), the isoelectronic NUN, and thorium (ThO, ThO 2+ ) and nobelium (NoO, NoO 2 ) oxides are studied.

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

Combining the GW formalism with the polarizable continuum model: A state-specific non-equilibrium approach

We have implemented the polarizable continuum model within the framework of the many-body Green’s function GW formalism for the calculation of electron addition and removal energies in solution. The present formalism includes both ground-state and non-equilibrium polarization effects. In addition, the polarization energies are state-specific, allowing to obtain the bath-induced renormalisation energy of all occupied and virtual energy levels. Our implementation is validated by comparisons with ΔSCF calculations performed at both the density functional theory and coupled-cluster single and double levels for solvated nucleobases. The present study opens the way to GW and Bethe-Salpeter calculations in disordered condensed phases of interest in organic optoelectronics, wet chemistry, and biology.

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

Sequential Bond Dissociation Energies of Th + (CO) x , x = 3–6: Guided Ion Beam Collision-Induced Dissociation and Quantum Computational Studies

Collision-induced dissociation (CID) of [Th,xC,xO] + , x = 3–6, with Xe is performed using a guided ion beam tandem mass spectrometer (GIBMS). Products are formed exclusively by the loss of CO ligands. Analyses of the kinetic energy-dependent CID product cross sections yield bond dissociation energies (BDEs) of (CO) x–1 Th + –CO at 0 K as 1.09 ± 0.05, 0.82 ± 0.07, 0.63 ± 0.05, and 0.70 ± 0.05 eV, respectively. Different structures of [Th,xC,xO] + were explored using various electronic structure methods, and BDEs for CO ligand loss from precursor [Th,xC,xO] + complexes were computed. Both experimental and theoretical results corroborate that the structures of [Th,xC,xO] + , x = 3–6, formed experimentally are homoleptic thorium cation carbonyl complexes, Th + (CO) x . The nonmonotonic trend in experimental BDEs is reproduced theoretically, although ambiguities in the spin states of the x = 4–6 complexes (doublet or quartet) remain. BDEs calculated at the coupled cluster with single, double, and perturbative triple excitations (CCSD(T))/cc-pVXZ//B3LYP/cc-PVXZ (X = T and Q) level and a complete basis set (CBS) extrapolation agree reasonably well with the experimental values for all complexes. Thorium oxide ketenylidene carbonyl cations, OTh + CCO(CO) y , y = 1–4, were calculated to be the most stable structures of [Th,xC,xO] + , x = 3–6, respectively; however, these are not observed in our experiment. Potential energy profiles (PEPs) having either quartet or doublet spin calculated at the B3LYP/cc-pVQZ level suggest that the failure to observe OTh + CCO(CO) y , y = 1–4, is the result of a barrier corresponding to the C–C bond formation, making the formation of OTh + CCO(CO) y inaccessible kinetically under the present experimental conditions.

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

Computational simulation and interpretation of the low-lying excited electronic states and electronic spectrum of thioanisole

Three singlet states, namely a closed-shell ground state and two excited states with 1 ππ* and 1 nσ* character, have been suggested to be responsible for the radiationless decay or photochemical reaction of photoexcited thioanisole. The correct interpretation of the electronic spectrum is critical for understanding the character of these low-lying excited states, but the experimental spectrum is yet to be fully interpreted. In the work reported here, we investigated the nature of those three states and a fourth singlet state of thioanisole using electronic structure calculations by multireference perturbation theory, by completely-renormalized equation-of-motion coupled cluster theory with single and double excitations and noniterative inclusion of connected triples (CR-EOM-CCSD(T)), and by linear-response time-dependent density functional theory (TDDFT). We clarified the assignment of the electronic spectrum by simulating it using a normal-mode sampling approach combined with TDDFT in the Tamm–Dancoff approximation (TDA). The understanding of the electronic states and of the accuracy of the electronic structure methods lays the foundation of our future work of constructing potential energy surfaces.

74 ATOMIC AND MOLECULAR PHYSICS↗