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 1,009 records · Page 56

Coupled Cluster Study of the Interactions of AnO 2 , AnO 2 + , and AnO 2 2+ (An = U, Np) with N 2 and CO

We report thermochemical and spectroscopic properties for actinyl complexes involving UO 2 2+/1+/0 and NpO 2 2+/1+/0 with N 2 and CO, together with the UO 2 -O 2 , UO 2 + -O 2 , and UO 2 + -NO complexes, have been studied for the first time using an accurate composite coupled cluster approach. Two general bonding motifs were investigated, end-on (η 1 ) and side-on (η 2 ) relative to the metal center of the actinyls. For end-on CO complexes, both C-coordinated (An-C) and O-coordinated (An-O) structures were examined, with the former always being lower in energy. All of the η 1 complexes were calculated to be stable, with dissociation energies ranging from 2 to 36 kcal/mol, except for that of UO 2 + -O 2 (the η 1 orientation for UO 2 + -NO was not amenable to single reference coupled cluster). In agreement with a previous study, the η 2 structure for UO 2 + -O 2 was calculated to be relatively strongly bound, by 22.3 kcal/mol in this work. The closely related NO complex, however, had a calculated dissociation energy of just 4.0 kcal/mol. The binding energy of O 2 to neutral UO 2 in a η 2 orientation was calculated to be very strong, 75.4 kcal/mol, and strongly resembled a UO 2 + (O 2 – ) complex at equilibrium. The N–N and C–O bonds were found to be somewhat activated for all the side-on (η 2 ) neutral An(IV) complexes, with stretching frequencies of N 2 or CO being red-shifted by as much as 480 cm –1 with a 0.06 Å bond length elongation. Dissociation energies for the η 1 complexes are strongly correlated with the extent of electron transfer from ligand to actinyl. The nature of bonding in the actinyl complexes is examined using natural resonance theory (NRT). The correlation between bonding motif and small molecule activation is in agreement with experiments in condensed phases.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mutual spin-phonon driving effects and phonon eigenvector renormalization in nickel (II) oxide

The physics of mutual interaction of phonon quasiparticles with electronic spin degrees of freedom, leading to unusual transport phenomena of spin and heat, has been a subject of continuing interests for decades. Despite its pivotal role in transport processes, the effect of spin-phonon coupling on the phonon system, especially acoustic phonon properties, has so far been elusive. By means of inelastic neutron scattering and first-principles calculations, anomalous scattering spectral intensity from acoustic phonons was identified in the exemplary collinear antiferromagnetic nickel (II) oxide, unveiling strong spin-lattice correlations that renormalize the polarization of acoustic phonon. In particular, a clear magnetic scattering signature of the measured neutron scattering intensity from acoustic phonons is demonstrated by its momentum transfer and temperature dependences. The anomalous scattering intensity is successfully modeled with a modified magneto-vibrational scattering cross-section, suggesting the presence of spin precession driven by phonon. The renormalization of phonon eigenvector is indicated by the observed “geometry-forbidden” neutron scattering intensity from transverse acoustic phonon. Importantly, the eigenvector renormalization cannot be explained by magnetostriction but instead, it could result from the coupling between phonon and local magnetization of ions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Quantitative study of the effect of grain boundary parameters on the slip system level Hall-Petch slope for basal slip system in Mg-4Al

Several theoretical studies have reported that the geometry and structure of grain boundaries in polycrystalline materials could impose a significant effect on the Hall-Petch slope. However, experimental observations are primarily limited by the ability of the techniques to accurately quantify the grain boundary strength and validate these theoretical models. Using high-resolution electron backscatter diffraction (HR-EBSD), the local stress tensor ahead of a slip band blocked by a grain boundary was quantified and coupled with a continuum dislocation pile-up model to assess the barrier strength of specific grain boundaries to specific slip systems, referred to as micro-Hall-Petch coefficient. For basal slip system in a deformed Mg-4Al alloy, the micro-Hall-Petch coefficient ($k^{basal}_{μ}$) varied significantly, from 0.054 to 0.184 MPa - m 1/2 for nine different grain boundaries. These results were correlated with geometric descriptors of the respective grain boundaries, with three-dimensional GB profile additionally measured via focused ion beam milling. It was found that the angle between the two slip plane traces on the grain boundary plane was the most sensitive parameter affecting $k^{basal}_{μ}$, followed by the angle between the slip directions. A functional form for calculation of $k^{basal}_{μ}$ depending on these two angles is proposed to augment crystal plasticity constitutive models with slip resistance dependent on some measure of the grain size. The method allows a new pathway to calibrate grain size strengthening parameters in crystal plasticity models, allowing further computational investigations of the interrelationship between texture, grain morphology, and the Hall Petch effect.

36 MATERIALS SCIENCE↗

Magnetically confined surface and bulk excitons in a layered antiferromagnet

The discovery of two-dimensional van der Waals magnets has greatly expanded our ability to create and control nanoscale quantum phases. A unique capability emerges when a two-dimensional magnet is also a semiconductor that features tightly bound excitons with large oscillator strengths that fundamentally determine the optical response and are tunable with magnetic fields. Here, in this study, we report a previously unidentified type of optical excitation-a magnetic surface exciton-enabled by the antiferromagnetic spin correlations that confine excitons to the surface of CrSBr. Magnetic surface excitons exhibit stronger Coulomb attraction, leading to a higher binding energy than excitons confined in bulk layers, and profoundly alter the optical response of few-layer crystals. Distinct magnetic confinement of surface and bulk excitons is established by layer- and temperature-dependent exciton reflection spectroscopy and corroborated by ab initio many-body perturbation theory calculations. By quenching interlayer excitonic interactions, the antiferromagnetic order of CrSBr strictly confines the bound electron-hole pairs within the same layer, regardless of the total number of layers. Our work unveils unique confined excitons in a layered antiferromagnet, highlighting magnetic interactions as a vital approach for nanoscale quantum confinement, from few layers to the bulk limit.

36 MATERIALS SCIENCE↗

All-electron ab initio Bethe-Salpeter equation approach to neutral excitations in molecules with numeric atom-centered orbitals

The Bethe-Salpeter equation (BSE) based on GW quasiparticle levels is a successful approach for calculating the optical gaps and spectra of solids and also for predicting the neutral excitations of small molecules. We here present an all-electron implementation of the GW+BSE formalism for molecules, using numeric atom-centered orbital (NAO) basis sets. We present benchmarks for low-lying excitation energies for a set of small organic molecules, denoted in the literature as “Thiel’s set.” Literature reference data based on Gaussian-type orbitals are reproduced to about one millielectron-volt precision for the molecular benchmark set, when using the same GW quasiparticle energies and basis sets as the input to the BSE calculations. For valence correlation consistent NAO basis sets, as well as for standard NAO basis sets for ground state density-functional theory with extended augmentation functions, we demonstrate excellent convergence of the predicted low-lying excitations to the complete basis set limit. A simple and affordable augmented NAO basis set denoted “tier2+aug2” is recommended as a particularly efficient formulation for production calculations. We finally demonstrate that the same convergence properties also apply to linear-response time-dependent density functional theory within the NAO formalism.

Liu, Chi↗

Characterization of rattling in relation to thermal conductivity: Ordered half-Heusler semiconductors

We report that the factors that affect the thermal conductivity of semiconductors is a topic of great scientific interest, especially in relation to thermoelectrics. Key developments have been the concept of the phonon-glass-electron-crystal (PGEC) and the related idea of rattling to achieve this. We use first-principles phonon and thermal conductivity calculations to explore the concept of rattling for stoichiometric-ordered half-Heusler compounds. These compounds can be regarded as filled zinc blende materials, and the filling atom could be viewed as a rattler if it is weakly bound. We use two simple metrics, one related to the frequency and the other to bond frustration and anharmonicity. We find that both measures correlate with thermal conductivity. Lastly, this suggests that both may be useful in screening materials for low thermal conductivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Chiral effective field theory calculation of neutrino reactions in warm neutron-rich matter

Neutrino scattering and absorption rates of relevance to supernovae and neutron star mergers are obtained from nuclear matter dynamical structure functions that encode many-body effects from nuclear mean fields and correlations. In this study, we employ nuclear interactions from chiral effective field theory to calculate the density, spin, isospin, and spin-isospin response functions of warm β-equilibrium nuclear matter. We include corrections to the single-particle energies in the mean field approximation as well as vertex corrections resummed in the random phase approximation (RPA), including both direct and exchange diagrams. We find that correlations included through the RPA redistribute the strength of the response to higher energy for neutrino absorption and lower energy for antineutrino absorption. This tends to suppress the absorption rate of electron neutrinos across all relevant energy scales. In contrast, the inclusion of RPA correlations enhances the electron antineutrino absorption rate at low energy and suppresses the rate at high energy. These effects are especially important at high-density and in the vicinity of the neutrino decoupling region. Implications for heavy element nucleosynthesis, electromagnetic signatures of compact object mergers, supernova dynamics, and neutrino detection from galactic supernovae are discussed briefly.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Deep Electric Field Predictions by Drift-Reduced Braginskii Theory with Plasma-Neutral Interactions Based on Experimental Images of Boundary Turbulence

We present 2-dimensional turbulent electric field calculations via physics-informed deep learning consistent with (i) drift-reduced Braginskii theory under the framework of an axisymmetric fusion plasma with purely toroidal field and (ii) experimental estimates of the fluctuating electron density and temperature on open field lines obtained from analysis of gas puff imaging of a discharge on the Alcator C-Mod tokamak. The inclusion of effects from the locally puffed atomic helium on particle and energy sources within the reduced plasma turbulence model are found to strengthen correlations between the electric field and electron pressure. The neutrals are also directly associated with broadening the distribution of turbulent field amplitudes and increasing ${\bf E \times B}$ shearing rates. As a result, this demonstrates a novel approach in plasma experiments by solving for nonlinear dynamics consistent with partial differential equations and data without encoding explicit boundary nor initial conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effect of ohmic, mass-transfer, and kinetic resistances on linear-sweep voltammetry in a cylindrical-pore electrode

A means is presented to account for the effect of ohmic, mass-transfer, and kinetic resistances on linear-sweep voltammograms by modeling a pore in a porous matrix as a cylindrical-pore electrode, and solving the mass and charge conservation equations in the context of this geometry for the simply redox reaction O + ne(-) yield R where both O and R are soluble species. Both analytical and numerical techniques are used to solve the governing equations. The calculated peak currents and potentials are correlated by empirical formulas to the measurable parameters: sweep rate, concentration of the redox species, diffusion coefficient, conductivity of the electrolyte, and pore dimensions. Using the correlations, a methodology is established for determining if the redox reaction kinetics are irreversible or reversible (Nernstian). If the reaction is irreversible, it is shown how the standard rate constant and the transfer coefficient may be extracted from linear-sweep voltammetry data, or, if the reaction is reversible, how the number of electrons transferred may be deduced.

Weidner, John W.↗

The computation of C-C and N-N bond dissociation energies for singly, doubly, and triply bonded systems

The bond dissociation energies (De) of C2H2, C2H4, C2H6, N2, N2H2, and N2H4 are studied at various levels of correlation treatment. The convergence of De with respect to the one-particle basis is studied at the single-reference modified coupled-pair-functional (MCPF)level. At all levels of correlation treatment, the errors in the bond dissociation energies increase with the degree of multiple bond character. The multireference configuration-interaction (MRCI) De values, corrected for an estimate of higher excitations, are in excellent agreement with those determined using the size-extensive averaged-coupled-pair-functional (ACPF) method. The full-valence complete-active-space self-consistent-field (CASSCF)/MRCI calculations are reproduced very well by MRCI calculations based on a CASSCF calculation that includes in the active space only those electrons involved in the C-C or N-N bonds.

Langhoff, Stephen R.↗

The determination of accurate dipole polarizabilities alpha and gamma for the noble gases

The static dipole polarizabilities alpha and gamma for the noble gases helium through xenon were determined using large flexible one-particle basis sets in conjunction with high-level treatments of electron correlation. The electron correlation methods include single and double excitation coupled-cluster theory (CCSD), an extension of CCSD that includes a perturbational estimate of connected triple excitations, CCSD(T), and second order perturbation theory (MP2). The computed alpha and gamma values are estimated to be accurate to within a few percent. Agreement with experimental data for the static hyperpolarizability gamma is good for neon and xenon, but for argon and krypton the differences are larger than the combined theoretical and experimental uncertainties. Based on our calculations, we suggest that the experimental value of gamma for argon is too low; adjusting this value would bring the experimental value of gamma for krypton into better agreement with our computed result. The MP2 values for the polarizabilities of neon, argon, krypton and zenon are in reasonabe agreement with the CCSD and CCSD(T) values, suggesting that this less expensive method may be useful in studies of polarizabilities for larger systems.

Rice, Julia E.↗

Parameter-free treatment of a layered correlated van der Waals magnet: CrPS 4

The electronic and magnetic structure of CrPS 4 , a two-dimensional (2D) magnetic semiconductor is examined by employing the SCAN meta-GGA density functional. We find the resulting magnetic moment and band gap are in excellent agreement with experiment. From the bulk magnetic configurations, we confirm the experimentally observed A-type antiferromagnetic (A-AFM) ordered ground state with a magnetic moment of 2.78 µ B per chromium atom and band gap of 1.34 eV. To gain insight into the evolution of the ground state with layers, the total energy of each magnetic configuration is calculated for a variety of thicknesses. Monolayer CrPS 4 is predicted to be a ferromagnetic insulator with a band gap of 1.37 eV, and A-AFM for bilayer and trilayer, with band gaps of 1.35 and 1.30 eV, respectively. The electronic structure is reported for the single-, two-, and three-layer and bulk CrPS 4 . Finally, we explore the optical properties of the 2D structure and report the dielectric tensor components and Kerr parameters for the monolayer.

2-dimensional systems↗

Static Subspace Approximation for Random Phase Approximation Correlation Energies: Applications to Materials for Catalysis and Electrochemistry

Modeling complex materials using high-fidelity, ab initio methods at low cost is a fundamental goal for quantum chemical software packages. The GW approximation and random phase approximation (RPA) provide a unified description of both electronic structure and total energies using the same physics in a many-body perturbative approach that can be more accurate than generalized-gradient density functional theory (DFT) methods. However, GW/RPA implementations have historically been limited to either specific materials classes or application toward small chemical systems. Here, the static subspace approximation allows for reduced cost full-frequency GW/RPA calculations and has previously been benchmarked thoroughly for GW calculations. Here, we describe our approach to including partial occupations of electronic orbitals in full-frequency GW and RPA calculations for the study of electrocatalysts. We benchmarked RPA total energy calculations using the subspace approximation across a diverse test suite of materials for a variety of computational parameters. The benchmarking quantifies the impact of different extrapolation procedures for representing the static polarizability at infinite screened cutoff, and shows that using screened cutoffs above 20-25 Ryd result in diminishing accuracy returns for predicting RPA total energies. Additionally, for moderately sized electrocatalytic models, 2-3 times fewer computational resources are used to compute RPA total energies by representing the static polarizability with 20-30% of the static subspace basis, with an error of approximately 0.01 eV or better in RPA adsorption energy calculations. Finally, we show that for these electrochemical models RPA can shift DFT adsorption energy shifts by up to 0.5 eV and that GW can frequently shift DFT eigenvalues of surface and adsorbate states by approximately 0.5-1 eV.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Spectrally sharp magnetic excitations above the critical temperature in a frustrated Weyl semimetal

The rare-earth α-pyrochlore iridates are a prospective class of conducting frustrated magnets where electronic correlations, large spin-orbit coupling, and geometrical frustration interplay, leading to a rich set of magnetic and electronic phases. Despite their intriguing properties, the magnetic order and excitations in this fundamental class of topological quantum materials remain poorly understood due to challenges in growing large single crystals and insufficient microscopic information on their temperature-dependent phases. Here, by combining state-of-the-art thin-film synthesis, resonant elastic and inelastic X-ray scattering, spin wave analysis, and dynamical spin susceptibility calculations, we unequivocally reveal the presence of spectrally sharp, gapped magnetic excitations in Y 2 Ir 2 O 7 that surprisingly persist well above the Néel transition temperature, signaling the presence of a quasi-universal regime connected to fluctuations on frustrated lattices. This finding implies the existence of a highly unusual cooperative paramagnetic (CP) phase above the ordering temperature and offers an explanation for the puzzling high-temperature magnetic behavior observed across the family of metallic pyrochlore crystals. Understanding such magnetic excitations at technologically relevant temperatures opens up possibilities for novel topological spintronic devices.

36 MATERIALS SCIENCE↗

Lunar ion flux and energy

The dynamics of the lunar ionosphere and resulting flux of lunar ions to the lunar surface are reviewed. In the Lunar rest frame, ions formed from the neutral lunar atmosphere are accelerated by the interplanetary electric and magnetic fields. The trajectories of heavier ions are primarily along the electric field; the ion flux is in a direction perpendicular to the solar wind flow, correlated to the orientation of the interplanetary magnetic field, and impacts the surface with energies of tens of electron volts to a few keV. Thus we predict a relatively energetic (compared to thermal energies), highly directional, lunar ion flux but with the possibility that light ions such as hydrogen can execute orbits that return them to portions of the lunar surface not directly exposed to the solar wind. The effects of surface electric and magnetic fields are discussed, as is the ion energy spectrum. We calculate the trapping of atmospheric ions which impact the surface and from this the density of neutral Ar 40 in the lunar atmosphere. We show that using the acceleration model, the lunar atmosphere total neutral number density can be calculated from ion detector data.

Manka, R. H.↗

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↗

La 2 ⁢O 3 ⁢Mn 2 ⁢Se 2 : A correlated insulating layered d-wave altermagnet

Altermagnets represent a new class of magnetic phases without net magnetization, invariant under a combination of rotation and time reversal. Unlike conventional collinear antiferromagnets (AFM), altermagnets could lead to new correlated states and important material properties deriving from their nonrelativistic spin-split band structure. Indeed, they serve as the magnetic analogue of unconventional superconductors and can yield spin-polarized electrical currents in the absence of external magnetic fields, making them promising candidates for next-generation spintronics. Here, we report altermagnetism in the correlated insulator, magnetically ordered tetragonal oxychalcogenide, L⁢a 2 ⁢O 3⁢ M⁢n 2 ⁢S⁢e 2 . Symmetry analysis reveals a 𝑑 𝑥 2 −𝑦 2 -wave-like spin-momentum locking arising from the M⁢n 2 ⁢O Lieb lattice, supported by density functional theory (DFT) calculations. Magnetic measurements confirm the AFM transition below ∼166K while neutron pair distribution function analysis reveals a 2D short-range magnetic order that persists above the Néel temperature. Single crystals are grown and characterized using x-ray diffraction, optical and electron microscopy, and micro-Raman spectroscopy to confirm the crystal structure, stoichiometry, and uniformity. Furthermore, our findings establish L⁢a 2⁢ O 3 ⁢M⁢n 2⁢ S⁢e 2 as a model altermagnetic system realized on a Lieb lattice.

36 MATERIALS SCIENCE↗

Small Al clusters. I - The effect of basis set and correlation on the geometry of small Al clusters

A detailed study is presented of the basis set requirements and effects of correlation on the geometry and structure of small Al(n) clusters n = 2, 4, and 13. An effective core potential (ECP) is developed from the Huzinaga basis which accurately reproduces the all-electron results. It is found that depolarization functions are very important in computing the bond length, and that the bond contraction obtained is about as large for a 13-atom cluster as for the Al(4) rhombus structure. With d functions on each center in Al(13), a bond distance shorter than the bulk is obtained, as expected. The inclusion of correlation in Al(4) is found to be less important than the addition of polarization functions for determining the bond length, but increases substantially the calculated atomization energy. These calibration calculations allow an accurate treatment for Al(n) clusters to be designed.

Bauschlicher, Charles W., Jr.↗