Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “electronic structure theory”

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 271 records · Page 15

Relaxation time approximations in PAOFLOW 2.0

Abstract Regardless of its success, the constant relaxation time approximation has limited validity. Temperature and energy dependent effects are important to match experimental trends even in simple situations. We present the implementation of relaxation time approximation models in the calculation of Boltzmann transport in 2.0 and apply those to model band-structures. In addition, using a self-consistent fitting of the model parameters to experimental conductivity data, we provide a flexible tool to extract scattering rates with high accuracy. We illustrate the approximations using simple models and then apply the method to GaAs, Si, $$\hbox {Mg}_3\hbox {Sb}_2$$ Mg 3 Sb 2 , and $$\hbox {CoSb}_3$$ CoSb 3 .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Construction of meta-GGA functionals through restoration of exact constraint adherence to regularized SCAN functionals

The strongly constrained and appropriately normed (SCAN) meta-GGA exchange–correlation functional is constructed as a chemical environment-determined interpolation between two separate energy densities: one describes single-orbital electron densities accurately and another describes slowly varying densities accurately. To conserve constraints known for the exact exchange–correlation functional, the derivatives of this interpolation vanish in the slowly varying limit. While theoretically convenient, this choice introduces numerical challenges that degrade the functional’s efficiency. We have recently reported a modification to the SCAN meta-GGA, termed restored-regularized-SCAN (r 2 SCAN) , that introduces two regularizations into SCAN, which improve its numerical performance at the expense of not recovering the fourth order term of the slowly varying density gradient expansion for exchange. Here, we show the derivation of a progression of density functional approximations [regularized SCAN (rSCAN), r++SCAN, r 2 SCAN, and r 4 SCAN] with increasing adherence to exact conditions while maintaining a smooth interpolation. The greater smoothness of r 2 SCAN seems to lead to better general accuracy than the additional exact constraint of SCAN or r 4 SCAN does.

, correlation energy↗

Phonons of metallic hydrogen with quantum Monte Carlo

Here, we describe a simple scheme to perform phonon calculations with quantum Monte Carlo (QMC) methods and demonstrate it on metallic hydrogen. Because of the energy and length scales of metallic hydrogen and the statistical noise inherent to QMC methods, the conventional manner of calculating force constants is prohibitively expensive. We show that our alternate approach is nearly 100 times more efficient in resolving the force constants needed to calculate the phonon spectrum in the harmonic approximation. This requires only the calculation of atomic forces, as in the conventional approach, and otherwise little or no programmatic modification.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Taming the virtual space for incremental full configuration interaction

Incremental full configuration interaction (iFCI) closely approximates the FCI limit with polynomial cost through a many-body expansion of the correlation energy, providing highly accurate total energies within a given basis set. To extend iFCI beyond previous basis set limitations, this work introduces a novel natural orbital (NO) screening approach, incremental NO full configuration interaction (iNO-FCI). By consideration of the importance of virtual orbital selection in the convergence of iFCI, iNO-FCI maximizes the consistency between orbitals selected for each correlated body. iNO-FCI employs a principle of cancellation of errors and ensures that the same set of virtual NOs is used for interdependent terms. Here, this strategy significantly reduces computational cost without compromising precision. Computational savings of up to 95% are demonstrated, allowing access to larger basis sets that were previously computationally prohibitive. iNO-FCI is herein introduced and benchmarked for several difficult test cases involving double-bond dissociation, biradical systems, conjugated π systems, and the spin gap of a Cu-based transition metal complex.

Correlation energy↗

A chemical dynamics study on the gas-phase formation of triplet and singlet C 5 H 2 carbenes

Significance Carbenes represent key reactive intermediates in molecular mass growth processes leading to carbonaceous nanostructures in the interstellar medium and in combustion systems. However, due to their short lifetimes and tendency for dimerization, carbenes represent one of the foremost obscured classes of reactive intermediates, with the preparation of carbenes in the gas phase remaining largely elusive. By merging molecular beams with electronic structure and quasi-classical trajectory calculations and exploiting triplet pentadiynylidene (HCCCCCH) and singlet ethynylcyclopropenylidene (c-C 5 H 2 ) carbene as benchmarks, we present a versatile protocol to unravel the dynamics leading to exotic carbenes under single-collision conditions. These elementary mechanisms are of significance to our understanding of the fundamental processes leading to unsaturated organic transients in extreme, hydrocarbon-rich environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

First-principles computations of the Stark shift of a defect-bound exciton: The case of the T center in silicon

The T center in silicon has recently drawn a lot of attention for its potential in quantum information science. The sensitivity of the zero-phonon line (ZPL) to electrical field was recently investigated by a combination of different experimental methods but there are still few first principles studies on the Stark shift of the T center. Dealing with the defect-bound exciton nature of the excited state is particularly challenging using density functional theory because of the large spatial delocalization associated with the wavefunction. Here, in this work, we tackle this issue by performing a convergence study over the supercell size. We obtain an exciton binding energy of 28.5 meV, in good agreement with experimental results. We then calculate the Stark shift through the dipole moment change of the ZPL transition of the T center using the modern theory of polarization formalism and find a modest linear coefficient of Δ⁢𝜇=0.79⁢D along X and Δ⁢𝜇=0.03⁢𝐷 along Y. We discuss our results in light of the recent experimental measurements of the Stark shift. Our analysis suggests that bound-exciton defects could be particularly sensitive to local field effect as a result of their large spatial extent.

color centers↗

Pauli potential formalism at finite temperature

At zero temperature, the Pauli potential—the functional derivative of the Pauli kinetic energy density functional—is the key to the accuracy of the orbital-free density functional theory (OFDFT) as it is supposed to capture all the effects associated with the Pauli exclusion principle. Here, we extend this concept to finite temperature by defining Pauli free energy and the modified Pauli free energy, both representing the natural generalizations of the Pauli term from zero- T to finite- T . We discuss their physical interpretation, the mathematical nuances, and the applicability, arguing that the modified Pauli potential should be used as an extension of the zero- T counterpart within the OF-DFT framework. Through analytical and numerical methods, we then analyze some of the exact properties concerning the modified Pauli free- and kinetic-energy terms and examine the temperature dependence of the modified Pauli potential.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Nuclear cusps and singularities in the nonadditive kinetic potential bifunctional from analytical inversion

The non-additive kinetic potential $v^{\text{NAD}}$ is a key quantity in density-functional theory (DFT) embedding methods, such as frozen density embedding theory and partition DFT. $v^{\text{NAD}}$ is a bi-functional of electron densities $\rho_{\rm B}$ and $\rho_{\rm tot} = \rho_{\rm A} + \rho_{\rm B}$. It can be evaluated using approximate kinetic-energy functionals, but accurate approximations are challenging. The behavior of $v^{\text{NAD}}$ in the vicinity of the nuclei has long been questioned, and singularities were seen in some approximate calculations. In this article, the existence of singularities in $v^{\text{NAD}}$ is analyzed analytically for various choices of $\rho_{\rm B}$ and $\rho_{\rm tot}$, using the nuclear cusp conditions for the density and Kohn-Sham potential. It is shown that no singularities arise from smoothly partitioned ground-state Kohn-Sham densities. We confirm this result by numerical calculations on diatomic test systems HeHe, HeLi$^+$, and H$_2$, using analytical inversion to obtain a numerically exact $v^{\rm NAD}$ for the local density approximation. We examine features of $v^{\rm NAD}$ which can be used for development and testing of approximations to $v^{\rm NAD}[\rho_{\rm B},\rho_{\rm tot}]$ and kinetic-energy functionals.

74 ATOMIC AND MOLECULAR PHYSICS↗

K 2 CoS 2 : A two-dimensional in-plane antiferromagnetic insulator

Recent discovery of two-dimensional (2D) magnetic materials has brought magnetism to the flatland and opened up exciting opportunities for the exploration of fundamental physics as well as novel device applications. Here, we predict a thermodynamically stable 2D magnetic material, ${\mathrm{K}}_{2}{\mathrm{CoS}}_{2}$, which retains its in-plane bulk antiferromagnetic (AFM) order down to the monolayer and bilayer limits. Magnetic moments $(2.5{{\mu}}_{B}/\mathrm{Co})$ are found to form a quasi-one-dimensional antiferromagnetically ordered chain of Co atoms. The nonmagnetic electronic spectrum of the monolayer film is found to host flatbands and van Hove singularities, which play a key role in stabilizing the magnetic ground state. Based on classical Monte Carlo simulations, we estimate the Neel temperature for the AFM monolayer to be ${\approx}15\mathrm{K}$. Our study demonstrates that ${\mathrm{K}}_{2}{\mathrm{CoS}}_{2}$ hosts a robust AFM state which persists from the monolayer limit to the bulk material.

2-dimensional systems↗

High-throughput search for magnetic topological materials using spin-orbit spillage, machine learning, and experiments

Magnetic topological insulators and semi-metals have a variety of properties that make them attractive for applications including spintronics and quantum computation. Here, we use systematic high-throughput density functional theory calculations to identify magnetic topological materials from the ≈ 40000 three-dimensional materials in the JARVIS-DFT database. First, we screen materials with net magnetic moment > 0.5 μB and spin-orbit spillage > 0.25, resulting in 25 insulating and 564 metallic candidates. The spillage acts as a signature of spin-orbit induced band-inversion. Then, we carry out calculations of Wannier charge centers, Chern numbers, anomalous Hall conductivities, surface bandstructures, and Fermi-surfaces to determine interesting topological characteristics of the screened compounds. We also train machine learning models for predicting the spillage, bandgaps, and magnetic moments of new compounds, to further accelerate the screening process. We experimentally synthesize and characterize a few candidate materials to support our theoretical predictions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Oxygen magnetic polarization, nodes in spin density, and zigzag spin order in oxides

Recent studies on Ba 2 CoO 4 (BCO) and SrRuO 3 (SRO) have unveiled a variety of intriguing phenomena, such as magnetic polarization on oxygens, unexpected nodes in the spin density profile along bonds, and unusual zigzag spin patterns in triangular lattices. Here, using simple model calculations supplemented by density functional theory we explain the presence of nodes based on the antibonding character of the dominant singly occupied molecular orbitals along the transition metal (TM) to oxygen bonds. Our simple model also allows us to explain the net polarization on oxygen as originated from the hybridization between atoms and mobility of the electrons with spins opposite to those of the closest TM atoms. Overall, our results are not limited to BCO and SRO, but they are generic and qualitatively predict the net polarization expected on any ligands, according to the spin order of the closest TM atoms and the number of intermediate ligand atoms. Finally, we propose that a robust easy-axis anisotropy would suppress the competing 120° antiferromagnetic order to stabilize the zigzag pattern order as ground state in a triangular lattice. Our generic predictions should be applicable to any other compound with characteristics similar to those of BCO and SRO.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Chemical control of the Rashba spin splitting size of α-GeTe(111) surface states by adjusting the potential at the topmost atomic layer

As control of the Rashba spin splitting size is highly desirable for spintronic devices, intensive studies have been performed to vary the splitting size by, for example, applying an electric field or designing novel heterostructures. However, direct observation of Rashba spin splitting size change via spectroscopic measurements has not been done so far. Here, we report results of angle-resolved photoemission studies on ferroelectric α-GeTe(111). We observe that the Rashba splitting size of α-GeTe(111) surface states is reduced upon dosing with potassium (K) which has a very low electron affinity. Based on density functional theory calculations, we find that the electric potential energy at the topmost atomic layer and the surface potential energy barrier is reduced upon K dosing. Furthermore, this change in the surface potential induces both delocalization of the surface states and reduction of the splitting size. We expect to increase the splitting size by dosing other elements or molecules with high electron affinity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

First-principles theory for Schottky barrier physics

We develop a first-principles theory for Schottky barrier physics. The Poisson equation is solved selfconsistently with the electrostatic charge density over the entire barrier using the density functional theory (DFT) electronic structure converged locally, allowing computation of a Schottky barrier entirely from DFT involving thousands of atomic layers in the semiconductor (SC). Here, the induced charge in the bulk consists of conduction and valence band charges from doping and band bending, as well as charge from the evanescent states in the gap of the SC. The Schottky barrier height (SBH) is determined when the induced charge density and the induced electrostatic potential reach self-consistency. Tests on the GaAs-graphene and Si/Al heterostructures yield SBH, width, along with depletion and inversion layers obtained self-consistently as functions of temperature and bulk doping.

36 MATERIALS SCIENCE↗

Fast and robust all-electron density functional theory calculations in solids using orthogonalized enriched finite elements

Here, we present a computationally efficient approach to perform systematically convergent real-space all-electron Kohn-Sham density functional theory calculations for solids using an enriched finite element (FE) basis. The enriched FE basis is constructed by augmenting the classical FE basis with atom-centered numerical basis functions, comprising of atomic solutions to the Kohn-Sham problem. Notably, to improve the conditioning, we orthogonalize the enrichment functions with respect to the classical FE basis, without sacrificing the locality of the resultant basis. In addition to improved conditioning, this orthogonalization procedure also renders the overlap matrix block diagonal, greatly simplifying its inversion. Subsequently, we use a Chebyshev polynomial based filtering technique to efficiently compute the occupied eigenspace in each self-consistent field iteration. We demonstrate the accuracy and efficiency of the proposed approach on periodic unit cells and supercells. The benchmark studies show a staggering 130× speedup of the orthogonalized enriched FE basis over the classical FE basis. We also present a comparison of the orthogonalized enriched FE basis with the linearized augmented plane-wave + local orbitals basis, both in terms of accuracy and efficiency. Notably, we demonstrate that the orthogonalized enriched FE basis can handle large system sizes of ~10 000 electrons. Finally, we observe good parallel scalability of our implementation with 92% efficiency at 22× speedup for a system with 620 electrons

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Sparse expansions of multicomponent oxide configuration energy using coherency and redundancy

We report that compressed sensing has become a widely accepted paradigm to construct high dimensional cluster expansion models used for statistical mechanical studies of atomic configuration in complex multicomponent crystalline materials. However, strict sampling requirements necessary to obtain minimal coherence measurements for compressed sensing to guarantee accurate estimation of model parameters are difficult and in some cases impossible to satisfy due to the inability of physical systems to access certain configurations. Nevertheless, the dependence of energy on atomic configuration can still be adequately learned without these strict requirements by using compressed sensing by way of coherent measurements using redundant function sets known as frames. We develop a particular frame constructed from the union of all occupancy-based cluster expansion basis sets. We illustrate how using this highly redundant frame yields sparse expansions of the configuration energy of complex oxide materials that are competitive and often surpass the prediction accuracy and sparsity of models obtained from standard cluster expansions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Electronic excitations of α - Fe 2 O 3 heteroepitaxial films measured by resonant inelastic x-ray scattering at the Fe L edge

Resonant inelastic x-ray scattering (RIXS) spectra of hematite (α-Fe 2 O 3 ) were measured at the Fe L 3 edge for heteroepitaxial thin films which were undoped and doped with 1% Ti, Sn, or Zn, in the energy-loss range in excess of 1 eV to study electronic transitions. The spectra were measured for several momentum transfers q, conducted at both low temperature (T=14 K) and room temperature. While we cannot rule out dispersive features possibly owing to propagating excitations, the coarse envelopes of the general spectra did not appreciably change shape with q, implying that the bulk of the observed L-edge RIXS intensity originates from (mostly) nondispersive ligand field excitations. Summing the RIXS spectra over q and comparing the results at T=14 K to those at T=300 K revealed pronounced temperature effects, including an intensity change and energy shift of the ≈1.4 eV peak, a broadband intensity increase of the 3-4 eV range, and higher energy features. The q-summed spectra and their temperature dependencies are virtually identical for nearly all of the samples with different dopants, save for the temperature dependence of the Ti-doped sample's spectrum, which we attribute to being affected by a large number of free charge carriers. Comparing with magnetization measurements for different temperatures and dopings likewise did not show a clear correlation between the RIXS spectra and the magnetic ordering states. To clarify the excited states, we performed spin multiplet calculations which were in excellent agreement with the RIXS spectra over a wide energy range and provide detailed electronic descriptions of the excited states. Furthermore, the implications of these findings to the photoconversion efficiency of hematite photoanodes is discussed.

36 MATERIALS SCIENCE↗

Spin-polarized imaging of strongly interacting fermions in the ferrimagnetic state of the Weyl candidate CeBi

CeBi has an intricate magnetic phase diagram whose fully polarized state has recently been suggested as a Weyl semimetal, though the role of f states in promoting strong interactions has remained elusive. Here we focus on the less-studied but also time-reversal symmetry-breaking ferrimagnetic phase of CeBi, where our density functional theory (DFT) calculations predict additional Weyl nodes near the Fermi level E F . We use spin-polarized scanning tunneling microscopy and spectroscopy to image the surface ferrimagnetic order on the itinerant Bi p states, indicating their orbital hybridization with localized Ce f states. We observe suppression of this spin-polarized signature at E F , coincident with a Fano line shape in the conductance spectra, suggesting the Bi p states partially Kondo screen the f magnetic moments, and this p – f hybridization causes strong Fermi-level band renormalization. Furthermore, the p-band flattening is supported by our quasiparticle interference measurements, which also show band splitting in agreement with DFT, painting a consistent picture of a strongly interacting magnetic Weyl semimetal.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗