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,171 records · Page 65

Automation of Active Space Selection for Multireference Methods via Machine Learning on Chemical Bond Dissociation

Predicting and understanding the chemical bond is one of the major challenges of computational quantum chemistry. Kohn–Sham density functional theory (KS-DFT) is the most common method, but approximate density functionals may not be able to describe systems where multiple electronic configurations are equally important. Multiconfigurational wave functions, on the other hand, can provide a detailed understanding of the electronic structures and chemical bonds of such systems. In the complete active space self-consistent field (CASSCF) method, one performs a full configuration interaction calculation in an active space consisting of active electrons and active orbitals. However, CASSCF and its variants require the selection of these active spaces. This choice is not black box; it requires significant experience and testing by the user, and thus active space methods are not considered particularly user-friendly and are employed only by a minority of quantum chemists. Our goal is to popularize these methods by making it easier to make good active space choices. Here, we present a machine learning protocol that performs an automated selection of active spaces for chemical bond dissociation calculations of main group diatomic molecules. The protocol shows high prediction performance for a given target system as long as a properly correlated system is chosen for training. Good active spaces are correctly predicted with a considerably better success rate than random guess (larger than 80% precision for most systems studied). Our automated machine learning protocol shows that a “black-box” mode is possible for facilitating and accelerating the large-scale calculations on multireference systems where single-reference methods such as KS-DFT cannot be applied.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Subnanometer Scale Mapping of Hydrogen Doping in Vanadium Dioxide

Hydrogen donor doping of correlated electron systems such as vanadium dioxide (VO 2 ) profoundly modifies the ground state properties. The electrical behavior of H x VO 2 is strongly dependent on the hydrogen concentration; hence, atomic scale control of the doping process is necessary. It is however a nontrivial problem to quantitatively probe the hydrogen distribution in a solid matrix. As hydrogen transfers its sole electron to the material, the ionization mechanism is suppressed. Here in this study, a methodology mapping the doping distribution at subnanometer length scale is demonstrated across a H x VO 2 thin film focusing on the oxygen-hydrogen bonds using electron energy loss spectroscopy (EELS) coupled with first-principles EELS calculations. The hydrogen distribution was revealed to be nonuniform along the growth direction and between different VO 2 grains, calling for intricate hydrogenation mechanisms. Our study points to a powerful approach to quantitatively map dopant distribution in quantum materials relevant to energy and information science

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Ultrafast photoinduced band splitting and carrier dynamics in chiral tellurium nanosheets

Trigonal tellurium (Te) is a chiral semiconductor that lacks both mirror and inversion symmetries, resulting in complex band structures with Weyl crossings and unique spin textures. Detailed time-resolved polarized reflectance spectroscopy is used to investigate its band structure and carrier dynamics. The polarized transient spectra reveal optical transitions between the uppermost spin-split H4 and H5 and the degenerate H6 valence bands (VB) and the lowest degenerate H6 conduction band (CB) as well as a higher energy transition at the Lpoint. Surprisingly, the degeneracy of the H6 CB (a proposed Weyl node) is lifted and the spin-split VB gap is reduced upon photoexcitation before relaxing to equilibrium as the carriers decay. Using ab initio density functional theory (DFT) calculations, we conclude that the dynamic band structure is caused by a photoinduced shear strain in the Te film that breaks the screw symmetry of the crystal. The band-edge anisotropy is also reflected in the hot carrier decay rate, which is a factor of two slower along the c-axis than perpendicular to it. The majority of photoexcited carriers near the band-edge are seen to recombine within 30 ps while higher lying transitions observed near 1.2 eV appear to have substantially longer lifetimes, potentially due to contributions of intervalley processes in the recombination rate. These new findings shed light on the strong correlation between photoinduced carriers and electronic structure in anisotropic crystals, which opens a potential pathway for designing novel Te-based devices that take advantage of the topological structures as well as strong spin-related properties.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Error-mitigated nonorthogonal quantum eigensolver via shadow tomography

We present a shadow-tomography-enhanced nonorthogonal quantum eigensolver (NOQE) for more efficient and accurate electronic structure calculations on near-term quantum devices. By integrating shadow tomography into the NOQE, the measurement cost scales linearly rather than quadratically with the number of reference states, while also reducing the required qubits and circuit depth by half. This approach enables extraction of all matrix elements via randomized measurements and classical postprocessing. We analyze its sample complexity and show that, for small systems, it remains constant in the high-precision regime, while for larger systems, it scales linearly with the system size. We further apply shadow-based error mitigation to suppress noise-induced bias without increasing quantum resources. Demonstrations on the hydrogen molecule in the strongly correlated regime achieve chemical accuracy under realistic noise, showing that our method is both resource-efficient and noise-resilient for practical quantum chemistry simulations in the near term.

quantum algorithms & computation↗

Thermal hybrid exchange-correlation density functional for improving the description of warm dense matter

Finite-temperature density functional theory (FT-DFT) has become a standard tool for first principles calculations of the properties of warm dense matter (WDM) relevant to high-energy density physics (HEDP) applications. In this work we introduce theoretical grounds of thermal hybrid exchange-correlation (XC) functionals within the generalized Mermin-Kohn-Sham scheme for an improved description of WDM. Building on the previously developed KSDT (Karasiev-Sjostrom-Dufty-Trickey) [Karasiev et al., Phys. Rev. Lett. 112, 076403 (2014)] local density approximation (LDA) and the KDT16 (Karasiev-Dufty-Trickey 2016) [Karasiev et al., Phys. Rev. Lett. 120, 076402 (2018)] generalized-gradient approximation (GGA) XC free-energy density functionals, we construct a new thermal hybrid XC functional, referred to here as KDT0. The KDT0 model at low temperature reduces to the popular ground-state PBE0 hybrid due to properties of the used KDT16 density functional approximation. Application to static calculations of electronic band gap and band structure at a wide range of temperatures for various systems of interest to HEDP show that KDT0 provides a significant improvement to the lower LDA and GGA rung XC functionals and to the ground-state PBE0 hybrid.

36 MATERIALS SCIENCE↗

Anion-Activated Bases and Nucleophiles Characterized by Photoelectron Spectroscopy

In this study, negative ion photoelectron spectra at 20 K along with ab initio (CCSD(T)) and M06-2X density function theory calculations are reported for a series of six basic and nucleophilic pyridine derivatives with an anionic substituent (i.e., 3- and 4-PyrBX 3 –, where X = F, 4-t-BuC 6 H 4 , 4- MeOC 6 H 4 , and 3,5-(MeO) 2 C 6 H 3 ). Vertical detachment energies spanning from 4.50 – 5.85 eV are observed and well reproduced by M06-2X/aug-cc-pVTZ and CCSD(T)/maug-cc-pVTZ computations. Surprisingly, the VDEs are found to correlate with the S N 2 reactivity of the PPh 4 + salts of the anions with 1-iodooctane. Spin densities of the neutral structures formed from the anionic species are also reported and indicate where the excess electron density in the anions resides and likely cation-interaction sites in their corresponding salts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

β decay of Ba 141

The β-decay strength function of nuclides produced in fission is important as it dictates the distribution of decay energy between electrons, neutrinos, and γ rays and so is critical for calculating decay heat in reactors and for estimating the reactor antineutrino spectrum. Several experimental techniques are available to determine this strength function, including electron spectroscopy, γ-ray calorimetry (TAGS spectroscopy), and detailed, high-resolution spectroscopy with modern large high-purity germanium arrays. This work investigates the decay of the well-known and strongly produced fission fragment 141 Ba. A beam of 141 Cs was implanted at the target position of the Gammasphere and the subsequent decay of the daughter 141 Ba was studied. Extensive decay spectroscopy was possible up to the decay Q value of 3.197(7) MeV, including a significant extension of the level scheme and detailed angular correlation measurements for all levels with greater than 0.25% β feeding. The distribution of the β-decay strength was then inferred and compared to previous calorimetric studies. In conclusion, the agreement was excellent and provides a benchmark for comparing strength function methods and data for a more detailed understanding of the structure of 141 La.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A simple interpretation of the Fe2(-) photoelectron spectrum

The photoelectron spectrum of Fe2(-) can be simply interpreted in terms of electron detachment from the mildly bonding 4s sigma sup * sub u orbital of a (4s sigma sub s)(sup 2) (4s sigma sup * sub u)(sup 2) (3d)(sup 13) anion. This interpretation implies a (4s sigma sub g)(sup 2) (4s sigma sup * sub u)(sup1) (3d)(sup13) configuration for the ground state of Fe2, correlating with one ground state (4s2 3d6) and one excited state (4s1 3d7) Fe atom. A comparison of the bond length and vibrational frequency of Fe2 to values for transition metal dimers containing single 4s-4s bonds is suggestive of 3d-3d bonding in this molecule. The results of preliminary full-valence configuration interaction calculations provide stron support for the proposed Fe2 and Fe2(-) configurations.

Leopold, Doreen G.↗

Strengthening model development and effects of low diffusing solutes to coarsening resistance in aluminum alloys

A modified Orowan strengthening model is proposed to account for finite rod-shaped precipitates with hemispherical caps in aluminum alloy systems. A combined computational and experimental approach is used to study the influences of anisotropic Orowan looping and solute-dislocation interaction on temperature-dependent yield strength. Here, the strengthening model is validated with a dataset containing 297 experimental precipitate geometries, chemistries, temperatures, and strength measurements, and achieves a strong predictive correlation of 0.8713 with experimentally measured yield strengths. Under conditions that are applicable to coarsening, constant particle volume fraction and aspect ratio, the model predicts that short rod precipitates provide far superior strengthening effects compared to plate precipitates. A cast Al-Si-Mg-Cu alloy with rod-shaped Q-phase (Al 3 Cu 2 Mg 9 Si 7 ) precipitates was developed with a novel chemistry exploring the use of low-diffusivity elements (Mn, Ni, V, Zr) to limit precipitate coarsening. The thermodynamic behavior of Mn, Ni, V, and Zr across the Q-phase interface is examined using transmission electron microscopy (TEM), first-principles density-functional theory (DFT) calculations. and atom-probe tomography (APT). DFT calculations utilizing TEM identified Q-phase/Al-matrix interfaces show that Mn, Ni, V, and Zr preferentially segregate to the Q-phase precipitate boundaries which suggests inhibition of precipitate coarsening and higher strengths after temperature exposure. Atom-probe tomography confirms solute atom partitioning/segregation at the Q-phase/Al-matrix interface, found in the modified commercial AS7GU alloy (A356 +0.5%Cu), which supports these observations.

36 MATERIALS SCIENCE↗

Quantum Spin Probe of Single Charge Dynamics

Electronic defects in semiconductors form the basis for emerging quantum technologies, but many defect centers are difficult to access at the single-particle level. A method for probing optically inactive spin defects would reveal semiconductor physics at the atomic scale and advance the study of new quantum systems. Here, we exploit the intrinsic correlation between the charge and spin states of defect centers to measure the charge populations and dynamics of single substitutional nitrogen spin defects in diamond. By probing their steady-state spin population, read out at the single-defect level with a nearby nitrogen vacancy center, we directly measure the defect ionization-corroborated by first-principles calculations-an effect we do not have access to with traditional coherence-based quantum sensing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Highly Accurate Quartic Force Field and Rovibrational Spectroscopic Constants for the Azirinyl Cation (c-C2NH2+) and Its Isomers

Azirinyl cation is an aromatic cyclic molecule that is isoelectronic with cyclopropenylidene, c-C3H2, and c-C3H3+. Cyclopropenylidene has been shown to be ubiquitous, existing in many different astrophysical environments. Given the similar chemistry between C and N, and the relative abundances between C and N in astrophysical environments, it is expected that there should be aromatic ringed molecules that incorporate N in the ring, but as yet, no such molecule has been identified. To address this issue, the present study uses high levels of electronic structure theory to compute a highly accurate quartic force field (QFF) for azirinyl cation and its two lowest lying isomers, the cyanomethyl and isocyanomethyl cations. The theoretical approach uses the singles and doubles coupled-cluster method that includes a perturbative correction for connected triple excitations, CCSD(T), together with extrapolation to the one-particle basis set limit and corrections for scalar relativity and core-correlation. The QFF is then used in a second order vibrational perturbation theory analysis (VPT2) to compute the fundamental vibrational frequencies and rovibrational spectroscopic constants for all three C2NH2+ isomers. The reliability of the VPT2 vibrational frequencies is tested by comparison to vibrational configuration interaction (VCI) calculations and excellent agreement is found between the two approaches. Fundamental vibrational frequencies and rovibrational spectroscopic constants for all singly substituted 13C, 15N, and D isotopologues are also reported. It is expected that the highly accurate spectroscopic data reported herein will be useful in the identification of these cations in high resolution experimental or astronomical observations.

Partha P Bera↗

Probing configurational disorder in ZnGeN 2 using cluster-based Monte Carlo

ZnGeN 2 is sought as a semiconductor with comparable lattice constant to GaN and tunable band gap for integration in optoelectronic devices. Configurational disorder on the cation sublattice of ZnGeN 2 can strongly modify the electronic structure compared to the ordered material, and both ordered and disordered forms of ZnGeN 2 are candidates for light-emitting diodes and other emitters. The nonisovalent character of the disordered species (Zn 2+ and Ge 4+ ) subjects the cation ordering to strong short-range order effects. To model these effects, we use Monte Carlo (MC) simulations utilizing a cluster expansion to approximate formation enthalpy. Representative disordered configurations in 1024-atom supercells are relaxed in density functional theory calculations. From the MC structures, we extract a short-range order parameter (the N-cation coordination motif), and two long-range order parameters (Bragg-Williams and stretching parameters), and examine their correlations. We perform a thermodynamic integration to determine the mixing entropy and free energy. ZnGeN 2 exhibits a first-order phase transition with pronounced discontinuities in enthalpy and entropy, as well as in the structural order parameters. We discuss the relationship between the effective temperature used in the MC simulation and the growth temperatures in experiment in relation to the crossover from the nonequilibrium to the equilibrium growth regime. This work expands on current models of site disorder in ZnGeN 2 and provides atomic structure models with a systematic variation of the degree of cation disorder.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Relativistic magnetospheric electrons: Lower ionospheric conductivity and long-term atmospheric variability

Long term observations of relativistic electrons in the earth's outer magnetosphere show a strong solar cycle dependence with a prominent intensity maximum during the approach to solar minimum. This population therefore closely corresponds to the presence of high speed solar wind streams emanating from solar coronal holes. Using a numerical code, the precipitating electron energy deposition in the earth's upper and middle atmosphere were calculated. Observed events (typically persisting several days) would have maximum effect in the 40 to 60 km altitude range with peak energy depositions greater than 110 keV/cu cm-s. It is suggested that this electron population could play an important long term role in modulating lower D region ionization and middle atmospheric ozone chemistry. Methods are described of observing middle atmospheric and lower ionospheric effects of the electrons including balloon, riometer, and space-based ozone sensor systems. A particularly promising approach may involve the monitoring of global Schumann resonance modes which are sensitive to global changes in the properties of the earth-ionosphere cavity. Present work indicates that Schumann resonance properties are moderately correlated with the flux of precipitating relativistic electrons thus offering the possibility of continuously monitoring this aspect of magnetosphere-atmosphere coupling.

Baker, D. N.↗

Equilibrium and Dynamics Properties of Poly(oxyethylene) Melts and Related Poly(alkylethers) from Simulations and Ab Initio Calculations

Molecular dynamics simulations of POE melts have been performed utilizing a potential force field parameterized to reproduce conformer energies and rotational energy barriers in dimethoxyethane as determined from ab initio electronic structure calculations. Chain conformations and dimensions of POE from the simulations were found to be in good agreement with predictions of a rotational isomeric state (RIS) model based upon the ab initio conformational. energies. The melt chains were found to be somewhat extended relative to chains at theta conditions. This effect will be discussed in light of neutron scattering experiments which indicate that POE chains are extended in the melt relative to theta solutions. The conformational characteristics of POE chains will also be compared with those of other poly (alkylethers), namely poly(oxymethylene), poly(oxytrimethylene) and poly(oxytetramethylene). Local conformational dynamics were found to be more rapid than in polymethylene. Calculated C-H vector correlation times were found to be in reasonable agreement with experimental values from C-13 NMR spin-lattice relaxation times. The influence of ionic salts on local conformations and dynamics will also be discussed.

Smith, Grant D.↗

Assessing Ground State Energy of Molecules and Energy Profile of the NH3 Capturing CO2 System Using the Quantum Computing Algorithms

Molecule size correlates with the number of electrons on electronic energies and strength of anharmonicity on vibrational properties, however, it is challenging to address using classical computing. In this study, variational quantum eigensolver (VQE) algorithm was implemented on a quantum simulator to quantify electronic and vibrational energies and reaction pathways of CO2 + NH3 = NH2COOH. The VQE-based Hartree-Fock-Embedding algorithm was adopted to benchmark electronic energies for a series of molecules (doi.org/10.1063/5.0188249) and quantify the reaction energy profile of the CO2 capture reaction (doi.org/10.1116/5.0137750). The generated reaction profile is in good agreement with the classical high-level Coupled-Cluster-Singles-and-Doubles (CCSD) results. The quantum computing algorithm also helps enhance the calculation of vibrational ground-state energies by considering the many-body coupling using the Vibrational Self-Consistent Field method, providing results for CO2 and NH3 molecules with accuracy comparable to the direct diagonalization method. Our approach indicates quantum computing can be applied to solve practical problems.

Lee, Yueh-Lin↗

Probing basal planes and edge sites in polygranular nuclear graphite by gas adsorption: Estimation of active surface area

Polygranular nuclear graphite, manufactured at 2800–3000 °C from a carbonized filler and binder mix, has high graphitization degree, crystalline structure close to perfect graphite, and about 20% porosity. The pore surfaces expose large regions of rough, defective prismatic edges of graphite crystallites which are the locus of graphite materials surface sites active for oxidation, chemisorption, and electron transfer. However, we show that high-resolution N 2 and Kr first monolayer adsorption on polygranular graphite (P/P 0 < 0.015) occurs in many ways like adsorption on graphitized carbon blacks. This proves the presence of energetically homogeneous basal planes domains in graphite porosity, which was not fully acknowledged before. Using classical analysis methods (Langmuir, Hill-de Boer, adsorption potential distribution) we quantity the basal plane area (BPA) of several polygranular graphite types and correlate it with their microstructure. Finally, we propose that gas adsorption is uniquely positioned to reliably estimate the active surface area of graphite by subtracting BPA from the BET total surface area. Direct estimation based on adsorption is preferable to indirect calculations based on microstructural information.

36 MATERIALS SCIENCE↗

Self-Consistent Field Methods for Excited States in Strong Magnetic Fields: A Comparison Between Energy - and Variance-based Approaches

An efficient implementation of geometrical derivatives at the Hartree–Fock (HF) and current-density functional theory (CDFT) levels is presented for the study of molecular structure in strong magnetic fields. The required integral derivatives are constructed using a hybrid McMurchie–Davidson and Rys quadrature approach, which combines the amenability of the former to the evaluation of derivative integrals with the efficiency of the latter for basis sets with high angular momentum. In addition to its application to evaluating derivatives of four-center integrals, this approach is also applied to gradients using the resolution-of-the-identity approximation, enabling efficient optimization of molecular structure for many-electron systems under a strong magnetic field. The CDFT contributions have been implemented for a wide range of density functionals up to and including the meta-GGA level with current-density dependent contributions and (range-separated) hybrids for the first time. Illustrative applications are presented to the OH and benzene molecules, revealing the rich and complex chemistry induced by the presence of an external magnetic field. Challenges for geometry optimization in strong fields are highlighted, along with the requirement for careful analysis of the resulting electronic structure at each stationary point. The importance of correlation effects is examined by comparison of results at the HF and CDFT levels. The present implementation of molecular gradients at the CDFT level provides a cost-effective approach to the study of molecular structure under strong magnetic fields, opening up many new possibilities for the study of chemistry in this regime.

Comparison↗