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 343 records · Page 19

P-Wave Electron-Hydrogen Scattering

A variational wave function incorporating short range correlations via Hylleraas type functions plus long-range polarization terms of the polarized orbital type but with smooth cut-off factors has been used to calculate P-wave phase shifts for electron-hydrogen scattering. This approach gives the direct r(exp -4) potential and a non-local optical potential which is definite. The resulting phase shifts have rigorous lower bounds and the convergence is much faster than those obtained without the modification of the target function. Final results will be presented at the conference.

Bhtia, Anand↗

Capturing 3D atomic defects and phonon localization at the 2D heterostructure interface

The three-dimensional (3D) local atomic structures and crystal defects at the interfaces of heterostructures control their electronic, magnetic, optical, catalytic, and topological quantum properties but have thus far eluded any direct experimental determination. Here, we use atomic electron tomography to determine the 3D local atomic positions at the interface of a MoS 2 -WSe 2 heterojunction with picometer precision and correlate 3D atomic defects with localized vibrational properties at the epitaxial interface. We observe point defects, bond distortion, and atomic-scale ripples and measure the full 3D strain tensor at the heterointerface. By using the experimental 3D atomic coordinates as direct input to first-principles calculations, we reveal new phonon modes localized at the interface, which are corroborated by spatially resolved electron energy-loss spectroscopy. We expect that this work will pave the way for correlating structure-property relationships of a wide range of heterostructure interfaces at the single-atom level.

42 ENGINEERING↗

Quantum Monte Carlo Calculations of Chemical Binding and Reactions

The auxiliary field quantum Monte Carlo method developed by the PIs has been shown to provide the most accurate description of strongly correlated electronic systems, from molecules to solids. Unlike other explicitly many‐body approaches, the quantum Monte Carlo method scales as a low order polynomial of systems size, similar to mean‐field methods such as density functional theory. However, the auxiliary field quantum Monte Carlo algorithm is significantly more expensive than traditional density functional calculations. This creates a bottleneck for applications to extended systems, such as large molecules and solids. One principal objective of this proposal was to develop new auxiliary field quantum Monte Carlo computational strategies to achieve improved scaling with system size, using downfolding and localization schemes, without sacrificing the predictive power of the calculations. A second goal is to extend the reach of auxiliary field quantum Monte Carlo to calculate excited states. This final report summarizes what has been achieved during the course the project toward these goals.

97 MATHEMATICS AND COMPUTING↗

Cross-correlation spectra in interacting quantum dot systems

Two-color spin-noise spectroscopy of interacting electron spins in singly charged semiconductor quantum dots provides information on the interquantum dot interactions. We investigate the spin cross-correlation function in a quantum dot ensemble employing a modified semiclassical approach. Spin-correlation functions are calculated using a Hamilton quaternion approach that maintains local quantum mechanical properties of the spins. This method takes into account the effects of the nuclear-electric quadrupolar interactions, the randomness of the coupling constants, and the variation of the electron g factor on the spin-noise power spectra. We demonstrate that the quantum dot ensemble can be mapped on an effective two-quantum dot problem and discuss how the characteristic length scale of the interdot interaction modifies the low-frequency cross-correlation spectrum. Here, we argue that details on the interaction strength distribution can be extracted from the cross-correlation spectrum when applying a longitudinal or a transversal external magnetic field.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

From Electronic Structure to Ion Transport: Photoelectron Spectroscopy and Molecular Dynamics Simulations Reveal the Role of Anions in Lithium Battery Electrolytes

Electrolyte anions are pivotal for lithium battery performance, yet their fundamental electronic structural properties are not well understood. In this work, we employ a combination of negative-ion photoelectron spectroscopy (NIPES), ab initio calculations, and molecular dynamics (MD) simulations to investigate the electronic structures of three representative electrolyte anions. This multiscale approach enables us to elucidate how their intrinsic electronic properties govern anion–solvent interactions in gas-phase clusters, as well as lithium-ion (Li + ) solvation structures and ion transport behavior in the condensed phase. NIPES reveals that difluoro(oxalato)borate (DFOB – ), bis(fluorosulfonyl)imide (FSI – ), and bis(oxalato)borate (BOB – ) all exhibit high electron binding energies, with vertical/adiabatic detachment energies increasing from DFOB – (6.09/5.70 eV) to FSI – (6.80/6.10 eV) to BOB – (6.82/6.40 eV), correlating with enhanced oxidation stability. Ab initio calculations reveal that DFOB – /FSI – –solvent complexes bind Li + ∼ 10 kcal/mol stronger than BOB – series, aligning with the strength of a Li + –anion model. DFOB – exhibits pronounced charge localization on both oxygen and fluorine atoms, enabling their involvement in Li + coordination. In contrast, fluorine atoms in FSI – are largely electron-depleted and remain excluded from direct Li + binding. MD simulations further demonstrate that LiDFOB and LiFSI systems exhibit Li + diffusion coefficients three and five times higher than those of LiBOB across four common solvents. Notably, LiFSI salt in acetonitrile (AN) exhibits the fastest Li + diffusion among 12 electrolyte systems, highlighting the synergistic effect of FSI – and AN in promoting ion mobility. In conclusion, these findings provide a molecular-level understanding of the critical roles of anion and its microsolvation in optimizing Li + diffusion dynamics, once again emphasizing the positioning of FSI – and DFOB – as prime candidates for next-generation electrolytes.

25 ENERGY STORAGE↗

Electronic structure and lattice dynamics of 1⁢𝑇−VSe 2 : Origin of the three-dimensional charge density wave

To characterize in detail the charge density wave (CDW) transition of 1⁢𝑇−VSe 2 , its electronic structure and lattice dynamics are comprehensively studied by means of x-ray diffraction, muon spectroscopy, angle resolved photoemission (ARPES), diffuse and inelastic x-ray scattering, and state-of-the-art first-principles density functional theory calculations. Resonant elastic x-ray scattering does not show any resonant enhancement at either V or Se, indicating that the CDW peak at the 𝐾 edges describes a purely structural modulation of the electronic ordering. ARPES experiments identify (i) a pseudogap at 𝑇 > 𝑇⁢ CDW , which leads to a depletion of the density of states in the ML-M'L' plane at 𝑇 < 𝑇⁢ CDW , and (ii) anomalies in the electronic dispersion reflecting a sizable impact of phonons on it. A diffuse scattering precursor, characteristic of soft phonons, is observed at room temperature (RT) and leads to the full collapse of the low-energy phonon (𝜔 1 ) with propagation vector (0.25 0 −0.3) r.l.u. Here, we show that the frequency and linewidth of this mode are anisotropic in momentum space, reflecting the momentum dependence of the electron-phonon interaction (EPI), hence demonstrating that the origin of the CDW is, to a much larger extent, due to the momentum dependent EPI with a small contribution from nesting. The pressure dependence of the 𝜔 1 soft mode remains nearly constant up to 13 GPa at RT, with only a modest softening before the transition to the high-pressure monoclinic C2/m phase. The wide set of experimental data is well captured by our state-of-the art first-principles anharmonic calculations with the inclusion of van der Waals corrections in the exchange-correlation functional. The comprehensive description of the electronic and dynamical properties of VSe 2 reported here adds important pieces of information to the understanding of the electronic modulations in the family of transition-metal dichalcogenides.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Exploiting a derivative discontinuity estimate for accurate G0W0 ionization potentials and electron affinities

Abstract The GW approximation has become an important tool for predicting charged excitations of isolated molecules and condensed systems. Its popularity can be attributed to many factors, including a favorable scaling and relatively good accuracy. In practical applications, the GW is often performed as a one-shot perturbation known as G 0 W 0 . Unfortunately, G 0 W 0 suffers from a strong starting point dependence and is often not as accurate as one would need. Self-consistent GW methodologies alleviate these problems but come with a marked increase in computational cost. In this manuscript, we propose the use of an estimate of the exchange-correlation derivative discontinuity to provide a remarkably good starting point for G 0 W 0 calculations, yielding ionization potentials and electron affinities with eigenvalue self-consistent GW quality at no additional cost. We assess the quality of the resulting methodology with the GW 100 benchmark set and compare its advantages over other similar methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Cautionary Tale: Failure of the Valence CASSCF to Describe the Hallmark of Hydrogen Bonding

Valence CASSCF (vCAS) calculations for the hydrogen-bonded (H2O)2, (HF)2, and HF–H2O dimers fail to predict the lengthening of the hydrogen donor bond and the corresponding red shift in the donor bond vibrational stretching frequency. Analysis of the active orbitals in the vCAS calculation reveals that one of the weakly occupied orbitals in each dimer is used to correlate the electrons in the acceptor lone pairs rather than those in the donor bonding orbitals. This has significant implications for subsequent treatments of the dynamical correlation that are based on the vCAS wavefunction.

Nguyen, Long H.↗

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

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

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

High-throughput determination of Hubbard $U$ and Hund $J$ values for transition metal oxides via the linear response formalism

DFT+U provides a convenient, cost-effective correction for the self-interaction error (SIE) that arises when describing correlated electronic states using conventional approximate density functional theory (DFT). The success of a DFT+U(+J) calculation hinges on the accurate determination of its Hubbard U and Hund J parameters, and the linear response (LR) methodology has proven to be computationally effective and accurate for calculating these parameters. This study provides a high-throughput computational analysis of the U and J values for transition metal d-electron states in a representative set of over 1000 magnetic transition metal oxides (TMOs), providing a frame of reference for researchers who use DFT+U to study transition metal oxides. In order to perform this high-throughput study, an ATOMATE workflow is developed for calculating U and J values automatically on massively parallel supercomputing architectures. Here, to demonstrate an application of this workflow, the spin-canting magnetic structure and unit cell parameters of the multiferroic olivine LiNiPO4 are calculated using the computed Hubbard U and Hund J values for Ni-d and O-p states, and are compared with experiment. Both the Ni-d U and J corrections have a strong effect on the Ni-moment canting angle. Additionally, including a O-pU value results in a significantly improved agreement between the computed lattice parameters and experiment

36 MATERIALS SCIENCE↗

Use of Density Functional Method to Study Molecular Vibrations

Use of a scale factor reduces the errors in computed vibrational frequencies (errors due to anharmonicities, basis set deficiencies, and approximate treatment of electron correlation). Additional errors come from modelling solids with small molecules. Nevertheless, we hope that the present calculations add to the understanding of surface IR measurements of the heat tiles.

Chong, Delano P.↗

Collisional Dissociation of CO: ab initio Potential Energy Surfaces and Quasiclassical Trajectory Rate Coefficients

We have generated accurate global potential energy surfaces for CO+Ar and CO+O that correlate with atom-diatom pairs in their ground electronic states based on extensive ab initio electronic structure calculations and used these potentials in quasi-classical trajectory nuclear dynamics calculations to predict the thermal dissociation rate coefficients over 5000- 35000 K. Our results are not compatible with the 20-45 year old experimental results. For CO + Ar we obtain fairly good agreement with the experimental rate coefficients of Appleton et al. (1970) and Mick and Roth (1993), but our computed rate coefficients exhibit a stronger temperature dependence. For CO + O our dissociation rate coefficient is in close agreement with the value from the Park model, which is an empirical adjustment of older experimental results. However, we find the rate coefficient for CO + O is only 1.5 to 3.3 times larger than CO + Ar over the temperature range of the shock tube experiments (8000-15,000 K). The previously accepted value for this rate coefficient ratio is 15, independent of temperature. We also computed the rate coefficient for the CO + O ex- change reaction which forms C + O2. We find this reaction is much faster than previously believed and is the dominant process in the removal of CO at temperatures up to 16,000 K. As a result, the dissociation of CO is accomplished in two steps (react to form C+O2 and then O2 dissociates) that are endothermic by 6.1 and 5.1 eV, instead of one step that requires 11.2 eV to break the CO bond.

exchange rate coefficient↗

Complex and Unusual Excited-State Relaxation Dynamics of 9,9′-Bifluorenylidene Revealed by Comprehensive Time-Resolved Spectroscopy and MRSF-TDDFT Calculations

Design of novel photochemical molecular motors often requires molecular building blocks that exhibit rather unusual photoactivity, for which conventional analyses of spectroscopic data can lead to conflicting interpretations. We here systematically investigated the excited-state relaxation dynamics of one such molecule, 9,9′-bifluorenylidene (BF), through comprehensive and complementary integration of ultrafast transient absorption (TA) and femtosecond stimulated Raman (FSRS) spectroscopies and first-principles mixed-reference spin-flip time-dependent density functional theory (MRSF-TDDFT). TA and FSRS identified two sequentially formed transients following photoexcitation. The decay kinetics of the two intermediates differ in response to excitation wavelengths and the viscosity/polarity of solvents. MRSF-TDDFT calculations reveal a direct, barrierless internal-conversion pathway from the bright Franck−Condon state to a dark S1 minimum, where the excited-state population is transiently trapped, accounting for the first transient species observed in spectroscopic experiments. Further tracking down along the PES with MRSF-TDDFT mapped out two nonradiative relaxation pathways via conical intersections that connect the dark S 1 state to three configurations in the ground-state manifolds, within which a ring structure with a C8−C8′ bond and the vibrationally excited ground-state BF were identified from spectroscopic and kinetic data. The complexity of relaxation kinetics was attributed to the flexible torsional and twisting motions about the C9−C9′ bridge bond enabled by the diradical character of the S 1 state. These findings clarify unusual photoactive relaxation dynamics stemming from a novel correlation between structural flexibility and shifting electronic characteristics, and they demonstrate the importance of integrating spectroscopic and advanced electronic structure calculation studies for judicious clarification of complex, competing relaxation pathways of excited states.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Nonempirical Range-Separated Hybrid Functional with Spatially Dependent Screened Exchange

Electronic structure calculations based on density functional theory (DFT) have successfully predicted numerous ground-state properties of a variety of molecules and materials. However, exchange and correlation functionals currently used in the literature, including semilocal and hybrid functionals, are often inaccurate to describe the electronic properties of heterogeneous solids, especially systems composed of building blocks with large dielectric mismatch. Here, in this study, we present a dielectric-dependent range-separated hybrid functional, screened-exchange range-separated hybrid (SE-RSH), for the investigation of heterogeneous materials. We define a spatially dependent fraction of exact exchange inspired by the static Coulomb-hole and screened-exchange (COHSEX) approximation used in many-body perturbation theory, and we show that the proposed functional accurately predicts the electronic structure of several nonmetallic interfaces, three- and two-dimensional, pristine, and defective solids and nanoparticles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Carrier tuning of Stoner ferromagnetism in ThCr 2 Si 2 -structure cobalt arsenides

CaCo 2–y As 2 is an unusual itinerant magnet with signatures of extreme magnetic frustration. The conditions for establishing magnetic order in such itinerant frustrated magnets, either by reducing frustration or increasing electronic correlations, is an open question. Here, we use results from inelastic neutron scattering and magnetic susceptibility measurements and density functional theory calculations to show that hole doping in Ca(Co 1–x Fe x ) 2–y As 2 suppresses magnetic order by quenching the magnetic moment while maintaining the same level of magnetic frustration. In this work, the suppression is due to tuning the Fermi energy away from a peak in the electronic density of states originating from a flat conduction band. This results in the complete elimination of the magnetic moment by $\textit{x}$ ≈ 0.25, providing a clear example of a Stoner-type transition.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Extending Orbital-Optimized Density Functional Theory to L-Edge XPS and Beyond: Spin–Orbit Coupling via Nonorthogonal Quasi-Degenerate Perturbation Theory

Quantum mechanical calculations of core electron binding energies (CEBEs) are relevant to interpreting X-ray photoelectron spectroscopy (XPS). Orbital-optimized density functional theory (OO-DFT) accurately predicts K-edge CEBEs but is challenged by the presence of significant spin–orbit coupling (SOC) at L- and higher edges involving inner-shell orbitals with nonzero angular momentum. Here, to extend OO-DFT to L-edges and higher, our method utilizes scalar-relativistic, spin-restricted open-shell OO-DFT to construct a minimal, quasi-degenerate basis of core-hole states corresponding to a chosen inner-shell (e.g., ionizing all six possible 2p spin orbitals). Nonorthogonal configuration interaction (NOCI) is then used to obtain the matrix elements of the full Hamiltonian including SOC in this quasi-degenerate model space of determinants. Using a screened 1-electron SOC operator parametrized with the Dirac-Coulomb-Breit (DCB) Hamiltonian results in doublet splitting (DS) values for third row elements that are nearly in quantitative agreement with experiment. The resulting NOCI eigenvalues are shifted by the average of the (scalar) OO-DFT CEBEs to yield CEBEs (split by SOC) corrected for dynamic correlation. Comparing calculations on gas phase molecules with experimental results establishes that NO-QDPT with the SCAN functional (NO-QDPT/SCAN), using the DCB screened 1-electron SOC operator is accurate to about 0.2 eV for L-edge CEBEs of molecules containing third row atoms. However, this NO-QDPT approach becomes less accurate for fourth-row elements starting in the middle of the 3d transition metal series, with errors increasing as atomic number increases.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

HfZr_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys hafnium-zirconium (Hf-Zr). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Hf and Zr. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions.

36 MATERIALS SCIENCE↗