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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

PyCCE: A python package for custer correlation expansion simulations of spin qubit dynamics.

PyCCE, an open-source Python library for simulating the dynamics of spin qubits in a spin bath, is presented using the cluster-correlation expansion (CCE) method. PyCCE includes modules to generate realistic spin baths, employing coupling parameters computed from first principles with electronic structure codes, and enables the user to run simulations with either the conventional or generalized CCE method. Three use cases of the Python library are illustrated: the calculation of the Hahn-echo coherence time of the nitrogen-vacancy in diamond; the calculation of the coherence time of the basal divacancy in silicon carbide at avoided crossings; and the calculation for magnetic field orientation-dependent dynamics of a shallow donor in silicon. The complete documentation, downloadable tutorials, and installation instructions are available at https://pycce.readthedocs.io/en/latest/.

coherence↗
DOE OSTI2021

Synthesis and Photophysical Properties of Light-Harvesting Gold Nanoclusters Fully Functionalized with Antenna Chromophores

The development of efficient light-harvesting systems is important to understand the key aspects of solar-energy conversion processes and to utilize them in various photonic applications. In this study, atomically well-defined gold nanoclusters are reported as a new platform to fabricate artificial light-harvesting systems. An efficient amide coupling method is developed to synthesize water-soluble Au 22 clusters fully protected with pyrene chromophores by taking advantage of their facile phase-transfer reaction. The synthesized Au 22 clusters with densely packed 18 pyrene chromophores (Au 22 –PyB 18 ) exhibit triple-emission in blue, green, and red wavelength regions arising respectively from pyrene monomer, pyrene excimer, and Au 22 emission, producing bright white light emission together. The photoluminescence of Au 22 is enhanced by more than tenfold, demonstrating that pyrenes at the periphery efficiently channel the absorbed energy to the luminescent Au 22 at the center. A combination of femtosecond transient absorption and anisotropy measurements of Au 22 –PyB 18 explicitly reveals three main decay components of 220 fs, 3.5 ps, and 160 ps that can be assigned to energy migration between pyrenes and energy transfer processes from pyrene monomer and excimer to the central Au 22 , respectively.

14 SOLAR ENERGY↗
DOE OSTI2022

Quadratic Unitary Coupled-Cluster Singles and Doubles Scheme: Efficient Implementation, Benchmark Study, and Formulation of an Extended Version

An efficient implementation of the quadratic unitary coupled-cluster singles and doubles (qUCCSD) scheme for calculations of electronic ground and excited states using an unrestricted molecular spin–orbital formulation and an efficient tensor contraction library is reported. The accuracy of the qUCCSD scheme and the efficiency of the present implementation are demonstrated using extensive benchmark calculations of excitation energies and an application to S o → S 1 vertical excitation energies for cis- and trans-4a,4b-dihydrotriphenylene. The qUCCSD scheme has been shown to provide improved excitation energies compared with the UCC3 scheme formulated based on perturbation theory. Here, a UCC truncation scheme that can provide excitation energies correct through the fourth order is also presented to further improve the accuracy of the qUCCSD scheme.

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

Relaxed active space: Fixing tailored-CC with high order coupled cluster. II

Due to the steep increase in computational cost with the inclusion of higher-connected cluster operators in coupled-cluster applications, it is usually not practical to use such methods for larger systems or basis sets without an active space partitioning. This study generates an active space subject to unambiguous statistical criteria to define a space whose size permits treatment at the CCSDT level. The automated scheme makes it unnecessary for the user to judge whether a chosen active space is sufficient to correctly solve the problem. Two demanding applications are presented: twisted ethylene and the transition states for the bicyclo[1,1,0]butane isomerization. As bi-radicals both systems require at least a CCSDT level of theory for quantitative results, for the geometries and energies.

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

Effective many-body interactions in reduced-dimensionality spaces through neural network models

Accurately describing properties of challenging problems in physical sciences often requires complex mathematical models that are unmanageable to tackle head on. Therefore, developing reduced-dimensionality representations that encapsulate complex correlation effects in many-body systems is crucial to advance the understanding of these complicated problems. However, a numerical evaluation of these predictive models can still be associated with a significant computational overhead. To address this challenge, in this paper we discuss a combined framework that integrates recent advances in the development of active-space representations of coupled cluster (CC) downfolded Hamiltonians with neural network approaches. The primary objective of this effort is to train neural networks to eliminate the computationally expensive steps required for evaluating hundreds or thousands of Hugenholtz diagrams, which correspond to multidimensional tensor contractions necessary for evaluating a many-body form of downfolded effective Hamiltonians. Using small molecular systems (the H 2 O and HF molecules) as examples, we demonstrate that training neural networks employing effective Hamiltonians for a few nuclear geometries of molecules can accurately interpolate or extrapolate their forms to other geometrical configurations characterized by different intensities of correlation effects. We also discuss differences between effective interactions that define CC downfolded Hamiltonians with those of bare Hamiltonians defined by Coulomb interactions in the active spaces. Published by the American Physical Society 2024

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

RT-EOM-CCSD Calculations of Inner and Outer Valence Ionization Energies and Spectral Functions

Photoelectron spectroscopy (PES) is a standard experimental method for material characterization, but its interpretation can be hampered by its reliance on standard materials. To facilitate the study of unknown systems, theoretical methods are desirable. Here we present a real-time equation-of-motion coupled cluster (RT-EOM-CC) approach for valence PES, extending our core-level development. Here, we demonstrate that RT-EOM-CC yields ionization energies and spectral functions in good agreement with experiment and CI-based methods, even for some more correlated cases.

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

Modeling Absolute Redox Potentials of Ferrocene in the Condensed Phase

Absolute thermodynamic quantities for critical chemical reactions are needed to determine the role of solvents and reactive environments in catalysis and electrocatalysis beyond the relative scales typically employed. In principle, theoretical methods can provide such quantification but are often hindered by the innate complexity of strong electron correlation and dynamic relaxation of solvent environments. Here, we present and validate a protocol for calculating the redox potentials of ferrocene/ferrocenium redox pair in the acetonitrile. Equation-of-motion ionization potential coupled-cluster single-double (EOM-IP-CCSD) and effective fragment potential (EFP) methods are used to characterize the adiabatic and vertical ionization potentials (IP) as well as the electron affinity processes. We benchmark molecular mechanics against the EFP model to show the differences in ferrocene electronic polarizability in two redox states. Our best estimate of the redox potential (4.94 eV) agrees well with the experimental value (4.93 eV). This demonstrated the ability of modern computational methods to predict absolute redox potentials quantitatively and, more critically, quantify the correlation of dynamic effects, which underlie their origin.

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

The Generality of the GUGA MRCI Approach in COLUMBUS for Treating Complex Quantum Chemistry

The core part of the program system COLUMBUS allows highly efficient calculations using variational multireference (MR) methods in the framework of configuration interaction with single and double excitations (MR-CISD) and averaged quadratic coupled-cluster calculations (MR-AQCC), based on uncontracted sets of configurations and the graphical unitary group approach (GUGA). The availability of analytic MR-CISD and MR-AQCC energy gradients and analytic nonadiabatic couplings for MR-CISD enables exciting applications including, e.g., investigations of π-conjugated biradicaloid compounds, calculations of multitudes of excited states, development of diabatization procedures, and furnishing the electronic structure information for on-the-fly surface nonadiabatic dynamics. With fully variational uncontracted spin-orbit MRCI, COLUMBUS provides a unique possibility of performing high-level calculations on compounds containing heavy atoms up to lanthanides and actinides. Crucial for carrying out all of these calculations effectively is the availability of an efficient parallel code for the CI step. Configuration spaces of several billion in size now can be treated quite routinely on standard parallel computer clusters. Emerging developments in COLUMBUS, including the all configuration mean energy multiconfiguration self-consistent field method and the graphically contracted function method, promise to allow practically unlimited configuration space dimensions. Spin density based on the GUGA approach, analytic spin-orbit energy gradients, possibilities for local electron correlation MR calculations, development of general interfaces for nonadiabatic dynamics, and MRCI linear vibronic coupling models conclude this overview.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Spherical coupled-cluster theory for open-shell nuclei

A microscopic description of nuclei is important to understand the nuclear shell model from fundamental principles. This is difficult to achieve for more than the lightest nuclei without an effective approximation scheme. Define and evaluate an approximation scheme that can be used to study nuclei that are described as two particles attached to a closed (sub-)shell nucleus. The equation-of-motion coupled-cluster formalism has been used to obtain ground- and excited-state energies. This method is based on the diagonalization of a non-Hermitian matrix obtained from a similarity transformation of the many-body nuclear Hamiltonian. A chiral interaction at the next-to-next-to-next-to leading order (N 3 LO) using a cutoff at 500 MeV was used. The ground-state energies of 6 Li and 6 He were in good agreement with a no-core shell-model calculation using the same interaction. Several excited states were also produced with overall good agreement. Only the $J^π = 3^+$ excited state in 6 Li showed a sizable deviation. The ground-state energies of 18 O, 18 F, and 18 Ne were converged but underbound compared to experiment. Moreover, the calculated spectra were converged and comparable to both experiment and shell-model studies in this region. Some excited states in 18 O were high or missing in the spectrum. It was also shown that the wave function for both ground and excited states separates into an intrinsic part and a Gaussian for the center-of-mass coordinate. Spurious center-of-mass excitations are clearly identified. Results are converged with respect to the size of the model space and the method can be used to describe nuclear states with simple structure. Especially the ground-state energies were very close to what has been achieved by exact diagonalization. To obtain a closer match with experimental data, effects of three-nucleon forces, the scattering continuum, as well as additional configurations in the coupled-cluster approximations are necessary.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2013

Local unitary transformation method toward practical electron correlation calculations with scalar relativistic effect in large-scale molecules

In order to perform practical electron correlation calculations, the local unitary transformation (LUT) scheme at the spin-free infinite-order Douglas–Kroll–Hess (IODKH) level [J. Seino and H. Nakai, J. Chem. Phys.136, 244102 (2012); J. Seino and H. Nakai, J. Chem. Phys.137, 144101 (2012)], which is based on the locality of relativistic effects, has been combined with the linear-scaling divide-and-conquer (DC)-based Hartree–Fock (HF) and electron correlation methods, such as the second-order Møller–Plesset (MP2) and the coupled cluster theories with single and double excitations (CCSD). Numerical applications in hydrogen halide molecules, (HX){sub n} (X = F, Cl, Br, and I), coinage metal chain systems, M{sub n} (M = Cu and Ag), and platinum-terminated polyynediyl chain, trans,trans-((p-CH{sub 3}C{sub 6}H{sub 4}){sub 3}P){sub 2}(C{sub 6}H{sub 5})Pt(C≡C){sub 4}Pt(C{sub 6}H{sub 5})((p-CH{sub 3}C{sub 6}H{sub 4}){sub 3}P){sub 2}, clarified that the present methods, namely DC-HF, MP2, and CCSD with the LUT-IODKH Hamiltonian, reproduce the results obtained using conventional methods with small computational costs. The combination of both LUT and DC techniques could be the first approach that achieves overall quasi-linear-scaling with a small prefactor for relativistic electron correlation calculations.

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

Coupled cluster benchmarks of large noncovalent complexes: The L7 dataset as well as DNA–ellipticine and buckycatcher–fullerene

In this work, benchmark binding energies for dispersion-bound complexes in the L7 dataset, the DNA–ellipticine intercalation complex, and the buckycatcher–C60 complex with 120 heavy atoms using a focal-point method based on the canonical form of second-order Møller−Plesset theory (MP2) and the domain based local pair natural orbital scheme for the coupled cluster with single, double, and perturbative triple excitations [CCSD(T)] extrapolated to the complete basis set (CBS) limit are reported. This work allows for increased confidence given the agreement with respect to values recently obtained using the local natural orbital CCSD(T) for L7 and the canonical CCSD(T)/CBS result for the coronene dimer (C2C2PD). Therefore, these results can be considered pushing the CCSD(T)/CBS binding benchmark to the hundred-atom scale. The disagreements between the two state-of-the-art methods, CCSD(T) and fixed-node diffusion Monte Carlo, are substantial with at least 2.0 (∼10%), 1.9 (∼5%), and 10.3 kcal/mol (∼25%) differences for C2C2PD in L7, DNA–ellipticine, and buckycatcher–C60, respectively. Such sizable discrepancy above “chemical accuracy” for large noncovalent complexes indicates how challenging it is to obtain benchmark binding interactions for systems beyond small molecules, although the three up-to-date density functionals, PBE0+D4, ωB97M-V, and B97M-V, agree better with CCSD(T) for these large systems. In addition to reporting these values, different basis sets and various CBS extrapolation parameters for Hartree–Fock and MP2 correlation energies were tested for the first time in large noncovalent complexes with the goal of providing some indications toward optimal cost effective routes to approach the CBS limit without substantial loss in quality.

Ballesteros, Francisco (ORCID:0000000325732307)↗
DOE OSTI2019

Nonempirical Interactions for the Nuclear Shell Model: An Update

The nuclear shell model has perhaps been the most signifcant conceptual and computational paradigm for the understanding of the structure of atomic nuclei. While the shell model has been used predominantly in a phenomenological context, there have been efforts stretching back more than half a century to derive shell model parameters based on a realistic interaction between nucleons. More recently, several ab initio many-body methods—in particular, many-body perturbation theory, the no-core shell model, the in-medium similarity renormalization group, and coupled-cluster theory—have developed the capability to provide effective shell model Hamiltonians. We offer an update on the status of these methods and investigate the connections between them and their potential strengths and weaknesses, with a particular focus on the in-medium similarity renormalization group approach. Three-body forces are demonstrated to be important for understanding the modifications needed in phenomenological treatments. We then review some applications of these methods to comparisons with recent experimental measurements, and conclude with some remaining challenges in ab initio shell model theory.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Coupled Cluster Green's function formulations based on the effective Hamiltonians

In this work, we demonstrate that the effective Hamiltonians obtained with the downfolding procedure based on double unitary coupled cluster (DUCC) ansatz can be used in the context of Green’s function coupled cluster (GFCC) formalism to calculate spectral functions of molecular systems. This combined approach (DUCC-GFCC) provides a significant reduction of numerical effort and good agreement with the corresponding all-orbital GFCC methods in energy windows that are consistent with the choice of active space. These features are demonstrated on the example of two benchmark systems: H 2 O and N 2 , where DUCC-GFCC calculations were performed for active spaces of various sizes.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2019

What Levels of Coupled Cluster Theory Are Appropriate for Transition Metal Systems? A Study Using Near-Exact Quantum Chemical Values for 3d Transition Metal Binary Compounds

Transition metal compounds are traditionally considered to be challenging for standard quantum chemistry approximations like coupled cluster (CC) theory, which are usually employed to validate lower level methods like density functional theory (DFT). To explore this issue, we present a database of bond dissociation energies (BDEs) for 74 spin states of 69 diatomic species containing a 3d transition metal atom and a main group element, in the moderately sized def2-SVP basis. The presented BDEs appear to have an (estimated) 3σ error less than 1 kJ/mol relative to the exact solutions to the nonrelativistic Born-Oppenheimer Hamiltonian. These benchmark values were used to assess the performance of a wide range of standard single reference CC models, as the results should be beneficial for understanding the limitations of these models for transition metal systems. We find that interactions between metals and monovalent ligands like hydride and fluoride are well described by CCSDT. Similarly, CCSDTQ appears to be adequate for bonds between metals and nominally divalent ligands like oxide and sulfide. However, interactions with polyvalent ligands like nitride and carbide are more challenging, with even CCSDTQ(P)Λ yielding errors on the scale of a few kJ/mol. We also find that many perturbative and iterative approximations to higher order terms either yield disappointing results or actually worsen the performance relative to the baseline low level CC method, indicating that complexity does not always guarantee accuracy.

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

Potential energy surface of the CO{sub 2}–N{sub 2} van der Waals complex

Four-dimensional potential energy surface (4D-PES) of the atmospherically relevant CO{sub 2}–N{sub 2} van der Waals complex is generated using the explicitly correlated coupled cluster with single, double, and perturbative triple excitation (CCSD(T)-F12) method in conjunction with the augmented correlation consistent triple zeta (aug-cc-pVTZ) basis set. This 4D-PES is mapped along the intermonomer coordinates. An analytic fit of this 4D-PES is performed. Our extensive computations confirm that the most stable form corresponds to a T-shape structure where the nitrogen molecule points towards the carbon atom of CO{sub 2}. In addition, we located a second isomer and two transition states in the ground state PES of CO{sub 2}–N{sub 2}. All of them lay below the CO{sub 2} + N{sub 2} dissociation limit. This 4D-PES is flat and strongly anisotropic along the intermonomer coordinates. This results in the possibility of the occurrence of large amplitude motions within the complex, such as the inversion of N{sub 2}, as suggested in the recent spectroscopic experiments. Finally, we show that the experimentally established deviations from the C{sub 2v} structure at equilibrium for the most stable isomer are due to the zero-point out-of-plane vibration correction.

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

27 AlNMR chemical shift of Al(OH)$^{-}_{4}$ calculated from first principles: Assessment of error cancellation in chemically distinct reference and target systems

Predicting accurate nuclear magnetic resonance chemical shieldings relies upon cancellation of different types of errors between the theoretically calculated shielding constant of the analyte of interest and the reference. Often, the intrinsic error in computed shieldings due to basis sets, approximations in the Hamiltonian, description of the wave function, and dynamic effects is nearly identical between the analyte and reference, yet if the electronic structure or sensitivity to local environment differs dramatically, this cannot be taken for granted. Detailed prior work has examined the octahedral trivalent cation Al(H 2 O)$^{3+}_{6}$, accounting for ab initio intrinsic errors. However, the use of this species as a reference for the chemically distinct tetrahedral anion Al(OH)$^{-}_{4}$ requires an understanding of how these errors cancel, in order to define the limits of accurately predicting 27 Al chemical shielding in Al(OH)$^{-}_{4}$. In this work, we estimate the absolute shielding of the 27 Al nucleus in Al(OH)$^{-}_{4}$ at the coupled cluster level (515.1± 5.3 ppm). Shielding sensitivity to the choice of method approximation and atomic basis sets used has been evaluated. Solvent and thermal effects are assessed through ensemble averaging techniques using ab-initio molecular dynamics. The contribution of each type of intrinsic error is assessed for the Al(H 2 O)$^{3+}_{6}$ and Al(OH)$^{-}_{4}$ ions, revealing significant differences that fundamentally hamper the ability to accurately calculate the 27 Al chemical shift of Al(OH)$^{-}_{4}$ from first principles.

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

Core spectroscopy of oxazole

We have measured, analyzed and simulated the ground state valence photoelectron (PES), X-ray absorption (XAS), X-ray photoelectron (XPS) as well as normal and resonant Auger-Meitner electron (AES) spectra of oxazole at the carbon, oxygen and nitrogen K-edge in order to understand its electronic structure. Experimental data is compared to theoretical calculations performed at the coupled cluster, restricted active space perturbation theory to second order and time-dependent density functional levels of theory. We demonstrate that both N and O K-edge XAS are sensitive to the amount of dynamical electron correlation included in the theoretical description, and that for a complete description of XPS, additional orbital correlation and orbital relaxation effects need to be considered. The normal AES are dominated by a singlet excitation channel and well described by theory. Furthermore, the resonant AES, however, are more complicated. While the participator decay channels, dominating at higher kinetic energies, are well described by coupled cluster theory, spectator channels can only be described satisfactorily using a method that combines restricted active space perturbation theory to second order for the bound part and a one-center approximation for the continuum.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

A Zonal Volt/VAR Control Mechanism for High PV Penetration Distribution Systems

This paper presents a zonal Volt/VAR control scheme for regulating voltage in unbalanced 3-phase distribution systems using inverter-based resources (IBR). First, the dependency between nodal voltage changes and IBR reactive power injections is derived via voltage sensitivity studies. Then, a fast-incremental clustering method is used to divide the distribution circuit into weakly-coupled zones based on correlations between nodal voltage sensitivities. The weak-coupling allows the voltage to be regulated independently within each zone using a rule-based voltage controller to dispatch IBR for voltage corrections. Simulation results on actual distribution feeder show that the proposed zone-based Volt/VAR control method maintains system voltages within their operational limits while reducing the runtime of the Volt/VAR controller from tens of second to a couple of milliseconds compared with centralized, optimization-based Volt/VAR control methods.

Alrushoud, Asmaa↗