Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Electron correlation calculations”

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 145 records · Page 8

Route to Chemical Accuracy for Computational Uranium Thermochemistry

Benchmark spinor-based relativistic coupled-cluster calculations for the ionization energies of the uranium atom, the uranium monoxide molecule (UO), and the uranium dioxide molecule (UO 2 ) and for the bond dissociation energies of UO and UO 2 are reported. The accuracy of these calculations in the treatments of relativistic, electron correlation, and basis-set effects is analyzed. Furthermore, the intrinsic convergence of the computed results and the favorable comparison with the experimental values demonstrate the unique applicability of the spinor representation of quantum-chemical methods to open-shell uranium-containing atomic and molecular species with uranium oxidation states ranging from U(0) to U(V).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Definitive Assessment of the Accuracy, Variationality, and Convergence of Relativistic Coupled Cluster and Density Matrix Renormalization Group in 100-Orbital Space

Accuracy, variationality, and convergence underpin the reliability of modern electronic structure methods, yet definitive benchmarks in the relativistic regime remain elusive due to the absence of numerically exact full configuration interaction (CI) references. Recent algorithmic advances in the CI framework, enabled by the small-tensor-product (STP) decomposition approach, have dramatically extended the tractable size of the configuration space, making numerically exact CI calculations feasible in large active spaces previously beyond reach. In this paper, we employ the recently developed STP-CI framework to perform large-scale numerically exact CI calculations and directly benchmark relativistic coupled cluster and density matrix renormalization group methods. Definitive benchmarking of approximate relativistic electronic structure methods is ensured through the application of the gap theorem, which provides rigorous error bounds on the CI reference and establishes a controlled standard for assessing accuracy, variationality, and convergence.

Chemical calculations↗

Precision calculation of the lowest 1S resonance in e-H scattering

The position and width of the lowest resonance in electron-hydrogen scattering have been calculated using a Hylleraas correlation function with up to 95 terms in the optical potential formalism. The results should be useful as calibration points for experimental electron scattering purposes. A formula relating the conventional (Breit-Wigner) width with the Feschbach formalism is derived.

Ho, Y. K.↗

Electronic Properties of UN and UN – from Photoelectron Spectroscopy and Correlated Molecular Orbital Theory

In this work, the results of calculations of the properties of the anion UN – including electron detachment are described, which further expand our knowledge of this diatomic molecule. High-level electronic structure calculations were conducted for the UN and UN – diatomic molecules and compared to photoelectron spectroscopy measurements. The low-lying Ω states were obtained using multireference CASPT2 including spin-orbit effects up to ~20,000 cm –1 . At the Feller–Peterson–Dixon (FPD) level, the adiabatic electron affinity (AEA) of UN is estimated to be 1.402 eV and the vertical detachment energy (VDE) is 1.423 eV. The assignment of the UN excited states shows good agreement with the experimental results with a VDE of 1.424 eV. An Ω = 4 ground state was obtained for UN – which is mainly associated with the 3 H ΛS state. Thermochemical calculations estimate a bond dissociation energy (BDE) for UN – (U – + N) of 665.9 kJ/mol, ~15% larger than that of UN and UN + . The NBO analysis reveals U–N triple bonds for the UN, UN – , and UN + species.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

ComDMFT v.2.0: Fully self-consistent ab initio GW+EDMFT for the electronic structure of correlated quantum materials

ComDMFT is a parallel computational package designed to study the electronic structure of correlated quantum materials from first principles. Our approach is based on the combination of first-principles methods and dynamical mean field theories. In version 2.0, we implemented fully-diagrammatic GW+EDMFT from first-principles self-consistently. In this approach, correlated electrons are treated within full GW+EDMFT and the rest are treated within full-GW, seamlessly. Further, this implementation enables the electronic structure calculation of quantum materials with weak, intermediate, and strong electron correlation without prior knowledge of the degree of electron correlation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Methane Adsorption on Heteroatom-Modified Maquettes of Porous Carbon Surfaces

Experimental and theoretical studies disagree on the energetics of methane adsorption on carbon materials. However, this information is critical for the rational design and optimization of the structure and composition of adsorbents for natural gas storage. The delicate nature of dispersion interactions, polarization of both the adsorbent and the adsorbate, interplay between H-bonding and tetrel bonding, and induced dipole/Coulomb interactions inherent to methane physisorption require computational treatment at the highest possible level of theory. In this study, we employed the smallest reasonable computational model, a maquette of porous carbon surfaces with a central site for substitution and methane binding. The most accurate predictions of methane adsorption energetics were achieved by electron-correlated molecular orbital theory CCSD(T) and hybrid density functional theory MN15 calculations employing a saturated, all-electron basis set. Here, the characteristic geometry of methane adsorption on a carbon surface (“lander approach”) arises due to bonding interactions of the adsorbent π-system with the proximal H–C bonds of methane, in addition to tetrel bonding between the antibonding orbital of the distal C–H bond and the central atom of the maquette (C, B, or N). The polarization of the electron density, structural deformations, and the comprehensive energetic analysis clearly indicate a ∼3 kJ mol –1 preference for methane binding on the N-substituted maquette. The B-substituted maquette showed a comparable or lower binding energy than the unsubstituted, pure C model, depending on the level of theory employed. The calculated thermodynamic results indicate a strategy for incorporating electron-enriched substitutions (e.g., N) into carbon materials as a way to increase methane storage capacity over electron-deficient (e.g., B) modifications. The thermochemical analysis was revised for establishing a conceptual agreement between the experimental isosteric heat of adsorption and the binding enthalpies from statistical thermodynamics principles.

03 NATURAL GAS↗

Investigation of $\gamma$-irradiated polyvinylidene fluoride and its acute toxicity

Polyvinylidene fluoride (PVDF) is an industrial polymer with many applications. The EPR spectrum of polyvinylidene fluoride (PVDF) irradiated by γ-rays at low temperature (-196°C) provides evidence for the presence of radicals formed by loss of a fluorine atom from PVDF, and polyenyl radicals. Polyenyl radicals dominate in the polymer irradiated at room temperature. Peroxide radicals formed as a result of oxidation of primary free radicals by atmospheric oxygen are observed in the EPR spectrum of a polymer irradiated in air, or when air is introduced to a polymer irradiated in a vacuum. Carbonization and oxidation of γ-irradiated PVDF were observed in the XPS spectrum. A model for polymer degradation based on correlated G3(MP2) molecular orbital theory electronic structure calculations is described in this report. The presence of HF in the products and a low yield of products with the composition CxFyHz indicate the absence of noticeable thermal destruction of the main polymer chain upon heating to 220°C of both the initial and γ-irradiated PVDF. Although any HF released during the PVDF dehydrofluorination would be toxic, a study of the acute toxicity of γ-irradiated PVDF showed that the polymer does not exhibit toxicity after a single intragastric route of administration to mice.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The scattering of low energy positrons by helium

Kohn's variational method is used to calculate the positron-helium scattering length and low energy S-wave phase shifts for a quite realistic Hylleraas type of helium function containing an electron-electron correlation term. The zero energy wavefunction is used to calculate the value of the annihilation rate parameter Z sub eff. All the results are significantly different from those for Drachman's helium model B, but are in better agreement with the available experimental data.

Humberston, J. W.↗

Langevin Dynamics modeling of gas-phase ion-ion recombination (Final Technical Report)

A self-consistent trajectory simulation approach to model MN reactions (Fig. 1) which incorporates the probability of electron transfer as a Monte Carlo operator (Fig. 2) was developed and published as Liu et al. J. Chem. Phys. 159, 114111 (2023). The electron transfer probability p ET estimated using the two-state Landau-Zener (LZ) theory was incorporated into classical trajectory simulations to elicit predictions of MN reaction cross-section σ (vacuum) or rate constant β (finite pressure). Electronic structure calculations with multireference configuration interaction (MRCI) and large correlation consistent basis sets were used to derive inputs to the LZ theory. The key advance of our trajectory simulation approach is the incorporation of electron transfer probability and the inclusion of the effect of ion-neutral interactions on MN using a Langevin representation of the effect of neutral gas on ions. For H + – H - and Li + – H(D) - pairs, our approach quantitatively agrees with measured speed-dependent cross-sections for up to ~10 5 m/s. For the ion pair Ne + – Cl - , our predictions of the MN rate constant at ~1 torr are a factor of ~2 – 3 higher than the experimentally measured value. Similarly, for Xe + – F - in the pressure range of ~20000 – 80000 Pa, our predictions of the MN rate constant are ~20% lower but are in excellent qualitative agreement with experimental data. The paradigm of using trajectory simulations to self-consistently model MN reactions is the basis for inclusion of additional non-classical, and static magnetic and electric field effects. Subsequent work, published as Roy et al. focused on modeling recombination rate constant for three ion pairs (rare gas Ar + cation and halide anions): Ar + – Cl - , Ar + – Br - , Ar + – I - , 2) considering spin-orbit couplings in the electronic structure calculations to obtain high-fidelity estimates of the electron transfer probability and incorporated within the classical trajectory simulations to elicit predictions. In addition to calculations of ion-ion recombination rate constants, a classical trajectory simulation technique (published as Roy et al. J. Chem. Phys. 162(9), 094104 (2023)) that uses quaternions to represent orientation of non-spherical particles (ions or aerosol particles) was developed to simulate the recombination of diatomic or more generally, polyatomic molecules. Finally, several other ion pairs such as Ne + – Cl - , Kr + – Cl - , were explored using the developed semi-classical trajectory simulations to understand various challenges in tackling electronic structure calculations. Using empirical approaches to parameterize the electron transfer radius, trajectory simulations were also used to probe the effect of ion number density on MN rate constant.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Electron-phonon coupling from GW perturbation theory: Practical workflow combining BerkeleyGW, ABINIT, and EPW

Here we present a workflow of practical calculations of electron-phonon (e-ph) coupling with many-electron correlation effects included using the GW perturbation theory (GWPT). This workflow combines BerkeleyGW, ABINIT, and EPW software packages to enable accurate e-ph calculations at the GW self-energy level, going beyond standard calculations based on density functional theory (DFT) and density-functional perturbation theory (DFPT). This workflow begins with DFT and DFPT calculations (ABINIT) as starting point, followed by GW and GWPT calculations (BerkeleyGW) for the quasiparticle band structures and e-ph matrix elements on coarse electron k- and phonon q-grids, which are then interpolated to finer grids through Wannier interpolation (EPW) for computations of various e-ph coupling determined physical quantities such as the electron self-energies or solutions of anisotropic Eliashberg equations, among others. A gauge-recovering symmetry unfolding technique is developed to reduce the computational cost of GWPT (as well as DFPT) while fulfilling the gauge consistency requirement for Wannier interpolation.

74 ATOMIC AND MOLECULAR PHYSICS↗

Optical conductivity and vibrational spectra of the narrow-gap semiconductor FeGa 3

Intermetallic narrow-gap semiconductors have been intensively explored due to their large thermoelectric power at low temperatures and a possible role of strong electronic correlations in their unusual thermodynamic and transport properties. Here we study the optical spectra and vibrational properties of FeGa 3 . The optical conductivity indicates that FeGa 3 has a direct band gap of ≈ 0.7 eV, consistent with density functional theory (DFT) calculations. Most importantly, we find a substantial spectral weight also below 0.4 eV, which is the energy of the indirect (charge) gap found in resistivity measurements and ab initio calculations. We find that the spectral weight below the gap decreases with increasing temperature, which indicates that it originates from the impurity states and not from the electronic correlations. Interestingly, we did not find any signatures of the impurity states in vibrational spectra. The infrared and Raman vibrational lines are narrow and weakly temperature dependent. The vibrational frequencies are in excellent agreement with our DFT calculations, implying a modest role of electronic correlations. Narrow Mössbauer spectral lines also indicate high crystallinity of the sample.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Electron Transfer Theory Elucidates the Hidden Role Played by Triethylamine and Triethanolamine during Photocatalysis

Triethylamine (TEA) and triethanolamine (TEOA) are renowned, in part, for their ability to reductively quench excited states by outer-sphere electron transfer with vast and still growing applications as sacrificial electron donors for photocatalytic systems. Upon amine oxidation, the resulting TEA •+ and TEOA •+ radical cations undergo proton transfer (or hydrogen atom transfer), resulting in the formation of a chemical reductant that has an α-carbon centered radical adjacent to the nitrogen center (TEA • and TEOA • ). In this contribution, we have electrochemically and spectroscopically characterized a set of electron acceptors which, upon accepting an electron, are a series of photocatalysts, [ReCl(R 1 R 2 -bpy)(CO) 3 ] •– , where R 1 and R 2 are electron-donating and electron-withdrawing groups in the 4,4′- and 5,5′-positions on the bipyridyl ligand. Here, we substantiated the formation of the electron donors, TE(O)A • , by spin trapping electron paramagnetic resonance spectroscopy, where TE(O)A • reacts with 2,4,6-tri-tert-butylnitrosobenzene to generate N-centered and O-centered radical adducts. Having established the chemical behaviors of the electron acceptors and donors individually, the electron transfer rate constants were determined across a 1.43 V range in driving force. The redox potential of TEA • was benchmarked to within ±80 mV on an absolute scale in V vs Fc + /Fc in CH 3 CN by using an empirical rate vs free-energy correlation, electron transfer theory, and density functional theory calculations. The equilibrium potentials for TEA • and TEOA • were determined to be −1.98 V and −1.76 V, respectively. Based on the kinetic and thermochemical analysis presented for TEA • and TEOA • , these transient radicals can be broadly considered strong homogeneous chemical reductants within the wider context of photoredox potentials. Thus, this work clarifies a frequently unnoticed secondary function for these sacrificial electron donors during photocatalysis and rationalizes the possibility of a one-photon/two-electron conversion process that is dependent on the free-energy exchange between TE(O)A • and photocatalysts.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Theoretical studies of the transition metal-carbonyl systems MCO and M(CO)2, M = Ti, Sc, and V

Large Gaussian basis sets and an electron correlation treatment are used in ab initio calculations on the transition metal-carbonyl systems MCO and N(CO)2 for M = Ti, Sc, and V. The results show that high-spin ground states are favored for the monocarbonyl molecules, and that for the dicarbonyl molecules there is a competition between high, intermediate, and low spin states which are very close in energy. Dissociation energies of 0.62 eV for Ti-CO and 1.02 for Ti(CO)2 are found, relative to the ground state Ti atomic asymptote and CO 1Sigma(+). The binding energy per carbonyl for the three metal atoms is shown to be significantly lower for the dicarbonyl than for the the monocarbonyl molecules.

Barnes, Leslie A.↗

A New Potential Energy Surface for N+O2: Is There an NOO Minimum?

We report a new calculation of the N+02 potential energy surface using complete active space self-consistent field internally contracted configuration interaction with the Dunning correlation consistent basis sets. The peroxy isomer of N02 is found to be a very shallow minimum separated from NO+O by a barrier of only 0.3 kcal/mol (excluding zero-point effects). The entrance channel barrier height is estimated to be 8.6 kcal/mol for ICCI+Q calculations correlating all but the Ols and N1s electrons with a cc-p VQZ basis set.

Walch, Stephen P.↗

A laser ultrasonics-based approach for rapid screening of high entropy alloys

This project demonstrates a laser ultrasonics-based characterization methodology for rapid metallic materials design and discovery via in situ determination of phases, microstructures and elastic properties with respect to temperature. Ultrasonic waves are strongly affected by material microstructure, and therefore serve as a facile means to probe elastic properties, phase content and their size distributions and volume fractions. In this study, a laser ultrasonic technique is used to systematically study the evolution of properties in a set of interrelated simple binary alloys and high entropy alloys (HEA). Phase transformations and microstructural changes inferred from the ultrasonic signals will be correlated with electron microscopy data and predictions using calculations of phase diagrams (CALPHAD). The non-contact and non-destructive ultrasonic testing approach developed in this study could overcome several limitations associated with current material characterization methods for materials discovery. It is expected that laser ultrasonics-based methodology developed in this work will be utilized to evaluate novel graded composition HEAs currently being developed at the Idaho National Laboratory (INL) using advanced manufacturing methods based on direct energy deposition and spark plasma sintering processes, and contribute to accelerate the discovery of HEAs.

36 MATERIALS SCIENCE↗

A laser ultrasonics-based approach for rapid screening of high entropy alloys

The primary objective of this seed project is to develop a laser ultrasonics-based characterization methodology for rapid metallic materials design and discovery via in situ determination of phases, microstructures and elastic properties with respect to temperature. Ultrasonic waves are strongly affected by material microstructure, and therefore, serve as a facile means to probe elastic properties, phase content and their size distributions and volume fractions. In this study, a laser ultrasonic technique will be used to systematically study the evolution of properties in a set of interrelated simple binary alloys and high entropy alloys (HEA). Phase transformations and microstructural changes inferred from the ultrasonic signals will be correlated with electron microscopy data and predictions using calculations of phase diagrams (CALPHAD). The non-contact and non-destructive ultrasonic testing approach developed in this study could overcome several limitations associated with current material characterization methods for materials discovery. It is expected that the products of the proposed work will be utilized to evaluate novel graded composition HEAs currently being developed at the Idaho National Laboratory (INL) using advanced manufacturing methods based on direct energy deposition and spark plasma sintering processes, and contribute to accelerate the discovery of HEAs.

36 MATERIALS SCIENCE↗

Nickel Enhances InPd-Catalyzed Nitrate Reduction Activity and N 2 Selectivity

Palladium–indium (PdIn) is a well-established bimetallic composition for reductively degrading nitrate anions, one of the most ubiquitous contaminants in the groundwater. However, the scarcity and the variable price of these rare-earth and platinum group critical metals may hinder their use for water treatment. Nickel (Ni), a nonprecious metal in the same element group as Pd, could partially replace and lower Pd usage if the resulting trimetallic composition is sufficiently catalytically active. Herein, we report the synthesis and nitrate reduction catalysis of activated carbon-supported “In-on-Pd-on-Ni” catalysts (InPdNi/AC). While bimetallic InPd/AC (0.05 wt % In, 1.3 wt % Pd) was expectedly active, trimetallic InPdNi/AC containing the same In amount, much less Pd (0.1 wt %), and 1 wt % Ni was >17 more active (k cat ≈ 20 vs 349 L min –1 g surface metal –1 ). X-ray photoelectron spectroscopy (XPS) and density functional theory (DFT) calculations showed that Pd gained electron density from Ni, correlating to the increased nitrate reduction activity. Ammonium byproduct selectivity for InPdNi/AC (18% at 50% nitrate conversion) was lower compared to that of InPd/AC (48%), suggestive of the higher surface coverage of NO or its greater reactivity with NO 2 – , which led to more N 2 . Accounting for the catalyst precursor, manufacturing costs, and spent metal recovery, we calculated that Ni incorporation lowered the net catalyst cost significantly (from $\$$1028/kg to $\$$170/kg). The trimetallic composition lowered, by ∼26 times, the catalyst cost of a stirred tank reactor sized to the same treatment capacity as that for the bimetallic case. In conclusion, the results demonstrate that the partial replacement of the precious metal with an earth-abundant one leads to a higher efficiency and lower cost denitrification catalyst, via a material strategy that should be beneficial for other clean-water catalytic systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Comparison of state-interaction and spinor-representation calculations of spin-orbit coupling within exact two-component coupled-cluster theories

Here a benchmark study of state-interaction and spinor-representation calculations of spin-orbit coupling using the exact two-component Hamiltonians with atomic mean-field integrals (the X2CAMF schemes) at the equation-of-motion coupled-cluster singles and doubles level is reported. We adopt a version of the X2CAMF scheme with spin-orbit integrals correct to first order for the state-interaction calculations and a version correct to infinite order for the spinor-representation calculations. The differences between the state-interaction calculations with minimum active spaces to account for (quasi-)degeneracy and the spinor-representation calculations thus correspond to higher-order spin-orbit contributions. This state-interaction approach with scalar-relativistic effects accurately included in the reference functions and in the spin-orbit integrals is shown to exhibit robust performance for elements across the periodic table. On the other hand, the more rigorous spinor representation shows more rapid convergence with respect to the number of correlated electrons and is the preferred choice for accurate calculations for heavy elements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗