Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density functional calculation”

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,639 records · Page 91
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

An automated cluster surface scanning method for exploring reaction paths on metal-cluster surfaces

Metal-cluster surfaces present a wide variety of unique coordination environments. This complexity makes it difficult to manually probe the surface reactivity of such clusters. Here, we present a simple and automated method to systematically discover reaction pathways on cluster surfaces, based on the automated cluster surface scanning (ACSS) technique for mapping out potential energy surfaces. We showcase our method on 55-atom icosahedral Cu and Ag clusters, where we determine the activation energies of four elementary steps common in heterogeneous catalysis – hydrogen recombination (H* + H* → H 2 * + *), oxygen recombination (O* + O* → O 2 * + *), water formation (OH* + H* → H 2 O* + *), and CO oxidation (CO* + O* → CO 2 * + *) – with density functional theory calculations (DFT-PBE + D3). We show that the ACSS method requires significantly less human effort than the established manually performed climbing-image nudged elastic band (MP + CI-NEB) technique and locates transition states with comparable accuracy (root-mean-squared error of 0.10 eV) and similar computational cost. Rigorous sampling of the potential energy surface with the ACSS method allows one to locate all lowest-energy reaction pathways obtained via the MP+CI-NEB approach, as well as alternative pathways that one may have missed with the MP+CI-NEB approach due to the many possible pathways available on these clusters. The accuracy and efficiency afforded by the ACSS method could enable high-throughput exploration of the diverse reactivity of metal clusters.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Ab initio property predictions of quinary solid solutions using small binary cells

The Set of Small Ordered Structures (SSOS) approach is an ab initio technique for modelling random solid solutions in which many small structures are averaged so that their correlation functions match those of a desired composition. SSOS has been shown to be effective in reducing the cost of density functional theory calculations relative to other well-known techniques such as cluster expansions and special quasirandom structures for modelling solid solutions. Here in this work, we demonstrate that SSOS’s can be constructed using cells with only a subset of elements while still accurately modelling multi-component systems. Specifically, we show that small binary cells can effectively model two quinary high entropy alloys – NbTaTiHfZr and MoNbTaVW – accurately capturing properties such as formation energy, lattice parameters, elastic constants, and root-mean-square atomic displacements. Overall, this insight is useful for those looking to construct databases of such small structures for predicting the properties of multi-component solid solutions, as it greatly decreases the number of structures that needs to be considered.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Quantum-Confined Lifshitz Transition on Weyl Semimetal T d -MoTe 2

Adsorption of alkali atoms onto material surfaces is widely utilized for controlling electronic properties and is particularly effective for two-dimensional materials. While tuning the chemical potential and band gap and creating quantum-confined states are well established for alkali adsorption on semiconductors, the effects on semimetallic systems remain largely elusive. Here, utilizing angle-resolved photoemission spectroscopy measurements and density functional theory calculations, we disclose the creation of two-dimensional electron gas and the quantum-confined Lifshitz transition at the surface of a Weyl semimetal T d -MoTe 2 by potassium adsorption. Electrons from potassium adatoms are shown to be transferred mainly to the lowest unoccupied band within the gapped part of the Brillouin zone, which, in turn, induces strong surface band bending and quantum confinement in the topmost layer. The quantum-confined topmost layer evolves from a semimetal to a strong metal with a Lifshitz transition departing substantially from the bulk band. The present finding and its underlying mechanism can be exploited for the creation of electronic heterojunctions in van der Waals semimetals.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Ferro and antiferro orbital ordering in Fe{sub 0.5}Mn{sub 0.5}V{sub 2}O{sub 4}

Using density functional theory calculations, we have investigated the orbital ordering in Fe{sub 0.5}Mn{sub 0.5}V{sub 2}O{sub 4} where Fe and V sites are orbitally active. Our first principles study within GGA+U and GGA+U+SO shows ferro-orbital ordering of d{sub x2−y2} orbital at all Fe sites, whereas A-type antiferro-orbital ordering at V sites where one 3d electron occupies d{sub xy} orbital at every V site and another electron occupies either 1/√2 (d{sub xz} + d{sub yz}) or 1/√2 (d{sub xz} - d{sub yz}) orbital alternatively along c axis. Insulating nature and the orbital ordering of this compound are found to be correlation driven while the effect of spin-orbit interaction on orbital ordering is not significant.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Data-Driven Prediction of Formation Mechanisms of Lithium Ethylene Monocarbonate with an Automated Reaction Network

Interfacial reactions are notoriously difficult to characterize, and robust prediction of the chemical evolution and associated functionality of the resulting surface film is one of the grand challenges of materials chemistry. The solid-electrolyte interphase (SEI), critical to Li-ion batteries (LIBs), exemplifies such a surface film, and despite decades of work, considerable controversy remains regarding the major components of the SEI as well as their formation mechanisms. Here we use a reaction network to investigate whether lithium ethylene monocarbonate (LEMC) or lithium ethylene dicarbonate (LEDC) is the major organic component of the LIB SEI. Our data-driven, automated methodology is based on a systematic generation of relevant species using a general fragmentation/recombination procedure which provides the basis for a vast thermodynamic reaction landscape, calculated with density functional theory. The shortest pathfinding algorithms are employed to explore the reaction landscape and obtain previously proposed formation mechanisms of LEMC as well as several new reaction pathways and intermediates. For example, we identify two novel LEMC formation mechanisms: one which involves LiH generation and another that involves breaking the (CH 2 )O-C(=O)OLi bond in LEDC. Most importantly, we find that all identified paths, which are also kinetically favorable under the explored conditions, require water as a reactant. This condition severely limits the amount of LEMC that can form, as compared with LEDC, a conclusion that has direct impact on the SEI formation in Li-ion energy storage systems. Finally, the data-driven framework presented here is generally applicable to any electrochemical system and expected to improve our understanding of surface passivation.

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

First-Principles Investigation of Adsorbate–Adsorbate Interactions on Ni(111), Ni(211), and Ni(100) Surfaces

Periodic density functional theory calculations were employed to investigate the coadsorption patterns between CO and key water gas shift reaction (WGSR) intermediates, i.e., H2O, H, OH, O, and COOH, on Ni(111), Ni(100), and Ni(211) single-crystal surfaces. It has been shown that although the nature of these adsorbate pair interactions are predominantly repulsive, as commonly assumed in the literature, the interaction pattern for each adsorbate pair can be complicated by facet lattice structures and potential intermolecular hydrogen bonds, which play an unsubtle role in influencing the WGSR redox and carboxyl pathways. This investigation can enrich our descriptions of the adsorbate lateral interactions for surface chemistry modeling and will also provide theoretical insights into tuning the properties of nanoscale nickel catalysts by manipulating the surface structure and morphology for hydrogen fuel production.

water-gas shift reaction↗
DOE OSTI2020

Beyond Volume Variation: Anisotropic and Protrusive Lithiation in Bismuth Nanowire

Materials storing energy via an alloying reaction are promising anode candidates in rechargeable lithium-ion batteries (LIBs) due to their much higher energy density than the current graphite anode. Until now, the volumetric expansion of such electrode particles during lithiation has been considered as solely responsible for cycling-induced structural failure. Here, we report different structural failure mechanisms using single-crystalline bismuth nanowires as the alloying-based anode. The Li–Bi alloying process exhibits a two-step transition, that is, Bi–Li 1 Bi and Li 1 Bi–Li 3 Bi. Interestingly, the Bi–Li 1 Bi phase transition occurs not only in the bulk Bi nanowire but also on the particle surface showing its characteristic behavior. The bulk alloying kinetics favors a Bi-(012)-facilitated anisotropic lithiation, whose mechanism and energetics are further studied using the density functional theory calculations. More importantly, the protrusion of Li 1 Bi nanograins as a result of anisotropic Li–Bi alloying is found to dominate the surface morphology of Bi particles. The growth kinetics of Li 1 Bi protrusions is understood atomically with the identification of two different controlling mechanisms, that is, the dislocation-assisted strain relaxation at the Bi/Li 1 Bi interface and the short-range migration of Bi supporting the off-Bi growth of Li 1 Bi. As loosely rooted to the bulk substrate and easily peeled off and detached into the electrolyte, these nanoscale protrusions developed during battery cycling are believed to be an important factor responsible for the capacity decay of such alloying-based anodes at the electrode level.

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

Growth of Single- and Bilayer ZnO on Au(111) and Interaction with Copper

The stoichiometric single- and bi-layer ZnO(0001) have been prepared by reactive deposition of Zn on Au(111) and studied in detail with X-ray photoelectron spectroscopy, scanning tunneling microscopy, and density functional theory calculations. Both single- and bi-layer ZnO(0001) adopt a planar, graphite-like structure similar to freestanding ZnO(0001) due to the weak van der Waals interactions dominating their adhesion with the Au(111) substrate. At higher temperature, the single-layer ZnO(0001) converts gradually to bi-layer ZnO(0001) due to the twice stronger interaction between two ZnO layers than the interfacial adhesion of ZnO with Au substrate. It is found that Cu atoms on the surface of bi-layer ZnO(0001) are mobile with a diffusion barrier of 0.31 eV, and likely to agglomerate and form nanosized particles at low coverages; while Cu atoms tend to penetrate a single layer of ZnO(0001) with a barrier of 0.10 eV, resulting in a Cu free surface.

36 MATERIALS SCIENCE↗
DOE OSTI2014

A first principles analysis of the effect of hydrogen concentration in hydrogenated amorphous silicon on the formation of strained Si-Si bonds and the optical and mobility gaps

In this paper, we use a model of hydrogenated amorphous silicon generated from molecular dynamics with density functional theory calculations to examine how the atomic geometry and the optical and mobility gaps are influenced by mild hydrogen oversaturation. The optical and mobility gaps show a volcano curve as the hydrogen content varies from undersaturation to mild oversaturation, with largest gaps obtained at the saturation hydrogen concentration. At the same time, mid-gap states associated with dangling bonds and strained Si-Si bonds disappear at saturation but reappear at mild oversaturation, which is consistent with the evolution of optical gap. The distribution of Si-Si bond distances provides the key to the change in electronic properties. In the undersaturation regime, the new electronic states in the gap arise from the presence of dangling bonds and strained Si-Si bonds, which are longer than the equilibrium Si-Si distance. Increasing hydrogen concentration up to saturation reduces the strained bonds and removes dangling bonds. In the case of mild oversaturation, the mid-gap states arise exclusively from an increase in the density of strained Si-Si bonds. Analysis of our structure shows that the extra hydrogen atoms form a bridge between neighbouring silicon atoms, thus increasing the Si-Si distance and increasing disorder in the sample.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Formic Acid Decomposition on Au catalysts: DFT, Microkinetic Modeling, and Reaction Kinetics Experiments

A combined theoretical and experimental approach is presented that uses a comprehensive mean-field microkinetic model, reaction kinetics experiments, and scanning transmission electron microscopy imaging to unravel the reaction mechanism and provide insights into the nature of active sites for formic acid (HCOOH) decomposition on Au/SiC catalysts. All input parameters for the microkinetic model are derived from periodic, self-consistent, generalized gradient approximation (GGA-PW91) density functional theory calculations on the Au(111), Au(100), and Au(211) surfaces and are subsequently adjusted to describe the experimental HCOOH decomposition rate and selectivity data. It is shown that the HCOOH decomposition follows the formate (HCOO) mediated path, with 100% selectivity toward the dehydrogenation products (CO21H2) under all reaction conditions. An analysis of the kinetic parameters suggests that an Au surface in which the coordination number of surface Au atoms is 4 may provide a better model for the active site of HCOOH decomposition on these specific supported Au catalysts.

active sites↗
DOE OSTI2025

Environmental controls on the kinetics of iron-sulfur cluster nucleation and nanoparticle formation

Anoxic, sulfidic conditions have been prevalent since the early Proterozoic and favor aqueous iron-sulfur (FeS aq ) clusters as a major fraction of the soluble, reduced iron and sulfur pool. FeS aq cluster formation and nucleation is driven by the high affinity between ferrous iron (Fe(II)) and sulfide (HS − ), ultimately yielding particles that precipitate as iron sulfide minerals. FeS aq clusters were recently shown to be bioavailable sources of iron and sulfur for a variety of anaerobes, yet little is known of the factors that influence the kinetics of their formation and nucleation. Here we apply computational and spectroscopic approaches to investigate the dynamics of FeS aq nucleation, cluster growth, precipitation, and redissolution as a function of Fe(II)/HS − concentration, temperature, and pH. Experiments were conducted under excess HS − to mimic euxinic conditions common to contemporary anaerobic aquatic ecosystems and those of the Proterozoic. Density functional theory calculations reveal the key role of water oxygen-iron interactions in stabilizing small FeS aq clusters and promoting solubility. Dynamic light scattering revealed a concentration-dependent increase in the kinetics of FeS aq nucleation and cluster aggregation. Increasing temperature promoted FeS aq cluster nucleation and aggregation while also enhancing dissolution. Alkaline pH also promoted FeS aq nucleation and cluster aggregation. At 25 °C, pH 7.0, and at reactant concentrations of 30 µM, FeS aq clusters < 10 nm in diameter remained in solution for > 2 h. These results underscore the importance of temperature, pH, and reactant concentration in the kinetics of FeS aq nucleation and cluster growth that, in turn, influence their bioavailability in anaerobic ecosystems.

Aquatic ecosystems↗
DOE OSTI2015

Thin film lubrication of hexadecane confined by iron and iron oxide surfaces: A crucial role of surface structure

A comparative analysis of thin film lubrication of hexadecane between different iron and its oxide surfaces has been carried out using classical molecular dynamic simulation. An ab initio force-field, COMPASS, was applied for n-hexadecane using explicit atom model. An effective potential derived from density functional theory calculation was utilized for the interfacial interaction between hexadecane and the tribo-surfaces. A quantitative surface parameterization was introduced to investigate the influence of surface properties on the structure, rheological properties, and tribological performance of the lubricant. The results show that although the wall-fluid attraction of hexadecane on pure iron surfaces is significantly stronger than its oxides, there is a considerable reduction of shear stress of confined n-hexadecane film between Fe(100) and Fe(110) surfaces compared with FeO(110), FeO(111), Fe{sub 2}O{sub 3}(001), and Fe{sub 2}O{sub 3}(012). It was found that, in thin film lubrication of hexadecane between smooth iron and iron oxide surfaces, the surface corrugation plays a role more important than the wall-fluid adhesion strength.

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

Coordination and activation of nitrous oxide by iron zeolites

Iron-containing zeolites are heterogeneous catalysts that exhibit remarkable activity in the selective oxidation of inert hydrocarbons and catalytic decomposition of nitrous oxide (N 2 O). The reduction of N 2 O is critical to both these functions, but experimental data tracking the iron active sites during N 2 O binding and activation are limited. Here, the N 2 O-ligated Fe(ii) active site in iron-exchanged zeolite beta is isolated and characterized by variable-temperature Mössbauer, diffuse reflectance UV-vis-NIR and Fourier transform infrared spectroscopy. N 2 O binds through the terminal nitrogen atom with substantial backbonding from the Fe(ii) centre at low temperature. At higher temperatures, the Fe–N 2 O interaction is weakened, facilitating isomerization to the O-bound form, which is competent in O-atom transfer. Density functional theory calculations show the geometric and electronic structure requirements for N 2 O binding and activation. A geometric distortion imposed by the zeolite lattice plays an important role in activating N 2 O. As a result, this highlights a mechanism for structural control over function in Fe-zeolite catalysts.

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

Gold-titania interface toughening and thermal conductance enhancement using an organophosphonate nanolayer

We demonstrate that a mercaptan-terminated organophosphonate nanolayer at gold-titania interfaces can give rise to two- to three-fold enhancement in the interfacial fracture toughness and thermal conductance. Electron spectroscopy reveals that interfacial delamination occurs at the metal-molecule interface near the gold-sulfur bonds, consistent with density functional theory calculations of bond energies. Qualitative correlation between interfacial fracture toughness and bond energies suggest that organophosphonate nanolayers are resilient to humidity-induced degradation. These results, and the versatility of organophosphonates as surface functionalization agents for technologically relevant materials, unlock uncharted avenues for molecular engineering of interfaces in materials and devices for a variety of applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Experimental and theoretical study of hole scattering in RF sputtered p-type Cu 2 O thin films

One of the key materials of interest for p-type oxide semiconductor thin film electronics is cuprous oxide (Cu 2 O), due to its relatively high hole mobility. Here, in this work, we use experiments, analytical models, and density functional theory calculations to study the scattering mechanisms that determine Hall mobility in two Cu 2 O samples. First, we examine a polycrystalline Cu 2 O thin film deposited by RF magnetron sputtering, and second, a single-crystalline Cu2O bulk substrate. Temperature-dependent Hall measurements indicate that neutral impurity and grain boundary scattering are dominant for the polycrystalline Cu 2 O thin film, while phonon scattering is dominant for single-crystalline Cu 2 O. Our first-principles calculations show that the room-temperature intrinsic hole mobility of Cu 2 O is 106 cm 2 V –1 s –1 , indicating the great promise of the material for p-type electronic devices. This intrinsic mobility is limited by phonon scattering, with the most dominant scattering modes having phonon energies of 88.4 and 17.1 meV. These results indicate that the key pathways to increase the hole mobility in Cu 2 O thin films are by reducing the impurity concentration and by increasing grain size. Our work thus sets the stage for the future development of high performance Cu 2 O-based p-type thin film transistors.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Single-Step Direct Laser Writing of Multimetal Oxygen Evolution Catalysts from Liquid Precursors

We investigate a laser direct-write method to synthesize and deposit metastable, mixed transition metal oxides and evaluate their performance as oxygen evolution reaction catalysts. This laser processing method enabled the rapid synthesis of diverse heterogeneous alloy and oxide catalysts directly from cost-effective solution precursors, including catalysts with a high density of nanocrystalline metal alloy inclusions within an amorphous oxide matrix. The nanoscale heterogenous structures of the synthesized catalysts were consistent with reactive force-field Monte Carlo calculations. By evaluating the impact of varying transition metal oxide composition ratios, we created a stable Fe 0.63 Co 0.19 Ni 0.18 O x /C catalyst with a Tafel slope of 38.23 mV dec -1 and overpotential of 247 mV, a performance similar to that of IrO 2 . Synthesized Fe 0.63 Co 0.19 Ni 0.18 O x /C and Fe 0.14 Co 0.46 Ni 0.40 O x /C catalysts were experimentally compared in terms of catalytic performance and structural characteristics to determine that higher iron content and a less crystalline structure in the secondary matrix decreases the charge transfer resistance and thus is beneficial for electrocatalytic activity. Furthermore, this conclusion is supported by density-functional theory calculations showing distorted active sites in ternary metal catalysts are key for lowering overpotentials for the oxygen evolution reaction.

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

Raman studies of hydrogen trapped in As 4 O 6 ·2H 2 at high pressure and low temperature

In this work, Raman spectroscopic measurements of the arsenolite–hydrogen inclusion compound As 4 O 6 ·2H 2 were performed in diamond anvil cells at high pressure and variable temperature down to 80 K. The experimental results were complemented by ab initio molecular dynamics simulations and phonon calculations. Observation of three hydrogen vibrons in As 4 O 6 ·2H 2 is reported in the entire temperature and pressure range studied (up to 24 GPa). While the experiments performed with protium and deuterium at variable temperatures allowed for the assignment of two vibrons as Q 1 (1) and Q 1 (0) transitions of ortho and para spin isomers of hydrogen trapped in the inclusion compound, the origin of the third vibron could not be unequivocally established. Low-temperature spectra revealed that the lowest-frequency vibron is actually composed of two overlapping bands of A g and T 2g symmetries dominated by H 2 stretching modes as predicted by our previous density functional theory calculations. We observed low-frequency modes of As 4 O 6 ·2H 2 vibrations dominated by H 2 “librations,” which were missed in a previous study. A low-temperature fine structure was observed for the J = 0 → 2 and J = 1 → 3 manifolds of hydrogen trapped in As 4 O 6 ·2H 2 , indicating the lifting of degeneracy due to an anisotropic environment. A non-spherical distribution was captured by molecular dynamics simulations, which revealed that the trajectory of H 2 molecules is skewed along the crystallographic $\langle$111$\rangle$ direction. Last but not least, low-temperature synchrotron powder x-ray diffraction measurements on As 4 O 6 ·2H 2 revealed that the bulk structure of the compound is preserved down to 5 K at 1.6 GPa.

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

Strain-modulated intercalated phases of Pb monolayer with dual periodicity in SiC(0001)-graphene interface

Intercalation of metal atoms at the SiC(0001)-graphene (Gr) interface can provide confined 2D metal layers with interesting electronic properties. The intercalated Pb monolayer (ML) has shown the coexistence of the Gr(10 x 10)-moiré and a stripe phase, which still lacks understanding. Using density functional theory calculation and thermal annealing with ab initio molecular dynamics as motivated by experiment, we have studied the formation energy of Gr/Pb/SiC(0001) for different Pb coverages. Near the coverage of a Pb(111)-like ML mimicking the (10 x 10)-moiré, we find a slightly more stable stripe structure, where one half of the structure has compressive strain with Pb occupying the Si-top sites and the other half has tensile strain with Pb off the Si-top sites. This stripe structure along the Gr zigzag direction has a periodicity of 2.3 nm across the [1$\overline{2}$10] direction agreeing with the previous observations using scanning tunneling microscopy. Analysis with electron density difference and density of states show the tensile region has a more metallic character than the compressive region, while both are dominated by charge transfer from Pb ML to SiC(0001). As a result, the small energy difference between the stripe and Pb(111)-like structures means the two phases are almost degenerate and can coexist, which explains the experimental observations.

36 MATERIALS SCIENCE↗