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

Quantum Chemistry Common Driver and Databases (QCDB) and Quantum Chemistry Engine (QCEngine): Automation and interoperability among computational chemistry programs

We report that community efforts in the computational molecular sciences (CMS) are evolving toward modular, open, and interoperable interfaces that work with existing community codes to provide more functionality and composability than could be achieved with a single program. The Quantum Chemistry Common Driver and Databases (QCDB) project provides such capability through an application programming interface (API) that facilitates interoperability across multiple quantum chemistry software packages. In tandem with the Molecular Sciences Software Institute and their Quantum Chemistry Archive ecosystem, the unique functionalities of several CMS programs are integrated, including CFOUR, GAMESS, NWChem, OpenMM, Psi4, Qcore, TeraChem, and Turbomole, to provide common computational functions, i.e., energy, gradient, and Hessian computations as well as molecular properties such as atomic charges and vibrational frequency analysis. Both standard users and power users benefit from adopting these APIs as they lower the language barrier of input styles and enable a standard layout of variables and data. These designs allow end-to-end interoperable programming of complex computations and provide best practices options by default.

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

Spin-state gaps and self-interaction-corrected density functional approximations: Octahedral Fe(II) complexes as case study

Accurate prediction of a spin-state energy difference is crucial for understanding the spin crossover phenomena and is very challenging for density functional approximations, especially for local and semi-local approximations due to delocalization errors. Here, we investigate the effect of the self-interaction error removal from the local spin density approximation (LSDA) and Perdew–Burke–Ernzerhof generalized gradient approximation on the spin-state gaps of Fe(II) complexes with various ligands using recently developed locally scaled self-interaction correction (LSIC) by Zope et al. [J. Chem. Phys. 151, 214108 (2019)]. The LSIC method is exact for one-electron density, recovers the uniform electron gas limit of the underlying functional, and approaches the well-known Perdew–Zunger self-interaction correction (PZSIC) as a particular case when the scaling factor is set to unity. Our results, when compared with reference diffusion Monte Carlo results, show that the PZSIC method significantly overestimates spin-state gaps favoring low spin states for all ligands and does not improve upon density functional approximations. The perturbative LSIC-LSDA using PZSIC densities significantly improves the gaps with a mean absolute error of 0.51 eV but slightly overcorrects for the stronger CO ligands. Finally, the quasi-self-consistent LSIC-LSDA, such as coupled-cluster single double and perturbative triple [CCSD(T)], gives a correct sign of spin-state gaps for all ligands with a mean absolute error of 0.56 eV, comparable to that of CCSD(T) (0.49 eV).

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

First Joint Oscillation Analysis of Super-Kamiokande Atmospheric and T2K Accelerator Neutrino Data

The Super-Kamiokande and T2K Collaborations present a joint measurement of neutrino oscillation parameters from their atmospheric and beam neutrino data. It uses a common interaction model for events overlapping in neutrino energy and correlated detector systematic uncertainties between the two datasets, which are found to be compatible. Using 3244.4 days of atmospheric data and a beam exposure of 19.7⁢(16.3) × 10 20 protons on target in (anti)neutrino mode, the analysis finds a 1.9⁢𝜎 exclusion of 𝐶⁢𝑃 conservation (defined as 𝐽 𝐶⁢𝑃 = 0) and a 1.2⁢𝜎 exclusion of the inverted mass ordering.

Neutrino detection↗
DOE OSTI2022

Does plasma jet sintering follow an Arrhenius-type expression?

Atmospheric pressure, ambient temperature plasma jets have become a promising candidate for material processing in parallel with developments in additive manufacturing. Recent work has shown that plasma jets can be used to sinter printed nanoparticles at temperatures much lower than typically required for conventional thermal sintering. In this report we conduct a mechanistic study on plasma jet sintering that correlates specific energy input with the electrical conductivity of printed silver films after sintering. Increasing the specific energy input accelerated the sintering process following an Arrhenius-like exponential trend across a large range of conditions, including both helium and argon plasma jets. Although an exponential relationship is also found with the plasma heated substrate temperature, independent studies indicate that heating is not the primary mechanism. These results suggest there is a general behavior that couples the plasma jet with the surface.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Enabling event-by-event precision in γ-ray cascades for neutron-induced reactions

Neutron-induced γ-ray spectra provide key inputs for modern active interrogation applications. A precise modeling of the nuclear reaction and subsequent emission of γ rays is challenging and often impossible due to limitations on evaluated data file formats and nuclear transport simulation codes. We present a framework that addresses these challenges by combining experimental data and reaction-model calculation outputs into an extended candidate version of the Generalized Nuclear Data Structure (GNDS) file, the successor format for the legacy Evaluated Nuclear Data File (ENDF-6). This proposed GNDS hierarchical format contains all the necessary ingredients for inline γ-ray cascade reproduction with event-by-event precision, including continuum–continuum and continuum–discrete transitions following neutron-capture and inelastic neutron scattering reactions. Cascade-event generation based on our approach demonstrates improved energy conservation on an event-by-event basis and permits the use of γ-γ coincidences in applications. This work offers, for the first time, a method to generate neutron-capture and inelastic neutron-scattering γ-ray cascades where energy conservation, correlations, and experimental primaries are fully accounted for.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

GPU Acceleration of Large-Scale Full-Frequency GW Calculations

Many-body perturbation theory is a powerful method to simulate electronic excitations in molecules and materials starting from the output of density functional theory calculations. By implementing the theory efficiently so as to run at scale on the latest leadership high-performance computing systems it is possible to extend the scope of GW calculations. Here, we present a GPU acceleration study of the full-frequency GW method as implemented in the WEST code. Excellent performance is achieved through the use of (i) optimized GPU libraries, e.g., cuFFT and cuBLAS, (ii) a hierarchical parallelization strategy that minimizes CPU-CPU, CPU-GPU, and GPU-GPU data transfer operations, (iii) nonblocking MPI communications that overlap with GPU computations, and (iv) mixed precision in selected portions of the code. A series of performance benchmarks has been carried out on leadership high-performance computing systems, showing a substantial speedup of the GPU-accelerated version of WEST with respect to its CPU version. Good strong and weak scaling is demonstrated using up to 25 920 GPUs. Finally, we showcase the capability of the GPU version of WEST for large-scale, full-frequency GW calculations of realistic systems, e.g., a nanostructure, an interface, and a defect, comprising up to 10 368 valence electrons.

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

Data-driven design of molecular nanomagnets

Three decades of research in molecular nanomagnets have raised their magnetic memories from liquid helium to liquid nitrogen temperature thanks to a wise choice of the magnetic ion and coordination environment. Still, serendipity and chemical intuition played a main role. In order to establish a powerful framework for statistically driven chemical design, here we collected chemical and physical data for lanthanide-based nanomagnets, catalogued over 1400 published experiments, developed an interactive dashboard (SIMDAVIS) to visualise the dataset, and applied inferential statistical analysis. Our analysis shows that the Arrhenius energy barrier correlates unexpectedly well with the magnetic memory. Furthermore, as both Orbach and Raman processes can be affected by vibronic coupling, chemical design of the coordination scheme may be used to reduce the relaxation rates. Indeed, only bis-phthalocyaninato sandwiches and metallocenes, with rigid ligands, consistently present magnetic memory up to high temperature. Analysing magnetostructural correlations, we offer promising strategies for improvement, in particular for the preparation of pentagonal bipyramids, where even softer complexes are protected against molecular vibrations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Tricyclononenes and tricyclononadienes as efficient monomers for controlled ROMP: understanding structure–propagation rate relationships and enabling facile post-polymerization modification

Grubbs 3rd-generation (G3) pre-catalyst-initiated ring-opening metathesis polymerization (ROMP) remains an indispensable tool in the polymer chemist's toolbox. Tricyclononenes (TCN) and tricyclononadienes (TCND) represent under-explored classes of monomers for ROMP that have the potential to both advance fundamental knowledge (e.g., structure-polymerization kinetics relationships) and serve as practical tools for the polymer chemist (e.g., post-polymerization functionalization). In this work, a library of TCN and TCND imides, monoesters, and diesters, along with their exo-norbornene counterparts, were synthesized to compare their behaviors in G3-initiated ROMP. Real-time 1 H NMR was used to study their polymerization kinetics; propagation rates ( k p ) were extracted for each monomer. To understand the relationships between monomer structure and ROMP propagation rates, density functional theory methods were used to calculate a variety of electronic and steric parameters for each monomer. While electronic parameters (e.g., HOMO energy levels) correlated positively with the measured k p values, steric parameters generally gave improved correlations, which indicates that monomer size and shape are better predictors for k p than electronic parameters for this data set. Furthermore, the TCND diester—which contains an electron-deficient cyclobutene that is resistant to ROMP—and its polymer p(TCND) are shown to be highly reactive toward DBU-catalyzed conjugate addition reactions with thiols, providing a protecting- and activating-group free strategy for post-polymerization modification.

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

Importance of cross reaction covariance data for user applications

The characterization of the uncertainty in radiation damage metrics presents many challenges. This paper examines the current approaches to characterizing radiation damage metrics such as hydrogen and helium gas production, material heating, trapped charge in microelectronics, and lattice displacement damage. Critical uncertainty aspects go beyond just the material cross sections and involve the consideration of energy-dependent cross reaction correlations, the recoil ion energy spectrum, and models used for the partitioning of the recoil ion energy into various forms of energy deposition. This paper starts with a review of terminology and then examines the current approaches in the characterization of uncertainty in radiation damage metrics for several applications. The major deficiencies in the uncertainty of the damage metric characterization are also identified.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Re-examining the silicon self-interstitial charge states and defect levels: A density functional theory and bounds analysis study

The self-interstitial atom (SIA) is one of two fundamental point defects in bulk Si. Isolated Si SIAs are extremely difficult to observe experimentally. Even at very low temperatures, they anneal before typical experiments can be performed. Given the challenges associated with experimental characterization, accurate theoretical calculations provide valuable information necessary to elucidate the properties of these defects. Previous studies have applied Kohn–Sham density functional theory (DFT) to the Si SIA, using either the local density approximation or the generalized gradient approximation to the exchange-correlation (XC) energy. The consensus of these studies indicates that a Si SIA may exist in five charge states ranging from -2 to +2 with the defect structure being dependent on the charge state. This study aims to re-examine the existence of these charge states in light of recently derived “approximate bounds” on the defect levels obtained from finite-size supercell calculations and new DFT calculations using both semi-local and hybrid XC approximations. We conclude that only the neutral and +2 charge states are directly supported by DFT as localized charge states of the Si SIA. Within the current accuracy of DFT, our results indicate that the +1 charge state likely consists of an electron in a conduction-band-like state that is coulombically bound to a +2 SIA. Furthermore, the -1 and -2 charge states likely consist of a neutral SIA with one and two additional electrons in the conduction band, respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Convergence acceleration of Monte Carlo many-body perturbation methods by using many control variates

Herein, the use of many control variates is proposed as a method to accelerate the second- and third-order Monte Carlo (MC) many-body perturbation (MC-MP2 and MC-MP3) calculations. A control variate is an exactly integrable function that is strongly correlated or anti-correlated with the target function to be integrated by the MC method. Evaluating both integrals as well as their covariances in the same MC run, one can effect a mutual cancellation of the statistical uncertainties and biases in the MC integrations, thereby accelerating its convergence considerably. Six and thirty-six control variates, whose integrals are known a priori, are generated for MC-MP2 and MC-MP3, respectively, by systematically replacing one or more two-electron-integral vertexes of certain configurations by zero-valued overlap-integral vertexes in their Goldstone diagrams. The variances and co- variances of these control variates are computed at a marginal cost, enhancing the overall efficiency of the MC-MP2 and MC-MP3 calculations by a factor of up to 14 and 20, respectively.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

A coupled cluster framework for electrons and phonons

Here, we describe a coupled cluster framework for coupled systems of electrons and harmonic phonons. Neutral and charged excitations are accessed via the equation-of-motion version of the theory. Benchmarks on the Hubbard–Holstein model allow us to assess the strengths and weaknesses of different coupled cluster approximations, which generally perform well for weak to moderate coupling. Finally, we report progress toward an implementation for ab initio calculations on solids and present some preliminary results on finite-size models of diamond with a linear electron–phonon coupling. We also report the implementation of electron–phonon coupling matrix elements from crystalline Gaussian type orbitals within the PySCF program package.

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

Spin-free formulation of the multireference driven similarity renormalization group: A benchmark study of first-row diatomic molecules and spin-crossover energetics

We report a spin-free formulation of the multireference (MR) driven similarity renormalization group (DSRG) based on the ensemble normal ordering of Mukherjee and Kutzelnigg (1997). This ensemble averages over all microstates of a given total spin quantum number, and therefore, it is invariant with respect to SU(2) transformations. As such, all equations may be reformulated in terms of spin-free quantities and they closely resemble those of spin-adapted closed-shell coupled cluster (CC) theory. The current implementation is used to assess the accuracy of various truncated MR-DSRG methods (perturbation theory up to third order and iterative methods with single and double excitations) in computing the constants of 33 first-row diatomic molecules. The accuracy trends for these first-row diatomics are consistent with our previous benchmark on a small subset of closed-shell diatomic molecules. We then present the first MR-DSRG application on transition-metal complexes by computing the spin splittings of the [Fe(H 2 O) 6 ] 2+ and [Fe(NH 3 ) 6 ] 2+ molecules. A focal point analysis (FPA) shows that third-order perturbative corrections are essential to achieve reasonably converged energetics. The FPA based on the linearized MR-DSRG theory with one- and two-body operators and up to a quintuple- ζ basis set predicts the spin splittings of [Fe(H 2 O) 6 ] 2+ and [Fe(NH 3 ) 6 ] 2+ to be –35.7 and –17.1 kcal mol –1 , respectively, showing good agreement with the results of local CC theory with singles, doubles, and perturbative triples.

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

Piecewise interaction picture density matrix quantum Monte Carlo

The density matrix quantum Monte Carlo (DMQMC) set of methods stochastically samples the exact N-body density matrix for interacting electrons at finite temperature. We introduce a simple modification to the interaction picture DMQMC (IP-DMQMC) method that overcomes the limitation of only sampling one inverse temperature point at a time, instead allowing for the sampling of a temperature range within a single calculation, thereby reducing the computational cost. At the target inverse temperature, instead of ending the simulation, we incorporate a change of picture away from the interaction picture. The resulting equations of motion have piecewise functions and use the interaction picture in the first phase of a simulation, followed by the application of the Bloch equation once the target inverse temperature is reached. We find that the performance of this method is similar to or better than the DMQMC and IP-DMQMC algorithms in a variety of molecular test systems.

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

Exchange correlation potentials from full configuration interaction in a Slater orbital basis

Ryabinkin–Kohut–Staroverov (RKS) theory builds a bridge between wave function theory and density functional theory by using quantities from the former to produce accurate exchange-correlation potentials needed by the latter. In this work, the RKS method is developed and tested alongside Slater atomic orbital basis functions for the first time. To evaluate this approach, full configuration interaction computations in the Slater orbital basis are employed to give quality input to RKS, allowing full correlation to be present along with correct nuclei cusps and asymptotic decay of the wavefunction. SlaterRKS is shown to be an efficient algorithm to arrive at exchange-correlation potentials without unphysical artifacts in moderately-sized basis sets. Furthermore, enforcement of the nuclear cusp conditions will be shown to be vital for the success of the Slater-basis RKS method. As a result, examples of weakly and strongly correlated molecular systems will demonstrate the main features of SlaterRKS.

74 ATOMIC AND MOLECULAR PHYSICS↗