SEARCH · Search NASA

Results for “Coupled cluster theory”

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 289 records · Page 16

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Solvothermal Synthesis of [Cr 7 S 8 (en) 8 Cl 2 ]Cl 3 ⋅ 2H 2 O with Magnetically Frustrated [Cr 7 S 8 ] 5+ Double‐Cubes**

Abstract A novel transition metal chalcohalide [Cr 7 S 8 (en) 8 Cl 2 ]Cl 3 ⋅ 2H 2 O, with [Cr 7 S 8 ] 5+ dicubane cationic clusters, has been synthesized by a low temperature solvothermal method, using dimethyl sulfoxide (DMSO) and ethylenediamine ( en ) solvents. Ethylenediamine ligand exhibits bi‐ and monodentate coordination modes; in the latter case ethylenediamine coordinates to Cr atoms of adjacent clusters, giving rise to a 2D polymeric structure. Although magnetic susceptibility shows no magnetic ordering down to 1.8 K, a highly negative Weiss constant, θ =−224(2) K, obtained from Curie‐Weiss fit of inverse susceptibility, suggests strong antiferromagnetic (AFM) interactions between S =3/2 Cr(III) centers. Due to the complexity of the system with (2 S +1) 7 =16384 microstates from seven Cr 3+ centers, a simplified model with only two exchange constants was used for simulations. Density‐functional theory (DFT) calculations yielded the two exchange constants to be J 1 =−21.4 cm −1 and J 2 =−30.2 cm −1 , confirming competing AFM coupling between the shared Cr 3+ center and the peripheral Cr 3+ ions of the dicubane cluster. The best simulation of the experimental data was obtained with J 1 =−20.0 cm −1 and J 2 =−21.0 cm −1 , in agreement with the slightly stronger AFM exchange within the triangles of the peripheral Cr 3+ ions as compared to the AFM exchange between the central and peripheral Cr 3+ ions. This compound is proposed as a synthon towards magnetically frustrated systems assembled by linking dicubane transition metal‐chalcogenide clusters into polymeric networks.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Near-Threshold Dipole Strength in 10 Be with Isoscalar Character

Isoscalar dipole transitions are a distinctive fingerprint of cluster structures. A 1 − resonance at 7.27(10) MeV, located just below the 𝛼-emission threshold, has been observed in the deuteron inelastic scattering reactions off 10 Be. The deformation lengths of the excited states in 10 Be below 9 MeV have been inferred from the differential cross sections using coupled channel calculations. This observed 1 − resonance has isoscalar characteristics and exhausts approximately 5%–15% of the isoscalar dipole energy-weighted sum rule, providing evidence for pronounced 𝛼 cluster structure in 10 Be. The Gamow coupled channel approach supports this interpretation and suggests the near-threshold effect might be playing an important role in this excitation energy domain. Here, the 𝛼 +𝛼 +𝑛 +𝑛 four-body calculation reproduces the observed enhanced dipole strength, implying that the four-body cluster structure is essential to describe the 1 − states in 10 Be.

Cluster models↗

DOE OSTI2023

Theoretical Kinetics Predictions for Reactions on the NH 2 O Potential Energy Surface

Recent modeling studies of ammonia oxidation, which are motivated by the prospective role of ammonia as a zero-carbon fuel, have indicated significant discrepancies among the existing literature mechanisms. In this study, high-level theoretical kinetics predictions have been obtained for reactions on the NH 2 O potential energy surface, including the NH 2 + O, HNO + H, and NH + OH reactions. These reactions have previously been highlighted as important reactions in NH 3 oxidation with high sensitivity and high uncertainty. The potential energy surface is explored with coupled cluster calculations, including large basis sets and high-level corrections to yield high-accuracy (~0.2 kcal/mol 2σ uncertainty) estimates of the stationary point energies. Variational transition state theory is used to predict the microcanonical rate constants, which are then incorporated in master equation treatments of the temperature- and pressure-dependent kinetics. For radical–radical channels, the microcanonical rates are obtained from variable reaction coordinate transition state theory implementing directly evaluated multireference electronic energies. The analysis yields predictions for the total rate constants as well as the branching ratios. We find that the NO + H 2 channel contributes 10% of the total NH 2 + O flux at combustion temperatures, although this channel is not included in modern NH 3 oxidation mechanisms. Modeling is used to illustrate the ramifications of these rate predictions on the kinetics of NH 3 oxidation and NO x formation. Finally, the present results for NH 2 + O are important for predicting the chain branching and formation of NO in the oxidation of NH 3 and thermal DeNO x .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Exactness of the normal-ordered two-body truncation of three-nucleon forces

Reference-state-based many-body methods start from Hamiltonians that are normal ordered with respect to the reference state. In low-energy nuclear physics applications, normal-ordered Hamiltonians consisting of two- and three-nucleon forces are usually truncated at the two-body rank with residual three-nucleon operators being discarded. Benchmark computations have shown that this truncation is accurate, but we lack an understanding about why it works. Here, we show that the normal-ordered two-body truncation is exact for zero-range three-body forces when nuclei are computed using the coupled cluster with singles and doubles method. As the nuclear three-nucleon force is short ranged and a three-body contact is a leading term in effective field theories of quantum chromodynamics, our result provides an analytical basis for the popular normal-ordered two-body approximation.

Rothman, Maxwell [Univ. of Tennessee, Knoxville, T↗

DOE OSTI2025

Review of radiation-induced defects in GaAs

Radiation-resilient optoelectronic materials are highly desired in space and nuclear applications, and understanding the relevant defects and electronic processes in those materials is crucial for both current and next-generation applications. GaAs is a prototypical semiconductor that serves as a foundational system for describing and understanding optoelectronic devices. In this review, both experimental studies and molecular dynamics (MD) simulations of irradiated GaAs are reviewed, with particular emphasis on the deep-level transient spectroscopy, irradiation-induced amorphization zones, and MD predictions of damage structures. The MD results are also compared to predictions of the non-ionized energy-loss model. Recent theoretical studies, in particular, density functional theory based calculations on the simple intrinsic defects and defect clusters in GaAs, are also reviewed. These defects have an important role in dictating the evolution of GaAs in radiation damage environments, as they impact the coupled dynamics of charge carriers. Finally, possible gaps and challenges toward the general understanding of defect evolution in GaAs are discussed.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Mechanical coupling in the nitrogenase complex

The enzyme nitrogenase reduces dinitrogen to ammonia utilizing electrons, protons, and energy obtained from the hydrolysis of ATP. Mo-dependent nitrogenase is a symmetric dimer, with each half comprising an ATP-dependent reductase, termed the Fe Protein, and a catalytic protein, known as the MoFe protein, which hosts the electron transfer P-cluster and the active-site metal cofactor (FeMo-co). A series of synchronized events for the electron transfer have been characterized experimentally, in which electron delivery is coupled to nucleotide hydrolysis and regulated by an intricate allosteric network. We report a graph theory analysis of the mechanical coupling in the nitrogenase complex as a key step to understanding the dynamics of allosteric regulation of nitrogen reduction. This analysis shows that regions near the active sites undergo large-scale, large-amplitude correlated motions that enable communications within each half and between the two halves of the complex. Computational predictions of mechanically regions were validated against an analysis of the solution phase dynamics of the nitrogenase complex via hydrogen-deuterium exchange. These regions include the P-loops and the switch regions in the Fe proteins, the loop containing the residue β-188Ser adjacent to the P-cluster in the MoFe protein, and the residues near the protein-protein interface. In particular, it is found that: (i) within each Fe protein, the switch regions I and II are coupled to the [4Fe-4S] cluster; (ii) within each half of the complex, the switch regions I and II are coupled to the loop containing β-188Ser; (iii) between the two halves of the complex, the regions near the nucleotide binding pockets of the two Fe proteins (in particular the P-loops, located over 130 Å apart) are also mechanically coupled. Notably, we found that residues next to the P-cluster (in particular the loop containing β-188Ser) are important for communication between the two halves.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2023

Excited-State Dynamics during Primary C–I Homolysis in Acetyl Iodide Revealed by Ultrafast Core-Level Spectroscopy

In typical carbonyl-containing molecules, bond dissociation events follow initial excitation to $nπ_{C=O}$$^*$ states. However, in acetyl iodide, the iodine atom gives rise to electronic states with mixed $nπ_{C=O}$$^*$ and $nπ_{C–I}$$^*$ character, leading to complex excited-state dynamics, ultimately resulting in dissociation. Using ultrafast extreme ultraviolet (XUV) transient absorption spectroscopy and quantum chemical calculations, we present an investigation of the primary photodissociation dynamics of acetyl iodide via time-resolved spectroscopy of core-to-valence transitions of the I atom after 266 nm excitation. The probed I 4d-to-valence transitions show features that evolve on sub-100-fs time scales, reporting on excited-state wavepacket evolution during dissociation. These features subsequently evolve to yield spectral signatures corresponding to free iodine atoms in their spin–orbit ground and excited states with a branching ratio of 1.1:1 following dissociation of the C–I bond. Calculations of the valence excitation spectrum via equation-of-motion coupled cluster with single and double substitutions (EOM-CCSD) show that initial excited states are of spin-mixed character. From the initially pumped spin-mixed state, we use a combination of time-dependent density functional theory (TDDFT)-driven nonadiabatic ab initio molecular dynamics and EOM-CCSD calculations of the N$_{4,5}$ edge to reveal a sharp inflection point in the transient XUV signal that corresponds to rapid C–I homolysis. Here, by examining the molecular orbitals involved in the core-level excitations at and around this inflection point, we are able to piece together a detailed picture of C–I bond photolysis in which d → σ* transitions give way to d → p excitations as the bond dissociates. We also report theoretical predictions of short-lived, weak 4d → 5d transitions in acetyl iodide, validated by weak bleaching in the experimental transient XUV spectra. This joint experimental–theoretical effort has thus unraveled the detailed electronic structure and dynamics of a strongly spin–orbit coupled system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

NESSi : The N on- E quilibrium S ystems S imulation package

The nonequilibrium dynamics of correlated many-particle systems is of interest in connection with pump–probe experiments on molecular systems and solids, as well as theoretical investigations of transport properties and relaxation processes. Nonequilibrium Green’s functions are a powerful tool to study interaction effects in quantum many-particle systems out of equilibrium, and to extract physically relevant information for the interpretation of experiments. Here, we present the open-source software package NESSi (The Non-Equilibrium Systems Simulation package) which allows to perform many-body dynamics simulations based on Green’s functions on the L-shaped Kadanoff–Baym contour. NESSi contains the library libcntr which implements tools for basic operations on these nonequilibrium Green’s functions, for constructing Feynman diagrams, and for the solution of integral and integro-differential equations involving contour Green’s functions. The library employs a discretization of the Kadanoff–Baym contour into time points and a high-order implementation of integration routines. The total integrated error scales up to $\mathcal{O}(N^{-7})$, which is important since the numerical effort increases at least cubically with the simulation time. A distributed-memory parallelization over reciprocal space allows large-scale simulations of lattice systems. We provide a collection of example programs ranging from dynamics in simple two-level systems to problems relevant in contemporary condensed matter physics, including Hubbard clusters and Hubbard or Holstein lattice models. The libcntr library is the basis of a follow-up software package for nonequilibrium dynamical mean-field theory calculations based on strong-coupling perturbative impurity solvers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Second Linear Response Theory and the Analytic Calculation of Excited-State Properties

We present a method based on second linear response time-dependent density func- tional theory (TDDFT) to calculate permanent and transition multipoles of excited states, which are required to compute excited-state absorption/emission spectra, multi- photon optical processes, among others. In previous work, we examined computations based on second linear response theory in which linear response TDDFT was employed twice. In contrast, the present methodology only requires information from a single linear response calculation to compute the excited-state properties. These are eval- uated analytically through various algebraic operations involving electron repulsion integrals and excitation vectors. The present derivation focuses on full many-body wave functions instead of single orbitals, as in our previous approach. Here, we test the proposed method by applying it to several diatomic and triatomic molecules. This shows that the computed excited-state dipoles are consistent with respect to reference equation-of-motion coupled cluster calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Conventional single-gap $s-$wave superconductivity and hidden peak effect in single crystals of Mo$_8$Ga$_{41}$ superconductor

London and Campbell penetration depths were measured in single crystals of the endohedral gallide cluster superconductor, Mo _8 8 Ga _41 41 . The full temperature range superfluid density, \rho_s(T) ρ s ( T ) , is consistent with the clean isotropic s- s − wave weak-coupling BCS theory without any signs of the second gap or strong coupling. The temperature dependence of the Campbell length is hysteretic between zero-field cooling (ZFC) and field-cooling (FC) protocols, indicating an anharmonic vortex pinning potential. The field dependence of the effective critical current density, j_{c}(H) j c ( H ) , reveals an unusual result. While in the ZFC protocol, j_{c}(H) j c ( H ) is monotonically suppressed by the magnetic field, it exhibits a profound “hidden” peak effect in the FC protocol, that is, without a vortex density gradient. We suggest a possible novel mechanism for such a peak effect, which involves both static and dynamic aspects of vortex pinning.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

The role of electron correlations in the electronic structure of putative Chern magnet TbMn6Sn6

Abstract A member of the RMn 6 Sn 6 rare-earth family materials, TbMn 6 Sn 6 , recently showed experimental signatures of the realization of a quantum-limit Chern magnet. In this work, we use quantum Monte Carlo (QMC) and density functional theory with Hubbard U (DFT + U ) calculations to examine the electronic structure of TbMn 6 Sn 6 . To do so, we optimize accurate, correlation-consistent pseudopotentials for Tb and Sn using coupled-cluster and configuration–interaction (CI) methods. We find that DFT + U and single-reference QMC calculations suffer from the same overestimation of the magnetic moments as meta-GGA and hybrid density functional approximations. Our findings point to the need for improved orbitals/wavefunctions for this class of materials, such as natural orbitals from CI, or for the inclusion of multi-reference effects that capture the static correlations for an accurate prediction of magnetic properties. DFT + U with Mn magnetic moments adjusted to the experiment predict the Dirac crossing in bulk to be close to the Fermi level, within ~120 meV, in agreement with the experiments. Our non-stoichiometric slab calculations show that the Dirac crossing approaches even closer to the Fermi level, suggesting the possible realization of Chern magnetism in this limit.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Ab initio benchmarks of neutrinoless double- β decay in light nuclei with a chiral Hamiltonian

Here, we report ab initio benchmark calculations of nuclear matrix elements (NMEs) for neutrinoless double-β (0νββ) decays in light nuclei with mass numbers ranging from A = 6 to A = 22. We use the transition operator derived from light-Majorana neutrino exchange and evaluate the NME with three different methods: two variants of in-medium similarity renormalization group (IMSRG) and importance-truncated no-core shell model (IT-NCSM). The same two plus three nucleon interaction from chiral effective field theory is employed, and both isospin-conserving (ΔT = 0) and isospin-changing (ΔT = 2) transitions are studied. We compare our resulting ground-state energies and NMEs to those of recent ab initio no-core shell model and coupled-cluster calculations, also with the same inputs. We show that the NMEs of ΔT = 0 transitions are in good agreement among all calculations, at the level of 10%. For ΔT = 2, relative deviations are more significant in some nuclei. The comparison with the exact IT-NCSM result allows us to analyze these cases in detail, and indicates the next steps toward improving the IMSRG-based approaches. The present study clearly demonstrates the power of consistent cross checks that are made possible by ab initio methodology. This capability is crucial for providing meaningful many-body uncertainties in the NMEs for the 0νββ decays in heavier candidate nuclei, where quasiexact benchmarks are not available.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

The Reaction Rates of Amidogen and Ammonia with Nitrous Oxide: Implications for Combustion Mechanisms

Pulsed laser photolysis experiments with laser-induced fluorescence detection of NH 2 set an upper limit to the rate constant for reaction with N 2 O of k < 1 × 10 −15 cm 3 molecule −1 s −1 at 513 K. Computations were based on geometries and anharmonic frequency analysis (B2PLYP-D3/cc-pVTZ) followed by coupled cluster calculations extrapolated to the infinite basis set limit, with corrections for core−valence electron correlation, scalar relativistic effects, and correlation up to CCSDT(Q). Species that showed multireference character were quantified with MRCI(7,7)+Q/cc-pVTZ theory. Rate constants were obtained for the dominant product channel H 2 NN + NO, along with HN 3 + OH, H 2 NO + N 2 , NNH + HNO, ON(NH) 2 and HNNH + NO. The last channel is slow even at 2500 K, contrary to an early empirical estimate and confirming recent suggestions. Modeling of literature experiments on oxidation of NH 3 by N 2 O shows that all channels are too slow to make a significant impact on the loss of N 2 O in ammonia flames. Similarly, the direct NH 3 + N 2 O reaction is found to be negligibly slow.

Ammonia↗

DOE OSTI2020

Performance of an atomic mean-field spin–orbit approach within exact two-component theory for perturbative treatment of spin–orbit coupling

The paper reports development of an atomic mean-field (AMF) spin–orbit approach within the spin-free exact two-component theory in its one-electron variant (SFX2C-1e), a new relativistic quantum-chemical method for treating spin-orbit coupling in molecules. The effective one-electron spin–orbit integral matrix in the four-component representation is assembled as a direct sum of one-centre spin–orbit integral matrices with the mean-field two-electron contributions evaluated using atomic SFX2C-1e Hartree–Fock density matrices. It is then transformed into two-component representation using analytic SFX2C-1e energy derivative formulation. The resulting two-component spin–orbit integral matrix is by design suitable for use in perturbative calculations of spin–orbit coupling, treating SFX2C-1e wavefunctions as unperturbed states. The accuracy of the present AMF approach has been demonstrated using benchmark calculations of spin–orbit splittings for representative diatomic radicals at the equation-of-motion coupled-cluster singles and doubles level. In conclusion, to demonstrate the applicability and accuracy of the present perturbative spin–orbit scheme in calculations of challenging heavy-element containing systems, a thorough computational investigation of six low-lying electronic states of ThO+ is reported.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2021

Complex magnetic properties associated with competing local and itinerant magnetism in ${\text {Pr}}_2 {\text {Co}}_{0.86} {\text {Si}}_{2.88}$

Ternary intermetallic compound ${\text {Pr}}_2 {\text {Co}}_{0.86} {\text {Si}}_{2.88}$ has been synthesized in single phase and characterized by x-ray diffraction, scanning electron microscopy with energy dispersive x-ray spectroscopy (SEM-EDX) analysis, magnetization, heat capacity, neutron diffraction and muon spin rotation/relaxation ( $\mu$ SR) measurements. The polycrystalline compound was synthesized in single phase by introducing necessary vacancies in Co/Si sites. Magnetic, heat capacity, and zero-field neutron diffraction studies reveal that the system undergoes magnetic transition below $\sim$ 4 K. Neutron diffraction measurement further reveals that the magnetic ordering is antiferromagnetic in nature with an weak ordered moment. The high temperature magnetic phase has been attributed to glassy in nature consisting of ferromagnetic clusters of itinerant (3 d ) Co moments as evident by the development of internal field in zero-field $\mu$ SR below 50 K. The density-functional theory (DFT) calculations suggest that the low temperature magnetic transition is associated with antiferromagnetic coupling between Pr 4 f and Co 3 d spins. Pr moments show spin fluctuation along with unconventional orbital moment quenching due to crystal field. The evolution of the symmetry and the crystalline electric field environment of Pr-ions are also studied and compared theoretically between the elemental Pr and when it is coupled with other elements such as Co. The localized moment of Pr 4 f and itinerant moment of Co 3 d compete with each other below $\sim$ 20 K resulting in an unusual temperature dependence of magnetic coercivity in the system.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Emergence of quasiparticles in a doped Mott insulator

How a Mott insulator develops into a weakly coupled metal upon doping is a central question to understanding various emergent correlated phenomena. To analyze this evolution and its connection to the high- T c cuprates, we study the single-particle spectrum for the doped Hubbard model using cluster perturbation theory on superclusters. Starting from extremely low doping, we identify a heavily renormalized quasiparticle dispersion that immediately develops across the Fermi level, and a weakening polaronic side band at higher binding energy. The quasiparticle spectral weight roughly grows at twice the rate of doping in the low doping regime, but this rate is halved at optimal doping. In the heavily doped regime, we find both strong electron-hole asymmetry and a persistent presence of Mott spectral features. Finally, we discuss the applicability of the single-band Hubbard model to describe the evolution of nodal spectra measured by angle-resolved photoemission spectroscopy (ARPES) on the single-layer cuprate La 2– x Sr x CuO 4 (0 ≤ x ≤ 0.15). This work benchmarks the predictive power of the Hubbard model for electronic properties of high- T c cuprates.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

A theoretical investigation of the hydrolysis of uranium hexafluoride: the initiation mechanism and vibrational spectroscopy

Depleted uranium hexafluoride (UF 6 ), a stockpiled byproduct of the nuclear fuel cycle, reacts readily with atmospheric humidity, but the mechanism is poorly understood. Here we compare several potential initiation steps at a consistent level of theory, generating underlying structures and vibrational modes using hybrid density functional theory (DFT) and computing relative energies of stationary points with double-hybrid (DH) DFT. A benchmark comparison is performed to assess the quality of DH-DFT data using reference energy differences obtained using a complete-basis-limit coupled-cluster (CC) composite method. The associated large-basis CC computations were enabled by a new general-purpose pseudopotential capability implemented as part of this work. Dispersion-corrected parameter-free DH-DFT methods, namely PBE0-DH-D3(BJ) and PBE-QIDH-D3(BJ), provided mean unsigned errors within chemical accuracy (1 kcal mol -1 ) for a set of barrier heights corresponding to the most energetically favorable initiation steps. The hydrolysis mechanism is found to proceed via intermolecular hydrogen transfer within van der Waals complexes involving UF 6 , UF 5 OH, and UOF 4 , in agreement with previous studies, followed by the formation of a previously unappreciated dihydroxide intermediate, UF 4 (OH) 2 . The dihydroxide is predicted to form under both kinetic and thermodynamic control, and, unlike the alternate pathway leading to the UO 2 F 2 monomer, its reaction energy is exothermic, in agreement with observation. Finally, harmonic and anharmonic vibrational simulations are performed to reinterpret literature infrared spectroscopy in light of this newly identified species.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Spinor-based coupled cluster thermochemistry: RgF and CnF 0/+ as compared to AuF and HgF 0/+

A relativistic coupled cluster approach that includes spin-orbit variationally in the eXact 2-component (X2C) approximation with the inclusion of the Gaunt interaction (X2Cg) was used to probe the thermochemistry and ground state spectroscopic constants of RgF, CnF, and CnF + . Utilizing large sequences of correlation consistent basis sets at the X2Cg-CCSD(T) level of theory, this work reports 0 K bond dissociation energies (BDEs) of AuF, RgF, HgF 0/+ , and CnF 0/+ , as well as ionization energies of Au, Rg, Hg, Cn, HgF, and CnF.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗