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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Vibrations in small Mg clusters

The equilibrium geometries, binding energies, and harmonic frequencies of Mg3 and Mg4 have been determined using large atomic natural orbital basis sets in conjunction with high levels of electron correlation. The correlation treatments comprise multireference configuration interaction (MRCI), singles and doubles coupled-cluster (CCSD) theory and the CCSD(T) extension that includes a perturbational estimate of connected triple excitations. As with previous studies of small Be clusters, the CCSD(T) method is found to reproduce the MRCI results with a remarkable degree of accuracy and at a fraction of the computational cost. Using the CCSD(T) method, full quartic force fields for Mg3 and Mg4 have been determined and anharmonic analyses have been performed using second-order perturbation theory. Vibration-rotation interaction constants and centrifugal distortion constants have also been determined.

Lee, Timothy J.↗
DOE OSTI2015

New XDM-corrected potential energy surfaces for Ar–NO(X{sup 2}Π): A comparison with CCSD(T) calculations and experiments

We report new potential energy surfaces for the ground state Ar–NO(X{sup 2}Π) van der Waals system calculated using the unrestricted Hartree-Fock (UHF) method with the addition of the Becke-Roussel correlation functional and exchange-hole dipole moment dispersion correction (XDM). We compare UHFBR-XDM surfaces and those previously reported by Alexander from coupled cluster CCSD(T) calculations [J. Chem. Phys. 111, 7426 (1999)]. The bound states of Ar–NO have been investigated with these new UHFBR-XDM surfaces, including relative energy-level spacing, adiabatic bender states and wave functions, and spectroscopic data. These results have been found to be in good agreement with calculations based on the CCSD(T) PESs. These new PESs are used to investigate the inelastic scattering of NO(X) by Ar. Full close-coupling integral cross sections at collision energies of 442 cm{sup −1}, 1774 cm{sup −1} and differential cross sections at collision energy of 530 cm{sup −1} were determined for transitions out of the lowest NO(X) rotational level (j = ω = 1/2,f). These cross sections are in good agreement with those calculated with CCSD(T) and accordingly in good agreement with the most recent initial and final state resolved experimental data. The UHFBR-XDM scheme yields high-quality potential surfaces with computational cost comparable to the Hartree-Fock method and our results may serve as a benchmark for application of this scheme to collisions between larger molecules.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2016

Towards quantifying the role of exact exchange in the prediction hydrogen bond spin-spin coupling constants involving fluorine

The effect of a fraction of Hartree-Fock exchange on the calculated spin-spin coupling constants involving fluorine through a hydrogen bond is analyzed in detail. Coupling constants calculated using wavefunction methods are revisited in order to get high-level calculations using the same basis set. Accurate MCSCF results are obtained using an additive approach. These constants and their contributions are used as a reference for density functional calculations. Within the density functional theory, the Hartree-Fock exchange functional is split in short- and long-range using a modified version of the Coulomb-attenuating method with the SLYP functional as well as with the original B3LYP. Results support the difficulties for calculating hydrogen bond coupling constants using density functional methods when fluorine nuclei are involved. Coupling constants are very sensitive to the Hartree-Fock exchange and it seems that, contrary to other properties, it is important to include this exchange for short-range interactions. Best functionals are tested in two different groups of complexes: those related with anionic clusters of type [F(HF){sub n}]{sup −} and those formed by difluoroacetylene and either one or two hydrogen fluoride molecules.

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

A comparative analysis of clustering algorithms: O{sub 2} migration in truncated hemoglobin I from transition networks

The ligand migration network for O{sub 2}–diffusion in truncated Hemoglobin N is analyzed based on three different clustering schemes. For coordinate-based clustering, the conventional k–means and the kinetics-based Markov Clustering (MCL) methods are employed, whereas the locally scaled diffusion map (LSDMap) method is a collective-variable-based approach. It is found that all three methods agree well in their geometrical definition of the most important docking site, and all experimentally known docking sites are recovered by all three methods. Also, for most of the states, their population coincides quite favourably, whereas the kinetics of and between the states differs. One of the major differences between k–means and MCL clustering on the one hand and LSDMap on the other is that the latter finds one large primary cluster containing the Xe1a, IS1, and ENT states. This is related to the fact that the motion within the state occurs on similar time scales, whereas structurally the state is found to be quite diverse. In agreement with previous explicit atomistic simulations, the Xe3 pocket is found to be a highly dynamical site which points to its potential role as a hub in the network. This is also highlighted in the fact that LSDMap cannot identify this state. First passage time distributions from MCL clusterings using a one- (ligand-position) and two-dimensional (ligand-position and protein-structure) descriptor suggest that ligand- and protein-motions are coupled. The benefits and drawbacks of the three methods are discussed in a comparative fashion and highlight that depending on the questions at hand the best-performing method for a particular data set may differ.

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

Nuclear charge radii of Na isotopes: Interplay of atomic and nuclear theory

The accuracy of atomic theory calculations limits the extraction of nuclear charge radii from isotope shift measurements of odd-proton nuclei. For Na isotopes, though precise spectroscopic measurements have existed for more than half a century, calculations by different methods offer a wide range of values. Here, we present accurate atomic calculations to reliably extract the Na charge radii. By combining experimental matter radii with nuclear coupled-cluster calculations based on nucleon-nucleon and three-nucleon forces, we constrain the parameters obtained from the atomic calculations. Furthermore, this study guides atomic theory and highlights the importance of using accurate atomic and nuclear computations in our understanding of the size of light nuclei.

20 ≤ A ≤ 38↗
DOE OSTI2019

Approaching coupled cluster accuracy with a general-purpose neural network potential through transfer learning

Computational modeling of chemical and biological systems at atomic resolution is a crucial tool in the chemist’s toolset. The use of computer simulations requires a balance between cost and accuracy: quantum-mechanical methods provide high accuracy but are computationally expensive and scale poorly to large systems, while classical force fields are cheap and scalable, but lack transferability to new systems. Machine learning can be used to achieve the best of both approaches. Here we train a general-purpose neural network potential (ANI-1ccx) that approaches CCSD(T)/CBS accuracy on benchmarks for reaction thermochemistry, isomerization, and drug-like molecular torsions. This is achieved by training a network to DFT data then using transfer learning techniques to retrain on a dataset of gold standard QM calculations (CCSD(T)/CBS) that optimally spans chemical space. The resulting potential is broadly applicable to materials science, biology, and chemistry, and billions of times faster than CCSD(T)/CBS calculations.

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

First-Principles Statistical Mechanics Study of Magnetic Fluctuations and Order–Disorder in the Spinel LiNi 0.5 Mn 1.5 O 4 Cathode

While significant magnetic interactions exist in lithium transition metal oxides, commonly used as Li-ion cathodes, the interplay between magnetic couplings, disorder, and redox processes remains poorly understood. In this work, we focus on the high-voltage spinel LiNi 0.5 Mn 1.5 O 4 (LNMO) cathode as a model system on which to apply a computational framework that uses first principles-based statistical mechanics methods to predict the finite temperature magnetic properties of materials and provide insights into the complex interplay between magnetic and chemical degrees of freedom. Density functional theory calculations on multiple distinct Ni–Mn orderings within the LNMO system, including the ordered ground-state structure (space group P4332), reveal a preference for a ferrimagnetic arrangement of the Ni and Mn sublattices due to strong antiferromagnetic superexchange interactions between neighboring Mn 4+ and Ni 2+ ions and ferromagnetic Mn–Mn and Ni–Ni couplings, as revealed by magnetic cluster expansions. These results are consistent with qualitative predictions using the Goodenough-Kanamori-Anderson rules. Simulations of the finite temperature magnetic properties of LNMO are conducted using Metropolis Monte Carlo. We find that a “semiclassical” Monte Carlo sampling method based on the Heisenberg Hamiltonian accurately predicts experimental magnetic transition temperatures observed in magnetometry measurements. This study highlights the importance of a robust computational toolkit that accurately captures the complex chemomagnetic interactions and predicts finite temperature magnetic behavior to help analyze experimental magnetic and magnetic resonance spectroscopy data acquired ex situ and operando.

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

Molecular structure determination: Equilibrium structure of pyrimidine (m-C 4 H 4 N 2 ) from rotational spectroscopy (reSE) and high-level ab initio calculation (re) agree within the uncertainty of experimental measurement

The pure rotational spectrum of pyrimidine (m-C 4 H 4 N 2 ), the meta-substituted dinitrogen analog of benzene, has been studied in the millimeter-wave region from 235 GHz to 360 GHz. The rotational spectrum of the ground vibrational state has been assigned and fit to yield accurate rotational and distortion constants. Over 1700 distinct transitions were identified for the normal isotopologue in its ground vibrational state and least-squares fit to a partial sextic S-reduced Hamiltonian. Transitions for all four singly substituted 13 C and 15 N isotopologues were observed at natural abundance and were likewise fit. Deuterium-enriched samples of pyrimidine were synthesized, giving access to all eleven possible deuterium-substituted isotopologues, ten of which were previously unreported. Experimental values of rotational constants and computed values of vibration–rotation interaction constants and electron-mass corrections were used to determine semi-experimental equilibrium structures (r e SE ) of pyrimidine. The r e SE structure obtained using coupled-cluster with single, double, and perturbative triple excitations [CCSD(T)] corrections shows exceptional agreement with the re structure computed at the CCSD(T)/cc-pCV5Z level (≤0.0002 Å in bond distance and ≤0.03° in bond angle). Of the various computational methods examined, CCSD(T)/cc-pCV5Z is the only method for which the computed value of each geometric parameter lies within the statistical experimental uncertainty (2σ) of the corresponding semi-experimental coordinate. The exceptionally high accuracy and precision of the structure determination is a consequence of the large number of isotopologues measured, the precision and extent of the experimental frequency measurements, and the sophisticated theoretical treatment of the effects of vibration–rotation coupling and electron mass. Taken together, these demanding experimental and computational studies establish the capabilities of modern structural analysis for a prototypical monocyclic aromatic compound.

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

From electronic structure to model application of key reactions for gasoline/alcohol combustion: Hydrogen-atom abstractions by $CH_3\dot{O}$ radical

H-atom abstraction by methoxy radical ($CH_3\dot{O}$) plays an important role in capturing the kinetics of reactions between gasoline components and alcohols. This study focuses on determining the reaction rates and thermodynamic properties of methoxy radical reactions with five gasoline fuel components: n-heptane, iso-octane, 1-hexene, cyclopentane and toluene. Electronic structure calculations were per-formed for all the stationary points with M06-2X/6 -311 ++ g(d,p) method. G3 composite method with atomization method is used for determining Δ f H 0 of all the closed shell and radical species, using which the necessary thermodynamic data of all the species was determined. Coupled cluster theory QCISD(T)/cc-pVXZ (where X = D and T) and Møller-Plesset perturbation theory MP2/cc-pVXZ (where X = D, T and Q) were used to calculate single point energies. Subsequently, rate constants for all hydrogen atom ab-straction channels have been performed using conventional transition state theory with unsymmetric tunneling corrections. A systematic comparison of rates for abstraction from different sites within the same species and same site from different species is done in order to get insights into this reaction class. Here, the computed thermodynamic properties and rate constants were incorporated into a recent gasoline mechanism to investigate the impact of the calculations performed in this work. A shift in predicted NTC (negative temperature coefficient) behavior and significant reduction in model reactivity is observed upon incorporating the rates calculated herein.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2018

Green’s Function Coupled-Cluster Approach: Simulating Photoelectron Spectra for Realistic Molecular Systems

In this paper we present an efficient implementation for the analytical energy- dependent coupled cluster Green’s function with singles and doubles (GFCCSD) approach with our first practice being computing spectral functions of real molecular systems. Due to its algebraic structure, the presented method is highly scalable, and is capable of computing spectral function for a given molecular system in any energy region. Several typical examples have been given to demonstrate its capability of computing spectral functions not only in the valence band, but also in the core-level energy region. Satellite peaks have been observed in the inner valence band and core-level energy region where many-body effect becomes significant and single particle picture of ionization often breaks down. The accuracy test has been carried out by extensively comparing our GFCCSD results with experimental results as well as the theoretical results from other methods. It turns out GFCCSD method is able to provide qualitative or semi-quantitative level of description of ionization processes in both the core and valence regimes. To significantly improve the GFCCSD results for the main states, larger basis set can usually be employed, whereas the improvement of the GFCCSD results for the satellite states needs higher order many-body terms to be included in the GFCC implementation.

electronic structure theory↗
DOE OSTI2025

Performance of Diffusion Monte Carlo Calculations for Predicting the Relative Energies of Quinoidal and Nonquinoidal Species

Coupled cluster singles and doubles with perturbative triples [CCSD(T)] and single determinant fixed-node diffusion Monte Carlo (SD-DMC) have emerged as two of the most useful methods for providing benchmark reaction and interaction energies of chemical systems without strong static correlation. The errors in DMC energies are dominated by an inexact description of the nodal surfaces for electron exchange. One of the main approaches to addressing the fixed-node error is to use multideterminant (MD) trial wave functions. We consider here the energy differences between pairs of related molecules with aromatic and quinoidal structures as well as between quinoidal isomers. Quinoidal systems tend to have some diradical character, leading one to anticipate that SD-DMC calculations may face challenges in accurately describing their energetics. The MD trial wave functions were generated from the complete active space calculations. A comparison is made with the predictions of well-converged CCSD(T) calculations.

basis sets↗
DOE OSTI2025

Mapping causal patterns in crystalline solids

The evolution of the atomic structures of the combinatorial library of Sm-substituted thin film BiFeO 3 along the phase transition boundary from the ferroelectric rhombohedral phase to the non-ferroelectric orthorhombic phase is explored using scanning transmission electron microscopy. Localized properties, including polarization, lattice parameter, and chemical composition, are parameterized from atomic-scale imaging, and their causal relationships are reconstructed using a linear non-Gaussian acyclic model. This approach is further extended to explore the spatial variability of the causal coupling using the sliding window transform method, which revealed that new causal relationships emerged at both the expected locations, such as domain walls and interfaces, and at additional regions forming clusters in the vicinity of the walls or spatially distributed features. While the exact physical origins of these relationships are unclear, they likely represent nanophase-separated regions in the morphotropic phase boundaries. Overall, we posit that an in-depth understanding of complex disordered materials away from thermodynamic equilibrium necessitates understanding not only the generative processes that can lead to observed microscopic states but also the causal links between multiple interacting subsystems.

Causal inference↗
NASA NTRS1994

Theoretical Study of the First Transition Row Oxides and Sulfides

The first transition row oxides and sulfides are studied using several different levels of theory. The calculations show the bonding mechanism in the sulfides and oxides to be very similar. For the oxides, accurate experimental data allow the theoretical methods to be calibrated. The same level of theory is used to study the sulfides where there is far less experimental information. For ScO through MnO and CuO the coupled cluster singles and doubles technique including a perturbational estimate of the unliked triple excitations [CCSD(T)] yields spectroscopic constants ((tau)e, (omega)e, and D0) in good agreement with experiment. The triple excitations are found to be very important in achieving this accuracy. For FeO to NiO, the single determinant self-consistent-field (SCF) approach yields pi orbitals that are localized on the metal or oxygen. This appears to cause problems for the single reference techniques; this is discussed in detail for NiO. The complete-active-space SCF/internally contracted averaged coupled pair functional approach (CASSCF/ICACPF) works well for FeO to NiO. The calculation of accurate dipole moments is found to be very difficult.

Bauschlicher, Charles W., Jt.↗
DOE OSTI2014

Method of increments for the halogen molecular crystals: Cl, Br, and I

Method of increments (MI) calculations reveal the n-body correlation contributions to binding in solid chlorine, bromine, and iodine. Secondary binding contributions as well as d-correlation energies are estimated and compared between each solid halogen. We illustrate that binding is entirely determined by two-body correlation effects, which account for >80% of the total correlation energy. One-body, three-body, and exchange contributions are repulsive. Using density-fitting (DF) local coupled-cluster singles, doubles, and perturbative triples for incremental calculations, we obtain excellent agreement with the experimental cohesive energies. MI results from DF local second-order Møller-Plesset perturbation (LMP2) yield considerably over-bound cohesive energies. Comparative calculations with density functional theory and periodic LMP2 method are also shown to be less accurate for the solid halogens.

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

Spin Reorientation in Antiferromagnetic Layered FePt 5 P

FePt 5 P, a substitutional variant of the anti-CeCoIn 5 structure type in the space group P4/mmm, was synthesized by a high-temperature solid-state method and structurally characterized by X-ray diffraction. FePt 5 P contains layers of FePt 12 clusters formed by magnetically active Fe and heavy Pt with strong spin-orbit coupling (SOC); the layers are separated by P atoms. The various Fe–Pt distances in FePt 12 clusters generate complex magnetic orders in FePt 5 P. According to temperature-dependent magnetic and specific heat measurements, FePt 5 P shows a stripe-type antiferromagnetic order at T N ≈ 90 K, which is also confirmed by resistivity measurements. Furthermore, a spin reorientation occurs at ~74 and ~68 K in and out of the ab plane based on the specific heat measurements. The temperature-dependent neutron powder diffraction patterns demonstrate the antiferromagnetic order in FePt 5 P, and the spins orientate up to 58.4° with respect to the c axis at 10 K. First-principles calculations of FePt 5 P show the band splitting at the Fermi level by strong SOC and the s–d hybridization between P and Fe/Pt electrons enhances the structural stability and affects the magnetic ordering.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Adaptive multiconfigurational wave functions

A method is suggested to build simple multiconfigurational wave functions specified uniquely by an energy cutoff Λ. These are constructed from a model space containing determinants with energy relative to that of the most stable determinant no greater than Λ. The resulting Λ-CI wave function is adaptive, being able to represent both single-reference and multireference electronic states. We also consider a more compact wave function parameterization (Λ+SD-CI), which is based on a small Λ-CI reference and adds a selection of all the singly and doubly excited determinants generated from it. We report two heuristic algorithms to build Λ-CI wave functions. The first is based on an approximate prescreening of the full configuration interaction space, while the second performs a breadth-first search coupled with pruning. The Λ-CI and Λ+SD-CI approaches are used to compute the dissociation curve of N{sub 2} and the potential energy curves for the first three singlet states of C{sub 2}. Special attention is paid to the issue of energy discontinuities caused by changes in the size of the Λ-CI wave function along the potential energy curve. This problem is shown to be solvable by smoothing the matrix elements of the Hamiltonian. Our last example, involving the Cu{sub 2}O{sub 2}{sup 2+} core, illustrates an alternative use of the Λ-CI method: as a tool to both estimate the multireference character of a wave function and to create a compact model space to be used in subsequent high-level multireference coupled cluster computations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2015

Exotic structure in light neutron-rich nuclei

In this presentation I discussed two subjects. One is the persistence of threefold symmetry in the ground state of {sup 12}C. Recently D{sub 3h} symmetry has been established in {sup 12}C, which reflects the geometric symmetry of the three α particles. Although the spin-orbit interaction plays a significant role and this interaction breaks the α clusters, we show that threefold symmetry of {sup 12}C is still there. We use AQCM approach and discuss that inclusion of spin-orbit interaction, which is absent in the conventional microscopic α cluster models, is possible keeping the threefold symmetry. The second subject is the appearance of rod shape in C isotopes, which has been investigated in the framework of the cranking covariant density functional theory. The relationship between the stability of such states and the spin and isospin degrees of freedom is discussed; adding valence neutrons and rotating the system. These two effects stabilize the rod shape, and in addition, their coherent effect has been found; the σ-orbits (parallel to the symmetry axis) of the valence neutrons, which enhances the rod shape, is lowered by the rotation of the system, and this σ-orbit pulls down the single particle energies of protons with linear configuration owing to the proton-neutron interaction effect.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗