Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “first principles 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,315 records · Page 73
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Ab initio study on fcc Pr with correlation matrix renormalization theory

Here, we studied fcc praseodymium (Pr) with the ab initio correlation matrix renormalization theory (CMRT) explicitly calculating the Coulomb interactions among basis orbitals and without using any adjustable parameters to work for strongly correlated electron systems. We calculated its total energy in a paramagnetic ground state and studied the role of the correlated 4f electrons in the system. Good agreement was obtained between CMRT and experiments in the pressure volume dependence of the fcc phase. We also compared the CMRT results against other theoretical methods including local density approximation + dynamical mean-field theory and showed consistent results among them. Moreover, we found the normalized local charge fluctuation of the 4f electrons as the function of volume exhibits a clear slope change at the volume collapse region, indicating a switch in their correlation nature there.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Electronic correlation induced expansion of Fermi pockets in δ-plutonium

Plutonium is a critically important material as the behavior of its 5f electrons stands midway between the metalliclike itinerant character of the light actinides and localized atomic-core-like character of the heavy actinides. The δ phase of plutonium (δ-Pu), whereas still itinerant, has a large coherent Kondo peak and strong electronic correlations coming from its near-localized character. Here, using sophisticated Gutzwiller wave function and dynamical mean-field theory correlated theories, we study the Fermi surface and associated mass renormalizations of δ-Pu together with calculations of the de Haas–van Alphen frequencies. We find a large (~200%) correlation induced volume expansion in both the hole and the electron pockets of the Fermi surface in addition to an intermediate mass enhancement. All of the correlated electron theories predict, approximately, the same hole pocket placement in the Brillouin zone, which is different from that obtained in conventional density-functional band-structure theory, whereas the electron pockets from all theories are in, roughly, the same place.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Electronic structure of van der Waals ferromagnet CrI 3 from self-consistent vertex corrected GW approaches

Electronic structure of the layered van der Waals ferromagnet CrI 3 is studied with self consistent diagrammatic approaches beyond the GW approximation. Considerable improvement in the calculated band gap as compared to the non-self-consistent G0W0 results has been found. Certain spectral features in the valence bands discovered recently by the angle resolved photoemission spectroscopy, are reproduced better when we use full frequency dependent self energy. Density functional theory and the quasiparticle self-consistent GW method which are based on the frequency-independent self energy are unable to resolve these features. The non-locality effects in the diagrams beyond the GW approximation are large for both polarizability and self energy. This finding can potentially have an impact on the development of methods like the GW+DMFT (dynamical mean field theory).

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Hybridization effect on the x-ray absorption spectra for actinide materials: Application to PuB 4

Studying the local moment and 5 f -electron occupations sheds insight into the electronic behavior in actinide materials. X-ray absorption spectroscopy (XAS) has been a powerful tool to reveal the valence electronic structure when assisted with theoretical calculations. However, the analysis currently taken in the community on the branching ratio of the XAS spectra generally does not account for the hybridization effects between local f orbitals and conduction states. In this paper, we discuss an approach which employs the density functional theory plus Gutzwiller rotationally invariant slave boson method to obtain a local Hamiltonian for the single-impurity Anderson model, and calculates the XAS spectra by the exact diagonalization (ED) method. A customized numerical routine was implemented for the ED XAS part of the calculation. By applying this technique to the recently discovered 5 f -electron topological Kondo insulator Pu B 4 , we determined the signature of 5 f -electronic correlation effects in the theoretical x-ray spectra. Furthermore, we found that the Pu 5 f - 6 d hybridization effect provides an extra channel to mix the j = 5 / 2 and 7 / 2 orbitals in the 5 f valence. As a consequence, the resulting electron occupation number and spin-orbit coupling strength deviate from the intermediate-coupling regime.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Robust antiferromagnetism in Y 2 Co 3

We report on a solution-growth based method to synthesise single crystals of Y 2 Co 3 and on its structural and magnetic properties. We find that Y 2 Co 3 crystallizes in the La 2 Ni 3 -type orthorhombic structure with space group Cmce (No. 64), with Co forming distorted kagome lattices. Y 2 Co 3 orders antiferromagnetically below TN=252 K. Magnetization measurements reveal that the moments are primarily aligned along the b axis with evidence for some canting. Band-structure calculations indicate that ferromagnetic and antiferromagnetic orders are nearly degenerate, at odds with experimental results. Magnetization measurements under pressure up to 1 GPa reveal that the Néel temperature decreases with the slope of -1.69 K/GPa. We observe a field-induced spin-flop transition in the magnetization measurements at 1.5 K and 21 T with magnetic field along the b direction. The magnetization is not saturated up to 35 T, indicating that the antiferromagnetic ordering in Y 2 Co 3 is quite robust, which is surprising for such a Co-rich intermetallic.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Electronic states dressed by an out-of-plane supermodulation in the quasi-two-dimensional kagome superconductor CsV 3 Sb 5

CsV 3 Sb 5 has attracted much recent attention as the first quasi-two-dimensional (2D) kagome superconductor. While the kagome layers are 2D in nature, increasing evidence has pointed to the importance of out-of-plane correlation in this material. However, it remains unclear whether such correlation can change the fundamental electronic structure of the quasi-2D system. Here, we reveal this missing piece of information, using angle-resolved photoemission spectroscopy, complemented by scanning tunneling microscope measurements. The three-dimensional electronic structures in the high-temperature state are revealed, which agree well with density-functional theory calculations. Electron energy bands are observed in the low-temperature state that exhibit additional periodicities along the out-of-plane momentum. Furthermore, these results reveal a direct response to the out-of-plane electronic supermodulation in the single-particle spectral function of CsV 3 Sb 5 , thus establishing an electronic platform to examine emergent phenomena beyond 2D limit in kagome superconductors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Elimination of the linearization error in APW/LAPW basis set: Dirac-Kohn-Sham equations

Here, a detailed account of the implementation of equations of the relativistic density functional theory (RDFT) using basis sets of APW/LAPW type with flexible extensions provided by local orbitals is given. Earlier discoveries of the importance of the high derivative local orbital (HDLO) extension of the APW/LAPW basis set for enhancing the accuracy of DFT calculations are confirmed using a fully relativistic approach and α – U as an example. High energy local orbitals (HELO's), however indispensable for GW calculations, are considerably less efficient in enhancing the accuracy of DFT applications. It is shown that a simplified approach to the relativistic effects, namely considering them only inside the muffin-tin (MT) spheres, produces basically identical results (as compared to fully relativistic approach) for the electronic free energy of the five materials considered in this work. By comparing the effect of the simplified approach on the electronic free energy with its effect on the electronic kinetic energy we conclude that the insensitivity of the free energy to the way we describe the relativistic effects in the interstitial region is related to the variational property of this quantity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Suppression of antiferromagnetic order and strong ferromagnetic spin fluctuations in Ca(Co 1-x Ni x ) 2-y As 2 single crystals

CaCo 2–y As 2 is a unique itinerant system having strong magnetic frustration. Here, we report the effect of electron doping on the physical properties resulting from Ni substitutions for Co. The single crystals of Ca(Co 1–x Ni x ) 2–y As 2 were characterized by single-crystal x-ray diffraction, energy-dispersive x-ray spectroscopy, magnetization M versus temperature T, magnetic field H, time t, and heat capacity C p (H, T) measurements. The A-type antiferromagnetic (AFM) transition temperature T N = 52 K for x = 0 decreases to 22 K with only 3% Ni substitution and is completely suppressed for x > 0.16. For 0.11 ≤ x ≤ 0.52 strong ferromagnetic (FM) fluctuations develop as revealed by magnetic susceptibility χ(T) = M(T) / H measurements. For x = 0.11 and 0.16 competing AFM and FM interactions result in a reentrant spin-glass behavior below T N , as evidenced by the observations of thermomagnetic hysteresis and magnetic relaxation. Enhanced FM fluctuations are also found for the x = 0.21 and 0.31 crystals, where χ c increases significantly at low T. A large χ anisotropy in these compositions where χ c is up to a factor of two larger than χ ab suggests that the FM spin fluctuations are quasi-1D in nature. Weak ferromagnetic contributions to the magnetization are found at T = 2 K for x = 0.11 –0.31. Heat-capacity Cp(T) measurements reveal the presence of FM quantum spin fluctuations for 0.11 ≤ x ≤ 0.52, where a logarithmic T dependence of Cp(T) / T is observed at low T. Here, the suppression of AFM order by the development of strong FM fluctuations in Ca(Co 1–x Ni x ) 2–y As 2 crystals suggests the presence of a FM quantum-critical point at x ≈ 0.20. Our density-functional theory (DFT) calculations confirm that FM fluctuations are enhanced by Ni substitutions for Co in CaCo 2–y As 2 . The Sommerfeld electronic heat-capacity coefficient is enhanced for x = 0, 0.21, and 0.42 by about a factor of two compared to DFT calculations of the density of states (DOS) at the Fermi energy, suggesting an enhancement of the DOS from electron-phonon and/or electron-electron interactions. The crystals with x > 0.52 do not exhibit FM spin fluctuations or magnetic order at T ≥ 1.8 K, which was found from the DFT calculations to arise from a Stoner transition. Superconductivity is not observed above 1.8 K for any of the compositions. Neutron-diffraction studies of crystals with x = 0.11 and 0.16 in the crossover regime (0.1 ≲ x ≲ 0.2) show no evidence of A-type ordering as observed in the parent compound with x = 0. Furthermore, no other common magnetic structures, such as ferromagnetic (FM), helical stacking of in-plane FM layers, or in-plane AFM structure, are found with an ordered moment greater than the uncertainty of 0.05 μ B per transition-metal atom.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Prediction of topological phases in metastable ferromagnetic $MPX_3$ monolayers

Density functional theory calculations are carried out to study the electronic and topological properties of MPX 3 (M = Mn, Fe, Co,Ni; X = S, Se) monolayers in the ferromagnetic (FM) metastable magnetic state. Here we find that FM MnPSe 3 monolayers host topological semimetal signatures that are gapped out when spin-orbit coupling is included. These findings are supported by explicit calculations of the Berry curvature and the Chern number. The choice of the Hubbard-U parameter to describe the d-electrons is thoroughly discussed, as well as the influence of using a hybrid-functional approach. The presence of band inversions and the associated topological features are found to be formalism-dependent. Nevertheless, routes to achieve the topological phase via the application of external biaxial strain are demonstrated. Within the hybrid-functional picture, topological band structures are recovered under a pressure of 15% (17 GPa). The present work provides a potential avenue for uncovering new topological phases in metastable ferromagnetic phases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Improved quasiparticle self-consistent electronic band structure and excitons in β – LiGaO 2

The band structure of β–LiGaO 2 is calculated using the quasiparticle self-consistent QSGWˆ method where the screened Coulomb interaction Wˆ is evaluated including electron-hole interaction ladder diagrams and G is the one-electron Green's function. Improved convergence compared to previous calculations leads to a significantly larger band gap of about 7.0 eV. However, exciton binding energies are found to be large and lead to an exciton gap of about 6.0 eV if also a zero-point-motion correction of about –0.4 eV is included. Furthermore, these results are in excellent agreement with recent experimental results on the onset of absorption. Besides the excitons observed thus far, the calculations indicate the existence of a Rydberg-like series of exciton excited states, which is however modified from the classical Wannier exciton model by the anisotropies of the material and the more complex mixing of Bloch states in the excitons resulting from the Bethe-Salpeter equation. The exciton fine structure and the exciton wave functions are visualized and analyzed in various ways.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Atomically Dispersed Single Ni Site Catalysts for Nitrogen Reduction toward Electrochemical Ammonia Synthesis Using N 2 and H 2 O

Ammonia (NH 3 ) electrosynthesis gains significant attention as NH 3 is essentially important for fertilizer production and fuel utilization. However, electrochemical nitrogen reduction reaction (NRR) remains a great challenge because of low activity and poor selectivity. Herein, a new class of atomically dispersed Ni site electrocatalyst is reported, which exhibits the optimal NH 3 yield of 115 µg cm -2 h -1 at –0.8 V versus reversible hydrogen electrode (RHE) under neutral conditions. High faradic efficiency of 21 ± 1.9% is achieved at -0.2 V versus RHE under alkaline conditions, although the ammonia yield is lower. The Ni sites are stabilized with nitrogen, which is verified by advanced X-ray absorption spectroscopy and electron microscopy. Density functional theory calculations provide insightful understanding on the possible structure of active sites, relevant reaction pathways, and confirm that the Ni-N 3 sites are responsible for the experimentally observed activity and selectivity. Extensive controls strongly suggest that the atomically dispersed NiN 3 site-rich catalyst provides more intrinsically active sites than those in N-doped carbon, instead of possible environmental contamination. Overall, this work further indicates that single-metal site catalysts with optimal nitrogen coordination is very promising for NRR and indeed improves the scaling relationship of transition metals.

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

Melting curves of atomic hydrogen and deuterium calculated using path-integral Monte Carlo

We calculate the melting line of atomic hydrogen and deuterium up to 900 GPa with path-integral Monte Carlo using a machine-learned interatomic potential. We improve upon previous simulations of melting by treating the electrons with reptation quantum Monte Carlo, and by performing solid and liquid simulations using isothermal-isobaric path-integral Monte Carlo. Here, the resulting melting line for atomic hydrogen is higher than previous estimates. There is a small but resolvable decrease in the melting temperature as pressure is increased, which can be attributed to quantum effects.

08 HYDROGEN↗
DOE OSTI2023

Time-dependent density-functional-theory calculations of the nonlocal electron stopping range for inertial confinement fusion applications

Nonlocal electron transport is important for understanding laser-target coupling for laser-direct-drive (LDD) inertial confinement fusion (ICF) simulations. Current models for the nonlocal electron mean free path in radiation-hydrodynamic codes are based on plasma-physics models developed decades ago; improvements are needed to accurately predict the electron conduction in LDD simulations of ICF target implosions. Here we utilized time-dependent density functional theory (TD-DFT) to calculate the electron stopping power (SP) in the so-called conduction-zone plasmas of polystyrene in a wide range of densities and temperatures relevant to LDD. Compared with the modified Lee-More model, the TD-DFT calculations indicated a lower SP and a higher stopping range for nonlocal electrons. We fit these electron SP calculations to obtain a global analytical model for the electron stopping range as a function of plasma conditions and the nonlocal electron kinetic energy. This model was implemented in the one-dimensional radiation-hydrodynamic code LILAC to perform simulations of LDD ICF implosions, which are further compared with simulations by the standard modified Lee-More model. In conclusion, results from these integrated simulations are discussed in terms of the implications of this TD-DFT-based mean-free-path model to ICF simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Hole Trapping at Acceptor Impurities and Alloying Elements in AlN

Hole localization is a root cause of difficulties in p-type doping of aluminum nitride. Herein, the stability of localized holes in AlN and their susceptibility to self-trapping within bulk material or in the presence of isovalent elements and acceptor impurities are calculated. It is found that self-trapped holes are metastable in bulk AlN, and also exhibit a very low barrier to detrapping. However, holes become trapped in the presence of acceptor impurities as well as isovalent elements such as B, In, and Sc, and may contribute to the broad luminescence observed in AlN alloys. It is also found that hole trapping contributes to the large ionization energies for acceptors in AlN. Most acceptors can also trap a second hole, becoming positively charged and further exacerbating difficulties in doping aluminum nitride p type. Furthermore, the stability of trapped holes is found to scale with atomic size mismatch to Al.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Correlation strength and orbital differentiation across the phase diagram of plutonium metal

We compare the trends on the strength of electronic correlations across the different phases of elemental Pu focusing on its site and orbital dependence, using a combination of density functional theory (DFT) and dynamical mean-field theory (DMFT) calculations within the vertex corrected one crossing approximation. We find that Pu-5f states are more correlated in δ -Pu, followed by some crystallographic sites in α and β phases. In addition, we observe that Pu-5f 5/2 and Pu-5f 7/2 orbital differentiation is a general feature of this material, as is site differentiation in the low-symmetry phases. The Pu-5f 5/2 states show Fermi liquid like behavior, whereas the Pu-5f 7/2 states remaining incoherent down to very low temperatures. We correlate the correlation strength in the different phases to their structure and the Pu-5 f occupancy of their crystallographic sites.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Hybrid interatomic potential for Sn

To design materials for extreme applications, it is important to understand and predict phase transitions and their influence on material properties under high pressures and temperatures. Atomistic modeling can be a useful tool to assess these behaviors. However, this can be difficult due to the lack of fidelity of the interatomic potentials in reproducing this high pressure and temperature extreme behavior. Here, in this work, a hybrid EAM-R—which is the combination of embedded atom method (EAM) and rapid artificial neural network potential—for Tin (Sn) is described which is capable of accurately modeling the complex sequence of phase transitions between different metallic polymorphs as a function of pressure. This hybrid approach ensures that a basic empirical potential like EAM is used as a lower energy bound. By using the final activation function, the neural network contribution to energy must be positive, assuring stability over the whole configuration space. This implementation has the capacity to reproduce density functional theory results at 6 orders of magnitude slower than a pair potential for molecular dynamics simulation, including elastic and plastic characteristics and relative energies of each phase. Using calculations of the Gibbs free energy, it is demonstrated that the potential precisely predicts the experimentally observed phase changes at temperatures and pressures across the whole phase diagram. At 10.2 GPa, the present potential predicts a first-order phase transition between body-centered tetragonal (BCT) β-Sn and another polymorph of BCT-Sn. This structure transforms into body-centered cubic near the experimentally reported value at 33 GPa. Thus, the Sn potential developed in this paper can be used to study complex deformation mechanisms under extreme conditions of high pressure and strain rates unlike existing potentials. Moreover, the framework developed in this paper can be extended for different material systems with complex phase diagrams.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Single pair of Weyl nodes in the spin-canted structure of EuCd 2 As 2

We report time reversal symmetry breaking Weyl semimetals are unique among Weyl materials in allowing the minimal number of Weyl points, thus offering the clearest signatures of the associated physics. Here we present neutron diffraction, density-functional theory, and transport measurement results which indicate that EuCd 2 As 2 , under ambient field, strain, and pressure, is such a material with a single pair of Weyl points. Our work reveals a magnetic structure (magnetic space group C2'/m') with Eu moments pointing along the [210] direction in plane and canted ~ 30° out of plane. Density functional theory calculations using this structure show that the observed canting drastically alters the relevant electronic bands, relative to the in-plane order, leading to a single set of well-defined Weyl points. Furthermore, we find the canting angle can tune the distance of the Weyl points above the Fermi level, with the smallest distance at low canting angles. Finally, transport measurements of the anomalous Hall effect and longitudinal magnetoresistance exhibit properties indicative of a chiral anomaly, thus supporting the neutron scattering and DFT results suggesting EuCd 2 As 2 is close to the ideal situation of the Weyl hydrogen atom.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Effect of localization on photoluminescence and zero-field splitting of silicon color centers

The study of defect centers in silicon has been recently reinvigorated by their potential applications in optical quantum information processing. A number of silicon defect centers emit single photons in the telecommunication O-band, making them promising building blocks for quantum networks between computing nodes. The two-carbon G-center, self-interstitial W-center, and spin-1/2 T-center are the most intensively studied silicon defect centers, yet despite this, there is no consensus on the precise configurations of defect atoms in these centers, and their electronic structures remain ambiguous. Here, in this work, we employ ab initio density functional theory to characterize these defect centers, providing insight into the relaxed structures, band structures, and photoluminescence spectra, which are compared to experimental results. Motivation is provided for how these properties are intimately related to the localization of electronic states in the defect centers. In particular, we present the calculation of the zero-field splitting for the excited triplet state of the G-center defect as the structure is linearly interpolated from the A-configuration to the B-configuration, showing a sudden increase in the magnitude of the D zz component of the zero-field-splitting tensor. By performing projections onto the local orbital states of the defect, we analyze this transition in terms of the symmetry and bonding character of the G-center defect, which sheds light on its potential application as a spin-photon interface.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗