Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Coupled cluster”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Energetic and Spectroscopic Properties of Astrophysically Relevant MgC 4 H Radicals Using High-Level Ab Initio Calculations

Considering the importance of magnesium-bearing hydrocarbon molecules (MgC n H; n = 2, 4, and 6) in the carbon-rich circumstellar envelopes (e.g., IRC+10216), a total of 28 constitutional isomers of MgC 4 H have been theoretically investigated using density functional theory (DFT) and coupled-cluster methods. The zero-point vibrational energy corrected relative energies at the ROCCSD(T)/cc-pCVTZ level of theory reveal that the linear isomer, 1-magnesapent-2,4-diyn-1-yl (1, 2 Σ + ), is the global minimum geometry on the MgC 4 H potential energy surface. The latter has been detected both in the laboratory and in the evolved carbon star, IRC+10216. The calculated spectroscopic data for 1 match well with the experimental observations (error ~ 0.78%) which validates our theoretical methodology. Plausible isomerization processes happening among different isomers are examined using DFT and coupled-cluster methods. CASPT2 calculations have been performed for a few isomers exhibiting multireference characteristics. The second most stable isomer, 1-ethynyl-1λ 3 -magnesacycloprop-2-ene-2,3-diyl (2, 2 A 1 , μ = 2.54 D), is 146 kJ mol –1 higher in energy than 1 and possibly the next promising candidate to be detected in the laboratory or in the interstellar medium in future.

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 OSTI2025

Slimmer Geminals For Accurate F12 Electronic Structure Models

The Slater-type F12 geminal length scales originally tuned for the second-order Mo̷ller-Plesset F12 method are too large for higher-order F12 methods formulated using the SP (diagonal fixed-coefficient spin-adapted) F12 ansatz. The new geminal parameters reported herein reduce the basis set incompleteness errors (BSIEs) of absolute coupled-cluster singles and doubles F12 correlation energies by a significant─and increase with the cardinal number of the basis─margin. The effect of geminal reoptimization is especially pronounced for the cc-pVXZ-F12 basis sets (specifically designed for use with F12 methods) relative to their conventional aug-cc-pVXZ counterparts. The BSIEs of relative energies are less affected, but substantial reductions can be obtained, especially for atomization energies and ionization potentials with the cc-pVXZ-F12 basis sets. The new geminal parameters are therefore recommended for all applications of high-order F12 methods, such as coupled-cluster F12 methods and transcorrelated F12 methods.

Powell, Samuel R. [Virginia Polytechnic Inst. and ↗
DOE OSTI2021

Shallow-circuit variational quantum eigensolver based on symmetry-inspired Hilbert space partitioning for quantum chemical calculations

Development of resource-friendly quantum algorithms remains highly desirable for noisy intermediate-scale quantum computing. Based on the variational quantum eigensolver (VQE) with unitary coupled-cluster Ansatz, we demonstrate that partitioning of the Hilbert space made possible by the point-group symmetry of the molecular systems greatly reduces the number of variational operators by confining the variational search within a subspace. In addition, we found that instead of including all subterms for each excitation operator, a single-term representation suffices to reach required accuracy for various molecules tested, resulting in an additional shortening of the quantum circuit by a factor of 4–8. With these strategies, VQE calculations on a noise-free quantum simulator achieve energies within a few meVs of those obtained with the full unitary coupled-cluster Ansatz with single and double excitations for the H 4 -square, H 4 -chain, and H 6 -hexagon molecules, while the number of cnot gates, a measure of the quantum-circuit depth, is reduced by a factor of as large as 35. Furthermore, we introduced an efficient “score” parameter to rank the excitation operators, so that the operators causing larger energy reduction can be applied first. Using the H 4 square and H 4 chain as examples, We demonstrated on noisy quantum simulators that the first few variational operators can bring the energy within the chemical accuracy, while additional operators do not improve the energy since the accumulative noise outweighs the gain from the expansion of the variational Ansatz.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1991

A parallel vectorized implementation of triple excitations in CCSD(T) - Application to the binding energies of the AlH3, AlH2F, AlHF2 and AlF3 dimers

An efficient method for various noniterative estimates of connected triple excitations in coupled-cluster theory is outlined and related to a similar expression occurring in Moller-Plesset perturbation theory. The method is highly vectorized and capable of utilizing multiple processors on a shared-memory machine, leading to computational rates in excess of one billion floating-point operations per second on four processors of a CRAY Y-MP. Using the new procedure, the binding energies of the D(2h) diborane-type dimers of AlH3, AlH2F, AlHF2, and AlF3 have been determined to be 32, 40, 20, and 47 kcal/mol, respectively. For Al2F6, the correlation procedure includes 232 molecular orbitals and over 1.5 x 10 to the 6th single and double coupled-cluster amplitudes, effectively accounting for over 2 x 10 to the 9th connected triple excitations.

Rendell, Alistair P.↗
DOE OSTI2020

How good are polarizable and flexible models for water: Insights from a many-body perspective

We present a systematic analysis of state-of-the-art polarizable and flexible water models from a many-body perspective, with a specific focus on their ability to represent the Born–Oppenheimer potential energy surface of water from the gas to the liquid phase. Using coupled cluster data in the completed basis set limit as a reference, we examine the accuracy of the polarizable models in reproducing individual many-body contributions to interaction energies and harmonic frequencies of water clusters and compare their performance with that of MB-pol, an explicit many-body model that has been shown to correctly predict the properties of water across the entire phase diagram. Based on these comparisons, we use MB-pol as a reference to analyze the ability of the polarizable models to reproduce the energy landscape of liquid water under ambient conditions. We find that, while correctly reproducing the energetics of minimum-energy structures, the polarizable models examined in this study suffer from inadequate representations of many-body effects for distorted configurations. To investigate the role played by geometry-dependent representations of 1-body charge distributions in reproducing coupled cluster data for both interaction and many-body energies, we introduce a simplified version of MB-pol that adopts fixed atomic charges and demonstrate that the new model retains the same accuracy as the original MB-pol model. Based on the analyses presented in this study, we believe that future developments of both polarizable and explicit many-body models should continue in parallel and would benefit from synergistic efforts aimed at integrating the best aspects of the two theoretical/computational frameworks.

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

Jordan–Wigner Transformation for the Description of Strong Correlation in Fermionic Systems

Seniority is a useful way of organizing Hilbert space for strongly correlated systems. The exact zero-seniority wave function, doubly occupied configuration interaction (DOCI), provides accurate results (given the right orbitals) for many strongly correlated electronic systems but has a combinatorial computational cost. In many cases, pair coupled cluster doubles provide a polynomial-cost approximation that closely reproduces the energies of DOCI, but it breaks down in some cases and, as shown herein, it does not provide particularly good density matrices. In this article, we demonstrate that by using the Jordan–Wigner transformation to turn the seniority zero problem back into a Fermionic one, we can provide mean-field variational results of DOCI quality for the Hubbard model and a few small molecular dissociation examples, with polynomial cost, both for the energies and for density matrices, all while being protected from collapse. This success is rooted in the proof we provide, showing that the Hartree–Fock wave function on the Jordan–Wigner-transformed Hamiltonian transforms back to variational coupled cluster doubles in the seniority zero representation, but restricted to have determinant rather than permanent amplitude coefficients, without compromising its overall accuracy.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS2011

Initiating Molecular Growth in the Interstellar Medium via Dimeric Complexes of Observed Ions and Molecules

A feasible initiation step for particle growth in the interstellar medium (ISM) is simulated by means of ab quantum chemistry methods. The systems studied are dimer ions formed by pairing nitrogen containing small molecules known to exist in the ISM with ions of unsaturated hydrocarbons or vice versa. Complexation energies, structures of ensuing complexes and electronic excitation spectra of the encounter complexes are estimated using various quantum chemistry methods. Moller-Plesset perturbation theory (MP2, Z-averaged perturbation theory (ZAP2), coupled cluster singles and doubles with perturbative triples corrections (CCSD(T)), and density functional theory (DFT) methods (B3LYP) were employed along with the correlation consistent cc-pVTZ and aug-cc-pVTZ basis sets. Two types of complexes are predicted. One type of complex has electrostatic binding with moderate (7-20 kcal per mol) binding energies, that are nonetheless significantly stronger than typical van der Waals interactions between molecules of this size. The other type of complex develops strong covalent bonds between the fragments. Cyclic isomers of the nitrogen containing complexes are produced very easily by ion-molecule reactions. Some of these complexes show intense ultraviolet visible spectra for electronic transitions with large oscillator strengths at the B3LYP, omegaB97, and equations of motion coupled cluster (EOM-CCSD) levels. The open shell nitrogen containing carbonaceous complexes especially exhibit a large oscillator strength electronic transition in the visible region of the electromagnetic spectrum.

Bera, Partha P.↗
DOE OSTI2020

Recent Developments in the General Atomic and Molecular Electronic Structure System

A discussion of many of the recently implemented features of GAMESS (General Atomic and Molecular Electronic Structure System) and LibCChem (the C++ CPU/GPU library associated with GAMESS) is presented. These features include fragmentation methods like the fragment molecular orbital, effective fragment potential and effective fragment molecular orbital methods, hybrid MPI/OpenMP approaches to Hartree-Fock and resolution of the identity second order perturbation theory. Many new coupled cluster theory methods have been implemented in GAMESS, as have multiple levels of density functional/tight binding theory. The role of accelerators, especially graphical processing units, is discussed in the context of the new features of LibCChem, as is the associated problem of power consumption as the power of computers increases dramatically. The process by which a complex program suite like GAMESS is maintained and developed is considered. Future developments are briefly summarized.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Thermofield Theory for Finite-Temperature Electronic Structure

Wave function methods have offered a robust, systematically improvable means to study ground-state properties in quantum many-body systems. Theories like coupled cluster and their derivatives provide highly accurate approximations to the energy landscape at a reasonable computational cost. Analogues of such methods to study thermal properties, though highly desirable, have been lacking because evaluating thermal properties involve a trace over the entire Hilbert space, which is a formidable task. Besides, excited-state theories are generally not as well studied as ground-state ones. In this mini-review, we present an overview of a finite-temperature wave function formalism based on thermofield dynamics to overcome these difficulties. Thermofield dynamics allows us to map the equilibrium thermal density matrix to a pure state, i.e., a single wave function, albeit in an expanded Hilbert space. Ensemble averages become expectation values over this so-called thermal state. Around this thermal state, we have developed a procedure to generalize ground-state wave function theories to finite temperatures. As explicit examples, we highlight formulations of mean-field, configuration interaction, and coupled cluster theories for thermal properties of Fermions in the grand-canonical ensemble. To assess the quality of these approximations, we also show benchmark studies for the one-dimensional Hubbard model, while comparing against exact results. We will see that the thermal methods perform similarly to their ground-state counterparts, while merely adding a prefactor to the asymptotic computational cost. Furthermore, they also inherit all the properties, good or bad, from the ground-state methods, signifying the robustness of our formalism and the scope for future development.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1995

The Heat of Formation of HNO

The HNO molecule is of interest in both combustion and atmospheric chemistry. For example, Guadagnini et al. have recently presented ab initio potential energy surfaces for the three lowest lying electronic states of HNO and then used these in examining several chemical reactions that take place in the combustion of nitrogen containing fuels and in the oxidation of atmospheric nitrogen. We have previously studied the ground state potential energy surface (i.e., stationary points along the HNO rev. reaction HON path), vibrational spectrum (using an accurate quartic force field), zero-point energy, and bonding of HNO using coupled-cluster ab initio methods. HNO is also very interesting because of the unique nature of its bonding characteristics. That is, the potential energy surface is very flat along the H-N bonding coordinate thereby giving unusual harmonic and fundamental vibrational frequencies, and the H-N bond energy is rather weak in comparison to other H-N bond energies. In fact, using experimental heats of formation for HO, H, and NO, the H- bond energy is computed to be only 49.9 kcal/ mol (298 K). However, ab initio calculations of isodesmic reaction energies involving HNO, FNO, ClNO, and several other molecules have shown that there is an inconsistency in the experimental heats of formation of the XNO (X=H, F, and Cl) species. Hence the motivation for this study was to determine a very accurate(DELTA)H(sup o)(sub f) value for HNO using state of-the-art ab initio methods. Based on many recent studies it is evident that the singles and doubles coupled-cluster method that includes a perturbational estimate of the effects of connected triple excitations, denoted CCSD(T), in conjunction with large one-particle basis sets should be reliable to better than +0.8 kcal/mol for this quantity. The computational methodology is described in the next section followed by our results and discussion. Conclusions are presented in the final section.

Lee, Timothy J.↗
NASA NTRS1995

The Heat of Formation of HNO

The HNO molecule is of interest in both combustion and atmospheric chemistry. For example, Guadagnini et al. have recently presented ab initio potential energy surfaces for the three lowest lying electronic states of HNO and then used these in examining several chemical reactions that take place in the combustion of nitrogen containing fuels and in the oxidation of atmospheric nitrogen. We have previously studied the ground state potential energy surface (i.e., stationary points along the HNO reversible reaction HON path), vibrational spectrum (using an accurate quartic force field), zero-point energy, and bonding of HNO using coupled-cluster ab initio methods. HNO is also very interesting because of the unique nature of its bonding characteristics. That is, the potential energy surface is very flat along the H-N bonding coordinate thereby giving unusual harmonic and fundamental vibrational frequencies, and the H-N bond energy is rather weak in comparison to other H-N bond energies. In fact, using experimental heats of formation for HNO, H, and NO, the H- bond energy is computed to be only 49.9 kcal/ mol (298 K). However, ab initio calculations of isodesmic reaction energies involving HNO, FNO, ClNO, and several other molecules have shown that there is an inconsistency in the experimental heats of formation of the XNO (X double bond H, F, and Cl) species. Hence the motivation for this study was to determine a very accurate (Delta)H(sub f, sup o) value for HNO using state-of-the-art ab initio methods. Based on many recent studies it is evident that the singles and doubles coupled-cluster method that includes a perturbational estimate of the effects of connected triple excitations, denoted CCSD(T), in conjunction with large one-particle basis sets should be reliable to better than +/- 0.8 kcal/mol for this quantity. The computational methodology is described in the next section followed by our results and discussion. Conclusions are presented in the final section.

Lee, Timothy J.↗
DOE OSTI2026

Microscopic Dynamics Controls Coupling and Cluster Formation in Brush Particle Solids

Thermodynamics-based models predict the structure of polymer-grafted nanoparticles (PGNs) as well as their assembly behavior based on geometric parameters such as particle size, degree of polymerization, and density of grafted chains. The role of microscopic polymer dynamics, such as the mobility of repeat units in the melt state, in the evolution of the structure and properties remains unknown. Brillouin light spectroscopy (BLS), due to its capability to concurrently discern the local and global elastic properties of PGN assemblies, enables the probing of microscopic processes, such as brush interdigitation, sensitive to the annealing of the assembly. For poly(methyl methacrylate) (PMMA)-grafted silica (SiO 2 ) PGNs in the dry powder state and annealed above the glass transition temperature, BLS revealed fully reversible local elasticity, indicative of limited interdigitation between adjacent PGNs. This contrasts with polystyrene (PS)−SiO 2 analogs that displayed ready (and irreversible) fusion of brush layers during annealing. The retardation of brush interdigitation in PMMA-grafted systems is surprising, given the similar thermomechanical properties of both polymers, and is rationalized as the consequence of higher friction between PMMA repeats compared to PS. Microscopic dynamics thus has a profound impact on the kinetic path of structure (and property) evolution and thus should be considered during the processing of PGNs into functional hybrid materials.

chemical structure↗
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 OSTI2020

TURBOMOLE: Modular program suite for ab initio quantum-chemical and condensed-matter simulations

TURBOMOLE is a collaborative, multi-national software development project aiming to provide highly efficient and stable computational tools for quantum chemical simulations of molecules, clusters, periodic systems, and solutions. The TURBOMOLE software suite is optimized for widely available, inexpensive, and resource-efficient hardware such as multi-core workstations and small computer clusters. TURBOMOLE specializes in electronic structure methods with outstanding accuracy–cost ratio, such as density functional theory including local hybrids and the random phase approximation (RPA), GW-Bethe–Salpeter methods, second-order Møller–Plesset theory, and explicitly correlated coupled-cluster methods. TURBOMOLE is based on Gaussian basis sets and has been pivotal for the development of many fast and low-scaling algorithms in the past three decades, such as integral-direct methods, fast multipole methods, the resolution-of-the-identity approximation, imaginary frequency integration, Laplace transform, and pair natural orbital methods. This review focuses on recent additions to TURBOMOLE’s functionality, including excited-state methods, RPA and Green’s function methods, relativistic approaches, high-order molecular properties, solvation effects, and periodic systems. A variety of illustrative applications along with accuracy and timing data are discussed. Moreover, available interfaces to users as well as other software are summarized. TURBOMOLE’s current licensing, distribution, and support model are discussed, and an overview of TURBOMOLE’s development workflow is provided. Challenges such as communication and outreach, software infrastructure, and funding are highlighted.

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

Excited states in variational Monte Carlo using a penalty method

Here, in this article, the authors present a technique using variational Monte Carlo to solve for excited states of electronic systems. This technique is based on enforcing orthogonality to lower energy states, which results in a simple variational principle for the excited states. Energy optimization is then used to solve for the excited states. This technique is applied to the well-characterized benzene molecule, in which ~10 000 parameters are optimized for the first 12 excited states. Agreement within ~0.2 eV is obtained with higher scaling coupled cluster methods; small disagreements with experiment are likely due to vibrational effects.

36 MATERIALS SCIENCE↗