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 343 records · Page 19

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Coupled cluster theory on modern heterogeneous supercomputers

This study examines the computational challenges in elucidating intricate chemical systems, particularly through ab-initio methodologies. This work highlights the Divide-Expand-Consolidate (DEC) approach for coupled cluster (CC) theory—a linear-scaling, massively parallel framework—as a viable solution. Detailed scrutiny of the DEC framework reveals its extensive applicability for large chemical systems, yet it also acknowledges inherent limitations. To mitigate these constraints, the cluster perturbation theory is presented as an effective remedy. Attention is then directed towards the CPS (D-3) model, explicitly derived from a CC singles parent and a doubles auxiliary excitation space, for computing excitation energies. The reviewed new algorithms for the CPS (D-3) method efficiently capitalize on multiple nodes and graphical processing units, expediting heavy tensor contractions. As a result, CPS (D-3) emerges as a scalable, rapid, and precise solution for computing molecular properties in large molecular systems, marking it an efficient contender to conventional CC models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Ab initio calculations of monopole sum rules: From finite nuclei to infinite nuclear matter

We compute moments of the isoscalar monopole response of 𝑁 = 𝑍 closed-shell nuclei based on chiral nucleon-nucleon plus three-nucleon interactions. We employ the random-phase approximation (RPA) and two ab initio many-body approaches, the in-medium similarity renormalization group (IMSRG) and coupled-cluster theory (CC). In the IMSRG framework, the moments are obtained as ground-state expectation values, whereas in the CC approach, they are evaluated through excited-state calculations. We find good agreement between the IMSRG and CC results across all nuclei studied. RPA provides a reasonable approximation to the correlated methods if the interaction is soft. From the calculated moments, we extract average energies of the monopole response, compute finite-nucleus incompressibilities, and estimate the incompressibility of symmetric nuclear matter by a fit to a leptodermous expansion. Our extrapolated values are lower than those obtained in nuclear-matter calculations with the same interactions, but the values are consistent with phenomenological ranges.

Bonaiti, Francesca [Michigan State Univ., East Lan↗

DOE OSTI2024

Continuous-variable quantum computation of the O(3) model in 1+1 dimensions

We formulate the $O(3)$ non-linear sigma model in $1+1$ dimensions as a limit of a three-component scalar field theory restricted to the unit sphere in the large squeezing limit. This allows us to describe the model in terms of the continuous variable (CV) approach to quantum computing. Here we construct the ground state and excited states using the coupled cluster ansatz and find excellent agreement with the exact diagonalization results for a small number of lattice sites. We then present the simulation protocol for the time evolution of the model using CV gates, estimate the discretization error, and present numerical results obtained from a photonic quantum simulator. We expect that the methods developed in this work will be useful for exploring interesting dynamics for a wide class of sigma models and gauge theories, as well as for simulating scattering events on quantum hardware in the coming decade.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Numerical simulations of noisy quantum circuits for computational chemistry

The opportunities afforded by near-term quantum computers to calculate the ground-state properties of small molecules depend on the structure of the computational ansatz as well as the errors induced by device noise. Here we investigate the behavior of these noisy quantum circuits using numerical simulations to estimate the accuracy and fidelity of the prepared quantum states relative to the ground truth obtained by conventional means. We implement several different types of ansatz circuits derived from unitary coupled cluster theory for the purposes of estimating the ground-state energy of sodium hydride using the variational quantum eigensolver algorithm. We show how relative error in the energy and the fidelity scale with the levels of gate-based noise, the internuclear configuration, the ansatz circuit depth, and the parameter optimization methods.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Morphological-informed Thermal Property Prediction in the Engineering Domain

Using first-principle, atomic scale methods to predict thermal carrier (electron and phonon) behavior at nanoscale yields highly resolved thermal properties. Density functional theory simulations can model myriad effects on the transport of phonons and electrons (e.g., phonon and electron coupling, electron correlation, defect presence, carrier scattering), and yield thermal properties for small geometric domains such as thermal and electrical conductivity, or heat capacity. While these thermal properties are accurate in atomistic and nanoscopic systems, they must be scaled up to the microstructural domain to be a useful predictor for experimental basis, as the microscale is where myriad changes and physical phenomenon occur within a material (e.g., grain boundaries, precipitate aggregations, interfaces, defect clusters). My talk will discuss the necessity of thermal property predictions for nuclear fuels applications and more. Additionally I will discuss the necessity in developing multi-physics, multi-scale methods for future applications in predictions for not only nuclear fuels, but other materials applications.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Block orbital-selective Mott insulators: A spin excitation analysis

We present a comprehensive study of the spin excitations—as measured by the dynamical spin structure factor S(q,ω)—of the so-called block-magnetic state of low-dimensional orbital-selective Mott insulators. We realize this state via both a multi-orbital Hubbard model and a generalized Kondo-Heisenberg Hamiltonian. Due to various competing energy scales present in the models, the system develops periodic ferromagnetic islands of various shapes and sizes, which are antiferromagnetically coupled. The 2×2 particular case was already found experimentally in the ladder material BaFe2Se3 that becomes superconducting under pressure. Here we discuss the electronic density as well as Hubbard and Hund coupling dependence of S(q,ω) using density matrix renormalization group method. Several interesting features were identified: (1) An acoustic (dispersive spin-wave) mode develops. (2) The spin-wave bandwidth establishes a new energy scale that is strongly dependent on the size of the magnetic island and becomes abnormally small for large clusters. (3) Optical (dispersionless spin excitation) modes are present for all block states studied here. In addition, a variety of phenomenological spin Hamiltonians have been investigated but none matches entirely our results that were obtained primarily at intermediate Hubbard U strengths. Our comprehensive analysis provides theoretical guidance and motivation to crystal growers to search for appropriate candidate materials to realize the block states, and to neutron scattering experimentalists to confirm the exotic dynamical magnetic properties unveiled here, with a rich mixture of acoustic and optical features.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2021

Numerical Simulations of Noisy Variational Quantum Eigensolver Ansatz Circuits

This is a case study of the variational quantum eigensolver (VQE) method using numerical simulations to test the influence of noise on the accuracy of the underlying circuit ansatz. We investigate a computational chemistry application of VQE to calculate the electronic ground state and its energy for Sodium Hydride (NaH), a prototypical two-electron problem. Using a one-parameter ansatz derived from unitary coupled cluster (UCC) theory, we simulate the effects of noise on the energy expectation value and variance with respect to the ansatz parameter. These numerical simulations provide insights into the accuracy of the prepared quantum state and the efficiency of the classical optimizer that iteratively refines the ansatz. We conduct a comparative study between analytical results derived for the UCC ansatz in the absence of noise and the noisy numerical simulation results obtained using an isotropic depolarizing noise model for each gate. We also compare the relative increase in noise on logically equivalent UCC ansatz circuits generated by randomized compiling. Notably, we observe that the intrinsic variance in the energy due to the simplicity of the ansatz itself compares with the noise induced by the bare circuit.

Gowrishankar, Meenambika↗

DOE OSTI2023

An Extension of the Athena++ Code Framework for Radiation-magnetohydrodynamics in General Relativity Using a Finite-solid-angle Discretization

We extend the general-relativistic magnetohydrodynamics (GRMHD) capabilities of Athena++ to incorporate radiation. The intensity field in each finite-volume cell is discretized in angle, with explicit transport in both space and angle properly accounting for the effects of gravity on null geodesics, and with matter and radiation coupled in a locally implicit fashion. Here we describe the numerical procedure in detail, verifying its correctness with a suite of tests. Motivated in particular by black hole accretion in the high-accretion-rate, thin-disk regime, we demonstrate the application of the method to this problem. With excellent scaling on flagship computing clusters, the port of the algorithm to the GPU-enabled AthenaK code now allows the simulation of many previously intractable radiation-GRMHD systems.

79 ASTRONOMY AND ASTROPHYSICS↗

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↗

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 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↗

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 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 OSTI2026

Computing nuclear response functions with time-dependent coupled-cluster theory

We compute nuclear response functions by solving the time-dependent 𝐴-body Schrödinger equation, recording the time-dependent transition moment and extracting spectral information via Fourier transforms. The solution of the time-dependent many-body problem accounts for correlations on top of the mean field by taking advantage of a time-dependent formulation of coupled-cluster theory. As a validation, we focus on electric dipole transitions in 4 He and 16 O and compare moments of the response function distribution to the results of an equivalent static framework, finding negligible discrepancies. We investigate how proton and neutron densities evolve in time, and we see the traditional picture of soft and giant dipole resonances as collective oscillations of protons and neutrons emerging from our calculations in 16 O and 24 O. Furthermore, this method also allows us to investigate the behavior of the nucleus in the presence of a strong electric field. In that regime, the behavior of the system becomes chaotic. Qualitatively, the spectral information obtained in this limit is in line with previous time-dependent mean-field results.

Ab initio calculations↗

DOE OSTI2026

Low-Temperature Non-Oxidative Coupling of Methane on Atomically Dispersed Titanium–Aluminum–Boron Nanopowder

Nonoxidative coupling of methane represents a long-standing challenge in heterogeneous catalysis, as it requires activation of the carbon–hydrogen (C–H) bond, controlled carbon–carbon (C–C) bond formation, and effective hydrogen management without relying on oxidants. Here, we report a low-temperature C–H activation and nonoxidative C–C coupling of methane over atomically dispersed titanium–aluminum–boron nanopowder (Ti–Al–B NP) utilizing a catalytic microreactor coupled to synchrotron single-photon photoionization reflectron time-of-flight mass spectrometry. The soft-ionization, in situ probing method detects the nascent reaction products and radical intermediates under operando conditions, including methyl radical, C2 hydrocarbons, and molecular hydrogen. Methane activation is initiated at 800 K, approximately 700 K below the gas-phase decomposition threshold, leading predominantly to ethylene formation with selectivity reaching up to 78% among the C–C coupled products. Electronic structure calculations on model Ti–Al–B clusters elucidate a cooperative catalytic mechanism in which titanium enables methane adsorption and C–H activation, boron acts as a reversible hydrogen reservoir, and aluminum stabilizes methylene intermediates, thereby facilitating selective C–C coupling and dehydrogenation. These findings establish a distinct catalyst architecture for nonoxidative methane coupling based on earth abundant elements alternative to expensive platinum and other noble metal-containing conventional catalysts and provide molecular-level design principles for controlling dehydrogenation and subsequent C–C bond formation in challenging light alkane conversions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗