Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “correlation energy”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Assessing MP2 frozen natural orbitals in relativistic correlated electronic structure calculations

The high computational scaling with the basis set size and the number of correlated electrons is a bottleneck limiting applications of coupled cluster algorithms, in particular for calculations based on two- or four-component relativistic Hamiltonians, which often employ uncontracted basis sets. This problem may be alleviated by replacing canonical Hartree–Fock virtual orbitals by natural orbitals (NOs). Here, in this paper, we describe the implementation of a module for generating NOs for correlated wavefunctions and, in particular, second order Møller–Plesset perturbation frozen natural orbitals (MP2FNOs) as a component of our novel implementation of relativistic coupled cluster theory for massively parallel architectures [Pototschnig et al. J. Chem. Theory Comput. 17, 5509, (2021)]. Our implementation can manipulate complex or quaternion density matrices, thus allowing for the generation of both Kramers-restricted and Kramers-unrestricted MP2FNOs. Furthermore, NOs are re-expressed in the parent atomic orbital (AO) basis, allowing for generating coupled cluster singles and doubles NOs in the AO basis for further analysis. By investigating the truncation errors of MP2FNOs for both the correlation energy and molecular properties—electric field gradients at the nuclei, electric dipole and quadrupole moments for hydrogen halides HX (X = F–Ts), and parity-violating energy differences for H 2 Z 2 (Z = O–Se)—we find MP2FNOs accelerate the convergence of the correlation energy in a roughly uniform manner across the Periodic Table. It is possible to obtain reliable estimates for both energies and the molecular properties considered with virtual molecular orbital spaces truncated to about half the size of the full spaces.

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

Ab initio calculation of atomic solid hydrogen phases based on Gutzwiller many-body wave functions

We apply two ab initio many-body methods based on Gutzwiller wave functions, i.e., correlation matrix renormalization theory (CMRT) and Gutzwiller conjugate gradient minimization (GCGM), to the study of crystalline phases of atomic hydrogen. Both methods avoid empirical Hubbard U parameters and are free from double-counting issues. CMRT employs a Gutzwiller-type approximation that enables efficient calculations, while GCGM goes beyond this approximation to achieve higher accuracy at higher computational cost. By benchmarking against available quantum Monte Carlo (QMC) results, we demonstrate that while both methods are more accurate than the widely used density-functional theory, GCGM systematically captures additional correlation energy missing in CMRT, leading to significantly improved total energy predictions. We also show that by including the correlation energy Ec from local density approximation in the CMRT calculation, CMRT + E c produces energy in better agreement with the QMC results in these hydrogen lattice systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1977

Diagrammatic perturbation theory - The ground state of the carbon monosulfide molecule

Diagrammatic many-body perturbation theory is employed in a study of the ground state of the carbon monosulfide molecule for bond lengths close to the equilibrium value. The calculations are complete through third order in the energy within the algebraic approximation. Two different zero-order Hamiltonians are considered, and all two-, three-, and four-body terms are determined for the corresponding perturbation expansions. Many-body effects are found to be very important. Pade approximants to the energy expansion are constructed, and upper bounds evaluated. Almost 53 percent of the estimated correlation energy is recovered. The variation of components of the correlation energy with nuclear separation is investigated. Spectroscopic constants are also calculated.

Wilson, S.↗
DOE OSTI2022

Analytic gradients for compressed multistate pair-density functional theory

Photochemical reactions often involve states that are closely coupled due to near degeneracies, for example by proximity to conical intersections. Therefore, a multistate method is used to accurately describe these states; for example, ordinary perturbation theory is replaced by quasidegenerate perturbation theory. Multiconfiguration pair-density functional theory (MC-PDFT) provides an efficient way to approximate the full dynamical correlation energy of strongly correlated systems, and we recently proposed compressed multistate pair-density functional theory (CMS-PDFT) to treat closely coupled states. In the present paper, we report the implementation of analytic gradients for CMS-PDFT in both OpenMolcas and PySCF, and here we illustrate the use of these gradients by applying the method to the excited states of formaldehyde and phenol.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Correlations between energy and $\gamma$-ray emission in 239 Pu($n, f$)

Here, we study $\gamma$-ray emission following $^{239}\mathrm{Pu}(n,\mathrm{f})$ over an incident neutron energy range of $2 < E_i < 40$ MeV. We present the first experimental evidence for positive correlations between the total angular momentum generated in fission and the excitation energy of the compound nucleus prior to fission. The $\gamma$-ray multiplicity increases linearly with incident energy below the 2\textsuperscript{nd}-chance fission threshold with a slope of $0.085 \pm 0.010$~MeV$^{-1}$. This linear trend appears to hold for the average excitation energy of the compound nucleus between $9 < \langle E_x \rangle < 19$ MeV. Most of the multiplicity increase comes from an enhancement around a $\gamma$-ray energy of 0.7 MeV, which we interpret as stretched quadrupole $\gamma$ rays that indicate an increase in total fission-fragment angular momentum with excitation energy.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Long-Range Azimuthal Correlation, Entanglement, and Bell Inequality Violation by Spinning Gluons at the Large Hadron Collider

We apply the recently developed concept of the nucleon energy–energy correlator (NEEC) for the gluon sector to investigate the long-range azimuthal angular correlations in proton–proton collisions at the Large Hadron Collider. The spinning gluon in these collisions will introduce substantial nonzero cos(2Φ) asymmetries in both Higgs boson and top quark pair productions, where Φ is the azimuthal angle between the forward and backward energy correlators in the NEEC observables. The genesis of the cos(2Φ) correlation lies in the intricate quantum entanglement. Owing to the substantial cos(2Φ) effect, the NEEC observable in Higgs boson and $t\bar{t}$ production emerges as a pivotal avenue for delving into quantum entanglement and scrutinizing the Bell inequality at high-energy colliders.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1996

The numerical simulation of multistage turbomachinery flows

The effect of the unsteady flow field in a multistage compressor on the time-averaged performance was assessed. The energy transport by the unsteady deterministic flow field was taken into account. The magnitude of the body force resulting from the aerodynamic loading on a blade row was compared to the gradient of the stress tensor associated with the unsteady time-resolved flow field generated by the blade row. The magnitude of the work performed by these forces was compared to the divergence of the energy correlations produced by the unsteady time-resolved flow field. The stress tensor and the energy correlations are non-negligible in the end wall regions. The results suggest that the turbulence is the primary source of flow mixing in the midspan region.

Adamczyk, John J.↗
DOE OSTI2024

Symmetry breaking and self-interaction correction in the chromium atom and dimer

Density functional approximations to the exchange–correlation energy can often identify strongly correlated systems and estimate their energetics through energy-minimizing symmetry-breaking. In particular, the binding energy curve of the strongly correlated chromium dimer is described qualitatively by the local spin density approximation (LSDA) and almost quantitatively by the Perdew–Burke–Ernzerhof generalized gradient approximation (PBE-GGA), where the symmetry breaking is antiferromagnetic for both. Here, we show that a full Perdew–Zunger self-interaction-correction (SIC) to LSDA seems to go too far by creating an unphysical symmetry-broken state, with effectively zero magnetic moment but non-zero spin density on each atom, which lies ~4 eV below the antiferromagnetic solution. A similar symmetry-breaking, observed in the atom, better corresponds to the 3d↑↑4s↑3d↓↓4s↓ configuration than to the standard 3d↑↑↑↑↑4s↑. For this new solution, the total energy of the dimer at its observed bond length is higher than that of the separated atoms. In conclusion, these results can be regarded as qualitative evidence that the SIC needs to be scaled down in many-electron regions.

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

Mixed Effect Modeling of Dose and Linear Energy Transfer Correlations With Brain Image Changes After Intensity Modulated Proton Therapy for Skull Base Head and Neck Cancer

Intensity modulated proton therapy (IMPT) could yield high linear energy transfer (LET) in critical structures and increased biological effect. For head and neck cancers at the skull base this could potentially result in radiation-associated brain image change (RAIC). The purpose of the current study was to investigate voxel-wise dose and LET correlations with RAIC after IMPT.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2022

Meta-GGA exchange-correlation free energy density functional to increase the accuracy of warm dense matter simulations

We discuss strategies for thermalization of the ground-state meta-generalized gradient approximation (meta-GGA) exchange-correlation (XC) functionals. A simple but accurate scheme is implemented via universal additive thermal correction to XC using a perturbative-like self-consistent approach. The additive correction with explicit temperature dependence is applied to the ground-state deorbitalized, strongly constrained and appropriately normed (SCAN-L) meta-GGA XC leading to thermal XC functional denoted here as T-SCAN-L. Thermal T-SCAN-L meta-GGA functional shows significant improvement in density functional theory calculation accuracy for warm dense matter by a factor of 3 to 10, achieving unprecedented accuracy of total pressure between a few tenths and 1% when compared to traditional XC functionals, as demonstrated by the comparison to pathintegral Monte Carlo simulations for helium equation of state. Furthermore, the T-SCAN-L calculations of dc conductivity of warm dense aluminum also give better agreement with experiments over other XC functionals such as PBE and SCAN-L.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Laplacian-level meta-generalized gradient approximation for solid and liquid metals

In this work, we derive and motivate a Laplacian-level, orbital-free meta-generalized-gradient approximation (LL-MGGA) for the exchange-correlation energy, targeting accurate ground-state properties of sp and sd metallic condensed matter, in which the density functional for the exchange-correlation energy is only weakly nonlocal due to perfect long-range screening. Our model for the orbital-free kinetic energy density restores the fourth-order gradient expansion for exchange to the r 2 SCAN meta-GGA [Furness et al., J. Phys. Chem. Lett. 11, 8208 (2020)], yielding a LL-MGGA we call OFR2. OFR2 matches the accuracy of SCAN for prediction of common lattice constants and improves the equilibrium properties of alkali metals, transition metals, and intermetallics that were degraded relative to the PBE GGA values by both SCAN and r 2 SCAN. We compare OFR2 to the r 2 SCAN -L LL-MGGA [Mejia-Rodriguez and Trickey, Phys. Rev. B 102, 121109(R) (2020)] and show that OFR2 tends to outperform r 2 SCAN -L for the equilibrium properties of solids, but r 2 SCAN -L much better describes the atomization energies of molecules than OFR2 does. For best accuracy in molecules and nonmetallic condensed matter, we continue to recommend SCAN and r 2 SCAN. Numerical performance is discussed in detail, and our paper provides an outlook to machine learning.

36 MATERIALS SCIENCE↗
DOE OSTI2022

CAFQA: A Classical Simulation Bootstrap for Variational Quantum Algorithms

Classical computing plays a critical role in the advancement of quantum frontiers in the NISQ era. In this spirit, this work uses classical simulation to bootstrap Variational Quantum Algorithms (VQAs). VQAs rely upon the iterative optimization of a parameterized unitary circuit (ansatz) with respect to an objective function. Since quantum machines are noisy and expensive resources, it is imperative to classically choose the VQA ansatz initial parameters to be as close to optimal as possible to improve VQA accuracy and accelerate their convergence on today’s devices. This work tackles the problem of finding a good ansatz initialization, by proposing CAFQA, a Clifford Ansatz For Quantum Accuracy. The CAFQA ansatz is a hardware-efficient circuit built with only Clifford gates. In this ansatz, the parameters for the tunable gates are chosen by searching efficiently through the Clifford parameter space via classical simulation. The resulting initial states always equal or outperform traditional classical initialization (e.g., Hartree-Fock), and enable high-accuracy VQA estimations. CAFQA is well-suited to classical computation because: a) Clifford-only quantum circuits can be exactly simulated classically in polynomial time, and b) the discrete Clifford space is searched efficiently via Bayesian Optimization. For the Variational Quantum Eigensolver (VQE) task of molecular ground state energy estimation (up to 18 qubits), CAFQA’s Clifford Ansatz achieves a mean accuracy of nearly 99% and recovers as much as 99.99% of the molecular correlation energy that is lost in Hartree-Fock initialization. CAFQA achieves mean accuracy improvements of 6.4x and 56.8x, over the state-of-the-art, on different metrics. Here, the scalability of the approach allows for preliminary ground state energy estimation of the challenging chromium dimer (Cr2) molecule. With CAFQA’s high-accuracy initialization, the convergence of VQAs is shown to accelerate by 2.5x, even for small molecules. Furthermore, preliminary exploration of allowing a limited number of non-Clifford (T) gates in the CAFQA framework, shows that as much as 99.9% of the correlation energy can be recovered at bond lengths for which Clifford-only CAFQA accuracy is relatively limited, while remaining classically simulable.

bayesian optimization↗
DOE OSTI2025

Pd–Methyl Bond Energy─Property Correlations, Noncorrelations, Machine Learning Models, and Application to Polymerization Catalysis

Metal–carbon bonds are a key intermediate in a variety of homogeneous organometallic transformations and often determine the critical thermodynamics and kinetics of catalytic processes. Surprisingly, the influence of different ligands on metal–carbon bond strengths has been largely overlooked. Here, in this study, we evaluated nearly 700 experimental Pd–methyl complexes by calculating their bond dissociation energies using density functional theory (DFT) and compared these bond strengths to several fundamental molecular properties, and this revealed several surprising correlations and noncorrelations. Most surprising was that several fundamental properties, such as the bond length, bond force constant, and bond electron density, have no correlation with bond strength, despite these correlations often holding for main-group compounds. We were indeed able to identify key ligand-dependent chemical features/descriptors that provided a highly accurate machine learning model and provided insight into the general factors that control the Pd–carbon bond strength, such as radical delocalization and nucleophilicity. Insights gained from the Pd–Me bond energy analysis were then applied to CO migratory insertion steps that are part of copolymerization reactions.

binding energy↗
DOE OSTI2022

First-principles wave-vector- and frequency-dependent exchange-correlation kernel for jellium at all densities

Here we propose a spatially and temporally nonlocal exchange correlation (XC) kernel for the spin-unpolarized fluid phase of ground-state jellium for use in time-dependent density functional and linear response calculations. The kernel is constructed to satisfy known properties of the exact XC kernel to accurately describe the correlation energies of bulk jellium and to satisfy frequency-moment sum rules at a wide range of bulk jellium densities, including those low densities that display strong correlation and symmetry breaking. These effects are easier to understand in the simple jellium model than in real systems. All exact constraints satisfied by the recent MCP07 kernel are maintained in the revised MCP07 (rMCP07) kernel, while others are added. The revision $f^{rMCP07}_{XC}$ (q, ω) differs from MCP07 only for nonzero frequencies ω. Only at densities much lower than those of real bulk metals is the frequency dependence of the kernel important for the correlation energy of jellium. As the wave vector q tends to zero, the kernel has a -4$πα(ω)/q^2$ divergence whose frequency-dependent ultranonlocality coefficient $α(ω)$ vanishes in jellium, and is predicted by rMCP07 to be extremely small for the real metals Al and Na.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Porting fragmentation methods to GPUs using an OpenMP API: Offloading the resolution-of-the-identity second-order Møller–Plesset perturbation method

Here, using an OpenMP Application Programming Interface, the resolution-of-the-identity second-order Møller–Plesset perturbation (RI-MP2) method has been off-loaded onto graphical processing units (GPUs), both as a standalone method in the GAMESS electronic structure program and as an electron correlation energy component in the effective fragment molecular orbital (EFMO) framework. First, a new scheme has been proposed to maximize data digestion on GPUs that subsequently linearizes data transfer from central processing units (CPUs) to GPUs. Second, the GAMESS Fortran code has been interfaced with GPU numerical libraries (e.g., NVIDIA cuBLAS and cuSOLVER) for efficient matrix operations (e.g., matrix multiplication, matrix decomposition, and matrix inversion). The standalone GPU RI-MP2 code shows an increasing speedup of up to 7.5× using one NVIDIA V100 GPU with one IBM 42-core P9 CPU for calculations on fullerenes of increasing size from 40 to 260 carbon atoms using the 6-31G(d)/cc-pVDZ-RI basis sets. A single Summit node with six V100s can compute the RI-MP2 correlation energy of a cluster of 175 water molecules using the correlation consistent basis sets cc-pVDZ/cc-pVDZ-RI containing 4375 atomic orbitals and 14 700 auxiliary basis functions in ~0.85 h. In the EFMO framework, the GPU RI-MP2 component shows near linear scaling for a large number of V100s when computing the energy of an 1800-atom mesoporous silica nanoparticle in a bath of 4000 water molecules. The parallel efficiencies of the GPU RI-MP2 component with 2304 and 4608 V100s are 98.0% and 96.1%, respectively.

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 OSTI2024

The rise and fall of stretched bond errors: Extending the analysis of Perdew–Zunger self-interaction corrections of reaction barrier heights beyond the LSDA

Incorporating self-interaction corrections (SIC) significantly improves chemical reaction barrier height predictions made using density functional theory methods. We present a detailed orbital-by-orbital analysis of these corrections for three semi-local density functional approximations (DFAs) situated on the three lowest rungs of Jacob’s ladder of approximations. The analysis is based on Fermi–Löwdin Orbital Self-Interaction Correction (FLOSIC) calculations performed at several steps along the reaction pathway from the reactants (R) to the transition state (TS) to the products (P) for four representative reactions selected from the BH76 benchmark set. For all three functionals, the major contribution to self-interaction corrections of the barrier heights can be traced to stretched bond orbitals that develop near the TS configuration. The magnitude of the ratio of the self-exchange–correlation energy to the self-Hartree energy (XC/H) for a given orbital is introduced as an indicator of one-electron self-interaction error. XC/H = 1.0 implies that an orbital’s self-exchange–correlation energy exactly cancels its self-Hartree energy and that the orbital, therefore, makes no contribution to the SIC in the FLOSIC scheme. For the practical DFAs studied here, XC/H spans a range of values. Here, the largest values are obtained for stretched or strongly lobed orbitals. We show that significant differences in XC/H for corresponding orbitals in the R, TS, and P configurations can be used to identify the major contributors to the SIC of barrier heights and reaction energies. Based on such comparisons, we suggest that barrier height predictions made using the strongly constrained and appropriately normed meta-generalized gradient approximation may have attained the best accuracy possible for a semi-local functional using the Perdew–Zunger SIC approach.

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

Geometrical picture of the electron–electron correlation at the large- D limit

In electronic structure calculations, the correlation energy is defined as the difference between the mean field and the exact solution of the non relativistic Schrödinger equation. Such an error in the different calculations is not directly observable as there is no simple quantum mechanical operator, apart from correlation functions, that correspond to such quantity. Here, we use the dimensional scaling approach, in which the electrons are localized at the large-dimensional scaled space, to describe a geometric picture of the electronic correlation. Both, the mean field, and the exact solutions at the large-D limit have distinct geometries. Thus, the difference might be used to describe the correlation effect. Moreover, correlations can be also described and quantified by the entanglement between the electrons, which is a strong correlation without a classical analog. Entanglement is directly observable and it is one of the most striking properties of quantum mechanics and bounded by the area law for local gapped Hamiltonians of interacting many-body systems. This study opens up the possibility of presenting a geometrical picture of the electron–electron correlations and might give a bound on the correlation energy. Furthermore, the results at the large-D limit and at D = 3 indicate the feasibility of using the geometrical picture to get a bound on the electron–electron correlations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗