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 1,315 records · Page 73

Experimental and theoretical study of the soot-forming tendencies of furans as potential biofuels

Sustainably-produced biofuels, including furans derived from lignocellulosic biomass, are being studied as potential next-generation fuels. Biofuels have unique structures that are absent in conventional fuels and they potentially emit fewer particulates than conventional fuels because of the presence of oxygen atoms. In this work, we focus on quantification and measurement of sooting tendency. A set of furans was selected and their sooting tendencies were characterized by Yield Sooting Index (YSI), which was obtained by measuring the line-of-sight spectral radiance (LSSR) of a coflow laminar diffusion flame doped with 1000 ppm of each test fuel. The measured YSIs are lower than traditional gasoline fuels and benzenoid aromatics (e.g., YSIbenzene = 100). We also performed numerical simulations with a detailed chemical kinetic mechanism to analyze the flame chemistry of 2-methylfuran (2MF). 4→5 H-atom migration in the furan ring causes 31.0% of the total decomposition and is the dominant fuel consumption pathway leading to soot. DFT calculations via the G4 composite method were carried out to obtain activation energies (Ea) and bond dissociation energies (BDE). The results showed that the 4→5 H-atom migration in the furan ring has the second lowest Ea and subsequently forms a 1,2-diene-substituted ketone derivative as a stable species. The ring-opening products of the 4→5 H-atom migration are propyne + a linear ketene for 2-n-alkylfurans, and 1-butyne + a methyl-substituted ketene for 2-t-butylfuran. These results can explain the trends of the measured YSIs of 2-alkylfurans.

09 BIOMASS FUELS↗

Assessing the Role of Interfacial and Metal Sites in Pt/TiO2-Catalyzed Acetic Acid Hydrodeoxygenation

Thermochemical conversion of biomass to produce drop-in quality biofuels typically involves hydrodeoxygenation (HDO) steps following catalytic fast pyrolysis (CFP) to remove excess oxygen and create a more-stable bio-oil product. HDO reactions are performed by co-feeding the CFP vapor-phase product and H2 gas over a bi-functional catalyst. Noble-metal catalysts supported on reducible metal oxides (e.g., Pt/TiO2) are active and selective toward these HDO reactions. Griffin and co-workers showed that Pt/TiO2 catalysts promote the desired deoxygenation steps for m-cresol HDO while mitigating undesired C-C bond-breaking steps that reduce the overall value/energy density of the biofuel. Such model-compound studies to inform the design of improved catalysts for HDO chemistry are necessary to improve the overall process economics/efficiencies for biofuels production. One important class of bio-derived compounds that has not been studied extensively with respect to HDO chemistry, particularly at the atomic level, is carboxylic acids. Past experimental work indicates that Pt/C and Pt/TiO2 catalysts are selective toward C-C and C-O bond-dissociation products for acetic acid HDO, respectively. To understand the role of Pt/TiO2 active sites in this observed change in selectivity and guide catalyst development, the work in this presentation focuses on modeling the role of Pt-metal and Pt-TiO2-interfacial sites in promoting key C-C bond-breaking, C-O bond-breaking, and (de)hydrogenation steps in acetic acid HDO. Density functional theory (DFT) calculations were performed using the Vienna Ab initio Simulation Package (VASP). The exchange-correlation functional was approximated by the Perdew-Burke-Ernzerhof functional. Dispersion interactions were captured using the D3 method. Projector augmented-wave potentials described electron-ion interactions, and electron wavefunctions were expanded via a planewave basis with an energy cutoff of 500 eV. Activation barriers were calculated using the climbing image nudged elastic band method. Results and Discussion: To discern the role of Pt-metal and Pt-TiO2 interface sites in promoting acetic acid HDO chemistry, Pt(111) slab and anatase TiO2(101)-supported Pt-nanowire (PtNW/OH-TiO2) surface models were constructed, respectively. Because H2 is co-fed in HDO reactions, the anatase support was terminated with OH groups. Interfacial vacancies have been shown to facilitate Ru/TiO2-catalyzed phenol HDO; thus, an interfacial model with an OH vacancy was also considered (PtNW/OHv-TiO2). Pt-TiO2-interface sites stabilize adsorption of all studied acetic acid HDO surface intermediates relative to terrace Pt-metal sites, particularly when an interfacial-OH vacancy is present. Oxygenated species prefer to bind at the OH vacancy through the O atom, suggesting a preference for C-O over C-C bond cleavage at these sites. This hypothesis is supported by net-negative and net-positive average shifts in the reaction energy and activation energy barriers for C-O and C-C bond-breaking steps, respectively, at Pt-TiO2-interface sites relative to Pt-metal sites. Using the calculated energetics, the predicted minimum-energy pathway was determined for each surface. Pt(111) is predicted to follow decarboxylation, producing undesired methane and carbon dioxide. Conversely, PtNW/OH-TiO2 and PtNW/OHv-TiO2 are both predicted to produce desired acetaldehyde and ethane. The interfacial vacancy may also play a key role in facilitating the first C-O bond-breaking step in acetic acid HDO, lowering the barrier by 0.6 eV relative to the defect-free interface model. These results demonstrate the critical role of the Pt-TiO2 interface in the shift in acetic acid HDO selectivity experimentally observed on Pt/C and Pt/TiO2 catalysts. The results herein demonstrate the important role of the Pt-TiO2 interface and interfacial oxygen vacancies in improving the carbon efficiency for HDO reactions in CFP upgrading.

BIOMASS FUELS,INORGANIC, ORGANIC, PHYSICAL, AND AN↗

Extraction of Ln(III) and An(III) from Alkaline Solutions by Sulfonamide Ligands and Analogs - 20216

As part of our effort to develop advanced separations involving combined extraction of cesium, strontium, and actinides from High-Level Waste (HLW) at the Savannah River Site, we are now reporting highly efficient extraction of Sm(III) from aqueous media of various alkalinity by dichloromethane solutions of several ligands, including tri-sulfonamides, o-mono-sulfonamido-phenols and o-phenylenediamine- derived disulfonamides. Distribution experiments and DFT calculations provide insights on the complexation of Samarium vs Americium by these ligands. Further optimization of these ligands for extraction into highly lipophilic solvents and application to HLW simulants are underway. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Ba₁₋ₓSrₓFeO₃₋ᵧ as an Improved Oxygen Carrier for Chemical Looping Air Separation: A Computational and Experimental Study

Chemical looping air separation (CLAS) is a promising method to generate pure carbon-dioxide from fuel combustion with a pure oxygen stream, which is produced through the capture, and targeted release, of oxygen from the atmosphere using a solid oxide carrier. The performance of this process depends on the redox characteristics of the oxide carrier. Using experimental oxygen-temperature-programmed desorption (TPD) and thermogravimetric analysis (TGA), the study shows that Ba₀․₇₅Sr₀․₂₅FeO₃ has improved oxygen storage capacity (OSC), oxidation, and reduction kinetics over pristine SrFeO₃ at temperatures ranging from 300-500 °C. The redox energetics computed by using first-principles density functional theory (DFT) calculations also depict the measured trend establishing it as an important descriptor of the measured performance. The Ba₁₋ₓSrₓFeO₃₋ᵧ depicts a Sr-substitution and oxygen stoichiometry induced structural phase transition from hexagonal at low temperatures to pseudo-cubic phase with higher OSC up to T=400 °C. Also, the mechanism leading to this structural phase transition is identified by performing electronic and vibrational structural calculations. This presentation was given on March 8, 2024 at the APS March Meeting in Minneapolis, MN.

Acharya, Shree Ram↗

Efficient CO2 and Natural Gas Refinery for Ethylene Production - Catalyst Development for Oxidative Coupling of Methane (OCM)

We designed a new electro-catalyst for oxidative coupling of methane (OCM) reaction, the structure of the Sc doped La0.3Sr0.7TiO3-d (LaSrTi) Sc doping to LaSrTi leading to lattice expansion and more distorted lattice, Investigated the effect of doping on the oxygen vacancy formation energy DFT calculations shows that Sc doping facilitate oxygen vacancy formation. Quantifying the surface basicity of the sample before and after Sc doping Sc doped sample exhibits high surface basicity.

catalysis↗

Using Density Functional Theory (DFT) for the Calculation of Atomization Energies

The calculation of atomization energies using density functional theory (DFT), using the B3LYP hybrid functional, is reported. The sensitivity of the atomization energy to basis set is studied and compared with the coupled cluster singles and doubles approach with a perturbational estimate of the triples (CCSD(T)). Merging the B3LYP results with the G2(MP2) approach is also considered. It is found that replacing the geometry optimization and calculation of the zero-point energy by the analogous quantities computed using the B3LYP approach reduces the maximum error in the G2(MP2) approach. In addition to the 55 G2 atomization energies, some results for transition metal containing systems will also be presented.

Bauschlicher, Charles W., Jr.↗

Mechanisms of Ethylene Epoxidation over Silver from Machine Learning-Accelerated First-Principles Modeling and Microkinetic Simulations

In this work we employed machine learning-augmented density functional theory (DFT) thermodynamic calculations to assess the stability of different AgO x structures under catalytic ethylene epoxidation reaction conditions. We found that there are multiple AgO x surface motifs that could co-exist under the relevant conditions. These included Ag surface oxides (e.g., AgO_p(4 × 4) and Ag 1.83 O) and atomic oxygen-covered Ag(111) surfaces. Furthermore, we employed DFT calculations to evaluate the energetics of different reaction mechanisms by which ethylene and oxygen can react on these surfaces. These studies revealed several energetically viable reaction pathways for ethylene epoxidation. Microkinetic modeling analyses, based on the DFT-calculated reaction pathways, showed that ethylene epoxidation can proceed on all surfaces and that multiple pathways, including those involving Langmuir–Hinshelwood and Eley–Rideal mechanisms, could be involved in selective and unselective reactions. The diversity of mechanisms that we discovered in the context of the relatively simple ethylene epoxidation reaction on Ag suggests that the richness and complexity of surface chemistry are most likely a rule rather than an exception in heterogeneous catalytic chemical transformations on metal surfaces and that the concept of a single or even a dominant mechanism and reaction intermediates might need to be revisited for many reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Boosting the activity of transition metal carbides towards methane activation by nanostructuring

The interaction of methane with pristine surfaces of bulk MoC and Mo 2 C is known to be weak. In contrast, a series of X-ray photoelectron spectroscopy (XPS) experiments, combined with thermal desorption mass spectroscopy (TDS), for MoC y (y = 0.5–1.3) nanoparticles supported on Au(111)—which is completely inert towards CH 4 —show that these systems adsorb and dissociate CH 4 at room temperature and low CH 4 partial pressure. This industrially-relevant finding has been further investigated with accurate density functional theory (DFT) based calculations on a variety of MoC y supported model systems. The DFT calculations reveal that the MoC y /Au(111) systems can feature low C–H bond scission energy barriers, smaller than the CH 4 adsorption energy. Our theoretical results for bulk surfaces of Mo 2 C and MoC show that a simple Brønsted–Evans–Polanyi (BEP) relationship holds for C–H bond scission on these systems. However, this is not the case for methane activation on the MoC y nanoparticles as a consequence of their unique electronic and chemical properties. Finally, the discovery that supported molybdenum carbide nanoparticles are able to activate methane at room temperature paves the road towards the design of a new family of active carbide catalysts for methane activation and valorisation, with important implications in climate change mitigation and carbon cycle closure.

36 MATERIALS SCIENCE↗

Using data-science approaches to unravel insights for enhanced transport of lithium ions in single-ion conducting polymer electrolyte

Solid polymer electrolytes have yet to achieve the an ionic conductivity > 1 mS/cm at room temperature for realistic applications. This target implies the need to reduce the effective energy barriers of ion transport in polymer electrolytes to around 20 kJ/mol. In this work, we combine information extracted from existing experimental results with theoretical calculations to provide insights into ion transport in single-ion conductors (SICs) with a focus on lithium ion SICs. Through the analysis of temperature-dependent ionic conductivity data obtained from the literature, we evaluate different methods of extracting energy barriers for lithium transport. The traditional Arrhenius fit to the temperature-dependent ionic conductivity data indicates that the Meyer-Neldel rule holds for SICs. However, the values of the fitting parameters remain unphysical. Our modified approach based on recent work (Macromolecules, 56, 15, 6051(2023)), which incorporates a fixed pre-exponential factor, reveals that the energy barriers exhibit temperature dependence over a wide range of temperatures. Using this approach, we identify a series of anions leading to the energy barriers less than 30 kJ/mol, which include trifluoromethane sulfonimide (TFSI), fluoromethane sulfonimide (FSI), and boron-based organic anions. In our efforts to design the next generation of anions, which can exhibit the energy barriers less than 20 kJ/mol, we focused on boron-containing SICs, and performed density functional theory (DFT) based calculations to connect the chemical structures via the binding energy of cation (lithium)-anion pairs with the experimentally derived effective energy barriers for ion transport. Not only have we identified a correlation between the binding energy and the energy barriers, but we also propose a strategy to design new boron-based anions by using the correlation. This combined approach involving experiments and theoretical calculations is capable of facilitating the identification of promising new anions, which can exhibit ionic conductivity $> 1$ mS/cm near room temperature, thereby expediting the development of novel superionic single-ion conducting polymer electrolytes. The published datasets include all the temperature-dependent ionic conductivity collected from the literature with literature DOIs, DFT calculated binding energies, and python scripts to analyze data, construct statistical models, and generate plots.

36 MATERIALS SCIENCE↗

OrbNet Denali: A machine learning potential for biological and organic chemistry with semi-empirical cost and DFT accuracy

We present OrbNet Denali, a machine learning model for an electronic structure that is designed as a drop-in replacement for groundstate density functional theory (DFT) energy calculations. The model is a message-passing graph neural network that uses symmetryadapted atomic orbital features from a low-cost quantum calculation to predict the energy of a molecule. OrbNet Denali is trained on a vast dataset of 2.3 × 10 6 DFT calculations on molecules and geometries. This dataset covers the most common elements in biochemistry and organic chemistry (H, Li, B, C, N, O, F, Na, Mg, Si, P, S, Cl, K, Ca, Br, and I) and charged molecules. OrbNet Denali is demonstrated on several well-established benchmark datasets, and we find that it provides accuracy that is on par with modern DFT methods while offering a speedup of up to three orders of magnitude. For the GMTKN55 benchmark set, OrbNet Denali achieves WTMAD-1 and WTMAD-2 scores of 7.19 and 9.84, on par with modern DFT functionals. For several GMTKN55 subsets, which contain chemical problems that are not present in the training set, OrbNet Denali produces a mean absolute error comparable to those of DFT methods. For the Hutchison conformer benchmark set, OrbNet Denali has a median correlation coefficient of R 2 = 0.90 compared to the reference DLPNO-CCSD(T) calculation and R 2 = 0.97 compared to the method used to generate the training data (ωB97X-D3/def2- TZVP), exceeding the performance of any other method with a similar cost. Similarly, the model reaches chemical accuracy for noncovalent interactions in the S66x10 dataset. For torsional profiles, OrbNet Denali reproduces the torsion profiles of ωB97X-D3/def2-TZVP with an average mean absolute error of 0.12 kcal/mol for the potential energy surfaces of the diverse fragments in the TorsionNet500 dataset.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Gapless spinons and a field-induced soliton gap in the hyperhoneycomb Cu oxalate framework compound [(C 2 H 5 ) 3 NH] 2 Cu 2 (C 2 O 4 ) 3

Here, we report a detailed study of the specific heat and magnetic susceptibility of single crystals of a spin-liquid candidate: the hyperhoneycomb Cu oxalate framework compound [(C 2 H 5 ) 3 NH] 2 Cu 2 (C 2 O 4 ) 3 . The specific heat shows no anomaly associated with a magnetic transition at low temperatures down to T ~ 180 mK in zero magnetic field. We observe a large linear-in-T contribution to the specific heat γT, γ = 98 (1) mJ / mol K 2 , at low temperatures, indicative of the presence of fermionic excitations despite the Mott insulating state. The low-T specific heat is strongly suppressed by applied magnetic fields H, which induce an energy gap, Δ(H), in the spin-excitation spectrum. We use the four-component relativistic density-functional theory (DFT) to calculate the magnetic interactions, including the Dzyaloshinskii-Moriya antisymmetric exchange, which causes an effective staggered field acting on one copper sublattice. The magnitude and field dependence of the field-induced gap, Δ (H) ∝ H 2/3 , are accurately predicted by the soliton mass calculated from the sine-Gordon model of weakly coupled antiferromagnetic Heisenberg chains with all parameters determined by our DFT calculations. Thus our experiment and calculations are entirely consistent with a model of [(C 2 H 5 ) 3 NH] 2 Cu 2 (C 2 O 4 ) 3 in which anisotropic magnetic exchange interactions due to Jahn-Teller distortion cause one copper sublattice to dimerize, leaving a second sublattice of weakly coupled antiferromagnetic chains. We also show that this model quantitatively accounts for the measured temperature-dependent magnetic susceptibility. Thus [(C 2 H 5 ) 3 NH] 2 Cu 2 (C 2 O 4 ) 3 is a canonical example of a one-dimensional spin-1/2 Heisenberg antiferromagnet and not a resonating-valence-bond quantum spin liquid, as previously proposed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Engineering metal–metal oxide surfaces for high-performance oxygen reduction on Ag–Mn electrocatalysts

Understanding fundamental material–property relationships in mixed-element catalyst systems is crucial to advancing the viability of renewable electrochemical energy technologies, an important part of creating a more sustainable future. Herein, we report our insight on the nature and dynamics of highly active silver–manganese oxide (Ag–MnOx) catalyst surfaces for the oxygen reduction reaction (ORR) via a combined experimental–theoretical approach. Experimentally, we synthesize well-mixed Ag–Mn co-deposited thin films that are measurably flat and smooth, despite Mn surface migration and oxidation upon air exposure and electrochemical measurements. Cyclic voltammetry in 0.1 M KOH demonstrates up to 10-fold specific activity enhancements over pure Ag at 0.8 V vs. RHE for Ag-rich films (70–95% Ag in bulk). To further understand the Ag–Mn system, separate samples were synthesized with small amounts of Mn sequentially deposited onto the surface of a pure Ag thin film (Mn@Ag), ranging from partial to full surface coverage (down to 0.3 nm Mn $cm^{-2}_{geo}$ ~ 0.2 μg Mn $cm^{-2}_{geo}$). These sequentially deposited Mn@Ag films show analogous performance to their co-deposited counterparts indicating similar enhanced active sites. With density functional theory (DFT), we calculate that this enhancement arises from the tuned d-band of these material surfaces owing to the optimal hybridization of the electronic structures in specific Ag and MnO x geometries. Together, electrochemical measurements, DFT calculations, X-ray absorption spectroscopy, and valence-band X-ray photoelectron spectroscopy suggest synergistic electronic interactions between Ag and MnOx yield enhanced oxygen adsorption, and thus ORR activity, with DFT highlighting the Ag–MnO x interface sites as the most enhanced. This work demonstrates how combined experimental–theoretical methods can help design electrocatalysts with enhanced electrocatalytic properties and understand the nature of complex mixed metal–metal oxide surfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dendrite Growth Morphology Modeling in Liquid and Solid Electrolytes

The main goal of this project is to develop a multi-scale modeling approach that connects micron-scale phase-field models and atomic-scale density functional theory (DFT)-based simulations via parameter- and relationship-passing in order to predict Li-metal dendrite morphology evolution, in both liquid and solid electrolytes. The key hypothesis of the DFT-informed phase-field multiscale modeling approach is that it can capture the electrochemical-mechanical driving forces and incorporate the roles of nano-meter-thin solid electrolyte interphase (SEI) in liquid electrolytes as well as of the microstructures of micro-meter-thick solid electrolytes (SEs) for all-solid-state batteries. In this project, we have formulated and implemented phase-field models to incorporate the electrochemical driving forces in liquid electrolytes and then incorporate mechanical driving forces to simulate dendrite growth in solid electrolytes with resolved microstructures. We have implemented two treatments for the SEI: an explicit model to include the microstructure of the SE or SEI in the phase field model and an implicit model to simulate the impact of nano-meter thick SEI in liquid electrolytes by varying the electrode/electrolyte interfacial properties. The key interfacial properties, including the electronic and ionic transport properties, the charge transfer reaction kinetics, and mechanical properties, were computed by DFT-based calculations. At the DFT-based model, one key advancement is to directly predict the charge transfer reaction kinetics at a complex Li/SEI/electrolyte interface by linking DFT with density functional tight binding (DFTB) calculations. As the main accomplishments, we have demonstrated two successful predictions in both solid electrolyte and liquid electrolyte based on this multiscale approach. The predicted intergranular Li dendrite growth in LLZO revealed the importance of trapped electrons at internal interfaces in the microstructure of LLZO. The predicted electroplating morphology of mossy Li and faceted Mg agreed well with experiments. The insights provided by the multiscale model and the model enabled electrolyte and SEI design will accelerate the development of Li-metal electrode for high energy density batteries, that meet DOE’s target on cell density (>350 Wh/kg) and cost below $100/kWhuse for EV applications.

25 ENERGY STORAGE↗

Time-dependent density-functional-theory calculations of the nonlocal electron stopping range for inertial confinement fusion applications

Nonlocal electron transport is important for understanding laser-target coupling for laser-direct-drive (LDD) inertial confinement fusion (ICF) simulations. Current models for the nonlocal electron mean free path in radiation-hydrodynamic codes are based on plasma-physics models developed decades ago; improvements are needed to accurately predict the electron conduction in LDD simulations of ICF target implosions. Here we utilized time-dependent density functional theory (TD-DFT) to calculate the electron stopping power (SP) in the so-called conduction-zone plasmas of polystyrene in a wide range of densities and temperatures relevant to LDD. Compared with the modified Lee-More model, the TD-DFT calculations indicated a lower SP and a higher stopping range for nonlocal electrons. We fit these electron SP calculations to obtain a global analytical model for the electron stopping range as a function of plasma conditions and the nonlocal electron kinetic energy. This model was implemented in the one-dimensional radiation-hydrodynamic code LILAC to perform simulations of LDD ICF implosions, which are further compared with simulations by the standard modified Lee-More model. In conclusion, results from these integrated simulations are discussed in terms of the implications of this TD-DFT-based mean-free-path model to ICF simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Modeling separation of lanthanides via heterogeneous ligand binding

Individual lanthanide elements have physical/electronic/magnetic properties that make each useful for specific applications. Several of the lanthanides cations (Ln 3+ ) naturally occur together in the same ores. They are notoriously difficult to separate from each other due to their chemical similarity. Predicting the Ln 3+ differential binding energies (ΔΔE) or free energies (ΔΔG) at different binding sites, which are key figures of merit for separation applications, will help design of materials with lanthanide selectivity. We apply ab initio molecular dynamics (AIMD) simulations and density functional theory (DFT) to calculate ΔΔG for Ln 3+ coordinated to ligands in water and embedded in metal–organic frameworks (MOFs), and ΔΔE for Ln 3+ bonded to functionalized silica surfaces, thus circumventing the need for the computational costly absolute binding (free) energies ΔG and ΔE. Perturbative AIMD simulations of water-inundated simulation cells are applied to examine the selectivity of ligands towards adjacent Ln 3+ in the periodic table. Static DFT calculations with a full Ln 3+ first coordination shell, while less rigorous, show that all ligands examined with net negative charges are more selective towards the heavier lanthanides than a charge-neutral coordination shell made up of water molecules. Amine groups are predicted to be poor ligands for lanthanide-binding. Finally, we also address cooperative ion binding, i.e., using different ligands in concert to enhance lanthanide selectivity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Valence transition and termination-dependent surface states in the topological Kondo semimetal YbPtBi

Here, the Yb-terminated and Bi-terminated (111) surface electronic structure of topological Kondo semimetal YbPtBi is investigated using both density-functional-theory (DFT) -based calculations and angle-resolved photoemission spectroscopy (ARPES). The cleavage plane is found to be between Yb-layers and Bi-layers in both experiment and theory, and the broken inversion symmetry ensures the Yb-terminated surface does not mix with the Bi-terminated surface. The ARPES results at the Bi-terminated surface are similar to those reported for YPtBi or LuPtBi, and they can be well explained using DFT calculations assuming trivalent Yb atoms. In contrast, at the Yb-terminated surface, a trivalent to divalent transition and reduced hybridization with conduction electrons for the topmost Yb atoms are observed, as a result of reduced bonding with Bi-atoms. In addition, the Fermi arc features induced by the triply degenerate-points, which are missing at the Bi-terminated surface, can be identified at the Yb-terminated surface. Our study demonstrates unambiguously the reconstruction of surface electronic structure due to the valence transition of the lanthanide element in a Kondo lattice system.

36 MATERIALS SCIENCE↗

A Multi-Scale Computational Platform for Predictive Modeling of Corrosion in Al-Steel Joints (Final Report)

The research team proposed to develop innovative multi-scale models to predict corrosion and the resulting mechanical performances in aluminum-steel joints. The methods of joining considered are resistance spot welding, self-piercing riveting, and rivet-welding, all suitable for mass production applications. The multi-scale models integrate high throughput first-principle calculations based on density functional theory (DFT), high throughput calculation of phase diagrams (CALPHAD) modeling, and finite element method (FEM) simulations. These models are to be validated through laboratory experiments. Furthermore, the models are available as open source so as to enable scientists and engineers in the community to adapt and contribute to the development and application. The approaches rely on the research team’s extensive experience on the prediction of properties of individual phases at finite temperatures and variable compositions through DFT calculations, and our broad expertise on dissimilar material joining and their corrosion. The proposed computational framework enables high throughput computations for improved predictions of corrosion and the associated mechanical performance in dissimilar material joints, resulting in significant reduction in computational time needed by the current state-of-the-art methods. With the participation of researchers from three universities, an auto manufacturer, two manufacturing technology/equipment suppliers, and a software developer/vendor, the interdisciplinary research team applies the technical development on both phase-based modeling and laboratory experiments into the automobile body joining processes for validation and technology demonstration. The global cost of corrosion was estimated at about 3.4% of the global GDP in 2013. By using available corrosion control practices, it is estimated a saving between 15-35% of the cost of corrosion. In the U.S., more than $276 billion is spent repairing corrosion damage. Prediction of the corrosion and its impact on performance of the dissimilar material joints is critical for reducing the massive number of the current corrosion-based recalls for automobiles. Thus, the project goal is to develop models to enable predictive maintenance and end-of-life planning of multi-metal joints with risk of corrosion under different conditions such as exposure to high temperatures in summer and salt solutions in winter, quantified through its pH. An academia-industry consortium led by the University of Michigan and including Pennsylvania State University, University of Illinois Urbana-Champaign, University of Georgia, General Motors Company, Livermore Software Technology Corporation, and Optimal Process Technologies, LLC. created multi-scale models for prediction of corrosion in aluminum-steel joint structures such of them used in vehicle subassemblies – chassis and transmission systems. Starting from the first principle calculations, the team developed mathematical and data-driven models to predict the metallic components, which are formed during joining of two metals, for example aluminum and steel - a lightweight multilateral system which is currently used in more than 60% car bodies. These models were used for simulating chemical reactions that are happening when the joining metallic components are exposed to high temperatures and different pH values. The team was able to predict how the corrosion installs on the metallic components and how they lead to a sudden failure of components in cars. Newly developed machine learning algorithms combining Science, Technology, Engineering and Math disciplines, advanced finite element simulation and experimental validations have been integrated in a platform for prediction of the corrosion evolution and prediction the failure of joints under mechanical loadings and fatigue. Moreover, based on machine learning and inverse analysis, the team proposed solutions for designing new metallic alloys less susceptible to corrosion when joining multi-material assembles. An average of 4% error compared with experiments was achieved for the most common joints that are used in vehicle subassemblies.

36 MATERIALS SCIENCE↗

Thermodynamic re-modeling of the Yb-Sb system aided by first-principles calculations

Here, the thermodynamic description of the Yb-Sb binary system is developed by means of the CALculations of PHAse Diagrams (CALPHAD) method by combining experimental data in the literature and predictions from first-principles calculations based on density functional theory (DFT) in the literature and the present work. Two pseudopotentials of Yb are compared in the present DFT-based calculations with 14 and 13 f-electrons frozen in the core, i.e., 5p 6 6s 2 and 5p 6 6s 2 5 d 1 electrons as valence electrons, termed Yb_2 and Yb_3, respectively. It is shown that the phonon spectrum of the YbSb phase calculated using the Yb_3 pseudopotential does not have imaginary phonon modes and is subsequently used to predict its temperature dependent thermodynamic properties by the DFT-based quasiharmonic phonon calculations. The present thermodynamic database includes the Yb 16 Sb 11 phase in addition to five intermetallic phases that were considered in previous modeling studies, i.e., YbSb 2 , YbSb, Yb 11 Sb 10 , Yb 4 Sb 3 , and Yb 5 Sb 3 . The high temperature orthorhombic structure of the Yb 5 Sb 3 phase is not considered in the present work as it was stabilized by hydrogen. The associate solution model is used to describe the short-range ordering behavior in the liquid phase. The calculations from the present thermodynamic model show good agreement with thermochemical and phase equilibrium data from both the present work and the literature.

36 MATERIALS SCIENCE↗