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

Anode power deposition in magnetoplasmadynamic thrusters

Results of anode heat-flux and anode fail measurements from a multimegawatt self-field quasi-steady magnetoplasmadynamic (MPD) thruster are presented. Measurements were obtained with argon and helium propellants for a variety of currents and mass flow rates. Anode heat flux was directly measured with thermocouples attached to the inner surface of a hollowed section. Anode falls were determined both from floating probes and through heat flux measurements. Comparison of data acquired through either method shows excellent agreement. Anode falls varied between 4-50 V with anode power fractions reaching 70 percent with helium at 150 kW, and 50 percent with argon at 1.9 MW. The anode fall was found to correlate well with electron Hall parameters calculated from triple Langmuir and magnetic probe data collected near the anode. Two possible explanations for this result are proposed: (1) the establishment of large electric fields at the anode to maintain current conduction across the strong magnetic fields; and (2) anomalous resistivity resulting from the onset of microturbulence in the plasma. To investigate the latter hypothesis, electric field, magnetic field, and current density profiles measured in the vicinity of the anode were incorporated into Ohm's law to estimate the electrical conductivity. Results of this analysis show a substantial deviation of the measured conductivity from that calculated with classical formulas. These results imply that anomalous effects are present in the plasma near the anode.

Gallimore, A. D.↗

Bond Dissociation Energies of the Actinide Halides AnX, An = Ac–Lr and X = F–I, Utilizing Relativistic Composite Coupled Cluster Approaches

Bond dissociation energies (BDEs) have been calculated for the set of actinide halides AnX with An=Ac, Pa, and Np-Lr and X=F-I. Two composite thermochemistry methods based on the Feller-Peterson-Dixon (FPD) approach have been utilized, one involving spinor-based relativistic CCSD(T) calculations where spin-orbit (SO) was included at the orbital level and another using scalar relativistic CCSD(T) with a posteriori SO contributions based on 2-component multireference configuration interaction calculations. The method that was chosen for a given actinide halide was based on which representation yielded the best single determinant reference determinant for the coupled cluster calculation. The spinor-based method was chosen for all cases except for AmX, CmX, and BkX. Both composite approaches included contributions accounting for basis set truncation, outer-core correlation, the Gaunt interaction, and QED. The scalar FPD results, as well as the spinor-based calculations for AcF, also included higher order electron correlation up through CCSDT(Q). In addition to BDEs, CCSD(T) equilibrium bond lengths, harmonic frequencies, and vibrational anharmonicity constants are reported for all species. Last, the FPD BDEs for the fluorides were used to confirm the trend across the actinide series previously predicted by Gibson using bonding models based atomic promotion energies that provide a single 6d electron for bonding. In particular the local minimum in the BDEs at AmF is confirmed in the present calculations. Furthermore, the BDEs for LrX are predicted to be slightly larger than those of AcX, making them the largest in the actinide halide series.

Actinides↗

Perspectives on UV and x-ray photoelectron spectroscopy

Atomic photoionization using electron spectroscopy was one of the scientific topics in which Professor David A. Shirley excelled and established a research legacy. We discuss selected highlights of atomic photoionization research done by the Shirley group beginning in the 1970s. Our main topics are electron correlation satellites, photoelectron angular distributions with tunable synchrotron radiation, and developments of theory. We also discuss examples of more recent research that extends the earlier work, including double K-shell photoionization, nondipole interactions, and the pervasive role of interchannel coupling. Furthermore, we conclude with examples of current experimental and theoretical research.

74 ATOMIC AND MOLECULAR PHYSICS↗

Random Phase Approximation Correlation Energy Using Real-Space Density Functional Perturbation Theory

We present a real-space method for computing the random phase approximation (RPA) correlation energy within Kohn–Sham density functional theory, leveraging the low-rank nature of the frequency-dependent density response operator. In particular, we employ a cubic-scaling formalism based on density functional perturbation theory that circumvents the calculation of the response function matrix, instead relying on the ability to compute its product with a vector through the solution of the associated Sternheimer linear systems. We develop a large-scale parallel implementation of this formalism using the subspace iteration method in conjunction with the spectral quadrature method while employing the Kronecker product-based method for the application of the Coulomb operator and the conjugate orthogonal conjugate gradient method for the solution of the linear systems. We demonstrate convergence with respect to key parameters and verify the method’s accuracy by comparing with plane-wave results. We show that the framework achieves good strong scaling to many thousands of processors, reducing the time to solution for a lithium hydride system with 128 electrons to around 150 s on 4608 processors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Electron correlation effects in attosecond photoionization of CO 2

A technique for measuring photoionization time delays with attosecond precision is combined with calculations of photoionization matrix elements to demonstrate how multielectron dynamics affect photoionization time delays in carbon dioxide. Electron correlation is observed to affect the time delays through two mechanisms: autoionization of molecular Rydberg states and accelerated escape from a continuum shape resonance.

74 ATOMIC AND MOLECULAR PHYSICS↗

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↗

Self-assembled Ag(111) nanostructures induced by Fermi surface nesting

Scanning tunneling microscopy measurements on Ag(111)/MoS 2 reveal atomically flat preferred, or “magic,” heights occurring at 6, 10, and 14 atomic layers. We note these results are consistent with Ag growth on a variety of semiconducting substrates and correlate with electronic energy savings in electronic structure calculations of freestanding Ag(111) films. Thus, under certain conditions, Ag will spontaneously form quantized structures independent of the substrate. To explain this, we have found Fermi surface nesting vectors in the bulk Ag band structure which account for these results and the fact Ag that is gapped along the surface normal. This model extends to a range of metallic systems which exhibit electronic confinement, epitaxial growth, and minimal strain. As with Au/MoS 2 , the Ag/MoS 2 system exhibits this behavior at unusually high temperatures so that these principles might be used for control over device features at the nanometer scale under standard fabrication conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Electronic correlation induced expansion of Fermi pockets in δ-plutonium

Plutonium is a critically important material as the behavior of its 5f electrons stands midway between the metalliclike itinerant character of the light actinides and localized atomic-core-like character of the heavy actinides. The δ phase of plutonium (δ-Pu), whereas still itinerant, has a large coherent Kondo peak and strong electronic correlations coming from its near-localized character. Here, using sophisticated Gutzwiller wave function and dynamical mean-field theory correlated theories, we study the Fermi surface and associated mass renormalizations of δ-Pu together with calculations of the de Haas–van Alphen frequencies. We find a large (~200%) correlation induced volume expansion in both the hole and the electron pockets of the Fermi surface in addition to an intermediate mass enhancement. All of the correlated electron theories predict, approximately, the same hole pocket placement in the Brillouin zone, which is different from that obtained in conventional density-functional band-structure theory, whereas the electron pockets from all theories are in, roughly, the same place.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

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↗

Photoemission and dynamical mean field theory study of electronic correlations in a t 2 g 5 metal SrRhO 3 thin film

Perovskite rhodates are characterized by intermediate strengths of both electronic correlation as well as spin-orbit coupling (SOC) and usually behave as moderately correlated metals. A recent publication on epitaxial SrRhO 3 thin films reported a bad-metallic behavior and suggested the occurrence of antiferromagnetism below 100 K. In this work, we have further studied this SrRhO 3 thin film by hard x-ray photoemission spectroscopy and found a very small density of states (DOS) at the Fermi level $(E_F)$, which is consistent with the previously reported bad-metallic behavior. However, the absence of DOS persists up to room temperature, which contradicts the explanation of antiferromagnetic transition at ~100 K. We also employed electronic structure calculations within the framework of density functional theory and dynamical mean-field theory. In contrast to the photoemission results, our calculations indicate metallic behaviors of both bulk SrRhO 3 and the SrRhO 3 thin film, and a stronger correlation effect was observed in the thin film than that in the bulk. The calculated uniform magnetic susceptibility is substantially larger in the thin film than that in the bulk. We also investigated the role of SOC and found only a moderate modulation of the band structure. Hence SOC is not expected to significantly affect the electronic correlation in SrRhO 3 . Extrinsic effects of finite-thickness effects beyond our calculations and localization effects may play important roles to induce the negligible spectral weight at $(E_F)$.

36 MATERIALS SCIENCE↗

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↗

Wannier–Koopmans method calculations for transition metal oxide band gaps

The widely used density functional theory (DFT) has a major drawback of underestimating the band gaps of materials. Wannier–Koopmans method (WKM) was recently developed for band gap calculations with accuracy on a par with more complicated methods. WKM has been tested for main group covalent semiconductors, alkali halides, 2D materials, and organic crystals. Here we apply the WKM to another interesting type of material system: the transition metal (TM) oxides. TM oxides can be classified as either with d 0 or d 10 closed shell occupancy or partially occupied open shell configuration, and the latter is known to be strongly correlated Mott insulators. We found that, while WKM provides adequate band gaps for the d0 and d10 TM oxides, it fails to provide correct band gaps for the group with partially occupied d states. This issue is also found in other mean-field approaches like the GW calculations. We believe that the problem comes from a strong interaction between the occupied and unoccupied d-state Wannier functions in a partially occupied d-state system. We also found that, for pseudopotential calculations including deep core levels, it is necessary to remove the electron densities of these deep core levels in the Hartree and exchange–correlation energy functional when calculating the WKM correction parameters for the d-state Wannier functions.

36 MATERIALS SCIENCE↗

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↗