Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DFT”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Theoretical assessments of CO 2 activation and hydrogenation pathways on transition-metal surfaces

Carbon dioxide (CO 2 ) hydrogenation on transition-metal active sites offers a promising carbon utilization route toward mitigating greenhouse gas emissions. C 1 products are often formed in parallel during CO 2 hydrogenation, prompting investigations into the intrinsic properties of transition metals that drive activity and product selectivity. Here, in this work, close-packed surfaces of a selection of transition-metal catalysts (Ni, Co, Rh, Ru, Pd, and Pt) were studied with density functional theory (DFT) calculations to understand their fundamental reactivities for CO 2 transformation reactions. Results indicate that CO 2 conversion proceeds through CO* formation and hydrogenation to form C 1 products (* denotes an adsorbed species). Ni, Co, Rh, and Ru favor CO/CH 4 formation, while Pd and Pt favor CO/CH 3 OH formation. The ability of a metal to dissociate C-O bonds drives selectivity between CH 4 and CH 3 OH, while competition between CO* desorption and surface hydrogenation describes CO selectivities. The C-O bond dissociation steps often impose the highest barrier along CH 4 formation reaction profiles, suggesting their kinetic relevance for CH 4 formation rates. The provided DFT-derived data sets detail a comprehensive reaction network of elementary steps relevant to C 1 chemistries, ultimately offering a benchmark for insights into design strategies for materials that exploit transition-metal active sites in carbon capture or utilization processes.

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

Product-promoted acetylene cyclotrimerization to benzene on Ag(111)

Acetylene cyclotrimerization to benzene on Ag(111) has recently been found to be 100% selective. However, it requires monolayer coverages, limiting its practical application due to the high pressure required to attain this surface coverage of acetylene near room temperature. Here, we show that co-adsorbed benzene, the reaction product, promotes cyclotrimerization at submonolayer acetylene coverages without compromising the 100% selectivity of Ag(111) via temperature-programmed desorption experiments. Scanning tunneling microscopy and density functional theory (DFT) calculations elucidate the origin of this effect whereby attractive acetylene–benzene interactions locally compress acetylene domains. DFT calculations further illustrate how co-adsorbed benzene reduces the transition-state energy for formation of the key C 4 intermediate, making cyclotrimerization competitive with desorption. Together, these findings reveal how product-reactant interactions can enhance C–C coupling.

Acetylene cyclotrimerization↗
DOE OSTI2022

Thermodynamic modeling of the Pd–Zn system with uncertainty quantification and its implication to tailor catalysts

Pd–Zn intermetallic catalysts show encouraging combinations of activity and selectivity on well-defined active site ensembles. Thermodynamic description of the Pd–Zn system, delineating phase boundaries and enumerating site occupancies within intermediate alloy phases, is essential to determining the ensembles of Pd–Zn atoms as a function of composition and temperature. Combining the present extensive first-principles calculations based on density functional theory (DFT) and available experimental data, the Pd–Zn system was remodeled using the CALculation of PHAse Diagrams (CALPHAD) approach. High throughput modeling tools with uncertainty quantification, i.e., ESPEI and PyCalphad, were incorporated in the phase analysis. Here, the site occupancies across the γ phase composition region were given special attention. A four-sublattice model was used for the γ phase owing to its four Wyckoff positions, i.e., the outer tetrahedral (OT) site 8c, the inner tetrahedral (IT) site 8c, the octahedral (OH) site 12e, and the cuboctahedral (CO) site 24g. The site fractions of Pd and Zn calculated from the present thermodynamic model show the occupancy preference of Pd in the OT and OH sublattices in agreement with experimental observations. The force constants obtained from DFT-based phonon calculations further supports the tendency of Pd occupying the OH sublattice compared with the IT and CO sublattice. The catalytic ensembles changing from Pd monomers (Pd 1 ) to trimers (Pd 3 ) on the surface of γ phase are attributed to the increase of Pd occupancy in the OH sublattice.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Evolution of PTCDA-derived seeds prior to graphene nanoribbon growth on Ge(001)

The synthesis of graphene nanoribbons (GNRs) can be realized via CH4 chemical vapor deposition (CVD) on substrates such as Ge(001) that promote highly anisotropic growth kinetics. Small polycyclic aromatic hydrocarbons (PAHs) are first sublimed onto Ge near room temperature to form graphene-like seeds that subsequently initiate GNR growth with CH4 exposure at 1173 K. The behaviors of PAHs in the time in between their sublimation onto Ge and GNR growth are unclear. Here, we study an archetypical PAH, perylene-3,4,9,10-tetracarboxyl acid dianhydride (PTCDA), on Ge(001) using both scanning tunneling microscopy (STM) and density functional theory (DFT) to characterize PAH conformation, surface diffusivity, and clustering, versus temperature. PTCDA becomes mobile above 673 K (consistent with a DFT calculated diffusion barrier of 2.01 eV). The mobile PTCDA molecules meet, cluster, and fuse at Ge step edges. These clusters have an apparent height of 0.4 nm, similar to small graphene islands, and grow laterally with increasing temperature – reaching 1.2 – 2.1 nm in width at 1173 K, consistent with the seed size extrapolated from CVD experiments. These results provide a plausible picture for how PTCDA forms seeds for anisotropic GNR CVD and show that PAHs with reduced surface diffusivity and inter-molecular reactivity are needed to enable more monodisperse PAH-seeded GNR synthesis.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Ab initio thermodynamics of Ni and Co incorporation in Mg hydroxide, carbonate, and hydroxycarbonate minerals

Ni and Co are critical elements needed for modern technologies, and a better understanding of the ability of Mg-based minerals to incorporate these elements would benefit strategy development for Ni and Co recovery from mafic and ultramafic deposits. Here, in this work, we performed density functional theory (DFT) calculations of Ni and Co incorporation in six potential products of the carbonation of mafic and ultramafic silicates: brucite (Mg(OH) 2 ), magnesite (MgCO 3 ), nesquehonite (MgCO 3 ⸱3H 2 O), lansfordite (MgCO 3 ⸱5H 2 O), artinite (Mg 2 CO 3 (OH) 2 ⸱3H 2 O), and hydromagnesite (Mg 5 (CO 3 ) 4 (OH) 2 ⸱4H 2 O). The DFT results were used in an ab initio thermodynamics framework to explore the pH 2 O–pCO 2 conditions at which the Mg-based minerals were predicted to be thermodynamically stable and to quantify the Gibbs free energy of Ni and Co substitution at Mg sites. Among the six Mg-based minerals, brucite and magnesite were predicted to have the lowest Ni and Co substitution free energy. An analysis of the effect of temperature indicated that, at low temperature (<100 K), brucite more readily accommodated Ni and Co, while, at higher temperature (>335 K), magnesite more favorably incorporated Ni and Co. Between 100 K and 335 K, Ni was predicted to preferentially substitute for Mg in brucite and Co for Mg in magnesite, thus leading to a driving force for separating Ni and Co in conditions where brucite and magnesite both form. Insights gained in this work could therefore help select experimental conditions that either promote or inhibit incorporation of these critical elements into Mg-based mineral phases.

Critical elements↗
DOE OSTI2023

Photoionization of C 60 at high energies

Calculations of the photoionization cross section of C60 including all valence and core subshells have been carried out for photon energies up to 1 keV using density functional theory (DFT) and time-dependent DFT (TDDFT) within the framework of a fully molecular model. The high-energy valence subshell cross sections behave rather differently than the results from model potentials where the cross section falls far too rapidly with energy. This unphysical behavior of the model calculations is traced to smearing out of the carbon nuclei in the model potentials which makes it difficult for momentum to be conserved at the higher energies. Comparison with 60 times carbon atomic cross section sheds light on similarities and differences. The high energy behavior of individual valence cross sections reflects the amount of C 2s/2p mixing in the molecular orbitals.

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

Surfactant inhibition mechanisms of carbonate mineral dissolution in shale

Surfactants are common additives to hydraulic fracturing and enhanced oil recovery (EOR) fluids, and are under consideration for amendment to supercritical carbon dioxide for geological carbon sequestration (GCS). The effect of a common anionic surfactant, internal olefin sulfonate (IOS), on mineral dissolution from shale into brine was evaluated. When added to brine at concentrations exceeding the critical micelle concentration (94 mg/L), IOS inhibited carbonate mineral dissolution in an Eagle Ford shale, as well as dissolution of optical quality calcite (the dominant carbonate in the shale). Laser profilometry images provide spatial resolution across > 3 orders of magnitude, and indicate that IOS addition to brine both enhances the formation of new etch pits in calcite, and impedes their further growth. Time-of-flight secondary ion mass spectrometry surface profiles show for the first time that IOS preferentially adsorbs at calcite pit edges versus flat calcite surfaces (i.e., terraces). Surface pressure calculations, sulfur K-edge near edge X-ray absorption fine structure (NEXAFS) spectroscopy results, and density functional theory (DFT) calculations support this observation; the DFT results indicate that the sulfonate head group of the IOS molecule binds strongly to the calcite step site as compared to the terrace site. The S K-edge NEXAFS results indicate that IOS adsorbed more to etched calcite surfaces compared to smooth calcite surfaces. Overall, the results indicate that weak adsorption on flat calcite surfaces (i.e., terraces) disrupts water structure and enhances mass transfer of dissolution, while strong adsorption on calcite pit edges displaces adsorbed water and inhibits further etch pit growth. This work provides the first direct evidence of preferential adsorption of IOS to etched calcite surfaces and links it to macroscopic dissolution kinetics. Finally, this work has implications for surfactant-containing fluids used in hydraulic fracturing, EOR and potentially GCS for subsurface injection into carbonate rich reservoirs.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Effects of cluster expansion on the locations of phase transition boundary as a first step to quantify uncertainty in first principles statistical mechanics framework

Predicting phase diagrams from first principle calculations eliminates the need of tedious experimental trials and errors. Fully automating first principle phase diagram calculations without any sort of human intervention has been a long daunting task and troubling scientific communities for decades. This grand problem remains not fully resolved, largely due to the vastly high-dimensional parameter space associated with density functional theory, cluster expansion, lattice Monte Carlo, and the substantial uncertainty propagating through a set of complex simulations. As a first step to tackle this grand problem, we reported a first demonstration of how sensitive phase boundary locations can be to various cluster expansion fittings and input DFT training data. To the best knowledge of the authors, this study reported the first ever attempt to quantify uncertainty in first principles statistical mechanics framework. In addition, a semi-automated phase transition detection algorithm has been devised in this paper to deal with the associated statistical errors and uncertainties from Monte Carlo method and its predecessor DFT calculations and cluster expansions. This algorithm has been applied in a classical cluster expansion Mg-Cd binary alloy system to detect phase transitions at various locations of phase diagram using different cluster expansions to demonstrate its predictive power and quantify uncertainties. The results suggested that using Chebyshev basis function shifted the transition locations toward the dilute solid solution phase and expanded the phase stability of concentrated ordered phases to wider composition range. The addition of perturbed defect configurations into training data set lowered the order-disorder transition temperature to be away from true transition temperature, suggesting the transition nature is indeed configurational disorder dominated rather than defect assisted. Finally, we found that the weighting has negligible effect on the transition locations, except for the case of using Chebyshev basis function to fit all non-weighted configurations that can be susceptible to Monte Carlo sampling hysteresis.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Theoretical predictive screening of noble-metal-containing M 3 AuC 2 (M = Ti, V, and Cr) MAX phases

Using first-principles density functional theory (DFT), the ground state physical properties of the newly synthesized noble-metal-containing Ti 3 AuC 2 MAX phase have been investigated. The effect of transition-element Ti replacement (V and Cr) on physical properties, including structural, elastic, electronic, thermal, and optical are presented. The optimized lattice parameters of Ti 3 AuC 2 MAX phase are in good agreement with the experimental values and decrease when we replace Ti with V and then Cr. The magnetic properties of M 3 AuC 2 are predicted by using generalized gradient approximation (GGA) incorporating onsite electron correlation parameter (Hubbard U parameter) within DFT. Out of three studied MAX phases, V 3 AuC 2 and Cr 3 AuC 2 phases exhibit magnetism. The thermodynamical stability of M 3 AuC 2 (M = Ti, V, and Cr) phases is discussed by the formation enthalpy with respect to their most competing phases. The electronic structure reveals that these phases have metallic nature and are electrically anisotropic. The M-element replacement has an effect on the bonding properties of M 3 AuC 2 . The bonding between M-C is in the order of Cr 3 AuC 2 > V 3 AuC 2 > Ti 3 AuC 2 according to their peak positions and heights of density of states in the occupied site, which is also confirmed from the charge density distribution. The elastic properties indicate that M 3 AuC 2 phases are ductile, machinable, less stiff, and better resistant to thermal shock and are elastically anisotropic. To understand the properties of M 3 AuC 2 for the extreme environment, we employed a quasi-harmonic Debye model at the pressure and temperature range of 0–50 GPa and 0–1600 K, respectively. Additionally, we evaluated thermal conductivity and melting temperature to further understand the potential of these MAX phases in coating applications at elevated temperature. In relation to optical properties, these MAX phases have a reasonable absorption coefficient in the visible as well as in the ultra-violet regions. The reflectivity of M 3 AuC 2 phases is polarization dependent and of 45 % against the visible region, which reveals its potential as coating material to minimize solar heating.

312 MAX phases↗
DOE OSTI2022

First-principles study of Al/Al 3 Ni interfaces

Al-Ni alloys have shown promise for high-temperature applications due to the strengthening of Al 3 Ni fibers resistant to coarsening and spheroidization up to 400°C. While the interface between Al and Al 3 Ni phases affects the coarsening rate of Al 3 Ni at elevated temperatures, its characteristics are largely unknown to date. Here, we have constructed various supercells to model this interface and performed a first-principles study based on density functional theory (DFT). We have considered three groups of Al/Al 3 Ni interfaces: experimentally reported orientation relationships from the solidification studies, crystallographically similar Fe-Fe 3 C pearlite interfaces, and the family of low-index (100) termination planes. We have analyzed the correlation between the DFT Al/Al 3 Ni interfacial energies and characteristic features, e.g., excess free volume and the number of broken bonds. We outline the further experimental and computational analysis required to improve the interface modeling of Al/Al 3 Ni.

36 MATERIALS SCIENCE↗
DOE OSTI2022

KSSOLV 2.0: An efficient MATLAB toolbox for solving the Kohn-Sham equations with plane-wave basis set

KSSOLV (Kohn-Sham Solver) is a MATLAB toolbox for performing Kohn-Sham density functional theory (DFT) calculations with a plane-wave basis set. KSSOLV 2.0 preserves the design features of the original KSSOLV software to allow users and developers to easily set up a problem and perform ground-state calculations as well as to prototype and test new algorithms. Furthermore, it includes new functionalities such as new iterative diagonalization algorithms, k-point sampling for electron band structures, geometry optimization and advanced algorithms for performing DFT calculations with local, semi-local, and hybrid exchange-correlation functionals. It can be used to study the electronic structures of both molecules and solids. So, we describe these new capabilities in this work through a few use cases. We also demonstrate the numerical accuracy and computational efficiency of KSSOLV on a variety of examples.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Electron-phonon coupling from GW perturbation theory: Practical workflow combining BerkeleyGW, ABINIT, and EPW

Here we present a workflow of practical calculations of electron-phonon (e-ph) coupling with many-electron correlation effects included using the GW perturbation theory (GWPT). This workflow combines BerkeleyGW, ABINIT, and EPW software packages to enable accurate e-ph calculations at the GW self-energy level, going beyond standard calculations based on density functional theory (DFT) and density-functional perturbation theory (DFPT). This workflow begins with DFT and DFPT calculations (ABINIT) as starting point, followed by GW and GWPT calculations (BerkeleyGW) for the quasiparticle band structures and e-ph matrix elements on coarse electron k- and phonon q-grids, which are then interpolated to finer grids through Wannier interpolation (EPW) for computations of various e-ph coupling determined physical quantities such as the electron self-energies or solutions of anisotropic Eliashberg equations, among others. A gauge-recovering symmetry unfolding technique is developed to reduce the computational cost of GWPT (as well as DFPT) while fulfilling the gauge consistency requirement for Wannier interpolation.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Radionuclide incorporation in negative thermal expansion α-Zr(WO 4 ) 2 : A density functional theory study

The incorporation of uranium, plutonium and technetium in the negative thermal expansion (NTE) α-Zr (WO 4 ) 2 has been investigated within the framework of density functional theory (DFT). It is found that the vacancy formation energies of the charged vacancies are overall larger than that of its counterpart neutral Frenkel defects and Schottky defects. DFT calculations suggest that U and Pu substitutions for the Zr site are preferred in α-Zr (WO 4 ) 2 . In case of Tc substitution, both Tc(IV) for the Zr site and Tc(VII) for the W site are considered under oxygen-poor and oxygen-rich conditions, while Tc(VII) substitution can be improved significantly by including Y 2 O 3 (charge compensation).

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

The role of nonmetallic ion substitution in perovskite LaCoO 3 for improved oxygen evolution reaction activity

Transition metal perovskite (ABO 3 ) is an emerging type of oxygen evolution reaction (OER) electrocatalyst that shows reasonably good activity and moderate stability. Although efforts have been made to improve perovskite’ OER performance by various element substitution at A/B-site, the influence of ion, particularly non-metallic ion, substitutions on the OER mechanism are rarely studied. More and more evidence has shown that the metal-center theory has failed to explain lots of OER-related phenomena. Therefore, it is urgent to understand how the cation and anion sites in perovskite determine OER performance. In this work, we used a Fe and P co-doped LaCoO 3 as a model system to explore the influence of substitution in perovskite by combining operando/ex-situ X-ray characterization and density functional theory (DFT). We observed enhanced OER catalytic activities in co-doped materials, which are attributed to the stronger transition-metal-oxygen-bonding-covalency (TMOBC). The detailed analyses by O K-edge XAS, electrochemical performance, and DFT suggest that the hybridization between O 2p and transition metal 3d e g orbitals could be a more credible descriptor of perovskite for OER, which is the combination of e g orbital theory and TMOBC theory. The finding in our work provides insights into the OER catalysis mechanism on metal oxides, which could guide new design of cost-effective oxide electrocatalysts.

25 ENERGY STORAGE↗
DOE OSTI2024

Core-to-core X-ray emission spectra from Wannier based multiplet ligand field theory

Recent advances using Density Functional Theory (DFT) to augment Multiplet Ligand Field Theory (MLFT) have led to ab-initio calculations of many formerly empirical parameters. This development makes MLFT more predictive instead of interpretive, thus improving its value for studies of highly correlated 3 d , 4 d , and ƒ-electron systems. Synchrotron time is always at a premium, and tools that provide predictive capabilities have clear value when it comes to planning studies. Here, in this work, we develop a DFT + MLFT based approach for core-to-core Kα x-ray emission spectra (XES) and evaluate its performance for a range of transition metal systems. We find good agreement between theory and experiment, as well as the ability to capture key spectral trends related to spin and oxidation state. We also discuss limitations of the model in the context of the remaining free parameters and suggest directions forward.

Ab initio theoretical treatment↗
DOE OSTI2022

High temperature stability and transport characteristics of hydrogen in alumina via multiscale computation

Here, the impact of hydrogen charge states on the stability and transport characteristics of hydrogen interstitials in alumina polymorphs is evaluated by multiscale computational methods including density functional theory (DFT), ab initio molecular dynamics (AIMD) and machine learned force fields. Thermodynamic calculations show that the protonic H i +1 interstitial is the most stable defect species for most values of the electronic bandgap in both and amorphous alumina (Al 2 O 3 ). Further, active learned Gaussian approximation potentials (GAP) were developed using AIMD data to study temperature dependent long time proton diffusion in alumina. Diffusivity calculations from GAP-MD simulations are found to be comparable with of the AIMD data, while being ~340 times faster and scalable to larger systems. Comparisons with diffusivity values for other interstitial charge states (H i 0 and H i -1 ) and published experimental literature indicate that H i +1 diffusion is the likely mechanism of hydrogen transport. A good agreement is obtained between H i +1 diffusivity calculated in α-Al 2 O 3 from DFT: 5.05 10 -3 exp(-0.81 eV/k B /T) cm 2 /s and reported experiment: 9.7X10 -4 exp(-0.83 eV/k B /T) cm 2 /s. Computationally and experimentally calculated energy barriers (0.81 and 0.83 eV respectively) only differ by 2.5%. Similarly, the pre-exponential diffusion coefficients only differ by 0.5 orders of magnitude. Moreover, the diffusivity of H i +1 in amorphous Al 2 O 3 in the 1000–2000 K range is calculated to be 2.53X10 -2 exp(-0.89 eV/k B /T), just one order of magnitude higher than the corresponding value in α-Al 2 O 3 . This suggests that local structural disorder does not significantly affect the energy landscape and diffusion behavior of H i +1 in Al 2 O 3 . Overall, these results show promise for the application of alumina polymorphs as hydrogen permeation barriers.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Effect of hydrogen on plasticity of $α$-Fe: A multi-scale assessment

A multi-scale study was carried out to quantify the effect of interstitial hydrogen concentration on plasticity in $α$-Fe. Here, in this work, the influence of hydrogen on the screw dislocation glide behavior was examined across several length-scales. The insights obtained were integrated to provide an accurate continuum description for the effect of hydrogen on the dislocation based plasticity in polycrystalline $α$-Fe. At the outset of this work, a new Fe—H interatomic potential was formulated that enhanced the atomistic estimation of the variation in dislocation glide behavior in presence of hydrogen. Next, the dislocation core reconstruction observed due to the addition of hydrogen using atomistic simulations was validated with the help of large-scale DFT calculations based on the DFT-FE framework. Several atomistic simulations were carried out to comprehensively quantify the effect of hydrogen on the non-Schmid behavior exhibited during the dislocation glide in $α$-Fe. Finally, crystal plasticity simulations were carried out to understand the effect of hydrogen on the meso-scale deformation behavior of polycrystalline $α$-Fe.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Relative phase stability of L1 2 and DO 22 /DO 23 structures in Al 3 Nb, Al 3 Zr and Al 3 V compounds

The relative stability of the different tri-aluminide (Al 3 M) phases in three binary systems (M = Zr, Nb and V) was assessed for their potential to form fine cubic L1 2 precipitates in additively manufactured alloys. Supersaturated thin films of Al-(8–30) at% M were sputtered and heat treated during in-situ x-ray diffraction (XRD) measurements to observe the temperature ranges of stability for each phase. As-sputtered films were then processed with laser tracks simulating additive manufacturing solidification conditions, and the formation of phases in the laser tracks was correlated with density functional theory (DFT) and nucleation rate calculations. We found that the metastable L1 2 structure is highly competitive with the stable DO 23 structure in the Al-Zr system, but much less stable than the DO 22 structure in the Al-Nb system, and both the DO 22 and Al 8 V 5 structure in the Al-V system. Furthermore, these experimental results were found to be in good agreement with the DFT and kinetic calculations, as we determined that the metastable L1 2 in Al-Zr only requires a small amount of undercooling to favor its nucleation over the stable DO 23 , suggesting additive manufacturing can be a viable pathway to develop Al-Zr alloys strengthened by a high volume fraction of L1 2 Al 3 Zr phase.

Perrin, Alice E. [Oak Ridge National Laboratory (O↗