Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “electronic structure method”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Electronic structures of HgTe and CdTe surfaces and HgTe/CdTe interfaces

A Green's-function method has been used to study the surface and interface electronic structures of the II-VI compounds HgTe and CdTe. Localized surface and resonance states near the cation-terminated (100) surface of CdTe and the anion-terminated surface of HgTe have been found for the ideal surfaces. The energies and strengths of these surface states are altered by surface perturbations. The bulk states near the surface are drastically modified by the creation of the surface, but the band gaps remain unchanged. Numerical evaluation of the local densities of states at the Gamma and J points shows that, at the (100) interface of HgTe/CdTe, the previously observed surface states are no longer present. However, in the interface region, bulk states of one material penetrate some distance into the other material.

Schick, J. T.↗
DOE OSTI2025

Equation-of-motion internally contracted multireference unitary coupled-cluster theory

The accurate computation of excited states remains a challenge in electronic structure theory, especially for systems with a ground state that requires a multireference treatment. In this work, we introduce a novel equation-of-motion (EOM) extension of the internally contracted multireference unitary coupled-cluster framework (ic-MRUCC), termed EOM-ic-MRUCC. EOM-ic-MRUCC follows the transform-then-diagonalize approach, in analogy to its non-unitary counterpart. By employing a projective approach to optimize the ground state, the method retains additive separability and proper scaling with system size. We show that excitation energies are size-intensive if the EOM operator satisfies the “killer” and the projective conditions. Furthermore, we propose to represent changes in the reference state upon electron excitation via projected many-body operators that span the active orbitals and show that the EOM equations formulated in this way are invariant with respect to active orbital rotations. We test the EOM-ic-MRUCC method truncated to single and double excitations by computing the potential energy curves for several excited states of a BeH2 model system, the HF molecule, and water undergoing symmetric dissociation. Across these systems, our method delivers accurate excitation energies and potential energy curves within 5 mE h (∼0.14 eV) from full configuration interaction. Here, we find that truncating the Baker–Campbell–Hausdorff series to fourfold commutators contributes negligible errors (on the order of 10 −5 E h or less), offering a practical route to highly accurate excited-state calculations with reduced computational overhead.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Extended Gutzwiller Approximation for Nonlocal Electron-Electron and Electron-Boson Correlations (I): The Theory

Understanding electron-electron and electron-photon correlations is central to uncovering the fundamental mechanisms governing material properties, particularly in systems where strong interactions give rise to emergent phenomena such as superconductivity, magnetism, and polaritonic effects. These correlations play a pivotal role in cavity quantum materials, where hybridized light-matter states enable quantum control over electronic properties. However, capturing both local and nonlocal correlations in these systems presents a significant theoretical challenge. In this work, we extend the Gutzwiller wavefunction method to include nonlocal electron-photon and electron-electron interactions, providing a unified framework to study the intricate interplay between these effects. Our approach accurately captures the long-range correlations induced by photon exchange, enabling the exploration of exotic quantum phases and the effects of cavity coupling on electronic structure. By benchmarking the method across coupling regimes, we reveal the critical role of nonlocal correlations in stabilizing phases, such as superconducting and insulating states, that are inaccessible through local interactions alone. This generalized Gutzwiller framework offers a versatile tool for understanding and designing materials that harness the transformative potential of strong light-matter coupling.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Light–matter interaction of a molecule in a dissipative cavity from first principles

Cavity-mediated light–matter coupling can dramatically alter opto-electronic and physico-chemical properties of a molecule. Ab initio theoretical predictions of these systems need to combine non-perturbative, many-body electronic structure theory-based methods with cavity quantum electrodynamics and theories of open-quantum systems. Here, we generalize quantum-electrodynamical density functional theory to account for dissipative dynamics of the cavity and describe coupled cavity–single molecule interactions in the weak-to-strong-coupling regimes. Specifically, to establish this generalized technique, we study excited-state dynamics and spectral responses of benzene and toluene under weak-to-strong light–matter coupling. By tuning the coupling, we achieve cavity-mediated energy transfer between electronically excited states. This generalized ab initio quantum-electrodynamical density functional theory treatment can be naturally extended to describe cavity-mediated interactions in arbitrary electromagnetic environments, accessing correlated light–matter observables and thereby closing the gap between electronic structure theory, quantum optics, and nanophotonics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Predicting electronic structures at any length scale with machine learning

The properties of electrons in matter are of fundamental importance. They give rise to virtually all material properties and determine the physics at play in objects ranging from semiconductor devices to the interior of giant gas planets. Modeling and simulation of such diverse applications rely primarily on density functional theory (DFT), which has become the principal method for predicting the electronic structure of matter. While DFT calculations have proven to be very useful, their computational scaling limits them to small systems. We have developed a machine learning framework for predicting the electronic structure on any length scale. It shows up to three orders of magnitude speedup on systems where DFT is tractable and, more importantly, enables predictions on scales where DFT calculations are infeasible. Our work demonstrates how machine learning circumvents a long-standing computational bottleneck and advances materials science to frontiers intractable with any current solutions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Adaptive finite differencing in high accuracy electronic structure calculations

Abstract A multi-order Adaptive Finite Differencing (AFD) method is developed for the kinetic energy operator in real-space, grid-based electronic structure codes. It uses atomic pseudo orbitals produced by the corresponding pseudopotential codes to optimize the standard finite difference (SFD) operators for improved precision. Results are presented for a variety of test systems and Bravais lattice types, including the well-known Δ test for 71 elements in the periodic table, the Mott insulator NiO, and borax decahydrate, which contains covalent, ionic, and hydrogen bonds. The tests show that an 8th-order AFD operator leads to the same average Δ value as that achieved by plane-wave codes and is typically far more accurate and has a much lower computational cost than a 12th-order SFD operator. The scalability of real-space electronic calculations is demonstrated for a 2016-atom NiO cell, for which the computational time decreases nearly linearly when scaled from 18 to 144 CPU-GPU nodes.

Briggs, E. L. (ORCID:0000000343983492)↗
DOE OSTI2023

Graph-based quantum response theory and shadow Born–Oppenheimer molecular dynamics

Graph-based linear scaling electronic structure theory for quantum-mechanical molecular dynamics simulations [A. M. N. Niklasson et al., J. Chem. Phys. 144, 234101 (2016)] is adapted to the most recent shadow potential formulations of extended Lagrangian Born–Oppenheimer molecular dynamics, including fractional molecular-orbital occupation numbers [A. M. N. Niklasson, J. Chem. Phys. 152, 104103 (2020) and A. M. N. Niklasson, Eur. Phys. J. B 94, 164 (2021)], which enables stable simulations of sensitive complex chemical systems with unsteady charge solutions. The proposed formulation includes a preconditioned Krylov subspace approximation for the integration of the extended electronic degrees of freedom, which requires quantum response calculations for electronic states with fractional occupation numbers. For the response calculations, we introduce a graph-based canonical quantum perturbation theory that can be performed with the same natural parallelism and linear scaling complexity as the graph-based electronic structure calculations for the unperturbed ground state. Further, the proposed techniques are particularly well-suited for semi-empirical electronic structure theory, and the methods are demonstrated using self-consistent charge density-functional tight-binding theory both for the acceleration of self-consistent field calculations and for quantum-mechanical molecular dynamics simulations. Graph-based techniques combined with the semi-empirical theory enable stable simulations of large, complex chemical systems, including tens-of-thousands of atoms.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Computational quantum chemistry of metal–organic frameworks

Metal–organic frameworks (MOFs) have premium exceptional properties for a variety of functions, such as gas separation and storage and catalysis. The large variety of possible inorganometallic nodes and organic linkers provide an almost unlimited number of combinations for assembling MOFs, which makes the experimental characterization and examination of all potentially useful combinations practically impossible. Furthermore, experimental studies of MOFs typically fall short in uncovering crucial details regarding their mechanisms of action or the molecular details responsible for their functional properties, such as the nature of adsorbate binding or the structures of transition states. Computational modeling has, therefore, become an efficient and important tool for strategizing the functionalization of MOFs and explicating the mechanisms of their functions. Furthermore, we review the computational methodologies used for computational studies of MOFs, especially Kohn–Sham density functional theory and combined quantum mechanical and molecular mechanical methods for calculating their structural, electronic, and magnetic properties, as well as for understanding the mechanisms of MOFs' applications to magetic devices, thermal conduction, gas adsorption, separation, storage, and sensing, thermal catalysis, photocatalysis, and electrocatalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Spatial locality of electronic correlations in LiFeAs

In this work, we address the question of the degree of spatial nonlocality of the self-energy in the iron-based superconductors, a subject which is receiving considerable attention. Using LiFeAs as a prototypical example, we extract the self-energy from angular-resolved photoemission spectroscopy data. We use two distinct electronic structure references: density functional theory in the local density approximation and linearized quasiparticle self-consistent GW (LQSGW). We find that with the LQSGW reference, spatially local dynamical correlations provide a consistent description of the experimental data, and account for some surprising aspects of the data such as the substantial out-of-plane dispersion of the electron Fermi surface having dominant xz/yz character. Hence, correlations effects can be separated into static nonlocal contributions well described by LQSGW and dynamical local contributions. Hall effect and resistivity data are shown to be consistent with this description.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Electronic structures on swollen hydrogels

This document describes a conformable substrate that includes a hydrogel having adhesion-promoting moieties, said adhesion-promoting moieties comprising one or more catechol groups. The conformable substrate includes an array of microelectrodes bonded to the hydrogel by the adhesion-promoting moieties via the one or more catechol groups. This document also describes a method for transfer printing of an electronic structure to a hydrogel. The method includes the steps of coating a donor substrate with a film of polyacrylic acid, crosslinking the film of polyacrylic acid in a solution comprising divalent ions, patterning a microelectrode array onto the crosslinked film of polyacrylic acid, laminating an adhesive hydrogel substrate onto the donor substrate coated by the crosslinked film of polyacrylic acid comprising the patterned microelectrode array, and separating the crosslinked film of polyacrylic acid from the donor substrate in a monovalent solution.

Bettinger, Christopher J.↗
NASA NTRS1993

Ballistic-Electron-Emission Microscopy Techniques for Nanometer-scale Characterization of Interfaces

Semiconductor interface properties are among the most important phenomena in materials science and technology. The study of metal/semiconductor Schottky barrier interfaces has been the primary focus of a large research and development community for decades. Throughout the long history of interface investigation, the study of interface defect electronic properties have been seriously hindered by the fundamental experimental difficulty of probing subsurface structures. A new method, Ballistic-Electron-Emission Microscopy (BEEM), has been developed which not only enables spectroscopic probing of subsurface interface properties, but also, provides nanometer-resolution imaging capabilities. BEEM employs Scanning Tunneling Microscopy (STM) and a unique spatially localized ballistic electron spectroscopy method...

Bell, L. D.↗
DOE OSTI2023

Solution-Phase Synthesis and Photoluminescence of Quaternary Chalcohalide Semiconductors

Mixed-metal chalcohalide semiconductors have emerged as promising candidates for photovoltaic applications. However, preparation of these multinary compounds using solution-phase techniques remains particularly challenging compared to traditional solid-state methods. To fully harness their potential, it is desirable to develop synthetic methods that enable control over both the phase purity and dimensionality of chalcohalides. In this work, we report the solution-phase synthesis of Pb 2 SbS 2 I 3 and Pb 2 BiS 2 I 3 quaternary chalcohalides using readily available precursors. Fine tuning of reaction parameters allows for the isolation of rod-like morphologies with tunable diameters and aspect ratios. The quaternary chalcohalides display photoluminescence as an ensemble as well as at the single particle level, as demonstrated using fluorescence microscopy. We further evaluate the relative stability and band gap of Pb 2 SbS 2 I 3 polymorphs and their coloring patterns using electronic structure calculations. The synthetic methods developed here will motivate the study of ever more complex chalcohalides and other multinary semiconductors for new technological applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Probing atomic-scale symmetry breaking by rotationally invariant machine learning of multidimensional electron scattering

Abstract The 4D scanning transmission electron microscopy (STEM) method maps the structure and functionality of solids on the atomic scale, yielding information-rich data sets describing the interatomic electric and magnetic fields, structural and electronic order parameters, and other symmetry breaking distortions. A critical bottleneck is the dearth of analytical tools that can reduce complex 4D-STEM data to physically relevant descriptors. We propose an approach for the systematic exploration of 4D-STEM data using rotationally invariant variational autoencoders (rrVAE), which disentangle the general rotation of the object from other latent representations. The implementation of purely rotational rrVAE is discussed as are applications to simulated data for graphene and zincblende structures. The rrVAE analysis of experimental 4D-STEM data of defects in graphene is illustrated and compared to the classical center-of-mass analysis. This approach is universal for probing symmetry-breaking phenomena in complex systems and can be implemented for a broad range of diffraction methods.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Analysis of diagonal G and subspace W approximations within fully self-consistent GW calculations for bulk semiconducting systems

Fully self-consistent GW (sc-GW) methods are now available to evaluate quasiparticle and spectral properties of various molecular and bulk systems. However, such techniques based on the full matrix of G and W are computationally demanding. Additionally, the routinely used single-shot GW approximation (G 0 W 0 ) has an undesirable dependency on the choice of initial exchange-correlation functional. In the literature, many so-called self-consistent GW methods are based on diagonal approximation of G and low-ranking approximation of W. It is thus worth checking how good such approximations are in comparison with the full matrix method. In this work, we consider AlAs, AlP, GaP, and ZnS as the prototype systems to perform sc-GW calculations by expressing the full G matrix using a plane-wave basis set. We compared our sc-GW results with the diagonal G and subspace W approximated sc-GW results (sc-GW-diagG and sc-GW-subW methods). In the sc-GW-diagG method, interacting G is expanded in the eigenvectors of noninteracting G such that only diagonal elements are retained, whereas the number of eigenmodes is truncated in sc-GW-subW calculations. A systematic analysis of the results obtained from the above techniques is presented. The differences in the quasiparticle band gap between the approximated and the full matrix sc-GW approaches are mostly less than 1.7%, which validates such widely adopted approximations, and also shows how such low-ranking approximation can be used to include higher-order terms such as the vertex correction.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Microstructurally Strained Pyrochlore–Perovskite Biphasic Electrocatalysts for the Oxygen Evolution Reaction

Efficiency of water splitting for hydrogen production is often limited by the sluggish kinetics of multiple electronic transfers required in the heterogeneous oxygen evolution reaction (OER). Catalyst design for reducing the high OER overpotential remains a major scientific challenge. Lattice-strain engineering, a method for tuning the electronic structure and surface geometric configuration of active sites, may greatly affect the interaction between adsorbates and catalytic surfaces for high activity and stability. Here, in this study, we present the synthesis of biphasic oxides of YPrSrRuMnO x , which consists of distinct phases of Y 2 Ru 2 O 7 pyrochlore and (Pr 0.7 Sr 0.3 )MnO 3 perovskite, and the development of a suitable analytical approach to study the strain–catalytic property relationship. Linear sweep voltammetry results reveal that the biphasic oxide exhibits approximately 3.1 times greater mass activity and 2.4 times larger turnover frequency (TOF) than single-phase Y 2 Ru 2 O 7 in the 0.1 M HClO 4 electrolyte. The biphasic catalyst is also about 3 times more stable than the single-phase oxide under acidic conditions. X-ray photoelectron spectroscopy, nitrogen isotherm, and electrochemical surface area analyses indicate that the oxidation state, specific surface area, and electrochemical surface area do not cause enough difference in the observed enhancement of OER performance. We examined the effects of microstrain on electrocatalysis, originating from lattice mismatch between different phases, using three different structural models. Specifically, we compared the Williamson–Hall method, standard stress–strain analysis, and Rietveld refinement in analyzing the structure–property relationship. Strain mapping using geometric phase analysis (GPA) further revealed significant microstrain and lattice dislocations localized near phase boundaries in the biphasic oxide, in contrast to the uniform strain in single-phase materials. The results reveal that the increased microstrain correlates well with the improved OER performance, as the biphasic oxide catalyst exhibits 2–3 times greater microstrain than Y 2 Ru 2 O 7 pyrochlore.

electrocatalysts↗
DOE OSTI2022

Electronic structure of α - RuCl 3 by fixed-node and fixed-phase diffusion Monte Carlo methods

Layered material α-RuCl 3 has caught wide attention due to its possible realization of Kitaev's spin liquid and its electronic structure that involves the interplay of electron-electron correlations and spin-orbit effects. Several DFT+U studies have suggested that both electron-electron correlations and spin-orbit effects are crucial for accurately describing the band gap. This work studies the importance of these two effects using fixed-node and fixed-phase diffusion Monte Carlo calculations both in spin-averaged and explicit spin-orbit formalisms. In the latter, the Slater-Jastrow trial function is constructed from two-component spin orbitals using our recent quantum Monte Carlo (QMC) developments and thoroughly tested effective core potentials. Our results show that the gap in the ideal crystal is already accurately described by the spin-averaged case, with the dominant role being played by the magnetic ground state with significant exchange and electron correlation effects. We find qualitative agreement between hybrid DFT, DFT+U, and QMC. In addition, QMC results agree very well with available experiments, and we identify the values of exact Fock exchange mixing that provide comparable gaps. Explicit spin-orbit QMC calculations reveal that the effect of spin-orbit coupling on the gap is minor, of the order of 0.2 eV, which corresponds to the strength of the spin orbit of the Ru atom.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Constrained nuclear–electronic orbital method for periodic density functional theory: Application to H 2 chemisorption on Si(001) surfaces

The nuclear–electronic orbital (NEO) method provides a powerful computational framework for incorporating nuclear quantum effects (NQE) in electronic structure calculations beyond the Born–Oppenheimer approximation. By incorporating additional constraints to the position operator on quantum particles like protons, the NEO method enables calculation of effective potential that accounts for NQE. Here, in this work, we present a new constrained NEO (cNEO) formulation for density functional theory (cNEO-DFT) calculations in the context of extended periodic systems. Using the nudged elastic band method, we discuss an application of the cNEO-DFT approach to studying the adsorption of a hydrogen molecule on the Si(001) surfaces. The calculation shows how NQE impacts the reaction energetics. The proton density changes are computed along the reaction pathways. This work demonstrates the capability of the new cNEO-DFT method to study a wide range of chemical processes, such as surface reactions where the quantum nature of light atoms like protons is non-negligible.

Chemical processes↗
DOE OSTI2023

Twisted bilayer graphene as topological heavy fermion: II. Analytical approximations of the model parameters

The recently-introduced topological heavy fermion (THF) model [1] of twisted bilayer graphene (TBG) aims to reconcile the quantum-dot-like electronic structure of the latter observed by scanning tunneling microscopy, with its electron delocalization seen in transport measurements. The THF model achieves this by coupling localized (heavy) fermions with anomalous conduction electrons. Originally, the parameters of the THF model were obtained numerically from the Bistritzer–Macdonald (BM) model of TBG [1]. In this work, we derive analytical expressions for the THF model parameters as a function of the twist angle, the ratio between the tunneling amplitudes at the AA and AB regions (w 0 /w 1 ), and the screening length of the interaction potential. By numerically computing the THF model parameters across an extensive experimentally-relevant parameter space, we show that the resulting approximations are remarkably good, i.e., within the 30% relative error for almost the entire parameter space. At the single-particle level, the THF model accurately captures the energy spectrum of the BM model over a large phase space of angles and tunneling amplitude ratios. When interactions are included, we also show that the THF description of TBG is good around the magic angle for realistic values of the tunneling amplitude ratios (0.6 ≤ w 0 /w 1 ≤ 1.0), for which the hybridization between the localized and conduction fermions γ is smaller than the onsite repulsion of the heavy fermions U 1 (i.e., |γ| < U 1 ).

Approximation methods↗