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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Explicitly correlated coupled-cluster theory with Brueckner orbitals

Brueckner orbitals are the optimal orbitals for use in F12 explicitly correlated coupled-cluster (CC) treatments. A novel approach, Brueckner coupled-cluster doubles with perturbative triples BCCD(T){sub (F12*)} is presented that avoids the expensive re-evaluation of F12 integrals throughout the orbital optimisation and includes a newly derived basis set correction to the Brueckner reference energy. The generalisation of F12 theory to arbitrary non-Hartree–Fock references and to Fock operators that include scalar relativistic effects is also presented. The performance of the new Brueckner F12 method is assessed for a test set of 50 open- and closed-shell reactions and for the ionisation potentials and electron affinities (EAs) of the first-row transition metal atoms. Benchmark basis set limit coupled-cluster singles, doubles and perturbative triples (CCSD(T)) and BCCD(T) values are reported for all energies in the test sets. BCCD(T){sub (F12*)} performs systematically better than CCSD(T)(F12*) for electron affinities where orbital relaxation effects are significant.

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

Defensive Islanding to Enhance the Resilience of Distribution Systems Against Cyber-Induced Failures

The extensive integration of communication, computation, and control technologies into cyber-physical power systems (CPPSs) has increased the vulnerabilities of CPPSs to cyberattacks. This calls for developing solutions that assess and reduce the impacts of cyber-induced failures on CPPSs. This paper proposes a defensive islanding strategy to isolate impacted parts of the CPPS and form self-sufficient islanded grids with an objective of minimum load curtailment. The defensive islanding aims to split a power system into smaller grids to improve its resilience against a potential extreme event. A clustering approach that leverages the hierarchical spectral clustering method is utilized for the optimal defensive islanding. The proposed approach captures the fragility behavior and loading conditions of power system components due to cyber-induced failures. A graphical-based coupling framework is used to map the impacts of cyber failures into operation of power system components. The proposed method is demonstrated on a modified 33-node distribution feeder system integrated with distributed energy resources. The amount of load curtailment and radiality constraints have been used to evaluate the performance of the proposed clustering strategies. The results show the capability of the proposed algorithm to create islands considering the cyber-induced failures for enhanced resilience.

cyber-induced failures↗
DOE OSTI2018

On the Coupling Time of the Heat-Bath Process for the Fortuin–Kasteleyn Random–Cluster Model

We consider the coupling from the past implementation of the random–cluster heat-bath process, and study its random running time, or coupling time. We focus on hypercubic lattices embedded on tori, in dimensions one to three, with cluster fugacity at least one. We make a number of conjectures regarding the asymptotic behaviour of the coupling time, motivated by rigorous results in one dimension and Monte Carlo simulations in dimensions two and three. Amongst our findings, we observe that, for generic parameter values, the distribution of the appropriately standardized coupling time converges to a Gumbel distribution, and that the standard deviation of the coupling time is asymptotic to an explicit universal constant multiple of the relaxation time. Perhaps surprisingly, we observe these results to hold both off criticality, where the coupling time closely mimics the coupon collector’s problem, and also at the critical point, provided the cluster fugacity is below the value at which the transition becomes discontinuous. Finally, we consider analogous questions for the single-spin Ising heat-bath process.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Classical Preoptimization Approach for ADAPT-VQE: Maximizing the Potential of High-Performance Computing Resources to Improve Quantum Simulation of Chemical Applications

The ADAPT-VQE algorithm is a promising method for generating a compact ansatz based on derivatives of the underlying cost function, and it yields accurate predictions of electronic energies for molecules. In this work, we report the implementation and performance of ADAPT-VQE with our recently developed sparse wave function circuit solver (SWCS) in terms of accuracy and efficiency for molecular systems with up to 52 spin orbitals. The SWCS can be tuned to balance computational cost and accuracy, which extends the application of ADAPT-VQE for molecular electronic structure calculations to larger basis sets and a larger number of qubits. Using this tunable feature of the SWCS, we propose an alternative optimization procedure for ADAPT-VQE to reduce the computational cost of the optimization. Furthermore, by preoptimizing a quantum simulation with a parametrized ansatz generated with ADAPT-VQE/SWCS, we aim to utilize the power of classical high-performance computing in order to minimize the work required on noisy intermediate-scale quantum hardware, which offers a promising path toward demonstrating quantum advantage for chemical applications.

ADAPT-VQE↗
DOE OSTI2013

Temperature dependences for the reactions of O{sub 2}{sup −} and O{sup −} with N and O atoms in a selected-ion flow tube instrument

Rate constants for the reactions of O{sub 2}{sup −} and O{sup −} with N and O atoms have been measured for the first time as a function of temperature from 173 to 500 K for O{sup −} reactions and 173 to 400 K for O{sub 2}{sup −} reactions. Room temperature rate constants for O{sub 2}{sup −} reacting with N and O are 3.1 × 10{sup −10} and 1.7 × 10{sup −10} cm{sup 3} s{sup −1}, respectively, and the corresponding O{sup −} rate constants are 1.7 × 10{sup −10} and 1.5 × 10{sup −10} cm{sup 3} s{sup −1}, in good agreement with previous values. Temperature dependences are about T{sup −1.7} for both O{sub 2}{sup −} reactions and T{sup −0.6} and T{sup −1.3} for the reactions of O{sup −} with N and O, respectively. Branching for the O{sub 2}{sup −} reaction with N is found to predominantly form O{sup −} (>85%) in contrast to previous measurements, which reported NO{sub 2}+e{sup −} as the main channel. Calculations point to the present results being correct. The potential energy surface for this reaction was calculated using density functional theory, coupled cluster with singles, doubles (triples), complete active space self-consistent field, and complete active space second-order perturbation methods and is found to be quite complex, with agreement between the calculated surface and the observed kinetic data only possible through the inclusion of dynamical correlation.

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

A phaseless auxiliary-field quantum Monte Carlo perspective on the uniform electron gas at finite temperatures: Issues, observations, and benchmark study

We investigate the viability of the phaseless finite-temperature auxiliary-field quantum Monte Carlo (ph-FT-AFQMC) method for ab initio systems using the uniform electron gas as a model. Through comparisons with exact results and FT coupled cluster theory, we find that ph-FT-AFQMC is sufficiently accurate at high to intermediate electronic densities. We show, both analytically and numerically, that the phaseless constraint at FT is fundamentally different from its zero-temperature counterpart (i.e., ph-ZT-AFQMC), and generally, one should not expect ph-FT-AFQMC to agree with ph-ZT-AFQMC in the low-temperature limit. With an efficient implementation, we are able to compare exchange-correlation energies to the existing results in the thermodynamic limit and find that the existing parameterizations are highly accurate. In particular, we found that ph-FT-AFQMC exchange-correlation energies are in better agreement with a known parameterization than is restricted path-integral MC in the regime of Θ ≤ 0.5 and rs ≤ 2, which highlights the strength of ph-FT-AFQMC.

Lee, Joonho (ORCID:0000000296671081)↗
DOE OSTI2023

The Generalized Green’s function Cluster Expansion: A Python package for simulating polarons

We present an efficient implementation of the Generalized Green’s function Cluster Expansion (GGCE), which is a new method for computing the ground-state properties and dynamics of polarons (single electrons coupled to lattice vibrations) in model electron-phonon systems. The GGCE works at arbitrary temperature and is well suited for a variety of electron-phonon couplings, including, but not limited to, site and bond Holstein and Peierls (Su-SchriefferHeeger) couplings, and couplings to multiple phonon modes with different energy scales and coupling strengths. Quick calculations can be performed efficiently on a laptop using solvers from NumPy and SciPy, or in parallel at scale using the PETSc sparse linear solver engine.

36 MATERIALS SCIENCE↗
NASA NTRS1990

Triple and quadruple excitation contributions to the binding in Be clusters - Calibration calculations on Be3

The contribution of connected triple and quadruple excitations to the binding in Be3 is investigated by comparing various coupled-cluster (CC) and truncated configuration-interaction (CI) treatments with multireference CI (MRCI) and full CI (FCI) calculations. The CC method with single and double excitations (CCSD) produces results that differ substantially from more elaborate treatments, but most extensions to CCSD that account approximately for connected triple excitations perform very well. In contrast, good agreement with CFI for Be2 can be achieved only with the highest level CC and MRCI methods.

Watts, John D.↗
DOE OSTI2020

An ab initio exciton model for singlet fission

Here we present an ab initio exciton model that extends the Frenkel exciton model and includes valence, charge-transfer, and multiexcitonic excited states. It serves as a general, parameter-free, yet computationally efficient and scalable approach for simulation of singlet fission processes in multichromophoric systems. A comparison with multiconfigurational methods confirms that our exciton model predicts consistent energies and couplings for the pentacene dimer and captures the correct physics. Calculations of larger pentacene clusters demonstrate the computational scalability of the exciton model and suggest that the mixing between local and charge-transfer excitations narrows the gap between singlet and multiexcitonic states. Local vibrations of pentacene molecules are found to facilitate singlet–multiexcitonic state-crossing and hence are important for understanding singlet fission. The exciton model developed in this work also sets the stage for further implementation of the nuclear gradients and nonadiabatic couplings needed for first principles nonadiabatic quantum molecular dynamics simulations of singlet fission.

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

Fourier imaging of non-linear structure formation

We perform a Fourier space decomposition of the dynamics of non-linear cosmological structure formation in ΛCDM models. From N -body simulations involving only cold dark matter we calculate 3-dimensional non-linear density, velocity divergence and vorticity Fourier realizations, and use these to calculate the fully non-linear mode coupling integrals in the corresponding fluid equations. Our approach allows for a reconstruction of the amount of mode coupling between any two wavenumbers as a function of redshift. With our Fourier decomposition method we identify the transfer of power from larger to smaller scales, the stable clustering regime, the scale where vorticity becomes important, and the suppression of the non-linear divergence power spectrum as compared to linear theory. Our results can be used to improve and calibrate semi-analytical structure formation models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Enhancing the accuracy of density functional tight binding models through ChIMES many-body interaction potentials

Semi-empirical quantum models such as Density Functional Tight Binding (DFTB) are attractive methods for obtaining quantum simulation data at longer time and length scales than possible with standard approaches. However, application of these models can require lengthy effort due to the lack of a systematic approach for their development. In this work, we discuss the use of the Chebyshev Interaction Model for Efficient Simulation (ChIMES) to create rapidly parameterized DFTB models, which exhibit strong transferability due to the inclusion of many-body interactions that might otherwise be inaccurate. We apply our modeling approach to silicon polymorphs and review previous work on titanium hydride. We also review the creation of a general purpose DFTB/ChIMES model for organic molecules and compounds that approaches hybrid functional and coupled cluster accuracy with two orders of magnitude fewer parameters than similar neural network approaches. In all cases, DFTB/ChIMES yields similar accuracy to the underlying quantum method with orders of magnitude improvement in computational cost. In conclusion, our developments provide a way to create computationally efficient and highly accurate simulations over varying extreme thermodynamic conditions, where physical and chemical properties can be difficult to interrogate directly, and there is historically a significant reliance on theoretical approaches for interpretation and validation of experimental results.

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

Hetero-site Double Core Ionization Energies with Sub-electronvolt Accuracy from Delta-Coupled-Cluster Calculations

Benchmark scalar-relativistic delta-coupled-cluster calculations of hetero-site double core ionization energies of small molecules containing second-row elements are reported. The present study has focused on the high-spin triplet components of two-site double core-ionized states, which are single reference in character and consistent with the use of standard coupledcluster methods. Contributions to computed double core ionization energies from electron-correlation and basis-set effects as well as corrections to the core-valence separation approximation have been analyzed. Furthermore, on the basis of systematic convergence of computational results with respect to these effects, delta-coupled-cluster calculations have been shown to be capable of providing accurate double core ionization energies with remaining errors estimated to be below 0.3 eV, and thus they are recommended for use to facilitate experimental studies of two-site double core-ionized states that are involved in X-ray pump/X-ray probe studies of electronic and molecular dynamics following inner shell ionization or excitation.

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

Effect of Coulomb Correlation on the Magnetic Properties of Mn Clusters

In spite of decades of research, a fundamental understanding of the unusual magnetic behavior of small Mn clusters remains a challenge. Experiments show that Mn 2 is antiferromagnetic while small clusters containing up to five Mn atoms are ferromagnetic with magnetic moments of 5 μB/atom and become ferrimagnetic as they grow further. Theoretical studies based on density functional theory (DFT), however, find Mn2 to be ferromagnetic, with ferrimagnetic order setting in at different sizes that depend upon the computational methods used. While quantum chemical techniques correctly account for the antiferromagnetic ground state of Mn 2 , they are computationally too demanding to treat larger clusters, making it difficult to understand the evolution of magnetism. These studies clearly point to the importance of correlation and the need to find ways to treat it effectively for larger clusters and nanostructures. In this report we show that the DFT+U method can be used to account for strong correlation. We determine the on-site Coulomb correlation, Hubbard U self-consistently by using the linear response theory and study its effect on the magnetic coupling of Mn clusters containing up to five atoms. With a calculated U value of 4.8 eV, we show that the ground state of Mn 2 is antiferromagnetic with a Mn–Mn distance of 3.34 Å, which agrees well with the electron spin resonance experiment. Equally important, we show that on-site Coulomb correlation also plays an important role in the evolution of magnetic coupling in larger clusters, as the results differ significantly from standard DFT calculations. We conclude that for a proper understanding of magnetism of Mn nanostructures (clusters, chains, and layers) one must take into account the effect of strong correlation.

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

Quantification of electron correlation for approximate quantum calculations

State-of-the-art many-body wave function techniques rely on heuristics to achieve high accuracy at an attainable computational cost to solve the many-body Schrödinger equation. By far, the most common property used to assess accuracy has been the total energy; however, total energies do not give a complete picture of electron correlation. In this work, we assess the von Neumann entropy of the one-particle reduced density matrix (1-RDM) to compare selected configuration interaction (CI), coupled cluster, variational Monte Carlo, and fixed-node diffusion Monte Carlo for benchmark hydrogen chains. A new algorithm, the circle reject method, is presented, which improves the efficiency of evaluating the von Neumann entropy using quantum Monte Carlo by several orders of magnitude. The von Neumann entropy of the 1-RDM and the eigenvalues of the 1-RDM are shown to distinguish between the dynamic correlation introduced by the Jastrow and the static correlation introduced by determinants with large weights, confirming some of the lore in the field concerning the difference between the selected CI and Slater–Jastrow wave functions.

Chemistry↗
DOE OSTI2017

Testing approximate predictions of displacements of cosmological dark matter halos

We present a test to quantify how well some approximate methods, designed to reproduce the mildly non-linear evolution of perturbations, are able to reproduce the clustering of DM halos once the grouping of particles into halos is defined and kept fixed. The following methods have been considered: Lagrangian Perturbation Theory (LPT) up to third order, Truncated LPT, Augmented LPT, MUSCLE and COLA. The test runs as follows: halos are defined by applying a friends-of-friends (FoF) halo finder to the output of an N-body simulation. The approximate methods are then applied to the same initial conditions of the simulation, producing for all particles displacements from their starting position and velocities. The position and velocity of each halo are computed by averaging over the particles that belong to that halo, according to the FoF halo finder. This procedure allows us to perform a well-posed test of how clustering of the matter density and halo density fields are recovered, without asking to the approximate method an accurate reconstruction of halos. We have considered the results at z =0,0.5,1, and we have analysed power spectrum in real and redshift space, object-by-object difference in position and velocity, density Probability Distribution Function (PDF) and its moments, phase difference of Fourier modes. We find that higher LPT orders are generally able to better reproduce the clustering of halos, while little or no improvement is found for the matter density field when going to 2LPT and 3LPT. Augmentation provides some improvement when coupled with 2LPT, while its effect is limited when coupled with 3LPT. Little improvement is brought by MUSCLE with respect to Augmentation. The more expensive particle-mesh code COLA outperforms all LPT methods, and this is true even for mesh sizes as large as the inter-particle distance. This test sets an upper limit on the ability of these methods to reproduce the clustering of halos, for the cases when these objects are reconstructed at the object-by-object level.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Analytic gradients for the single-reference driven similarity renormalization group second-order perturbation theory

Herein we derive and implement analytic energy gradients for the single-reference driven similarity renormalization group second-order perturbation theory (DSRG-PT2). The resulting equations possess an asymptotic scaling that is identical to that of the second-order Møller–Plesset perturbation theory (MP2), indicating that the exponential regularizer in the DSRG equations does not introduce formal difficulties in the gradient theory. We apply the DSRG-PT2 method to optimizing the geometries of 15 small molecules. The equilibrium bond lengths computed with DSRG-PT2 are found similar to those of MP2, yielding a mean absolute error of 0.0033 Å and a standard deviation of 0.0045 Å when compared with coupled cluster with singles, doubles, and perturbative triples.

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

Chemical accuracy from quantum Monte Carlo for the benzene dimer

We report an accurate study of interactions between benzene molecules using variational quantum Monte Carlo (VMC) and diffusion quantum Monte Carlo (DMC) methods. We compare these results with density functional theory using different van der Waals functionals. In our quantum Monte Carlo (QMC) calculations, we use accurate correlated trial wave functions including three-body Jastrow factors and backflow transformations. We consider two benzene molecules in the parallel displaced geometry, and find that by highly optimizing the wave function and introducing more dynamical correlation into the wave function, we compute the weak chemical binding energy between aromatic rings accurately. We find optimal VMC and DMC binding energies of −2.3(4) and −2.7(3) kcal/mol, respectively. The best estimate of the coupled-cluster theory through perturbative triplets/complete basis set limit is −2.65(2) kcal/mol [Miliordos et al., J. Phys. Chem. A 118, 7568 (2014)]. Our results indicate that QMC methods give chemical accuracy for weakly bound van der Waals molecular interactions, comparable to results from the best quantum chemistry methods.

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

Theoretical prediction of nuclear magnetic shieldings and indirect spin-spin coupling constants in 1,1-, cis-, and trans-1,2-difluoroethylenes

A theoretical prediction of nuclear magnetic shieldings and indirect spin-spin coupling constants in 1,1-, cis- and trans-1,2-difluoroethylenes is reported. The results obtained using density functional theory (DFT) combined with large basis sets and gauge-independent atomic orbital calculations were critically compared with experiment and conventional, higher level correlated electronic structure methods. Accurate structural, vibrational, and NMR parameters of difluoroethylenes were obtained using several density functionals combined with dedicated basis sets. B3LYP/6-311++G(3df,2pd) optimized structures of difluoroethylenes closely reproduced experimental geometries and earlier reported benchmark coupled cluster results, while BLYP/6-311++G(3df,2pd) produced accurate harmonic vibrational frequencies. The most accurate vibrations were obtained using B3LYP/6-311++G(3df,2pd) with correction for anharmonicity. Becke half and half (BHandH) density functional predicted more accurate {sup 19}F isotropic shieldings and van Voorhis and Scuseria's τ-dependent gradient-corrected correlation functional yielded better carbon shieldings than B3LYP. A surprisingly good performance of Hartree-Fock (HF) method in predicting nuclear shieldings in these molecules was observed. Inclusion of zero-point vibrational correction markedly improved agreement with experiment for nuclear shieldings calculated by HF, MP2, CCSD, and CCSD(T) methods but worsened the DFT results. The threefold improvement in accuracy when predicting {sup 2}J(FF) in 1,1-difluoroethylene for BHandH density functional compared to B3LYP was observed (the deviations from experiment were −46 vs. −115 Hz)

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗