Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Materials Modeling, Density Functional Theory, Computational Materials Science”

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

A prediction of the thermodynamic, thermophysical, and mechanical properties of CrTaO 4 from first principles

It is reported that the self-forming CrTaO 4 oxide scale can protect refractory high-entropy alloys from oxidation, superior to Cr 2 O 3 . In this paper, the phase stability, mechanical, and thermal properties of three polymorphous phases of CrTaO 4 are systematically investigated from first-principles density functional theory calculations. The mechanical properties predicted using the strain–energy methods indicated that all three phases are mechanically stable. The temperature dependence of elastic constants and polycrystalline moduli of three phases demonstrated the thermal softening as temperature increase. The Helmholtz-free energies as a function of volume and temperature are derived from phonon dispersions within the quasi-harmonic approximation at six strained volumes. The calculated apparent bulk coefficients of thermal expansion of these three phases are evaluated, the highest value approximately 13.4× 10 –6 K –1 within a temperature range of 500–2000 K for the rutile I4 1 md phase. The lattice thermal conductivity calculated by the Debye–Callaway model suggested that the rutile type I4 1 md phase has the lowest value of approximately 2.1 W/m/K at 1800 K. The other two phases, C2/m and P2/c, exhibit higher values due to relatively lower Grüneisen parameters and larger phonon velocities. The melting point of CrTaO 4 is predicted to be between 1975 and 2449 K using ab initio molecular dynamics simulations. Finally, this work provides a comprehensive theoretical understanding of the thermodynamic, mechanical, and thermal properties for the new material CrTaO 4 and serves as an example of a viable computational design strategy for improved oxidation resistance of refractory alloys at high temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Mechanical and Thermal Properties of Low-Density Al 20+x Cr 20-x Mo 20-y Ti 20 V 20+y Alloys

Refractory high-entropy alloys (RHEAs) Al 20+x Cr 20-x Mo 20-y Ti 20 V 20+y ((x, y) = (0, 0), (0, 10), and (10, 15)) were computationally studied to obtain a low density and a better mechanical property. The density functional theory (DFT) method was employed to compute the structural and mechanical properties of the alloys, based on a large unit cell model of randomly distributed elements. Debye–Grüneisen theory was used to study the thermal properties of Al 20+x Cr 20-x Mo 20-y Ti 20 V 20+y . The phase diagram calculation shows that all three RHEAs have a single body-centered cubic (BCC) structure at high temperatures ranging from 1000 K to 2000 K. The RHEA Al 30 Cr 10 Mo 5 Ti 20 V 35 has shown a low density of 5.16 g/cm 3 and a hardness of 5.56 GPa. The studied RHEAs could be potential candidates for high-temperature application materials where high hardness, ductility, and low density are required.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Effect of applied potential on metal surfaces: Surface energy, Wulff shape and charge distribution

Here we use grand canonical density functional theory to predict the surface energies, Wulff shapes, charge distributions and catalytically active sites of different metal surfaces under electrochemical conditions. We propose a method for computing surface energies from grand canonical density functional theory (GC-DFT) calculations of periodic slab models and use it to compute the surface energies of the facets of Pt, Cu, and Ag crystals to predict their Wulff shapes under electrochemical conditions. GC-DFT predicts that, for the pure metals studied, solvation only slightly affects the Wulff shape while applied potentials considerably affect the surface energies and corresponding Wulff shapes. We used Bader charge analysis of GC-DFT computed electron densities to investigate the effect of applied potential on the distribution of electron density over the atoms of the surfaces of Pt, Cu, Ag, and the 75–25 Ag-Pt and Au-Ni alloys. This analysis shows that, under an applied potential, the electron density is unevenly distributed over the surface atoms and that the charges of atoms more exposed to solvent are more sensitive to bias. Our results show that the most sensitive atom to bias can be used to identify the most favorable adsorption site and thus, the active sites of electrochemical reactions, which is computationally less demanding than calculating the adsorption energies on all possible adsorption sites.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Orthorhombic Structures as Inorganic Halide Perovskite Models for High-Throughput DFT Investigations

Although high-symmetry Pm-3m space group cubic models are computationally efficient for high-throughput density functional theory (DFT) calculations of inorganic ternary (ABX 3 ) halide perovskites (HPs), they frequently predict band gaps (E g ) that disagree with experiment. Conversely, while low-symmetry cubic polymorphous networks (PN) comprised of 160 to 320 atoms incur significantly greater computational cost, they predict E g ’s that are more closely aligned with experiment. Here, in this study, we compare the DFT total energies and E g ’s predicted by four high-symmetry structure models (Pnma orthorhombic, R3m trigonal, P4/mbm tetragonal, and Pm-3m cubic) to cubic PNs for 5 experimentally characterized ternary HPs and find that the orthorhombic model computes Eg’s with the smallest MAD of 0.23 eV relative to the PNs. Pair distribution functions and DFT-computed total energies show that octahedral tilting, which is present in the 20-atom orthorhombic and 160-atom cubic PN models but not in the 5-atom cubic models, stabilizes all 5 compositions in our benchmarking set. We also find that imposing PN constraints when generating and optimizing these orthorhombic structures by fixing the unit cell lattice vectors and displacing the atoms prior to ionic relaxation with DFT, which we call the orthorhombic surrogate model (OSM), lowers the MAD of E g predictions to 0.09 eV. Our OSM predicts the PN band gaps of an additional 95 theoretical inorganic ternary HPs with MAD of 0.08 eV, supporting its usage in high-throughput DFT investigations to closely estimate PN band gaps with much less computational expense.

14 SOLAR ENERGY↗
DOE OSTI2021

Tuning magnetic antiskyrmion stability in tetragonal inverse Heusler alloys

The identification of materials supporting complex, tunable magnetic order at ambient temperatures is foundational to the development of new magnetic device architectures. In this work, we report the design of Mn 2 XY tetragonal inverse Heusler alloys that are capable of hosting magnetic antiskyrmions whose stability is sensitive to elastic strain. Here We first construct a universal magnetic Hamiltonian capturing the short- and long-range magnetic order which can be expected in these materials. This model reveals critical combinations of magnetic interactions that are necessary to approach a magnetic phase boundary, where the magnetic structure is highly susceptible to small perturbations such as elastic strain. We then computationally search for quaternary Mn 2 (X 1 , X 2 ) Y alloys where these critical interactions may be realized and which are likely to be synthesizable in the inverse Heusler structure. We identify the Mn 2 Pt 1–z X z Ga family of materials with X = Au, Ir, Ni as an ideal system for accessing all possible magnetic phases, with several critical compositions where magnetic phase transitions may be actuated mechanically.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Computational modeling of grain boundary segregation: A review

Nearly all metals, alloys, ceramics, and their associated composites are polycrystalline in nature, with grain boundaries that separate well-defined crystalline regions that influence materials properties. In all but the most pure elemental systems, intentional solutes or impurities are present and can segregate to, or less commonly away from, the grain boundaries, in turn influencing boundary behavior, their stability, and associated materials properties. In some cases, grain-boundary segregation can also trigger “phase-like” structural transitions that dramatically alter the essential nature of the boundary. With the development of advanced electron microscopy techniques, researchers can directly observe grain-boundary structures and segregation with atomic precision. Despite such spatial resolution, the underlying mechanisms governing grain-boundary segregation remain difficult to characterize. As a result, computational modeling techniques such as density functional theory, molecular dynamics, mesoscale phase-field, continuum defect theory, and others are important complementary tools to experimental observations for studying grain-boundary segregation behavior. In conclusion, these computational methods offer the ability to explore the underlying formation mechanisms of grain-boundary segregation, elucidate complex segregation behavior, and provide insights into solutions to effectively controlling microstructure.

36 MATERIALS SCIENCE↗
DOE OSTI2020

w19_OMEC Scientific Highlight [Slides]

Scientific Achievement: Revealed the effect of graphite substrate on H binding / hydrogen evolution reaction (HER) activity of organic molecular electrocatalysts (OMECs); and demonstrated the effect of molecular modification on calculated HER activity descriptor. Significance and Impact: Understanding of HER processes is important for discovering low-cost OMECs for hydrogen production and energy application; design of OMEC structures remains challenging due to the unexplored redox properties of organic molecules; and density functional theory (DFT) has been proposed to provide needed insights into HER activity for improved materials design. Research Details: DFT calculations on the binding energy of HER intermediates are used for predicting reduction potentials and and calculated H adsorption energy as activity descriptor; and LANL Institutional Computing (IC) resources are essential for successful project execution given the size of structures involved in modeling molecule-support interaction.

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

A Comparison of Exact and Model Exchange–Correlation Potentials for Molecules

Accurate exchange–correlation (XC) potentials for three-dimensional systems–via solution of the inverse density functional theory (DFT) problem–are now available to test the quality of DFT approximations. Herein, the exact XC potential for seven molecules–dihydrogen at four different bond-lengths, lithium hydride, water, and ortho-benzyne–are computed from full configuration interaction reference densities. These are compared to model XC potentials from nonlocal (B3LYP, HSE06, SCAN0, and M08-HX) and semilocal/local (SCAN, PBE, and PW92) XC functionals. Whereas for most systems, relative errors in the ground-state densities are O(10 –3 –10 –2 ) , the model XC potentials have much higher errors of O(10 –1 –10 0 ) . Among the model XC functionals, SCAN0 offers the best agreement with the exact XC potential, underlining the significance of satisfying exact conditions as well as including nonlocal effects in XC functionals. This work indicates that tests against the exact XC potential will provide a promising new direction for building more accurate XC functionals for DFT.

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

Development of neural network force fields for corrosion studies

To fully understand the chemistry and physics of corrosion, novel methods of simulation must be developed. One approach is designing machine learning (ML) algorithms integrated with density functional theory to develop adaptive force fields to gain insight into corrosion behavior namely at the surface of metal oxides. Current methods of modeling corrosion are slow due to the computational cost of resolving both reaction mechanics and mass transport processes. Machine learning methods can be implemented to obtain structure-activity relationships at both the molecular and bulk scale while still retaining the accuracy of density functional theory (DFT) and significantly decreasing the time needed for simulations of complex chemical processes in the various environments of corrosion. Multiscale models are needed for corrosion studies to fully understand its processes not only at the atomic length scale (chemical bonding, energies, and forces), but also at the nano and meso length scales (solid-state physics and material science processes). Current methods of study include DFT, molecular dynamics, and Monte Carlo. The limitation of DFT is that only a small number of atoms or molecules can be simulated at that level of theory. Density functional theory is used to study the electronic structure of atoms and molecules, and calculate the force component of each atom. However, these calculations are limited to about 1000 atoms. Custom periodic boundary conditions (PBC) can be used to describe the various environments and defects that affect the atomic forces to produce a large data set from which a training set can be derived. Machine learning can be utilized to overcome the barrier of modeling macroscopic and multi-scale processes from ab initio calculations through the development of adaptive force fields. Local environments determine the atomic forces of a given system, therefore adaptive force fields must be created to produce reliable quantum mechanical calculations. This can be achieved by developing a learning algorithm that uses the mapped atomic forces or fingerprint as an input to produce energies and magnetic moments as output. A systematic approach was used to begin to build a data set in order to accurately describe the atomic forces in various environments. In Figure 4 below, a simple PBC cell of Fe{sub 2}O{sub 3} was first optimized. A surface optimization was performed next, followed by a hydroxylated surface optimization. Once this calculation has converged, the adsorption of halide species to the hydroxylated surface will be investigated. TensorFlow is an open source platform for machine learning developed by Google. Using a high level application program interface (API) such as Keras allows for building and training ML models easily in a number of different environments and languages. For this project, a neural network was developed within Anaconda in Python. Future Work: Further development of reference data set; Refining neural network and learning algorithm; Fingerprinting atomic environment to enable mapping of atomic force components; Choosing appropriate training set from reference data; Learning from training set and enabling non-linear mapping of training set fingerprints and the atomic forces; Estimation of uncertainty to identify ranges of outside applicability; Testing and analysis of molecular dynamic simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Institute for Catalysis in Energy Processes (ICEP) (Final Report)

The Institute for Catalysis in Energy Processes (ICEP) was a multi-PI program located at the Northwestern University Center for Catalysis and Surface Science from 2006-2024. Over several renewal cycles and several organizing themes, ICEP addressed fundamental questions in catalysis science. In turn, the scientific questions addressed were directly relevant to the efficient use of the nation’s resources to produce fuels and commodity chemicals, harness alternate energy sources for chemical reactions such as light or electricity, reduce emissions and waste, and minimize the impact of our use of plastics. Selective oxidation, (oxidative) dehydrogenation, deNOx, CO 2 reduction, hydrogen release, polymer decomposition / depolymerization, and many other reactions of intense interest to the US Department of Energy were studied in the center. Major ICEP strengths were in catalyst synthesis, in measurements that elucidated catalyst properties, in reaction mechanisms and kinetics, and in predictive theory and modeling. ICEP researchers were often the inventors or developers of materials and methods that were brought to bear on the Institute’s catalytic systems. Key characterization tools advanced through this project included resonance Raman and related techniques, sum frequency generation, and X-ray standing wave spectroscopy. Center members relied extensively on computational tools such as density functional theory and microkinetic modeling to fully investigate catalyst structures and catalytic mechanisms. ICEP researchers were early developers of atomic layer deposition (ALD) for catalyst synthesis and were early pioneers of metal organic frameworks for catalysis. ICEP innovations are found in areas of intense research such as single-atom catalysts, catalytic deconstruction of plastics, and catalytic MOFs, to name a few. Research areas initiated by ICEP indirectly led to DOE EFRCs, startup companies, and other achievements. Over its 18 years of existence, ICEP involved two different PIs (Stair and Notestein) and 24 senior investigators across several academic disciplines at Northwestern University. There were tight collaborations with Argonne National Laboratory, including specialized equipment and experiments permanently located there. The center supported ~400 person-years of effort by graduate students and postdocs, either directly through DOE support or indirectly by leveraging independent support provided to the trainees, e.g. through the NSF GRFP or internal fellowships. The project resulted in 328 manuscripts and numerous conference presentations.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Spectroscopic Probe Molecule Selection Using Quantum Theory, First-Principles Calculations, and Machine Learning

Probe molecule vibrational spectra have a long history of being used to characterize materials including metals, oxides, metal-organic frameworks, and even human proteins. Furthermore, recent advances in machine learning have enabled computationally generated spectra to aid in detailed characterization of complex surfaces with probe molecules. Despite widespread use of probe molecules, the science of probe molecule selection is underdeveloped. Here, we develop physical concepts, including orbital interaction energy and the energy overlap integral, to explain and predict the ability of probe molecules to discriminate structural descriptors. We resolve the crystal orbital overlap population (COOP) to specific molecular orbitals and quantify their bonding character, which directly influences vibrational frequencies. Using only a single adsorbate calculation from density function theory (DFT), we compute the interaction energy of individual adsorbate molecular orbitals with adsorption site atomic orbitals across many different sites. Combining the molecular orbital resolved COOP and changes in orbital interaction energy enables probe molecule selection for improved discrimination of various sites. We demonstrate these concepts by comparing the predicted effectiveness of carbon monoxide (CO), nitric oxide (NO), and ethylene (C 2 H 4 ) to probe Pt adsorption sites. Finally, using a previously developed machine learning framework, we show that models trained on hundreds of thousand C 2 H 4 spectra, computed from DFT, which regress surface binding-type and generalized coordination number (GCN), outperform those trained using CO and NO spectra. Lastly, a python package, pDOS_overlap, for implementing the electron density based analysis on any combination of adsorbates and materials, is also made available.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2023

Predicting electronic structures at any length scale with machine learning

The properties of electrons in matter are of fundamental importance. They give rise to virtually all material properties and determine the physics at play in objects ranging from semiconductor devices to the interior of giant gas planets. Modeling and simulation of such diverse applications rely primarily on density functional theory (DFT), which has become the principal method for predicting the electronic structure of matter. While DFT calculations have proven to be very useful, their computational scaling limits them to small systems. We have developed a machine learning framework for predicting the electronic structure on any length scale. It shows up to three orders of magnitude speedup on systems where DFT is tractable and, more importantly, enables predictions on scales where DFT calculations are infeasible. Our work demonstrates how machine learning circumvents a long-standing computational bottleneck and advances materials science to frontiers intractable with any current solutions.

36 MATERIALS SCIENCE↗
DOE OSTI2020

B 2 H 6 splitting on catalytic surfaces and role of BH 3 towards hydrogen spillover

A fundamental understanding of the spillover mechanism is an open and challenging problem and plays an important role in catalysis. In particular, bond-exchange spillover mechanism is considered to be effective for reversible storage and release of hydrogen at near ambient conditions. For this, three critical steps are needed: finding the right support (acceptor), the right catalyst to split H 2 , and ensuring that once H 2 is split, the H atoms can migrate on the surface with the help of secondary catalysts and eventually hydrogenate the entire material. In this paper we address these challenges using density functional theory. Here we show that BH 3 , a secondary catalyst, can be produced by symmetrically splitting its stable precursor, B 2 H 6 , on doped metal-free surfaces such as graphene and h-BN as well as on MOF5. In addition, to reduce computational cost, we develop structural descriptor and predictive model equation to effectively screen potential BH 3 binding sites. Symmetrical splitting of B 2 H 6 on different types of materials can address the hydrogen spillover challenge, making efficient storage of hydrogen possible.

08 HYDROGEN↗
DOE OSTI2022

Te/CdTe and Al/CdTe Interfacial Energy Band Alignment by Atomistic Modeling

A synergistic approach that incorporates first-principles atomistic modeling with numerical device simulations is used to systematically evaluate the role of heterointerfaces within metal-chalcogenide-based photovoltaic technologies. Two interfaces involving either a tellurium back contact or aluminum back electrode combined with a cadmium telluride absorber layer within cadmium-telluride-based thin-film solar cells are investigated on an atomic scale to determine the mechanisms contributing to variations in device performance. Electronic structures and predicted charge transport behavior with respect to cadmium and tellurium termination of the absorber layer are studied along the polar oriented CdTe{111} facets. The computational methodology reveals a noticeable contrast between the Schottky barrier forming Al/CdTe interface versus the Type I Te/CdTe heterojunction. Greater band bending features are exhibited by the cadmium termination as opposed to the tellurium termination for each interface case. Subsequent device modeling suggests that 3.6% higher photovoltaic conversion efficiency is achievable for the cadmium termination relative to the tellurium termination of the Te/CdTe interface. Based strictly on an idealistic representation, both interface models show the importance of atomic-scale interfacial properties for cadmium telluride solar cell device performance with their bulk properties being validated in comparison to published experimental data. Here, the synergistic approach offers a suitable method to analyze solar cell interfaces through a predictive computational framework for the engineering and optimization of metal-chalcogenide-based thin-film photovoltaic technologies.

Green’s function↗
DOE OSTI2022

Computational Modeling of Battery Materials

This chapter provides a broad overview of the various computational modeling techniques used to gain fundamental insights into coupled electrochemical processes that occur in battery materials at electronic-to-mesoscopic scales. Furthermore, representative successes of these techniques in modeling electrodes, electrolytes, and electrode-electrolyte interfaces are highlighted to establish the current state-of-the-art in the field.

25 ENERGY STORAGE↗
DOE OSTI2020

Thermochemical Data Fusion Using Graph Representation Learning

Large databases are required for “Big Data” applications in catalysis and materials science. Thermochemical databases can be created by combining data from various sources and by correcting low-fidelity datasets to higher accuracy with minimal computation. To achieve this “data fusion”, thermochemical quantities of interest, calculated at various levels of density functional theory (DFT), need to be mapped to the same, high levels of theory. In this work, a graph theoretical, statistical framework is proposed for such tasks. Subgraph frequencies are shown to provide a natural representation for learning these fusion maps. The maps are linear and are learnt with automated descriptor selection. Using a dataset of as few as ~1% from the QM9 database of 133,885 molecules, these models can predict multiple thermochemical quantities at a higher level of theory with an accuracy of 1 kcal/mol. Here, the method is explainable, generalizable, and provides a diagnostic tool for outlier identification

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

Transferable predictions of energetic and structural properties for refractory solid solution alloys across chemical compositions

We present a data-efficient approach to train graph neural networks (GNNs) on density functional theory (DFT) data for accurate and transferable predictions of energetic and structural properties of refractory solid solution alloys in the niobium-tantalum-vanadium (Nb-Ta-V) chemical space. We start by training the GNN model only on DFT data that describes refractory binary alloys niobium-tantalum (Nb-Ta), niobium-vanadium (Nb-V), and tantalum-vanadium (Ta-V) to predict formation enthalpy and root mean squared displacement. Once trained, the GNN predictions are tested on DFT data describing refractory ternary alloys Nb-Ta-V. While, unsurprisingly, direct transferability from binary to ternary is not sufficiently accurate, augmenting the training with only 1% of the available ternary data (uniformly distributed across the entire range of chemical compositions) improves significantly the quality of the GNN predictions. For comparison, we assess the transferability in the opposite direction by training GNN models on ternary Nb-Ta-V data and making predictions on binaries Nb-Ta, Nb-V, and Ta-V, which exhibits notably higher predictive errors. The proposed methodology, which favors transferability from lower-component to higher-component alloys, offers an efficient path towards avoiding the curse of dimensionality incurred when collecting DFT data for discovery and design of multi-component disordered alloys.

Density functional theory calculations↗
DOE OSTI2023

Dataset of theoretical multinary perovskite oxides

Perovskite oxides (ternary chemical formula ABO 3 ) are a diverse class of materials with applications including heterogeneous catalysis, solid-oxide fuel cells, thermochemical conversion, and oxygen transport membranes. However, their multicomponent (chemical formula $A_xA^{'}_{1-x}B_yB^{'}_{1-y}O_3$) chemical space is underexplored due to the immense number of possible compositions. To expand the number of computed $A_xA^{'}_{1-x}B_yB^{'}_{1-y}O_3$ compounds we report a dataset of 66,516 theoretical multinary oxides, 59,708 of which are perovskites. First, 69,407 $A_{0.5}A^{'}_{0.5}B_{0.5}B^{'}_{0.5}O_3$ compositions were generated in the a - b + a - Glazer tilting mode using the computationally-inexpensive Structure Prediction and Diagnostic Software (SPuDS) program. Next, we optimized these structures with density functional theory (DFT) using parameters compatible with the Materials Project (MP) database. Our dataset contains these optimized structures and their formation (ΔH f ) and decomposition enthalpies (ΔH d ) computed relative to MP tabulated elemental references and competing phases, respectively. This dataset can be mined, used to train machine learning models, and rapidly and systematically expanded by optimizing more SPuDS-generated $A_{0.5}A^{'}_{0.5}B_{0.5}B^{'}_{0.5}O_3$ perovskite structures using MP-compatible DFT calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗