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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2022

Watershed zonation through hillslope clustering for tractably quantifying above- and below-ground watershed heterogeneity and functions

Abstract. In this study, we develop a watershed zonation approach for characterizing watershed organization and functions in a tractable manner by integrating multiple spatial data layers. We hypothesize that (1) a hillslope is an appropriate unit for capturing the watershed-scale heterogeneity of key bedrock-through-canopy properties and for quantifying the co-variability of these properties representing coupled ecohydrological and biogeochemical interactions, (2) remote sensing data layers and clustering methods can be used to identify watershed hillslope zones having the unique distributions of these properties relative to neighboring parcels, and (3) property suites associated with the identified zones can be used to understand zone-based functions, such as response to early snowmelt or drought and solute exports to the river. We demonstrate this concept using unsupervised clustering methods that synthesize airborne remote sensing data (lidar, hyperspectral, and electromagnetic surveys) along with satellite and streamflow data collected in the East River Watershed, Crested Butte, Colorado, USA. Results show that (1) we can define the scale of hillslopes at which the hillslope-averaged metrics can capture the majority of the overall variability in key properties (such as elevation, net potential annual radiation, and peak snow-water equivalent – SWE), (2) elevation and aspect are independent controls on plant and snow signatures, (3) near-surface bedrock electrical resistivity (top 20 m) and geological structures are significantly correlated with surface topography and plan species distribution, and (4) K-means, hierarchical clustering, and Gaussian mixture clustering methods generate similar zonation patterns across the watershed. Using independently collected data, we show that the identified zones provide information about zone-based watershed functions, including foresummer drought sensitivity and river nitrogen exports. The approach is expected to be applicable to other sites and generally useful for guiding the selection of hillslope-experiment locations and informing model parameterization.

58 GEOSCIENCES↗
DOE OSTI2015

Investigation of the CH{sub 3}Cl + CN{sup −} reaction in water: Multilevel quantum mechanics/molecular mechanics study

The CH{sub 3}Cl + CN{sup −} reaction in water was studied using a multilevel quantum mechanics/molecular mechanics (MM) method with the multilevels, electrostatic potential, density functional theory (DFT) and coupled-cluster single double triple (CCSD(T)), for the solute region. The detailed, back-side attack S{sub N}2 reaction mechanism was mapped along the reaction pathway. The potentials of mean force were calculated under both the DFT and CCSD(T) levels for the reaction region. The CCSD(T)/MM level of theory presents a free energy activation barrier height at 20.3 kcal/mol, which agrees very well with the experiment value at 21.6 kcal/mol. The results show that the aqueous solution has a dominant role in shaping the potential of mean force. The solvation effect and the polarization effect together increase the activation barrier height by ∼11.4 kcal/mol: the solvation effect plays a major role by providing about 75% of the contribution, while polarization effect only contributes 25% to the activation barrier height. Our calculated potential of mean force under the CCSD(T)/MM also has a good agreement with the one estimated using data from previous gas-phase studies.

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

Communication: Multireference equation of motion coupled cluster: A transform and diagonalize approach to electronic structure

The novel multireference equation-of-motion coupled-cluster (MREOM-CC) approaches provide versatile and accurate access to a large number of electronic states. The methods proceed by a sequence of many-body similarity transformations and a subsequent diagonalization of the transformed Hamiltonian over a compact subspace. The transformed Hamiltonian is a connected entity and preserves spin- and spatial symmetry properties of the original Hamiltonian, but is no longer Hermitean. The final diagonalization spaces are defined in terms of a complete active space (CAS) and limited excitations (1h, 1p, 2h, …) out of the CAS. The methods are invariant to rotations of orbitals within their respective subspaces (inactive, active, external). Applications to first row transition metal atoms (Cr, Mn, and Fe) are presented yielding results for up to 524 electronic states (for Cr) with an rms error compared to experiment of about 0.05 eV. The accuracy of the MREOM family of methods is closely related to its favorable extensivity properties as illustrated by calculations on the O{sub 2}–O{sub 2} dimer. The computational costs of the transformation steps in MREOM are comparable to those of closed-shell Coupled Cluster Singles and Doubles (CCSD) approach.

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

Highly Efficient and Scalable Compound Decomposition of Two-Electron Integral Tensor and Its Application in Coupled Cluster Calculations

The representation and storage of two-electron integral tensors are vital in large- scale applications of accurate electronic structure methods. Low-rank representation and efficient storage strategy of integral tensors can significantly reduce the numerical overhead and consequently time-to-solution of these methods. In this paper, by combining pivoted incomplete Cholesky decomposition (CD) with a follow-up truncated singular vector decomposition (SVD), we develop a decomposition strategy to approximately represent the two-electron integral tensor in terms of low-rank vectors. A systematic benchmark test on a series of 1-D, 2-D, and 3-D carbon-hydrogen systems demonstrates high efficiency and scalability of the compound two-step decomposition of the two-electron integral tensor in our implementation. For the size of atomic basis set N_b ranging from ~ 100 up to ~ 2, 000, the observed numerical scaling of our implementation shows O(N_b^{2.5~3}) versus O(N_b^{3~4}) of single CD in most of other implementations. More importantly, this decomposition strategy can significantly reduce the storage requirement of the atomic-orbital (AO) two-electron integral tensor from O(N_b^4) to O(N_b^2 log_{10}(N_b)) with moderate decomposition thresholds. The accuracy tests have been performed using ground- and excited-state formulations of coupled- cluster formalism employing single and double excitations (CCSD) on several bench- mark systems including the C_{60} molecule described by nearly 1,400 basis functions. The results show that the decomposition thresholds can be generally set to 10^{-4} to 10^{-3} to give acceptable compromise between efficiency and accuracy.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Few-nucleon and many-nucleon systems with semilocal coordinate-space regularized chiral nucleon-nucleon forces

We employ a variety of ab initio methods, including Faddeev-Yakubovsky equations, no-core configuration interaction approach, coupled-cluster theory, and in-medium similarity renormalization group, to perform a comprehensive analysis of the nucleon-deuteron elastic and breakup reactions and selected properties of light and medium-mass nuclei up to 48 C a using the recently constructed semilocal coordinate-space regularized chiral nucleon-nucleon potentials. We compare the results with those based on selected phenomenological and chiral EFT two-nucleon potentials, discuss the convergence pattern of the chiral expansion, and estimate the achievable theoretical accuracy at various chiral orders using an approach to quantify truncation errors of the chiral expansion without relying on cutoff variation. We also address the robustness of this method and explore alternative ways to estimate the theoretical uncertainty from the truncation of the chiral expansion.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Contaminant source identification using semi-supervised machine learning

Identification of the original groundwater types present in geochemical mixtures observed in an aquifer is a challenging but very important task. Frequently, some of the groundwater types are related to different infiltration and/or contamination sources associated with various geochemical signatures and origins. The characterization of groundwater mixing processes typically requires solving complex inverse models representing groundwater flow and geochemical transport in the aquifer, where the inverse analysis accounts for available site data. Usually, the model is calibrated against the available data characterizing the spatial and temporal distribution of the observed geochemical types. Numerous different geochemical constituents and processes may need to be simulated in these models which further complicates the analyses. In this paper, we propose a new contaminant source identification approach that performs decomposition of the observation mixtures based on Non-negative Matrix Factorization (NMF) method for Blind Source Separation (BSS), coupled with a custom semi-supervised clustering algorithm. Our methodology, called NMFk, is capable of identifying (a) the unknown number of groundwater types and (b) the original geochemical concentration of the contaminant sources from measured geochemical mixtures with unknown mixing ratios without any additional site information. NMFk is tested on synthetic and real-world site data. Finally, the NMFk algorithm works with geochemical data represented in the form of concentrations, ratios (of two constituents; for example, isotope ratios), and delta notations (standard normalized stable isotope ratios).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Analytic Gradients for Equation-of-Motion Coupled Cluster with Single, Double, and Perturbative Triple Excitations

Understanding the process of molecular photoexcitation is crucial in various fields, including drug development, materials science, photovoltaics, and more. The electronic vertical excitation energy is a critical property, for example in determining the singlet-triplet gap of chromophores. However, a full understanding of excited-state processes requires additional explorations of the excited-state potential energy surface and electronic properties, which is greatly aided by the availability of analytic energy gradients. Owing to its robust high accuracy over a wide range of chemical problems, equation-of-motion coupled-cluster with single and double excitations (EOM-CCSD) is a powerful method for predicting excited state properties, and the implementation of analytic gradients of many EOM-CCSD (excitation energies, ionization potentials, electron attachment energies, etc.) along with numerous successful applications high- lights the flexibility of the method. In specific cases where a higher level of accuracy is needed or in more complex electronic structures, the inclusion of triple excitations becomes essential, for example, in the EOM-CCSD* approach of Saeh and Stanton. In this work, we derive and implement for the first time the analytic gradients of EOMEE-CCSD*, which also provides a template for analytic gradients of related ex- cited state methods with perturbative triple excitations. Here, the capabilities of analytic EOMEE-CCSD* gradients are illustrated by several representative examples.

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

Modeling Oxidative Dehydrogenation of Propane with Supported Vanadia Catalysts Using Multireference Methods

The oxidative dehydrogenation of propane over supported vanadium oxide catalysts poses significant computational challenges due to complex electronic structure changes along the reaction coordinate, driven primarily by changes in the oxidation states of vanadium. To address these challenges, we systematically test quantum chemical methods, including multireference (MR) approaches, domain-based local pair natural orbital coupled cluster theory (DLPNO-CCSD(T)), and density functional theory (DFT). The initial C–H bond-breaking transition state requires MR treatment due to its multireference character, while subsequent steps permit efficient single-reference calculations. For the rate-limiting C–H activation step mediated by the vanadyl moiety, complete 1 active space second-order perturbation theory (CASPT2) yields an apparent activation barrier (E app 600K ) of 138 kJ/mol, consistent with experimental values (134 ± 4 kJ/mol; Gruene et al. Catal. Today 2010, 157, 137). In contrast, DLPNO-CCSD(T) overestimates this barrier (198 kJ/mol), whereas DFT predictions span 125–150 kJ/mol, depending on the functional. Our multireference investigation of this transition metal oxide-catalyzed process demonstrates that an active space that incorporates the C–H σ and V=O σ/π bonding orbitals, oxygen lone pairs, and their antibonding counterparts adequately captures electronic structure changes along the chemical transformation. Furthermore, these findings provide a general strategy for active space selection in transition metal oxide-catalyzed C/O–H bond activation reactions. The reference dataset from this work, which includes MR calculations with manually selected active spaces for all intermediates and transition states in the propane ODH reaction network, will serve as a benchmark for automating active space selection in similar systems.

Catalysts↗
DOE OSTI2014

Communication: A combined periodic density functional and incremental wave-function-based approach for the dispersion-accounting time-resolved dynamics of {sup 4}He nanodroplets on surfaces: {sup 4}He/graphene

In this work we propose a general strategy to calculate accurate He–surface interaction potentials. It extends the dispersionless density functional approach recently developed by Pernal et al. [Phys. Rev. Lett. 103, 263201 (2009)] to adsorbate-surface interactions by including periodic boundary conditions. We also introduce a scheme to parametrize the dispersion interaction by calculating two- and three-body dispersion terms at coupled cluster singles and doubles and perturbative triples (CCSD(T)) level via the method of increments [H. Stoll, J. Chem. Phys. 97, 8449 (1992)]. The performance of the composite approach is tested on {sup 4}He/graphene by determining the energies of the low-lying selective adsorption states, finding an excellent agreement with the best available theoretical data. Second, the capability of the approach to describe dispersionless correlation effects realistically is used to extract dispersion effects in time-dependent density functional simulations on the collision of {sup 4}He droplets with a single graphene sheet. It is found that dispersion effects play a key role in the fast spreading of the {sup 4}He nanodroplet, the evaporation-like process of helium atoms, and the formation of solid-like helium structures. These characteristics are expected to be quite general and highly relevant to explain experimental measurements with the newly developed helium droplet mediated deposition technique.

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

Spectral function for 4 He using the Chebyshev expansion in coupled-cluster theory

Here, we compute spectral function for 4 He by combining coupled-cluster theory with an expansion of integral transforms into Chebyshev polynomials. Our method allows us to estimate the uncertainty of spectral reconstruction. The properties of the Chebyshev polynomials make the procedure numerically stable and considerably lower in memory usage than the typically employed Lanczos algorithm. We benchmark our predictions with other calculations in the literature and with electron-scattering data in the quasi-elastic peak. The spectral function formalism allows one to extend ab initio lepton-nucleus cross sections into the relativistic regime. This makes it a promising tool for modeling this process at higher-energy transfers. The results we present open the door for studies of heavier nuclei, important for the neutrino oscillation programs.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1980

Neutrinos and the age of the universe

The age of the universe should be calculable by independent methods with similar results. Previous calculations using nucleochronometers, globular clusters and dynamical measurements coupled with Friedmann models and nucleosynthesis constraints have given different values of the age. A consistent age is reported, whose implications for the constituent mass density are very interesting and are affected by the existence of a third neutrino flavor, and by allowing the possibility that neutrinos may have a non-zero rest mass.

Symbalisty, E. M. D.↗
NASA NTRS1994

The Molecular Structure of cis-FONO

The molecular structure of cis-FONO has been determined with the CCSD(T) correlation method using an spdf quality basis set. In agreement with previous coupled-cluster calculations but in disagreement with density functional theory, cis-FONO is found to exhibit normal bond distances. The quadratic and cubic force fields of cis-FONO have also been determined in order to evaluate the effect of vibrational averaging on the molecular geometry. Vibrational averaging is found to increase bond distances, as expected, but it does not affect the qualitative nature of the bonding. The CCSD(T)/spdf harmonic frequencies of cis-FONO support our previous assertion that a band observed at 1200 /cm is a combination band (upsilon(sub 3) + upsilon(sub 4)), and not a fundamental.

Lee, Timothy J.↗
DOE OSTI2022

Fathoming the anisotropic magnetoelasticity and magnetocaloric effect in GdNi

Intermetallic GdNi adopts a CrB type of crystal structure (space group Cmcm), and it orders ferromagnetically via a second-order phase transition at 70 K, exhibiting unusually strong spontaneous striction along the three independent crystallographic axes in the ferromagnetically ordered state. Here we introduce a microscopic model to describe anisotropic changes of lattice parameters and elastic contribution to magnetocaloric effect of GdNi. In the model, results of density functional theory (DFT) calculations are used as inputs into a Hamiltonian that includes elastic energy of an anisotropic crystal lattice, exchange interactions, and Zeeman effect. The magnetic and elastic Hamiltonians are coupled through an anisotropic Bean-Rodbell model of magnetoelastic interactions. This coupling gives rise to anisotropic changes in the lattice parameters observed experimentally, and the model reveals good to reasonable agreements between the current theoretical results and earlier experimental data, thus validating the model within the limits of assumptions made. We also show that DFT calculations with 4f electrons of Gd treated as core electrons lead to a more adequate estimate of elastic constants of GdNi in comparison with the LDA + U method where 4f electrons are treated as valence electrons.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Properties of the water to boron nitride interaction: From zero to two dimensions with benchmark accuracy

Molecular adsorption on surfaces plays an important part in catalysis, corrosion, desalination, and various other processes that are relevant to industry and in nature. As a complement to experiments, accurate adsorption energies can be obtained using various sophisticated electronic structure methods that can now be applied to periodic systems. The adsorption energy of water on boron nitride substrates, going from zero to 2-dimensional periodicity, is particularly interesting as it calls for an accurate treatment of polarizable electrostatics and dispersion interactions, as well as posing a practical challenge to experiments and electronic structure methods. Here, we present reference adsorption energies, static polarizabilities, and dynamic polarizabilities, for water on BN substrates of varying size and dimension. Adsorption energies are computed with coupled cluster theory, fixed-node quantum Monte Carlo (FNQMC), the random phase approximation, and second order Møller-Plesset theory. These wavefunction based correlated methods are found to agree in molecular as well as periodic systems. The best estimate of the water/h-BN adsorption energy is -107 ± 7 meV from FNQMC. In addition, the water adsorption energy on the BN substrates could be expected to grow monotonically with the size of the substrate due to increased dispersion interactions, but interestingly, this is not the case here. This peculiar finding is explained using the static polarizabilities and molecular dispersion coefficients of the systems, as computed from time-dependent density functional theory (DFT). Dynamic as well as static polarizabilities are found to be highly anisotropic in these systems. Finally, the many-body dispersion method in DFT emerges as a particularly useful estimation of finite size effects for other expensive, many-body wavefunction based methods.

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

Snapshots of Proton Accommodation at a Microscopic Water Surface: Understanding the Vibrational Spectral Signatures of the Charge Defect in Cryogenically Cooled H+(H2O)n=2 – 28 Clusters

In this Article, we review the role of gas-phase, size-selected protonated water clusters, H+(H2O)n, in the analysis of the microscopic mechanics responsible for the behavior of the excess proton in bulk water. We extend upon previous studies of the smaller, two-dimensional sheet-like structures to larger (n≥10) assemblies with three-dimensional cage morphologies which better mimic the bulk environment. Indeed, clusters in which a complete second solvation shell forms around a surface-embedded hydronium ion yield vibrational spectra where the signatures of the proton defect display strikingly similar positions and breadth to those observed in dilute acids. We investigate effects of the local structure and intermolecular interactions on the large red shifts observed in the proton vibrational signature upon cluster growth using various theoretical methods. We show that, in addition to sizeable anharmonic couplings, the position of the excess proton vibration can be traced to large increases in the electric field exerted on the embedded hydronium ion upon formation of the first and second solvation shells. MAJ acknowledges support from the U.S. Department of Energy under Grant No. DE-FG02- 06ER15800 as well as the facilities and staff of the Yale University Faculty of Arts and Sciences High Performance Computing Center, and by the National Science Foundation under Grant No. CNS 08-21132 that partially funded acquisition of the facilities. SMK and SSX acknowledge support from the US Department of Energy, Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences and Biosciences. Pacific Northwest National Laboratory (PNNL) is a multiprogram national laboratory operated for DOE by Battelle. This research used resources of the National Energy Research Scientific Computing Center, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231.

Fournier, Joseph A.↗