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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 829 records · Page 46

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Spin Decoherence in VOPc@graphene Nanoribbon Complexes

Carbon nanoribbons or nanographene qubit arrays can facilitate quantum-to-quantum transduction between light, charge, and spin, making them an excellent testbed for fundamental science in quantum coherent systems and for the construction of higher-level qubit circuits. Here, in this work, we study spin decoherence due to coupling with a surrounding nuclear spin bath of an electronic molecular spin of a vanadyl phthalocyanine (VOPc) molecule integrated onto an armchair-edged graphene nanoribbon (GNR). Density functional theory (DFT) is used to obtain ground-state atomic configurations. Decay of spin coherence in Hahn echo experiments is then simulated using the cluster correlation expansion method with a spin Hamiltonian involving hyperfine and electric field gradient tensors calculated from DFT. We find that the decoherence time T 2 is anisotropic with respect to magnetic field orientation and determined only by the hydrogen nuclear spins on both VOPc and GNR. Large electron spin echo envelope modulation (ESEEM) due to nitrogen and vanadium nuclear spins is present at specific field ranges and can be completely suppressed by tuning the magnetic field. The relation between these field ranges and the hyperfine interactions is analyzed. The effects of interactions with the nuclear quadrupole moments are also studied, validating the applicability and limitations of the spin Hamiltonian when they are disregarded.

Hamiltonians↗

DOE OSTI2020

Ab initio multishell valence-space Hamiltonians and the island of inversion

In the shell-model framework, valence-space Hamiltonians connecting multiple major-oscillator shells are of key interest for investigating the physics of neutron-rich nuclei, which have been the subject of intense experimental activity for decades. Here we present an extension of the ab initio valence-space in-medium similarity renormalization group, which allows the derivation of such Hamiltonians nonperturbatively. Starting from initial two- and three-nucleon forces from chiral effective field theory, we then calculate properties of nuclei in the important island-of-inversion region above oxygen, so far unexplored with ab initio methods. In this work, our results in the neon and magnesium isotopes indicate the importance of neutron excitation from the $\textit{sd}$ to $\textit{pf}$ shells and ground states dominated by intruder configurations around $\textit{N}$ = 20, consistent with the conclusions from phenomenological studies. We also benchmark the excitation spectrum of 16 O with coupled-cluster theory, finding generally good agreement, and discuss implications for ground-state energies and charge radii in oxygen and calcium isotopes. Finally we outline the proper procedure for treating the longstanding issue of center-of-mass contamination, and show that with a particular choice of valence space, these spurious states can be removed successfully.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1992

Computational studies of metal-metal and metal-ligand interactions

The geometric structure of Cr(CO)6 is optimized at the modified coupled-pair functional (MCPF), single and double excitation coupled-cluster (CCSD) and CCSD(T) levels of theory (including a perturbational estimate for connected triple excitations), and the force constants for the totally symmetric representation are determined. The geometry of Cr(CO)5 is partially optimized at the MCPF, CCSD and CCSD(T) levels of theory. Comparison with experimental data shows that the CCSD(T) method gives the best results for the structures and force constants, and that remaining errors are probably due to deficiencies in the one-particle basis sets used for CO. A detailed comparison of the properties of free CO is therefore given, at both the MCPF and CCSD/CCSD(T) levels of treatment, using a variety of basis sets. With very large one-particle basis sets, the SSCD(T) method gives excellent results for the bond distance, dipole moment and harmonic frequency of free CO. The total binding energies of Cr(CO)6 and Cr(CO)5 are also determined at the MCPF, CCSD and CCSD(T) levels of theory. The CCSD(T) method gives a much larger total binding energy than either the MCPF or CCSD methods. An analysis of the basis set superposition error (BSSE) at the MCPF level of treatment points out limitations in the one-particle basis used here and in a previous study. Calculations using larger basis sets reduced the BSSE, but the total binding energy of Cr(CO)6 is still significantly smaller than the experimental value, although the first CO bond dissociation energy of Cr(CO)6 is well described. An investigation of 3s3p correlation reveals only a small effect. The remaining discrepancy between the experimental and theoretical total binding energy of Cr(CO)6 is probably due to limitations in the one-particle basis, rather than limitations in the correlation treatment. In particular an additional d function and an f function on each C and O are needed to obtain quantitative results. This is underscored by the fact that even using a very large primitive se (1042 primitive functions contracted to 300 basis functions), the superposition error for the total binding energy of Cr(CO)6 is 22 kcal/mol at the MCPF level of treatment.

Barnes, Leslie A.↗

NASA NTRS1992

The structure and energetics of Cr(CO)6 and Cr(CO)5

The geometric structure of Cr(CO)6 is optimized at the modified coupled pair functional (MCPF), single and double excitation coupled-cluster (CCSD) and CCSD(T) levels of theory (including a perturbational estimate for connected triple excitations), and the force constants for the totally symmetric representation are determined. The geometry of Cr(CO)5 is partially optimized at the MCPF, CCSD, and CCSD(T) levels of theory. Comparison with experimental data shows that the CCSD(T) method gives the best results for the structures and force constants, and that remaining errors are probably due to deficiencies in the one-particle basis sets used for CO. The total binding energies of Cr(CO)6 and Cr(CO)5 are also determined at the MCPF, CCSD, and CCSD(T) levels of theory. The CCSD(T) method gives a much larger total binding energy than either the MCPF or CCSD methods. An analysis of the basis set superposition error (BSSE) at the MCPF level of treatment points out limitations in the one-particle basis used. Calculations using larger basis sets reduce the BSSE, but the total binding energy of Cr(CO)6 is still significantly smaller than the experimental value, although the first CO bond dissociation energy of Cr(CO)6 is well described. An investigation of 3s3p correlation reveals only a small effect. In the largest basis set, the total CO binding energy of Cr(CO)6 is estimated to be 140 kcal/mol at the CCSD(T) level of theory, or about 86 percent of the experimental value. The remaining discrepancy between the experimental and theoretical value is probably due to limitations in the one-particle basis, rather than limitations in the correlation treatment. In particular an additional d function and an f function on each C and O are needed to obtain quantitative results. This is underscored by the fact that even using a very large primitive set (1042 primitive functions contracted to 300 basis functions), the superposition error for the total binding energy of Cr(CO)6 is 22 kcal/mol at the MCPF level of treatment.

Barnes, Leslie A.↗

DOE OSTI2024

Characterization of ferredoxins involved in electron transfer pathways for nitrogen fixation implicates differences in electronic structure in tuning 2[4Fe 4S] Fd activity

Ferredoxins (Fds) are small proteins which shuttle electrons to pathways like biological nitrogen fixation. Physical properties tune the reactivity of Fds with different pathways, but knowledge on how these properties can be manipulated to engineer new electron transfer pathways is lacking. Recently, we showed that an evolved strain of Rhodopseudomonas palustris uses a new electron transfer pathway for nitrogen fixation. This pathway involves a variant of the primary Fd of nitrogen fixation in R. palustris, Fer1, in which threonine at position 11 is substituted for isoleucine (Fer1 T11I ). To understand why this substitution in Fer1 enables more efficient electron transfer, we used in vivo and in vitro methods to characterize Fer1 and Fer1 T11I . Electrochemical characterization revealed both Fer1 and Fer1 T11I have similar redox transitions (–480 mV and – 550 mV), indicating the reduction potential was unaffected despite the proximity of T11 to an iron-sulfur (Fe—S) cluster of Fer1. Additionally, disruption of hydrogen bonding around an Fe—S cluster in Fer1 by substituting threonine with alanine (T11A) or valine (T11V) did not increase nitrogenase activity, indicating that disruption of hydrogen bonding does not explain the difference in activity observed for Fer1 T11I . Electron paramagnetic resonance spectroscopy studies revealed key differences in the electronic structure of Fer1 and Fer1 T11I , which indicate changes to the high spin states and/or spin-spin coupling between the Fe—S clusters of Fer1. Finally, our data implicates these electronic structure differences in facilitating electron flow and sets a foundation for further investigations to understand the connection between these properties and intermolecular electron transfer.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2014

Revisiting Photoemission and Inverse Photoemission Spectra of Nickel Oxide from First Principles: Implications for Solar Energy Conversion

We use two different ab initio quantum mechanics methods, complete active space self-consistent field theory applied to electrostatically embedded clusters and periodic many-body G₀W₀ calculations, to reanalyze the states formed in nickel(II) oxide upon electron addition and ionization. In agreement with interpretations of earlier measurements, we find that the valence and conduction band edges consist of oxygen and nickel states, respectively. However, contrary to conventional wisdom, we find that the oxygen states of the valence band edge are localized whereas the nickel states at the conduction band edge are delocalized. We argue that these characteristics may lead to low electron-hole recombination and relatively efficient electron transport, which, coupled with band gap engineering, could produce higher solar energy conversion efficiency compared to that of other transition-metal oxides. Both methods find a photoemission/inverse-photoemission gap of 3.6-3.9 eV, in good agreement with the experimental range, lending credence to our analysis of the electronic structure of NiO.

Environmental Molecular Sciences Laboratory↗

DOE OSTI2016

Comparing simulations of AGN feedback

We perform adaptive mesh refinement (AMR) and smoothed particle hydrodynamics (SPH) cosmological zoom simulations of a region around a forming galaxy cluster, comparing the ability of the methods to handle successively more complex baryonic physics. In the simplest, non-radiative case, the two methods are in good agreement with each other, but the SPH simulations generate central cores with slightly lower entropies and virial shocks at slightly larger radii, consistent with what has been seen in previous studies. The inclusion of radiative cooling, star formation, and stellar feedback leads to much larger differences between the two methods. Most dramatically, at z=5, rapid cooling in the AMR case moves the accretion shock to well within the virial radius, while this shock remains near the virial radius in the SPH case, due to excess heating, coupled with poorer capturing of the shock width. On the other hand, the addition of feedback from active galactic nuclei (AGNs) to the simulations results in much better agreement between the methods. For our AGN model, both simulations display halo gas entropies of 100 keV cm{sup 2}, similar decrements in the star formation rate, and a drop in the halo baryon content of roughly 30%. This is consistent with the AGN growth being self-regulated, regardless of the numerical method. However, the simulations with AGN feedback continue to differ in aspects that are not self-regulated, such that in SPH a larger volume of gas is impacted by feedback, and the cluster still has a lower entropy central core.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Computational Screening of Cationic Zeolites for n -Butane/Methane Separations Using Quantitatively Accurate First-Principles-Derived Force Fields

We developed force fields for linear alkanes in Na and Ca-exchanged zeolites based on periodic first-principles calculations with coupled cluster corrections. These force fields were validated by comparing the simulated adsorption properties with extensive experimental data. As an example of using these force fields, we screened silica and aluminosilicate zeolites for separation of n-butane/methane mixtures in several different separation processes. A key step in this screening was the use of an interpolation method to create accurate zeolite structures that take account of the lattice constant changes for different Si/Al ratios for a given zeolite topology. We focused on zeolite structures that have known synthesis routes, so the materials selected from our calculations can be tested experimentally. Several promising materials that show good separation performance were selected for more detailed simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Towards a Description of the Ground State Properties of LiHm Complexes and Reaction Energetics of Hydrogenation of Small Li Clusters using VQE-based Quantum Chemical Calculations

This presentation describes the utilization of the variational quantum eigensolver (VQE) using a unitary coupled cluster with singles and doubles (UCCSD) ansatz to simulate a series of LiHn molecules (n=1-3), including their singly charged ions. The performance is considered using a simulator by systematically increasing system size, providing data to compare against future VQE calculations. The predicted energy profiles are in semi-quantitative agreement with highly accurate CCSD/FCI methods.

quantum computing↗

NASA NTRS2001

Enhancement of Tropical Land Cover Mapping with Wavelet-Based Fusion and Unsupervised Clustering of SAR and Landsat Image Data

The characterization and the mapping of land cover/land use of forest areas, such as the Central African rainforest, is a very complex task. This complexity is mainly due to the extent of such areas and, as a consequence, to the lack of full and continuous cloud-free coverage of those large regions by one single remote sensing instrument, In order to provide improved vegetation maps of Central Africa and to develop forest monitoring techniques for applications at the local and regional scales, we propose to utilize multi-sensor remote sensing observations coupled with in-situ data. Fusion and clustering of multi-sensor data are the first steps towards the development of such a forest monitoring system. In this paper, we will describe some preliminary experiments involving the fusion of SAR and Landsat image data of the Lope Reserve in Gabon. Similarly to previous fusion studies, our fusion method is wavelet-based. The fusion provides a new image data set which contains more detailed texture features and preserves the large homogeneous regions that are observed by the Thematic Mapper sensor. The fusion step is followed by unsupervised clustering and provides a vegetation map of the area.

LeMoigne, Jacqueline↗

DOE OSTI2022

Analysis of spin frustration in an Fe III 7 cluster using a combination of computational, experimental, and magnetostructural correlation methods

The synthesis, structure, and magnetic properties are reported for [Fe 7 O 3 (O 2 C t Bu) 9 (mda) 3 (H 2 O) 3 ] ( 1 ), where mdaH 2 is N -methyldiethanolamine. 1 was prepared from the reaction of [Fe 3 O(O 2 C t Bu) 6 (H 2 O) 3 ](NO 3 ) with mdaH 2 in a 1:~3 ratio in MeCN. The core of 1 consists of a central octahedral Fe III ion held within a non-planar Fe 6 loop by three μ 3 -O 2- and three μ 2 -RO - arms from the three mda 2- chelates. Variable-temperature dc and ac magnetic susceptibility studies revealed dominant antiferromagnetic coupling, leading to a ground state spin of S = 5 / 2 . The ground state was confirmed by a fit of magnetization data collected in the 0.1–7.0 T and 1.8–10.0 K ranges. The four Fe 2 pairwise exchange parameters ( J 1 - J 4 ) were estimated by independent methods: theoretical calculations using either broken symmetry energy differences (-46.3, -16.2, -3.9, and - 28.1 cm -1 , respectively) or Green’s function approximation methods (-41.4, -14.8, -13.2, and - 24.7 cm -1 ), and a magnetostructural correlation (MSC) previously developed for high nuclearity Fe III /O complexes (-39.5, -13.8, -6.7, and - 23.5 cm -1 ). Additionally, the J 1 - J 4 obtained from the MSC and theoretical methods were used with the program PHI to both simulate χ M T vs T as well as to serve as reasonable input values to fit the experimental data (-41.0, -11.4, -5.0, and - 27.3 cm -1 ). Analysis of the J ij led to identification of the spin frustration effects operative and the resultant spin vector alignments at each Fe III ion, thus allowing for the rationalization of the experimental ground state.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2018

The effect of Mn/Ni on thermodynamic properties of critical nucleus in Fe-Cu-Mn (Ni) ternary alloys

The aging- or radiation-induced hardening of Cu/Mn/Ni precipitates in Fe alloys is one of property degradation mechanisms in structural materials in nuclear reactors. Experiments show that aging or radiation leads the formation of Cu-rich precipitates, and the addition of Mn or Ni elements enhances the precipitation kinetics. In this study, the phase-field model coupled with the constrained string method have been applied to investigate the thermodynamic properties of critical nuclei such as the minimum energy path of Cu/Mn/Ni precipitation in Fe-Cu-Mn and Fe-Cu-Ni ternary alloys. The chemical free energies used in the model are taken from CALPHAD. The simulation results show that the formation of Cu/Mn/Ni clusters needs to overcome an energy barrier, and the precipitate has a Core-Shell structure. The thermodynamic properties of the critical nucleus are influenced by temperature and Cu/Mn/Ni overall concentrations, which are in accordance with the simulation results as well as the experimental observations.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2021

Charge radii of exotic potassium isotopes challenge nuclear theory and the magic character of N = 32

Nuclear charge radii are sensitive probes of different aspects of the nucleon–nucleon interaction and the bulk properties of nuclear matter, providing a stringent test and challenge for nuclear theory. Experimental evidence suggested a new magic neutron number at N = 32 in the calcium region, whereas the unexpectedly large increases in the charge radii open new questions about the evolution of nuclear size in neutron-rich systems. By combining the collinear resonance ionization spectroscopy method with β-decay detection, we were able to extend charge radii measurements of potassium isotopes beyond N = 32. Here we provide a charge radius measurement of 52 K. It does not show a signature of magic behaviour at N = 32 in potassium. The results are interpreted with two state-of-the-art nuclear theories. The coupled cluster theory reproduces the odd–even variations in charge radii but not the notable increase beyond N = 28. This rise is well captured by Fayans nuclear density functional theory, which, however, overestimates the odd–even staggering effect in charge radii. These findings highlight our limited understanding of the nuclear size of neutron-rich systems, and expose problems that are present in some of the best current models of nuclear theory.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Lowest-energy broad 𝛼-cluster resonances in 19 F

There is strong astrophysical interest in the structure of 19 F states near the α-decay threshold, as their properties are important for understanding the development of α clustering in the 20 Ne region. The emergence of clustered states, and more generally of states that couple strongly to reaction channels near their decay thresholds, is also a topic of current interest in theoretical nuclear physics. In this work, we determine the parameters of broad low-spin states in 19 F near the α-decay threshold and present a theoretical study of these states. The analysis is restricted to ℓ = 0 and 1 resonances in the α + 15 N system close to the α threshold in 19 F. Excitation functions for 15 N(α,α) elastic scattering were measured using the Thick Target Inverse Kinematics (TTIK) method, and these new data, together with older high–energy-resolution measurements, were analyzed within the R-matrix framework. The nuclear structure of 19 F was calculated using configuration-interaction methods with a recently developed effective interaction Hamiltonian. As a result, we identify a sequence of α-cluster resonances in 19 F and map the distribution of clustering strength, which is relevant for astrophysical reaction modeling and for the theoretical understanding of many-body dynamics and the emergence of clustering in loosely bound or unstable nuclei. Furthermore, the work advances theoretical insight into the origins of clustering and highlights open questions for future theoretical and experimental studies.

19F↗

DOE OSTI2015

Octahedral tilt transitions in the relaxor ferroelectric Na{sub 1/2}Bi{sub 1/2}TiO{sub 3}

The kinetics of octahedral tilt transitions in the lead-free relaxor material sodium bismuth titanate Na{sub 1/2}Bi{sub 1/2}TiO{sub 3} (NBT) is investigated by electronic structure calculations within density functional theory. Energy barriers for transitions between tetragonal, rhombohedral and orthorhombic tilts in cation configurations with [001]- and [111]-order on the A-sites are determined by nudged elastic band calculations. By tilting entire layers of octahedra simultaneously we find that the activation energy is lower for structures with 001-order compared to such with 111-order. The energetic coupling between differently tilted layers is, however, negligibly small. By introducing a single octahedral defect we create local tilt disorder and find that the deformation energy of the neighboring octahedra is less in a rhombohedral than in a tetragonal structure. By successively increasing the size of clusters of orthorhombic defects in a rhombohedral matrix with 001-order, we determine a critical cluster size of about 40 Å . Thus groups of about ten octahedra can be considered as nuclei for polar nanoregions, which are the cause of the experimentally observed relaxor behavior of NBT. - Graphical abstract: Nine orthorhombic oxygen octahedral tilt defects in a rhombohedral tilt configuration. - Highlights: • Chemical order influences energy barriers of octahedral tilt transitions. • The octahedral deformation energy is lower in rhombohedral phases. • Tilt defect clusters are more likely in rhombohedral structures. • Tilt defect clusters can act as nuclei for polar nanoregions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2016

Spectroscopic properties of oxygen vacancies in LaAlO 3

Oxygen vacancies in LaAlO 3 (LAO) play an important role in the formation of the two-dimensional electron gas observed at the LaAlO 3 /SrTiO 3 interface and affect the performance of MOSFETs using LAO as a gate dielectric. However, their spectroscopic properties are still poorly understood, which hampers their experimental identification. Here we predict the absorption spectra and ESR parameters of oxygen vacancies in LAO using periodic and embedded cluster methods and density functional theory (DFT). The structure, charge distribution, and spectroscopic properties of the neutral (V 0 O) and charged (V + O and V 2+ O ) oxygen vacancies in cubic and rhombohedral LaAlO 3 are investigated. The highest intensity optical transitions [calculated using time-dependent DFT (TDDFT)], from the oxygen vacancy states to the conduction-band states have onsets at 3.5 and 4.2 eV for V 0 O and 3.6 eV for V + O in rhombohedral LAO and 3.3 and 4.0 eV for V 0 O and 3.4 eV for V + O in cubic LAO, respectively. Also reported are the isotropic g value (2.004026) and hyperfine coupling constants of V + O , which are compared to the experimental data obtained using electron spin resonance (ESR) spectroscopy, and accurately predict both the position and the width (3 mT) of its ESR signature. These results may further facilitate the experimental identification of oxygen vacancies in LAO and help to establish their role at the LAO/STO interfaces and in nanodevices using LAO.

36 MATERIALS SCIENCE↗

DOE OSTI2019

Fully Quantum Embedding with Density Functional Theory for Full Configuration Interaction Quantum Monte Carlo

We here develop a fully quantum embedded version of initiator full configuration interaction quantum Monte Carlo (i–FCIQMC) and apply it to study an ionic bond (lithium hydride, LiH) and a covalent bond (hydrogen flouride, HF) physisorbed to a benzene molecule. The embedding is performed using a recently developed Huzinaga projection operator approach, which affords good synergy with i–FCIQMC by minimizing the number of orbitals in the calculation. When considering the dissociation energy of these bonds into closed-shell ionic fragments, we find that i–FCIQMC embedded in density functional theory (i–FCIQMC-in-DFT) delivers comparable accuracy with coupled cluster singles and doubles with perturbative triples embedded in density functional theory (CCSD(T)-in-DFT). In treating the bond dissociation energy curve of HF, i–FCIQMC-in-DFT has improved accuracy over CCSD(T)-in-DFT due to the presence of strong correlation. Lastly, we discuss the implications of the new i–FCIQMC-in-DFT method as applied to bond breaking in catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2005

Testing the Archivas Cluster (Arc) for Ozone Monitoring Instrument (OMI) Scientific Data Storage

The Ozone Monitoring Instrument (OMI) launched on NASA's Aura Spacecraft, the third of the major platforms of the EOS program on July 15,2004. In addition to the long term archive and distribution of the data from OM1 through the Goddard Earth Science Distributed Active Archive Center (GESDAAC), we are evaluating other archive mechanisms that can archive the data in a more immediately available method where it can be used for futher data production and analysis. In 2004, Archivas, Inc. was selected by NASA s Small Business Innovative Research (SBIR) program for the development of their Archivas Cluster (ArC) product. Arc is an online disk based system utilizing self-management and automation on a Linux cluster. Its goal is to produce a low cost solution coupled with the ease of management. The OM1 project is an application partner of the SBIR program, and has deployed a small cluster (5TB) based on the beta Archwas software. We performed extensive testing of the unit using production OM1 data since launch. In 2005, Archivas, Inc. was funded in SBIR Phase II for further development, which will include testing scalability with the deployment of a larger (35TB) cluster at Goddard. We plan to include Arc in the OM1 Team Leader Computing Facility (TLCF) hosting OM1 data for direct access and analysis by the OMI Science Team. This presentation will include a brief technical description of the Archivas Cluster, a summary of the SBIR Phase I beta testing results, and an overview of the OMI ground data processing architecture including its interaction with the Phase II Archivas Cluster and hosting of OMI data for the scientists.

Tilmes, Curt↗