Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DFT 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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

First principles calculations in support of Pu aging: calculating the effects of lattice imperfections on thermodynamics

Plutonium (Pu) has, in theory, well defined crystal structures: its atoms are arranged in regular spatial patterns. But Pu is radioactive, and as its nuclei decay those regular spatial patterns are interrupted. The interruptions are lattice imperfections, which are known to affect how materials respond to their environment. To adequately model Pu, we need to know which lattice imperfections are present, how they interact with each other, and how they affect the material’s response to its environment. The work presented here aims to use density functional theory (DFT) calculations to begin to answer the latter, in particular, how individual lattice imperfections affect measurable effects including thermal expansion (the change in volume in response to a change in temperature), heat capacity (the amount of thermal energy needed to change a material’s temperature), and elastic moduli (a material’s resistance to applied stresses). Pu poses many computational challenges. The $\textit{f}$ electrons require special attention. Of all the elements, Pu has the largest number of electrons that must be included in the calculations. Thermal effects require calculating the phonons (the lattice vibrations), which for systems containing lattice imperfections demand large, complex computational cells - but computational resources limit the size and complexity. With careful restructuring of how the calculations are performed, all these challenges have been met to enable calculations that provide insight into how lattice imperfections affect Pu’s response to its environment. Reported here are the computational challenges and the advances developed to meet them, along with the first results showing the strong effect that one prototype lattice imperfection (an interstitial Pu atom in a delta-phase Pu lattice) has on Pu’s response to its environment.

36 MATERIALS SCIENCE↗
DOE OSTI2023

High reactivity of mesoporous CeO 2 to dissociate chemical warfare agent sarin

We perform spectroscopic surface characterization on mesoporous ceria (CeO 2 ) under exposure of sarin gas, also known as GB. We find mesoporous CeO 2 displays extraordinary reactivity towards GB dissociation under high vacuum. The spectroscopic results combined with density functional theory (DFT) calculations suggest that the cleavage of the P–F and P–OR bonds are both possible pathways for GB to dissociate on the CeO 2 surfaces. The modeling reveals differences in dissociation pathways of GB on the crystallographic (111)- and (110)-surfaces of ceria. Here we also report some dissociation of GB on mesoporous CeO 2 under ambient conditions. The finding that mesoporous CeO 2 can degrade GB under mild conditions is encouraging and implies the strong potential of ceria-based materials to be applied in combating chemical warfare agents related to the GB series.

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

A Quantum Mechanical Study of Quartz (101) Interfacial Boundaries

Density functional theory (DFT) calculations were performed on periodic and molecular models to explore the energetics of bonding modes at quartz–quartz interfaces, common grain contacts in sandstones. Four interface types were modeled: H 2 O-mediated, silanol H-bond mediated, O–O peroxide bonds, and siloxane bonds. Each type of interface may exist in the subsurface of the Earth, depending upon water activity, temperature, and pressure, as interfacial structures are stable within DFT molecular dynamics simulations up to 473 K. The results predict interfacial energies SiOSi > SiOOSi > SiOH-6H 2 O > 2SiOH, as expected based on similar bond strengths in other systems. The proposed SiOOSi interface is not generally accounted for in chemomechanical models but could be important in sandstones at depth in the Earth. The implications for intergrain fracturing and H2 generation are discussed.

Hydrogen↗
DOE OSTI2020

A comprehensive diffusion mobility database comprising 23 elements for magnesium alloys

We report that reliable experimental diffusion coefficients of 10 key alloying elements in Mg obtained by the present authors together with experimental data in the literature enabled us to perform a systematic test of the reliability of diffusion coefficients obtained from DFT calculations. The computed activation energy values were found to be quite accurate (mostly within 0.2 eV) but the computed pre-factors were less reliable. Such insights allowed us to develop a practical and yet robust strategy to perform diffusion mobility assessments by adopting the computed activation energy while fitting only the pre-factor when available experimental data are limited to a narrow temperature range. The overall good agreement between the DFT data and experimental data also gave us the confidence to employ the computed data for those that were still missing or inaccessible from experimental measurements. A systematic assessment of both the measured and computed diffusion data in hcp Mg was performed using the above holistic approach to yield the most comprehensive open Mg mobility database to date, comprising 23 elements (Mg, Ag, Al, Be, Ca, Cd, Ce, Cu, Fe, Ga, Gd, In, La, Li, Mn, Nd, Ni, Pu, Sb, Sn, U, Y, Zn). This more reliable mobility database will contribute to future development of advanced Mg alloys. The holistic approach developed in this study will be very beneficial to the future establishment of reliable mobility databases for other alloy systems as well.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Magnetic states of the quasi-one-dimensional iron chalcogenide Ba 2 FeS 3

Quasi-one-dimensional iron-based ladders and chains, with the 3d iron electronic density n=6, are attracting considerable attention. Recently, a new iron chain system Ba 2 FeS 3 , also with n=6, was prepared under high-pressure and high-temperature conditions. Here the magnetic and electronic phase diagrams are theoretically studied for this quasi-one-dimensional compound. Based on first-principles calculations, a strongly anisotropic one-dimensional electronic band behavior near the Fermi level was observed. In addition, a three-orbital electronic Hubbard model for this chain was constructed. Introducing the Hubbard and Hund couplings and studying the model via the density matrix renormalization group (DMRG) method, we studied the ground-state phase diagram. A robust staggered ↑-↓-↑-↓ AFM region was unveiled in the chain direction, consistent with our density functional theory (DFT) calculations. Furthermore, at intermediate Hubbard U coupling strengths, this system was found to display an orbital selective Mott phase (OSMP) with one localized orbital and two itinerant metallic orbitals. At very large U/W (W=bandwidth), the system displays Mott insulator characteristics, with two orbitals half-filled and one doubly occupied. Our results for high pressure Ba 2 FeS 3 provide guidance to experimentalists and theorists working on this one-dimensional iron chalcogenide chain material.

1-dimensional systems↗
DOE OSTI2024

Molten Flux Synthesis of Plutonium (IV) Silicates K 2 PuSi 2 O 7 and Rb 2 PuSi 6 O 15

Few plutonium containing silicates have been reported in the literature and only one other tetravalent plutonium silicate structure has been reported using the flux growth method. With the goal of understanding actinide incorporation into silicates, we describe herein two new plutonium (IV) silicates, Rb 2 PuSi 6 O 15 and K 2 PuSi 2 O 7 , synthesized using a mixed alkali chloride/alkali fluoride flux growth. Both structures crystallize in the C2/c space group with lattice parameters a=24.2980(10) Å, b=7.1466(3) Å, c=17.2025(8) Å and $\blacksquare$96.6560(10)° for Rb 2 PuSi 6 O 15 , and a=9.9478(2) Å, b=5.59820 (10) Å, c=13.1848(2) Å, and $\blacksquare$105.7400(10)° for K 2 PuSi 2 O 7 . Comparisons to the synthesis of other reported alkali actinide silicates (Ce, U, and Th) are made. Raman spectra and DFT calculations of formation enthalpies and vibrational modes provide confirmation of structure refinements and enhanced product characterization for comparison with other materials.

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

The Intermetallic Reactivity Database: Compiling Chemical Pressure and Electronic Metrics toward Materials Design and Discovery

Here, the advent of high-throughput Density Functional Theory (DFT) calculations has supported the creation of large databases containing the quantitative output necessary for constructing theoretical phase diagrams and predicting physical properties. In this Article, we present a complementary resource, the Intermetallic Reactivity Database (IRD), focused on the chemical bonding features of solid-state structures and indicators of potential structural transformations. Each IRD entry augments common features, such as band structures and density of states (DOS) distributions, with chemically motivated information including DFT-Chemical Pressure (CP) schemes and visualizable representations of the atomic charges. Together, these data types enable the rationalization and prediction of potential structural phenomena encountered in intermetallic chemistry, as we illustrate with four examples: the origins of the Y 2 Ni 2 Mg structure in terms of CP features of its parent structures, the anticipation of intergrowth phases from the net atomic CPs collected in Al-containing binary phases, the correlation between trends in the CP schemes of CaCu 5 -type phases and experimentally observed structural variations, and finally, the development of theoretical methodology with the testing of a streamlined method generating DFT-CP schemes. Altogether, these examples highlight how the IRD supports the creation of models of structural chemistry that extend beyond the bounds of its entries.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Scaling of Transition State Vibrational Frequencies and Application of d -Band Theory to the Brønsted–Evans–Polanyi Relationship on Surfaces

Semiempirical energy relations provide a means of estimating thermodynamic properties. Specifically, linear scaling relationships (LSRs) and Brønsted–Evans–Polanyi (BEP) relationships correlate adsorption energies between adsorbates across surfaces and reaction energies with activation barriers, respectively. Although vibrational scaling relations (VSRs) exist between adsorbates at identical sites, scaling between vibrational frequencies of adsorbed local minima and transition states is lacking. Here, we present density functional theory calculations for AH X (A = C, N, O) diffusions on transition metal surfaces and reveal linear scaling between frequencies of local minima and the transition state between those minima. Using d-band theory and linear muffin tin orbital theory (LMTO), we derive the slopes of these transition state vibrational scaling relations (TSVSRs) and, in so doing, provide a rigorous theory extending the original BEP relations developed for solution chemistry to surface chemistry. Furthermore, with a single reference DFT calculation, we predict the slopes and quantify uncertainty in the predictions.

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

An automatically curated first-principles database of ferroelectrics

Ferroelectric materials have technological applications in information storage and electronic devices. The ferroelectric polar phase can be controlled with external fields, chemical substitution and size-effects in bulk and ultrathin film form, providing a platform for future technologies and for exploratory research. In this work, we integrate spin-polarized density functional theory (DFT) calculations, crystal structure databases, symmetry tools, workflow software, and a custom analysis toolkit to build a library of known, previously-proposed, and newly-proposed ferroelectric materials. With our automated workflow, we screen over 67,000 candidate materials from the Materials Project database to generate a dataset of 255 ferroelectric candidates, and propose 126 new ferroelectric materials. We benchmark our results against experimental data and previous first-principles results. The data provided includes atomic structures, output files, and DFT values of band gaps, energies, and the spontaneous polarization for each ferroelectric candidate. Finally, we contribute our workflow and analysis code to the open-source python packages atomate and pymatgen so others can conduct analogous symmetry driven searches for ferroelectrics and related phenomena.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Correlating binding energies of adsorbed CO and H on model surfaces with CO/H 2 selectivity from co-electrolysis of CO 2 and H 2 O over copper–palladium bimetallic catalysts

Binding energies of adsorbed CO and H are key descriptors governing the activity and selectivity of the co-electrolysis of CO 2 and H 2 O to produce syngas with desired CO/H 2 ratios. Palladium hydride (PdH), which forms in situ at negative overpotentials, has been identified as the active Pd phase for CO 2 reduction to syngas. Herein, binding energies of CO and H are determined using temperature programmed desorption (TPD) of CO and H 2 from Pd(111), PdH/Pd(111), and Cu/PdH/Pd(111) under ultra-high vacuum (UHV) conditions. TPD results reveal that desorption of H 2 from subsurface PdH occurs at 460 K, while desorption from surface PdH is more facile at 320 K. CO desorption temperatures shift 20 K lower on PdH/Pd(111) compared to on Pd(111). The presence of 0.7 ML Cu further increases the desorption temperature of H 2 by 30 K while simultaneously reducing CO desorption temperatures by 70 K. Density functional theory (DFT) calculations show that CO adsorption onto Pd sites is hindered on the 0.7 ML Cu/PdH/Pd(111) surface while the kinetic barrier for H 2 desorption is increased. The trends in the binding energies of CO and H on model surfaces are consistent with electrochemical measurements of CuPd powder catalysts in a membrane electrode assembly (MEA), where H 2 evolution is reduced while CO production is enhanced compared to unmodified Pd catalysts. Overall, the results from model surface studies (TPD and DFT) provide a prediction and explanation for the activity and CO/H 2 ratios observed in electrochemical experiments. This study also demonstrates that CuPd is a promising catalyst with reduced Pd-loading to produce CO-rich syngas.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Low-Temperature Activation and Coupling of Methane on MgO Nanostructures Embedded in Cu 2 O/Cu(111)

Here, the efficient conversion of methane into valuable hydrocarbons such as ethane and ethylene at relatively low temperatures without deactivation issues is crucial for advancing sustainable energy solutions. Herein, AP-XPS and STM studies show that MgO nanostructures (0.2-0.5 nm wide, 0.4-0.6 Å high) embedded in a Cu 2 O/Cu(111) substrate activate methane at room temperature, mainly dissociating it into CH x (x = 2 or 3) and H adatoms, with minimal conversion to C adatoms. These MgO nanostructures in contact with Cu 2 O/Cu(111) exhibit unique reactivity, enabling C-C coupling into ethane and ethylene at 500 K, a significantly lower temperature than that required for bulk MgO catalysts (>700 K), with negligible carbon deposition and no deactivation. DFT calculations corroborate these experimental findings. The CH 4,gas → *CH 3 +*H reaction is a downhill process on MgO/Cu 2 O/Cu(111) surfaces. The activation of methane is facilitated by an electron transfer from copper to MgO and the existence of Mg and O atoms with a low coordination number in the oxide nanostructures. The formation of O-CH 3 and O-H bonds overcomes the energy necessary for the cleavage of a C-H bond in methane. DFT studies reveal that smaller Mg 2 O 2 model clusters provide stronger binding and lower activation barriers for C-H dissociation in CH 4 , while larger Mg 3 O 3 clusters promote C-C coupling due to weaker *CH 3 binding. All these results emphasize the importance of size when optimizing the catalytic performance of MgO nanostructures in the selective conversion of methane.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Final Technical Report for "Cloud-based Low-Scaling Quantum Chemistry Simulations for Materials"

The objective of phase I was to develop the basic infrastructure for performing efficient hybrid DFT calculations on solid materials with a turnaround time amenable to the use in large-scale data-driven approaches. To fulfill the goal we have developed a novel algorithm that significantly reduces the cost of exchange matrix formation. The algorithm does so by making use of three ingredients (a) interpolative decomposition of the electron integrals (b) robust pseudospectral method and (c) occ-RI approach. Our published work demonstrates that the algorithm is orders of magnitude faster than any other hybrid-DFT method and is on the order of only 3-4 times slower than pure DFT. All steps of the algorithm developed during phase I can be accelerated to the point that, for large enough systems, the diagonalization of the Fock matrix will become the most expensive step. To treat such large systems we have created an interface to the ASCR-funded PEXSI library. The PEXSI method is used to compute the density matrix directly from the Fock matrix in a manner that preserves the sparsity of the local representation.

Shiozaki, Toru↗
DOE OSTI2025

Grain boundary zirconia-modified garnet solid-state electrolyte

Here, we report a method for promoting electrochemical stability in garnet Li 6.4 La 3 Zr 1.4 Ta 0.6 O 12 solid-state electrolyte based on a composite two-phase oxide–oxide microstructure. Grain boundary precipitation of the controlled distribution of amorphous zirconium oxide microparticles is achieved through the addition of reactive tantalum carbide. During ambient-atmosphere sintering, the carbide decomposes through an in situ reaction, the ‘extra’ Ta substituting for Zr within the Li 6.4 La 3 Zr 1.4 Ta 0.6 O 12 lattice. Density functional theory (DFT) calculations identify a thermodynamically favourable reaction path and show how substituting Ta 5+ at Zr 4+ sites affects the crystal structure as well as bulk ionic and electronic conductivities. Quantitative stereology highlights that zirconia also acts as a sintering aid, reducing compact porosity. Cryogenic focused-ion-beam scanning electron microscopy and fractography analysis of cycled solid-state electrolytes illustrates that near-universally observed intergranular Li-metal dendrite propagation is suppressed by the two-phase microstructure, favouring transgranular dendrites instead. Importantly, DFT demonstrates that compared with the Li 6.4 La 3 Zr 1.4 Ta 0.6 O 12 surface, the zirconium oxide surface per se is less electronically conductive and does not trap excess electrons to reduce Li ions. This is a key reason for the substantial improvement in the electrochemical properties over the single-phase baseline.

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 OSTI2020

C–H Bond Activation via U(II) in the Reduction of Heteroleptic Bis(trimethylsilyl)amide U(III) Complexes

Reduction of (C 5 Me 5 ) 2 U III (NR 2 ) and (C 5 Me 5 )U III (NR 2 ) 2 (R = SiMe 3 ) with potassium graphite in the presence of 2.2.2-cryptand (crypt) generates dark solutions that have UV-visible spectra consistent with time-dependent density functional theory (TDDFT) calculations on the U(II) products, [(C 5 Me 5 ) 2 U II (NR 2 )] - and [(C 5 Me 5 )U II (NR 2 ) 2 ] - . However, the solutions quickly change color and form the U(III) C–H bond activation products [K(crypt)][(C 5 Me 5 ) 2 U III (CH 2 SiMe 2 NSiMe 3 –κC,κN)], 1, and [K(crypt)][(C 5 Me 5 )U III (NR 2 )(CH 2 SiMe 2 NSiMe 3 –κC,κN)], 2, that were identified by X-ray crystallography. DFT calculations on the putative [(C 5 Me 5 ) 2 U II (NR 2 )] - and [(C 5 Me 5 )U II (NR 2 ) 2 ] - complexes revealed 5f 3 6d 1 ground state electron configurations as previously found in isolable [(C 5 H 4 SiMe 3 ) 3 U II ] - , which indicated that these low symmetry heteroleptic complexes are reasonable precursors for new U(II) complexes.

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

Magnetization reversal driven by electron localization-delocalization crossover in the inverse spinel Co 2 VO 4

Neutron diffraction, magnetization, and muon spin relaxation measurements, supplemented by density functional theory (DFT) calculations are employed to unravel temperature-driven magnetization reversal in inverse spinel Co 2 VO 4 . All measurements show a second-order magnetic phase transition at T C = 168 K to a collinear ferrimagnetic phase. Neutron diffraction measurements reveal two antiparallel ferromagnetic (FM) sublattices, belonging to magnetic ions on two distinct crystal lattice sites, where the relative balance between the two sublattices determine the net FM moment in the unit cell. As the evolution of the ordered moment with temperature differs between the two sublattices, the net magnetic moment reaches a maximum at T NC = 138 K and reverses its sign at T MR = 65 K. The DFT results suggest that the underlying microscopic mechanism for the reversal is a delocalization of the unfilled 3d -shell electrons on one sublattice just below T C , followed by a gradual localization as the temperature is lowered. Furthermore, this delocalized-localized crossover is supported by muon spectroscopy results, as strong T 1 relaxation observed below T C indicates fluctuating internal fields.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Solution-Based, Anion-Doping of Li 4 Ti 5 O 12 Nanoflowers for Lithium-Ion Battery Applications

Solution–based, anionic doping represents a convenient strategy with which to improve upon the conductivity of candidate anode materials such as Li 4 Ti 5 O 12 (LTO). As such, novel synthetic hydrothermally–inspired protocols have primarily been devised herein, aimed at the large–scale production of unique halogen–doped, micron–scale, three–dimensional, hierarchical LTO flower–like motifs. Although fluorine (F) doping has been explored, the use of chlorine (Cl) dopants is the primary focus here. Several experimental variables, such as dopant amount, lithium hydroxide concentration, and titanium butoxide purity, were probed and perfected. Furthermore, the Cl doping process did not damage the intrinsic LTO morphology. The analysis, based on interpreting a compilation of SEM, XRD, XPS, and TEM–EDS results, was used to determine an optimized dopant concentration of Cl. Electrochemical tests demonstrated an increased capacity via cycling of 12 % for a Cl–doped sample as compared with pristine LTO. Moreover, the Cl–doped LTO sample described in this study exhibited the highest discharge capacity yet reported at an observed rate of 2C for this material at 143mAh g –1 . Overall, these data suggest that the Cl dopant likely enhances not only the ion transport capabilities, but also the overall electrical conductivity of our as–prepared structures. Furthermore, to help explain these favorable findings, theoretical DFT calculations were used to postulate that the electronic conductivity and Li diffusion were likely improved by the presence of increased Ti 3+ ion concentration coupled with widening of the Li migration channel.

25 ENERGY STORAGE↗