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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Methane Activation on Metal-Doped (111) and (100) Ceria Surfaces with Charge-Compensating Oxygen Vacancies

Reducing the temperature for methane activation is an important objective that would benefit many technological applications. In this work we explore the possibility to achieve this goal using single atom catalysts (SACs) on ceria surfaces. In particular, we focus on Ag SACs, which have recently been suggested to be promising catalysts for both the oxidation of molecular hydrogen and CH 4 activation. The rate limiting step of methane activation on ceria is the breaking of the first C-H bond. Using first principles calculations, we investigate this reaction on CeO 2 (111) and (100) surfaces, which are both frequently exposed on ceria nanoparticles. For undoped ceria, our results show that the dissociation of CH 4 to CH 3 + H has a lower barrier on the (001) surface than on the (111) one. On both surfaces, the presence of Ag enhances the binding of CH 4 to the surface and reduces the activation energy for its dissociation. Interestingly, the formation of methanol via the Mars-van Krevelen mechanism has a slightly lower energy barrier than the dissociation of CH 4 to CH 3 + H on Ag-doped CeO 2 (111), suggesting that methanol is a likely product of methane activation on this surface. Here, a novel aspect of our work relative to previous studies on CH 4 activation over transition metal doped ceria is that we mainly focus on stable surface structures where each Ag dopant substituting a Ce atom forms a complex with a surface oxygen vacancy (VO). Our results show that these Ag-VO complexes play a critical role and accounting for their presence is essential for a proper description of the surface reactivity.

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

RbB 3 Si 3 : An Alkali Metal Borosilicide that is Metastable and Superconducting at 1 atm

The stability, electronic structure and potential superconductivity in AB 3 Si 3 (A = Na, K, Rb, and Cs) compounds that assume the clathrate-based sodalite structure whose frameworks consists of covalent B-Si bonds is investigated via first-principles calculations. This structure type has recently been predicted in a number of high temperature superconducting hydrides, but these are only stable under megabar pressures. Herein, we predict a novel superconducting phase, RbB 3 Si 3 , that could be synthesized under pressures that are a factor of ten smaller, ~10 GPa, and quenched to atmospheric conditions. Electron-phonon coupling calculations predict that RbB 3 Si 3 possesses a superconducting critical temperature, T c , of 14 K at 1 atm. Here, the dynamic stability of RbB 3 Si 3 and CsB 3 Si 3 at ambient pressure can be explained by considering the chemical pressure exerted on the B-Si framework that is caused by the size effect of the alkali metal atom.

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

Intercalation of Lithium into Graphite: Insights from First-Principles Simulations

Understanding ion intercalation at electrode–electrolyte interfaces is key to the development of energy storage and water desalination. Here, we investigate Li + kinetics at a prototypical interface between graphite anodes and an organic electrolyte, and we elucidate key factors that determine ion transport, using first-principles methodology coupling ab initio molecular dynamics simulations with a solvation model. We show that surface chemical composition significantly influences the kinetics of ion intercalation from the liquid into graphite. We find that this is partly related to the ion desolvation process, which varies notably for different graphite surface chemical terminations. In addition, interfacial polarization is found to play an important role in determining energy barriers for ion transfer. We also discuss the impact of electrode potentials, which is often neglected in conventional first-principles calculations despite being a key factor in device configurations. Our study provides insights into the coupling of electronic and ionic effects of interfacial chemistry on ion transport at complex electrode–electrolyte interfaces.

25 ENERGY STORAGE↗
DOE OSTI2020

Electronic Structure of Double-Layer Epitaxial Graphene on SiC(0001) Modified by Gd Intercalation

In this work, we systematically study the effects of Gd adsorption and intercalation on the electronic band structure of double-layer epitaxial graphene on Si-terminated SiC(0001) by first-principles calculations. We show that Gd adsorption and intercalation exhibit strong effects on the coupling between the graphene layers and between the buffer layer and substrate. Different adsorption/intercalation geometries can result in very different electron band structures. The number of Dirac cones and the positions of the Dirac cones relative to the Fermi level can be effectively manipulated through controlling the Gd adsorption/intercalation geometries. Our calculations provide useful insights to guide the experimental design of graphene-based materials with desirable functionalities for applications.

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

Mechanism of Metal Intercalation under Graphene through Small Vacancy Defects

Metal intercalation under graphene has attracted extensive experimental and theoretical research because of its capability to manipulate the electronic structure and properties of graphene. However, the pathways and mechanisms of intercalation are still not well understood. Here, we systematically investigate the intercalation process of metal atoms through graphene vacancies using first-principles calculations. We show that the energy barrier for metal atom penetration through the vacancies in graphene is small as long as the size of the vacancy is larger than a mono-vacancy. However, metal atoms are strongly bonded to the vacancy so that the detachment energy of a metal atom from the vacancy is extremely high. This inhibits the diffusion of the metal atom into the gallery beneath the surface to complete the intercalation process. On the other hand, our calculation results show that the detachment energy of a metal atom from a metal dimer at small vacancy defects is significantly reduced, making intercalation much easier. Therefore, the key step limiting the intercalation process is the detachment of the metal atoms from vacancy defects. This finding from our study provides useful insight into the defect-assisted intercalation mechanism.

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

Correlation between Spin and Orbital Dynamics during Laser-Induced Femtosecond Demagnetization

Spin and orbital angular momenta are two intrinsic properties of an electron and are responsible for the physics of a solid. How the spin and orbital evolve with respect to each other on several hundred femtoseconds is largely unknown, but it is at the center of laser-induced ultrafast demagnetization. In this paper, we introduce a concept of the spin–orbital correlation diagram, where spin angular momentum is plotted against orbital angular momentum, much like the position-velocity phase diagram in classical mechanics. We use four sets of highly accurate time-resolved X-ray magnetic circular dichroism data to construct four correlation diagrams for iron and cobalt. To our surprise, a pattern emerges. The trace on the correlation diagram for iron is an arc, and at the end of demagnetization, it has a pronounced cusp. The correlation diagram for cobalt is different and appears more linear but with kinks. We carry out first-principles calculations with two different methods: time-dependent density functional theory (TDDFT) and time-dependent Liouville density functional theory. These two methods agree that the experimental findings for both Fe and Co are not due to experimental errors. It is the spin–orbit coupling that correlates the spin dynamics to the orbital dynamics. Microscopically, Fe and Co have different orbital occupations, which leads to distinctive correlation diagrams. We believe that this correlation diagram presents a useful tool to better understand spin and orbital dynamics on an ultrafast time scale. A brief discussion on the magnetic anisotropy energy is also provided.

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

Computational Investigation of the Catalytic Hydrodeoxygenation of Propanoic Acid over a Cu(111) Surface

Cu-based alloy catalysts have recently been investigated experimentally for the hydrodeoxygenation (HDO) of biomass-derived organic acids. Here, the HDO of propanoic acid (PAc) has been studied over Cu(111) by mean-field microkinetic modeling based on parameters obtained from first-principles calculations. Models were developed for the gas- and liquid-phase HDO in condensed water and 1,4-dioxane. In agreement with experimental observations, the gas-phase PAc conversion rate is low at 573 K and increases in liquid water by 1 order of magnitude. In all reaction environments, the decarboxylation mechanism is dominant at low hydrogen partial pressures less than 0.1 bar, and the C–COO bond dissociation is the rate-controlling elementary step. This observation contrasts with the rate-controlling step identified over most group VIII metal surfaces, which is the C–OH bond dissociation in the decarbonylation mechanism. At high hydrogen (H 2 ) partial pressures greater than 10 bar, the HDO of PAc produces propionaldehyde that can readsorb and further react through decarbonylation to produce C 2 alkane products, which is conceptually different from the low H 2 partial pressure scenario. At high H 2 partial pressures, the initial hydrogenation at the carbonyl carbon of PAc becomes the rate-controlling elementary step.

Catalysts↗
DOE OSTI2021

Quantifying Temperature Dependence of Electronic Band Gaps and Optical Properties in SnO 2 and SnO via First-Principles Simulations

Tin metal oxides SnO x ( x = 1, 2) have gained interest as gas-sensing materials. For their applications as high-temperature sensors, however, a better understanding of their temperature dependence sensing responses is needed. Here, we comparatively quantify the temperature-dependent electronic band gaps and optical properties of SnO 2 and SnO using first-principles calculations. Without considering the temperature effect, SnO 2 and SnO are predicted to have direct and indirect band gaps of 2.18 and 1.66 eV, respectively, at the PBE + U -GGA level that we employed. The temperature effect on the electronic and optical properties is captured by taking account of the electron–phonon interaction. Band gap renormalization with temperature is calculated via the Allen–Heine–Cardona theory. For both oxides, we find a monotonic decrease in the electronic band gap such that renormalization at zero point (0 K) is ~–0.17 and ~–0.52 eV at 1000 K. These results are also analyzed by employing an analytical equation that helps characterize the band gap shift with temperature. In addition, the optical properties at finite temperatures are simulated using the frozen-phonon method that combines electron–phonon coupling with the momentum matrix. As temperature increases, the optical property spectra are smoothed because of the smearing effect, which diminishes optical constants at shorter wavelengths. Our results are of interest for high-temperature functional materials in applications of optical detection.

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

Straintronic Effect on Phonon-Mediated Superconductivity of Nb 2 CT 2 (T = O, S, Se, or Te) MXenes

Here, the electronic structures, phonon dispersions, and electron–phonon coupling of Nb 2 CT 2 (T = O, S, Se, or Te) MXenes were investigated via first-principles calculations. Different models of Nb 2 CT 2 were constructed, and the results show that the low-energy models of Nb 2 CT 2 are intrinsic phonon-mediated superconductors. Of the four Nb 2 CT 2 MXenes, Nb 2 CO 2 MXene exhibits the largest superconducting critical temperature ($T_c$) of 14.43 K. The existence of soft modes induced by Kohn anomalies and the contribution of Nb atoms to the Fermi level lead to strong electron–phonon coupling (λ = 0.92) in Nb 2 CO 2 MXene. The $T_c$ of Nb 2 CO 2 is further enhanced by biaxial tensile strain and reaches up to 18.28 K under 4% tensile strain. The predicted $T_c$ of Nb 2 CS 2 is 4.5 K, which is comparable with experimental data. These findings will further stimulate the search for superconducting MXenes.

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

Effect of Multiple Oxygen Vacancies on the Optical and Thermodynamic Properties of La 0.75 Sr 0.25 Co 0.25 Fe 0.75 O 3-δ Perovskite

We report the oxygen vacancy concentration of perovskite changes with the work condition in applications such as oxygen carriers and oxygen sensors. To explore the effect of multiple oxygen vacancies on the optical and thermodynamic properties, we perform first-principles calculations to investigate the formation of multiple oxygen vacancies in the paramagnetic La 0.75 Sr 0.25 Co 0.25 Fe 0.75 O 3-δ . Our results show that the higher VO concentration in La 0.75 Sr 0.25 Co 0.25 Fe 0.75 O 3-δ occurs at the conditions of higher temperatures and lower oxygen partial pressures. The electrons liberated upon oxygen vacancy formation are partially localized on the nearby Fe/Co sites and partially delocalized through the lattice, with the result that the unoccupied Fe/Co states near the Fermi level decrease and the Fermi level shows an upward shift. From the results of optical absorption, the peaks from 250 to 600 nm are weakened and blue-shifted with increase of oxygen vacancy concentration, due to the decrease of empty Fe states and an upward shift of the Fermi level. The obvious change of optical absorption with oxygen vacancy concentration indicates that La 0.75 Sr 0.25 Co 0.25 Fe 0.75 O 3-δ could be a good candidate for optical oxygen sensors in the visible spectrum.

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

Unraveling 1-Hexene Hydrogenation over Dilute Pd-in-Au Alloys

Dilute Pd-in-Au alloys are valuable materials for selectively hydrogenating alkynes and isomerizing alkenes. By diluting Pd in a Au host, the selectivity toward semihydrogenated alkene isomers can be significantly enhanced and the unfavorable overhydrogenation to alkanes is suppressed. However, a detailed mechanistic study on the origin of the enhanced alkene selectivity over dilute alloy catalysts is still missing. Here, we combine experiment and theory to unravel the reaction mechanism, identifying rate-limiting and selectivity-controlling steps in 1-hexene hydrogenation over dilute Pd-in-Au catalysts. Using isotope-exchange hydrogenation experiments, we show that 1-hexene and hydrogen over a bimetallic Pd4Au96 in silica catalyst preferentially form 1-hexene isomers, (trans and cis) 2- and 3-hexene and only a small amounts of hexane. The reaction is consistent with a Horiuti–Polanyi mechanism, similar to a monometallic Pd nanoparticle catalyst. Computation of the free-energy profiles for 1-hexene hydrogenation and isomerization over a single Pd atom in a Au surface using first principles calculations indicated that the isomerization of 1-hexene to 2-hexene is energetically favorable due to the relatively large barrier for H2 dissociation preventing hydrogenation to n-hexane. Microkinetic modeling established that H2 dissociation on the single-atom Pd sites and H spillover from these sites onto the Au host are rate-limiting and key in steering the selectivity of dilute Pd-in-Au alloys toward the hexene isomers. Furthermore, the mechanistic insights from this study contribute to the rational design of optimized dilute alloy catalysts for selective alkene isomerization.

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

Computational Screening of Cationic Zeolites for n -Butane/Methane Separations Using Quantitatively Accurate First-Principles-Derived Force Fields

We developed force fields for linear alkanes in Na and Ca-exchanged zeolites based on periodic first-principles calculations with coupled cluster corrections. These force fields were validated by comparing the simulated adsorption properties with extensive experimental data. As an example of using these force fields, we screened silica and aluminosilicate zeolites for separation of n-butane/methane mixtures in several different separation processes. A key step in this screening was the use of an interpolation method to create accurate zeolite structures that take account of the lattice constant changes for different Si/Al ratios for a given zeolite topology. We focused on zeolite structures that have known synthesis routes, so the materials selected from our calculations can be tested experimentally. Several promising materials that show good separation performance were selected for more detailed simulations.

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

THD-C Sheet: A Novel Nonbenzenoid Carbon Allotrope with Tetra-, Hexa-, and Dodeca-Membered Rings

Here, we propose a novel two-dimensional carbon-based structure with tetra-, hexa-, and dodeca-membered rings, which we refer to by the abbreviated name, THD-C. The structure presents a mixture of sp–sp 2 hybridization and can potentially be synthesized by the topological assembly of 4-ethynyldiphenylacetylene molecules. By employing first-principles calculations, the stability and ease of synthesis of the sheet are investigated and compared with various C-allotropes. We predict its metallic behavior and excellent kinetic and dynamic stability. Due to the crystal structure of the sheet, a strong mechanical anisotropy is observed. The effects of functionalization on the electronic properties of the material are also studied, and different semiconducting systems are obtained. The potential of THD-C for energy storage in metal-based batteries, hydrogen storage, and catalysis is also investigated, and we find a superior performance in comparison to graphite and other allotropes. The quantum confinement effect is investigated by constructing nanoribbons and nanotubes of various sizes. For ribbons, we find that tailor-made electronic and magnetic properties can be obtained and explored in potential spintronic devices. Additionally, we observe that nanotubes are conducting irrespective of their chirality and can potentially be used for capture, storage, and separation of industrially relevant small gas molecules.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Theoretical Investigations of Anharmonic Effects and Phonon Transport in the Cubic Phase of Crystalline Perovskite CsPbCl 3

The role of anharmonic effects on lattice dynamics and thermal transport has been investigated in the cubic phase of the CsPbCl 3 perovskite. Limitations of the harmonic approximation, which lead to phonon instabilities in the Brillouin zone, were addressed based on self-consistent phonon theory in combination with first-principles calculations and by incorporation of frequency renormalization effects from bubble and loop diagrams. This theoretical approach demonstrates significant improvements compared to results obtained using a self-consistent phonon dynamic matrix. Both the cubic-to-tetragonal phase transition temperature as well as the predicted lattice thermal conductivity values show good agreement with selected sets of available experimental data. In conclusion, the accuracy of the predicted thermal conductivity data is also tested against systems significantly larger than those allowed by quantum calculations, by using molecular dynamic simulations with machine-learning force fields.

Lattices↗
DOE OSTI2025

Stabilization of Preternatural Barium Oxidation States as an Unexpected Byproduct of β-Decay: Discovery of a New Halide Semiconductor Alloy

137 Cs has a wide range of roles in the nuclear industry. The solid material, safely encapsulated in CsCl as 137 CsCl, is stored as fission product waste from nuclear power production and legacy waste from nuclear weapons production; it has also served as a radiation source in food and sewage irradiators as well as medical devices. However, because of the solubility of the chloride salt and the relatively high specific activity of 137 Cs, damaged or broken capsules can lead to severe radiological accidents. Safe capsule design and material recycling are complicated by the unclear structural evolution during β-decay, which remains ambiguous due to the differing oxidation states of Cs (1+) and Ba (2+). Here, in this study, we use first-principles calculations to investigate the evolving structure–property relationships of Cs 1–x Ba x Cl during β-decay. Despite the well-established 2+ formal oxidation state of alkali-earth metals, we find that Ba 1+ can be stabilized in the form of a mixed-valence alloy at low concentrations. Specifically, we identify three regimes for the β-decay of 137 Cs into CsCl: Ba-doped CsCl (Ba ≤ 14%), wherein Ba has the expected 2+ oxidation state; Cs–Ba–Cl alloys, where Ba has a mix of the usual Ba 2+ and highly unusual Ba 1+ oxidation state in the form of a quasi-disordered mixed-valence alloy (Ba = 25%); and phase separation into a CsCl + BaCl 2 + Ba (m) mechanical mixture, where Ba reverts to its expected 2+ oxidation state (Ba > 25%). Surprisingly, the Cs 0.75 Ba 0.25 Cl mixed-valence alloy is a narrow indirect band gap semiconductor (1.05 eV) despite the insulating nature of both CsCl and BaCl 2 . It also exhibits strongly excitonic polarized optical properties, has glass-like ultralow thermal conductivity (directional average of 0.21 W/mK at 300 K), and shows greater resistance to deformation under both tensile and volumetric strengths compared with the original CsCl structure (e.g., shear and Young’s modulus of 9.04 and 31.62 GPa, respectively). These findings imply that transmutation of 137 Cs leads to highly unusual chemical bonding that stabilizes Ba 1+ in local regions of the quasi-disordered Cs 0.25 Ba 0.75 Cl, resulting in anomalous physical properties. Moreover, this discovery provides valuable insight for safe nuclear waste capsule design, which can aid in preventing environmental or human exposure to radioactive materials.

Fuhr, Addis S. [Oak Ridge National Laboratory (ORN↗
DOE OSTI2020

Ion Pairing and Redissociaton in Low-Permittivity Electrolytes for Multivalent Battery Applications

Detailed speciation of electrolytes as a function of chemical system and concentration provides the foundation for understanding bulk transport as well as possible decomposition mechanisms. In particular, multivalent electrolytes have shown a strong coupling between anodic stability and solvation structure. Furthermore, solvents that are found to exhibit reasonable stability against alkaline-earth metals generally exhibit low permittivity, which typically increases the complexity of the electrolyte species. To improve our understanding of ionic population and associated transport in these important classes of electrolytes, the speciation of Mg(TFSI)2 in monoglyme and diglyme systems is studied via a multiscale thermodynamic model using first-principles calculations for ion association and molecular dynamics simulations for dielectric properties. The results are then compared to Raman and dielectric relaxation spectroscopies, which independently confirm the modeling insights. We find that the significant presence of free ions in the low-permittivity glymes in the concentration range from 0.02 to 0.6 M is well-explained by the low-permittivity redissociation hypothesis. Here, salt speciation is largely dictated by long-range electrostatics, which includes permittivity increases due to polar contact ion pairs. The present results suggest that other low-permittivity multivalent electrolytes may also reach high conductivities as a result of redissociation.

25 ENERGY STORAGE↗
DOE OSTI2020

Inverse Design of Ultralow Lattice Thermal Conductivity Materials via Materials Database Screening of Lone Pair Cation Coordination Environment

The presence of lone pair (LP) electrons is strongly associated with the disruption of lattice heat transport, which is a critical component of strategies to achieve efficient thermoelectric energy conversion. By exploiting an empirical relationship between lattice thermal conductivity, κ L , and the bond angles of pnictogen group LP cation coordination environments, we develop an inverse design strategy based on a materials database screening to identify chalcogenide materials with ultralow κ L for thermoelectrics. Screening the ~635 000 real and hypothetical inorganic crystals of the Open Quantum Materials Database based on the constituent elements, nominal electron counting, LP cation coordination environment, and synthesizability, we identify 189 compounds expected to exhibit ultralow κ L . As a validation, we explicitly compute the lattice dynamical properties of two of the compounds (Cu 2 AgBiPbS 4 and MnTl 2 As 2 S 5 ) using first-principles calculations and successfully find both achieve ultralow κ L values at room temperature of ~0.3-0.4 W/(m·K) corresponding to the amorphous limit. Furthermore, our data-driven approach provides promising candidates for thermoelectric materials and opens new avenues for the design of phononic properties of materials.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Work Function Engineering of 2D Materials: The Role of Polar Edge Reconstructions

2D materials have attracted tremendous interest as functional materials because of their diverse and tunable properties, especially at their edges. A material’s work function is a critical parameter in many electronic devices; however, a fundamental understanding and a path toward large alterations of the work function in 2D materials still remain elusive. Here in this paper, we report the first evidence for anisotropy of the work function in 2D MoS 2 from first-principles calculations. We also demonstrate large work-function tunability (in the range of 3.45–6.29 eV) choosing the 2H phase of MoS 2 as a model system by sampling various edge configurations. We furthermore reveal the origin of this work function anisotropy and tunability by extending the existing work function relation to the local dipole moment at surfaces of 3D materials to those at edges in 2D materials. We then use machine-learning approaches to correlate work function with edge structures. These results pave the way for intrinsic edge engineering for electronic and catalytic applications.

36 MATERIALS SCIENCE↗