Search NASA⌕ Search

SEARCH · Search NASA

Results for “density functional theory (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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Data-Driven and Multiscale Modeling of DNA-Templated Dye Aggregates

Dye aggregates are of interest for excitonic applications, including biomedical imaging, organic photovoltaics, and quantum information systems. Dyes with large transition dipole moments (μ) are necessary to optimize coupling within dye aggregates. Extinction coefficients (ε) can be used to determine the μ of dyes, and so dyes with a large ε (>150,000 M –1 cm –1 ) should be engineered or identified. However, dye properties leading to a large ε are not fully understood, and low-throughput methods of dye screening, such as experimental measurements or density functional theory (DFT) calculations, can be time-consuming. In order to screen large datasets of molecules for desirable properties (i.e., large ε and μ), a computational workflow was established using machine learning (ML), DFT, time-dependent (TD-) DFT, and molecular dynamics (MD). ML models were developed through training and validation on a dataset of 8802 dyes using structural features. A Classifier was developed with an accuracy of 97% and a Regressor was constructed with an R2 of above 0.9, comparing between experiment and ML prediction. Using the Regressor, the ε values of over 18,000 dyes were predicted. The top 100 dyes were further screened using DFT and TD-DFT to identify 15 dyes with a μ relative to a reference dye, pentamethine indocyanine dye Cy5. Two benchmark MD simulations were performed on Cy5 and Cy5.5 dimers, and it was found that MD could accurately capture experimental results. The results of this study exhibit that our computational workflow for identifying dyes with a large μ for excitonic applications is effective and can be used as a tool to develop new dyes for excitonic applications.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

Machine Learning-Assisted Recovery of Delicate Kinetic Information from Transient Reactor Experiments

Identifying active sites and their roles in chemical reaction steps remains a vital challenge in heterogeneous catalysis. Transient experiments offer a unique way to probe active sites and distinguish subtle kinetic features. Although physics-based analysis methods may be well-developed, they can be highly susceptible to experimental noise, and smoothing methods may erase or even distort important features; a smooth curve is not always the best curve. We demonstrate a new workflow for the direct interpretation of intrinsic kinetic information from exit flux curves measured in transient reactor experiments. This workflow contains three artificial neural networks (ANNs), including a noise reducer, a concentration predictor, and a rate predictor to analyze experimental data, followed by the virtual TAP (VTAP) physics-based reactor model and density functional theory (DFT) calculations of adsorption energies on specific sites. We use this workflow to analyze the data from experiments titrating Pt/Al 2 O 3 and Pt/SiO 2 catalysts with carbon monoxide (CO) in the temporal analysis of products (TAP) reactor. Our workflow separates the time-evolving chemical reaction and mass transfer information contained in the TAP pulse response. The existence of strong- and weak-binding sites on the Pt/Al 2 O 3 catalyst is observed in the catalyst titration experiment in the transient reactor. The structures of the strong- and weak-binding sites are then identified by using DFT calculations. We find that the Pt/SiO 2 catalyst has only strong-binding sites, which aligns with the inactive support effect of SiO 2 . We demonstrate how machine learning methods provide unique insights with high-resolution data analysis that cannot be achieved by using state-of-the-art physics-based methods.

Adsorption↗
DOE OSTI2021

Flat-Band-Enabled Triplet Excitonic Insulator in a Diatomic Kagome Lattice

The excitonic insulator (EI) state is a strongly correlated many-body ground state, arising from an instability in the band structure of narrow-gap semiconductors towards exciton formation. Here we show that the flat valence and conduction bands of a Yin-Yang Kagome lattice, as exemplified in a superatomic graphene lattice, conspire to enable an interesting state of triplet EI, based on first-principles Density Functional Theory (DFT) calculations combined with many-body GW and Bethe-Salpeter Equation (BSE). As an intrinsic property of topological flat bands, highly localized electron and hole wavefunctions significantly reduce the screening and enhance the exchange interaction, leading to an unusually high triplet exciton binding energy (~1.2 eV) exceeding the GW band gap by ~0.2 eV and a huge singlet-triplet splitting of ~0.4 eV. Here, the flat-bands-enabled triplet EI state also points to the possibility of complete population inversion between the two topological flat bands for the realization of excited quantum spin Hall effect. Our findings enrich once again the intriguing physics of flat bands, which has drawn broad interest.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Predicting Partial Atomic Charges in Metal–Organic Frameworks: An Extension to Ionic MOFs

Molecular simulation is an invaluable tool to predict and understand the usage of metal–organic frameworks (MOFs) for gas storage and separation applications. Accurate partial atomic charges, commonly obtained from density functional theory (DFT) calculations, are often required to model the electrostatic interactions between the MOF and adsorbates, especially when the adsorbates have dipole or quadrupole moments, such as water and CO 2 . Machine learning (ML) models have been previously employed to predict partial charges and avoid the computational cost associated with DFT calculations. However, previous ML models suffer from small training data sets, which limit their scope of application. In this work, we introduce two novel machine learning models, PACMOF2-neutral and PACMOF2-ionic, aimed at predicting the density-derived electrostatic and chemical (DDEC6) partial atomic charges for both neutral and ionic MOFs. These models not only yield DFT-level accuracy at a fraction of the computational cost but also demonstrate a remarkable improvement in prediction of adsorption, as validated with grand canonical Monte Carlo simulations. Furthermore, the robustness and fast computational time of the PACMOF2 models, along with their transferability to other porous materials such as covalent organic frameworks and zeolites, underscores their potential in high-throughput screening of MOFs for diverse applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A first principles analysis of potential-dependent structural evolution of active sites in Fe-N-C catalysts

Fe-N-C (iron–nitrogen–carbon) electrocatalysts have emerged as potential alternatives to precious metal-based materials for the oxygen reduction reaction (ORR). However, the structure of these materials under electrochemical conditions is not well understood, and their poor stability in acidic environments poses a formidable challenge for successful adoption in commercial fuel cells. To provide molecular-level insights into these complex phenomena, we combine periodic density functional theory (DFT) calculations, exhaustive treatment of coadsorption effects for ORR reaction intermediates, including O and OH, and comprehensive analysis of solvation stabilization effects to construct voltage-dependent ab initio thermodynamic phase diagrams that describe the in situ structure of the active sites. These structures are further linked to activity and stability descriptors that can be compared with experimental parameters such as the half-wave potential for ORR and the onset potential for carbon corrosion and CO 2 evolution. The results indicate that pyridinic Fe sites at zigzag carbon edges, as well as other edge sites, exhibit high activity for ORR compared to sites in the bulk. However, edges neighboring the active sites are prone to instability via overoxidation and consequent site loss. The results suggest that it could be beneficial to synthesize Fe-N-C catalysts with small sizes and large perimeter edge lengths to enhance ORR activity, while voltage fluctuations should be limited during fuel cell operation to prevent carbon corrosion of overoxidized edges.

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

Data-driven search for promising intercalating ions and layered materials for metal-ion batteries

The rise in demand for lithium-ion batteries has led to a large-scale search for electrode materials and intercalating ion species to meet the demands of next-generation energy technologies. Recent efforts largely focus on searching for cathodes that can accommodate large amounts of intercalating ions, but similar work on anodes is relatively limited. This study utilizes machine learning methods to find alternative two-dimensional (2D) materials and intercalating ions beyond Li for metal-ion batteries with high-power efficiencies. The approach first uses density functional theory (DFT) calculations to estimate the theoretical capacities and voltages of various metal ions on 2D materials. The DFT-generated data also provide insights into the local structural accommodation upon ion intercalation on various 2D materials. Significant changes to the lattice can result in irreversible changes to the bonding environments in the anode material, resulting in poor cycling stability. Next, this study develops a binding energy and structural accommodation-based classification model to screen anode materials for next-generation batteries. The classification model selects intercalating ions and 2D material pairs suitable for batteries based on the calculated voltage and volumetric changes in the 2D material upon intercalation. Finally, this study builds a regression model to accurately predict the binding energies of the various intercalating ions on 2D materials. The approach highlights the importance of different elemental and structural features for classification and regression tasks. In conclusion, the insights gained from this study on the role of involved features, such as electronegativities of the constituent ions and the presence of unfilled electronic levels, will help to streamline further studies towards the search for future layered battery materials.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Equation of state for Hf, Ta, W, Re, Os, Ir, Pt, and Au to multi-terapascal pressures from density-functional theory

We present the zero-temperature equation of state (pressure dependence of compression) and phase stability predictions for the 5d-transition metals obtained from all-electron density-functional theory (DFT) calculations. The results compare favorably with experiments but extend beyond current experimental capabilities to 10 TPa. Our study reveals phase changes that are explained from the calculated electronic structure. The cubic face-centered and body-centered structures (fcc and bcc), together with two-, three-, and four-layered hexagonal structures, play major roles under compression. The results’ dependence on the electron exchange and correlation in the DFT approach is investigated, and it is shown that the impact of the choice, while significant at lower pressures, diminishes in the terapascal regime. We further illustrate that the normal parabolic trends in atomic volume and bulk modulus with atomic number, due to the occupation of bonding and anti-bonding 5d states, break down at TPa pressures, suggesting drastically different chemical bonding at these extreme conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Unraveling the CO Oxidation Mechanism over Highly Dispersed Pt Single Atom on Anatase TiO 2 (101)

Catalysts with noble metals deposited as single atoms on metal oxide supports have recently been studied extensively due to their maximized metal utilization and potential for performing difficult chemical conversions owing to their unique electronic properties. Understanding of the reaction mechanisms on supported single-metal atoms is still limited but is highly important for designing more efficient catalysts. In this study, we report the complexity of the CO oxidation reaction mechanism on Pt single atoms supported on anatase TiO 2 (Pt SA /a-TiO 2 ) by coupling density functional theory (DFT) calculations and microkinetic analysis with kinetic measurements, in situ/operando infrared, and X-ray absorption spectroscopies. Starting from the adsorbed Pt SA occupying an O vacancy induced by reductive pretreatment, we show that CO oxidation follows a complex mechanism consisting of initiation steps to reorganize the active site and multibranch reactive cycles, with the Pt SA /a-TiO 2 catalyst not returning to its initial configuration. The initiation step consists of CO and O 2 adsorption healing the O vacancy, followed by CO oxidation using gas-phase CO to form Pt(CO). The reactive cycle alternates O 2 adsorption and dissociation to oxidize the catalyst to Pt(O)(O)(CO) and branching pathways of competing Langmuir–Hinshelwood (LH)- or Eley–Rideal (ER)-type CO oxidation steps to reduce it again to Pt(CO). In situ/operando infrared experiments, including cryogenic CO adsorption and isotopic CO exchange, confirm the combined involvement of strongly adsorbed CO and gas-phase CO in an Eley–Rideal step along the reaction cycle. Microkinetic modeling shows that Pt single atoms are present in a mixture of Pt(CO), Pt(CO)(O 2 ), Pt(O)(CO)(O 2 ), and Pt(CO)(CO 3 ) structures as the main intermediates during steady-state CO oxidation, all having the C–O vibrational stretch close to the experimentally observed value of 2115 cm –1 . Microkinetic modeling also shows that the fractional orders of CO and O 2 measured experimentally originate from multiple steps with a high degree of rate control and not from a simple competitive adsorption. The results demonstrate the complex reaction pathways that even CO oxidation on a simple single-atom system can follow, providing mechanistic insights for designing efficient Pt-based single-atom catalysts. Here, we further show that microkinetic modeling results are sensitive to changes in energies of intermediate and transition states within errors of density functional theory, which can ultimately lead to incorrect conclusions regarding the reaction pathways and most abundant reaction intermediates if not accounted for by experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Odd-Parity Magnetism Driven by Antiferromagnetic Exchange

Realizing odd-parity, time-reversal-preserving, nonrelativistic spin splitting is a central goal for spintronics applications. We propose a group-theory-based microscopic framework to induce odd-parity spin splitting from coplanar antiferromagnetic (AFM) states without spin-orbit coupling (SOC). We develop phenomenological models for 421 conventional period-doubling AFM systems in nonsymmorphic space groups and construct minimal microscopic models for 119 of these. We find that these AFM states can attain three possible competing ground states. These ground states all break symmetries in addition to those broken by the usual AFM order. Specifically, they give rise to either odd-parity spin-splitting, nematic order, or scalar odd-parity order related to multiferroicity. Our microscopic theories reveal that the odd-parity spin-splitting energy scale is generically large and further reveal that the scalar odd-parity order gives a nonzero Berry curvature dipole without SOC. We identify 67 materials in the Magndata database for which our theory applies. We provide density-functional theory (DFT) calculations on Fe-based materials that reveal an ℎ-wave spin splitting consistent with our symmetry arguments and apply our microscopic model to determine the nonrelativistic Edelstein response for CeNiAsO.

Antiferromagnets↗
DOE OSTI2025

Pathways of Photocatalytic Oxidation of Formic Acid on Dry and Hydrated Anatase TiO 2 Surfaces

The photocatalytic oxidation of formic acid (FA), which is one of the most abundant volatile organic compounds, is a promising air remediation technology inspired by nature. However, the detailed mechanism of this photocatalytic reaction on the surface of TiO 2 , a typical photocatalyst, is not yet well-understood. In this work, we present a computational mechanistic study of the thermal vs photocatalytic oxidation of FA on dry and hydrated anatase TiO 2 (101) surfaces, based on periodic hybrid density functional theory (DFT) calculations, in which the photooxidation is treated as an excited-state process in a constrained triplet spin state. We first compare the adsorption modes of FA on the anatase (101) surface in the ground and excited states, followed by identification of the corresponding reaction intermediates that lead to the formation of CO 2 . We unveil the pivotal role of photogenerated holes localized at surface under-coordinated oxygen sites in mediating the C−H bond cleavage, thereby promoting CO 2 formation through a highly stable intermediate and an exergonic reaction step. Further investigation of the effect of coadsorbed water molecules shows that hydrogen bonding with water stabilizes FA in a monodentate configuration. This is favored over the unreactive bidentate structure that is the most stable under dry conditions, thus providing insight into the experimentally observed increase of the reaction rate in the presence of water.

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

Sub-Nanometer Nanoclusters of Copper Atop Single-Atom Copper Moieties toward Electrochemical CO 2 Hydrogenation to Methane

The electrochemical CO 2 reduction (eCO 2 R) offers a compelling route for converting CO 2 into value-added fuels and chemicals. Among CO 2 -derived products, methane (CH 4 ) occupies a distinct position, serving both as a key intermediate for emerging cascade electro-oxidation to oxygenates and as a strategically important extraterrestrial fuel that can be generated in situ from off-planet CO 2 resources. Although Cu-based catalysts capable of selectively producing CH 4 have been reported, they seldom sustain high selectivity at practically relevant current densities. Here, we created a single-step co-pyrolysis strategy toward generating and anchoring Cu sub-nanometer clusters (Cu SNC ) atop Cu-N x single-atom (SA) motifs embedded within N-doped carbon (NC), with controllable nanostructures through tuning of the synthesis parameters. Complementary spectroscopic analyses and density functional theory (DFT) calculations help reveal a structure−activity correlation that could guide the catalyst design. The Cu SNC @NC sample synthesized at 550 °C pyrolysis temperature (best described and modeled as Cu 3 -CuN 4 domains) represents the most effective combination of cluster size, metal-nitrogen coordination, and adsorption energetics needed to selectively promote CH 4 generation versus other eCO 2 R products. Incorporating pulsed electrolysis and hydrophobicity-modulated transport tuning at the triple-phase boundary (TPB) further enhanced CH 4 production achieving a partial CH 4 current density of ∼321 mA cm −2 , 53% Faradaic efficiency (FECH 4 ), and less than 4% combined FE for other eCO 2 R products, simplifying downstream CH 4 purification or upgrading. This work establishes generalizable principles for controlling Cu cluster atomicity and metal−nitrogen coordination, both of which are recognized determinants of CH 4 -efficient eCO 2 R.

CH4 production↗
DOE OSTI2021

From Intermolecular Interaction Energies and Observable Shifts to Component Contributions and Back Again: A Tale of Variational Energy Decomposition Analysis

Quantum chemistry in the form of density functional theory (DFT) calculations is a powerful numerical experiment for predicting intermolecular interaction energies. However, no chemical insight is gained in this way beyond predictions of observables. Energy decomposition analysis (EDA) can quantitatively bridge this gap by providing values for the chemical drivers of the interactions, such as permanent electrostatics, Pauli repulsion, dispersion, and charge transfer. These energetic contributions are identified by performing DFT calculations with constraints that disable components of the interaction. This review describes the second-generation version of the absolutely localized molecular orbital EDA (ALMO-EDA-II). The effects of different physical contributions on changes in observables such as structure and vibrational frequencies upon complex formation are characterized via the adiabatic EDA. Example applications include red- versus blue-shifting hydrogen bonds; the bonding and frequency shifts of CO, N 2 , and BF bound to a [Ru(II)(NH 3 ) 5 ] 2 + moiety; and the nature of the strongly bound complexes between pyridine and the benzene and naphthalene radical cations. Additionally, the use of ALMO-EDA-II to benchmark and guide the development of advanced force fields for molecular simulation is illustrated with the recent, very promising, MB-UCB potential.

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

Effect of Zn Addition on Phase Evolution in AlCrFeCoNiZn High–Entropy Alloy

The addition of Zn to AlCrFeCoNi high-entropy alloy (HEA) poses intriguing questions as to how it would affect phase evolution. Herein, the phase evolution in AlCrFeCoNiZn is studied using a combination of experimental techniques (X-ray diffraction, scanning electron microscopy, energy-dispersive spectroscopy, and differential scanning calorimetry) and computational (density-functional theory [DFT], calculation of phase diagrams, and machine-learning) methods. Mechanically alloyed and spark-plasma-sintered AlCrFeCoNiZn assumes a metastable single-phase, body-centered-cubic (BCC) structure that undergoes diffusion-controlled phase separation upon subsequent heat treatment to form separate (Al, Cr)-rich, (Fe, Co)-rich, and (Zn, Ni)-rich phases. The formation of (Al, Cr)-rich phase, not reported previously in AlCrFeCoNi-based HEAs, is attributed to strong clustering tendency of Cr–Zn and Cr–Ni pairs, combined with the strong ordering of Zn–Ni pair, driving out Cr that in turn combines with Al to form a (Al, Cr)-rich phase. In the DFT results, the formation of thermodynamically stable L1 2 phase is shown wherein Cr–Fe–Zn [Al–Ni-Co] preferably occupy1a (000) [3c (0 ½ ½)] positions. Furthermore, the sluggish diffusional transformation to L1 2 phase from BCC precursors is attributed to the small stacking-fault energy of AlCrFeCoNiZn. The equilibrated HEA exhibits a high microhardness of 8.24 GPa with an elastic modulus of 184 GPa.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Experimental and Theoretical Study of the Electronic Structures of Lanthanide Indium Perovskites LnInO 3

Ternary lanthanide indium oxides LnInO 3 (Ln = La, Pr, Nd, Sm) were synthesized by high-temperature solid-state reaction and characterized by Xray powder diffraction. Rietveld refinement of the powder patterns showed the LnInO 3 materials to be orthorhombic perovskites belonging to the space group Pnma, based on almost-regular InO 6 octahedra and highly distorted LnO 12 polyhedra. Experimental structural data were compared with results from density functional theory (DFT) calculations employing a hybrid Hamiltonian. Valence region X-ray photoelectron and K-shell X-ray emission and absorption spectra of the LnInO 3 compounds were simulated with the aid of the DFT calculations. Photoionization of lanthanide 4f orbitals gives rise to a complex final-state multiplet structure in the valence region for the 4f n compounds PrInO 3 , NdInO 3 , and SmInO 3 , and the overall photoemission spectral profiles were shown to be a superposition of final-state 4f n-1 terms onto the cross-section weighted partial densities of states from the other orbitals. The occupied 4f states are stabilized in moving across the series Pr-Nd-Sm. Band gaps were measured using diffuse reflectance spectroscopy. These results demonstrated that the band gap of LaInO 3 is 4.32 eV, in agreement with DFT calculations. This is significantly larger than a band gap of 2.2 eV first proposed in 1967 and based on the idea that In 4d states lie above the top of the O 2p valence band. However, both DFT and X-ray spectroscopy show that In 4d is a shallow core level located well below the bottom of the valence band. Band gaps greater than 4 eV were observed for NdInO 3 and SmInO 3 , but a lower gap of 3.6 eV for PrInO 3 was shown to arise from the occupied Pr 4f states lying above the main O 2p valence band.

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

Formic Acid: A Hydrogen-Bonding Cocatalyst for Formate Decomposition

Hydrogen bonding accelerates many catalytic reactions by orienting intermediates, stabilizing transition states, and even opening reaction pathways. However, most mechanistic studies regarding the decomposition of formic acid (FA), a promising hydrogen storage material, neglect hydrogen-bonding interactions even though FA is a strong hydrogen-bond donor and acceptor. Here, we probe the formation of bimolecular hydrogen-bonded complexes between FA and formate (FA–HCOO complexes) adsorbed on metal surfaces and how these complexes affect HCOO* decomposition. Using first-principles density functional theory (DFT) calculations on 12 close-packed (111)/(0001) and 8 open (100) surfaces of 12 transition metals—Ag, Au, Co, Cu, Ir, Ni, Os, Re, Pd, Pt, Rh, and Ru, we—show that FA–HCOO complexes are generally thermodynamically stable, even at elevated temperatures and pressures. We then illustrate that these complexes produce infrared spectroscopic signatures consistent with as yet unassigned experimental peaks. We last demonstrate that by stabilizing the dangling bond of monodentate HCOO*, these complexes significantly lower the barriers for rotation of HCOO* from a bidentate to a monodentate configuration, the rate-limiting step for HCOO* decomposition on many surfaces. FA thus acts as a cocatalyst for HCOO* decomposition. Our results may guide the community toward improved catalysts for reactions involving HCOO* such as FA decomposition, methanol steam reforming, and the water gas shift reaction. More broadly, our work highlights the ability of hydrogen bonding to modify the adsorbed structures of intermediates and lower the barriers for their reaction on heterogeneous catalysts. Lastly, this phenomenon can be relevant for other reactions involving ammonia, alcohols, and carboxylic acids.

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

Establishing the Role of Metal, Interface, and Vacancy Sites in Pt/TiO 2 -Catalyzed Acetic Acid Hydrodeoxygenation

Catalytic hydrodeoxygenation (HDO) following catalytic fast pyrolysis (CFP) offers an approach to convert the vapor-phase product of biomass pyrolysis to a stable bio-oil product by reducing the oxygen content. Fundamental insights into the HDO of carboxylic acids, which are a corrosive and acidic CFP product, on promising catalyst materials, such as Pt/TiO 2 , are needed to inform the design of multifunctional HDO catalysts with improved carbon efficiency. In this contribution, density functional theory (DFT) calculations were used to assess the role of Pt-metal and Pt-TiO 2 -interface sites on acetic acid HDO (AA-HDO), and to determine the effect of interfacial oxygen vacancies at the Pt-TiO 2 interface, by calculating the reaction energetics for key AA-HDO surface intermediates and elementary steps on each site type. Pt-metal sites, modeled via Pt(111), preferred to form undesired decarboxylation products (CH 4 and CO 2 ), whereas Pt-TiO 2 -interface sites, modeled via an anatase-supported Pt nanowire, favored the formation of desired deoxygenation products (acetaldehyde and ethane). Interfacial-vacancy sites lowered the activation energy barrier for the first C-O bond-scission step in AA-HDO, predicted to be the rate-limiting step for AA-HDO at the Pt-TiO 2 interface in the absence of a vacancy. These atomistic insights reveal the importance of metal-metal oxide interface sites in AA-HDO selectivity and can be used to inform the rational design of improved HDO catalysts.

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

Enhancing high-fidelity neural network potentials through low-fidelity sampling

The efficacy of neural network potentials (NNPs) critically depends on the quality of the configurational datasets used for training. Prior research using empirical potentials has shown that well-selected liquid–solid transitional configurations of a metallic system can be translated to other metallic systems. This study demonstrates that such validated configurations can be relabeled using density functional theory (DFT) calculations, thereby enhancing the development of high-fidelity NNPs. Training strategies and sampling approaches are efficiently assessed using empirical potentials and subsequently relabeled via DFT in a highly parallelized fashion for high-fidelity NNP training. Our results reveal that relying solely on energy and force for NNP training is inadequate to prevent overfitting, highlighting the necessity of incorporating stress terms into the loss functions. To optimize training involving force and stress terms, we propose employing transfer learning to fine-tune the weights, ensuring that the potential surface is smooth for these quantities composed of energy derivatives. This approach markedly improves the accuracy of elastic constants derived from simulations in both empirical potential-based NNPs and relabeled DFT-based NNPs. Overall, this study offers significant insights into leveraging empirical potentials to expedite the development of reliable and robust NNPs at the DFT level.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Dynamic Molecular Investigation of the Solid-Electrolyte Interphase of an Anode-Free Lithium Metal Battery Using In Situ Liquid SIMS and Cryo-TEM

Solid electrolyte interphase (SEI) has been widely perceived to play a critical role in the stable cycling of rechargeable batteries. However, associated with the fragile and air-sensitive nature of the SEI layer, delineation of the formation process and the nature of SEI remains a big challenge. Here, we use in situ liquid time-of-flight secondary ion mass spectroscopy (TOF-SIMS), cryo- transmission electron microscope (TEM) and density functional theory (DFT) calculation to delineate molecular process on the formation of SEI layer under the dynamic operating condition. We discover that the onset potential for SEI layer formation and the thickness of the SEI show dependence on the solvation shell structure. Using LiCoO 2 as a cathode and Cu film as an anode, the SEI is noticed to start to form at around 2.0 V and reach its final thickness (irreversible part, ~ 40-50 nm) at about 3.0 V in the 1 M LiPF 6 –EC/DMC electrolyte, while for the case of 1 M LiFSI–DME, the SEI starts to form at around 1.5 V and reaches its final thickness (~ 20 nm) at about 2.0 V. The in situ TOF-SIMS clearly indicates the outer SEI layer formation and dissipation upon charging and discharging, implying a continued evolution of electrolyte structure with extended cycling. In conclusion, the present work establishes a direct correlation between the molecular signature of SEI layer with solvation feature of electrolytes in lithium batteries, providing insights for tailoring SEI layer toward improved electrochemical properties of lithium batteries.

25 ENERGY STORAGE↗