SEARCH · Engineering Papers

Results for “computational molecular design”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 667 records · Page 37

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Enhanced Electrocatalytic and Selective CO 2 -to-CO Reduction by a Rhenium(I) Complex Bearing 6,6′-Substituted 2,2′-Bipyridines

The electrochemical reduction of CO 2 (CO 2 RR) into value-added chemicals offers a promising route toward a circular carbon economy and reduced reliance on fossil fuels. A detailed understanding of the structural and electronic factors governing the performance of molecular CO 2 RR electrocatalysts is essential for the design of efficient, tunable systems. Here, in this study, we report a series of rhenium(I) complexes, fac-[Re I (6,6′-(R) 2 -bpy)(CO) 3 Cl] (bpy = 2,2′-bipyridine; R = mesityl (mes), 2,4,6-triisopropylphenyl (trip), or isophthalic acid (phth)) and evaluate their electrocatalytic activity. Among these, fac-[Re I (6,6′-(mes) 2 -bpy)(CO) 3 Cl] exhibited the highest performance, enabling selective CO 2 -to-CO conversion for 1 hour with Faradaic efficiency (FE) > 97%, representing an unprecedented activity level for a Re-bpy catalysts. Single-crystal X-ray diffraction and density functional theory (DFT) calculations indicated that favorable CO 2 binding could be promoted by the tilting of the 6,6′-(mes)2-bpy ligand (from the Re-CO coordination plane), providing mechanistic insight into the observed enhancement. The study consequently demonstrates a rational correlation between the CO 2 electrocatalytic performance of Re-bpy catalysts and their structural variations, as derived from X-ray data and corroborated by computational modeling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Lattice thermal transport in two-dimensional alloys and fractal heterostructures

Abstract Engineering thermal transport in two dimensional materials, alloys and heterostructures is critical for the design of next-generation flexible optoelectronic and energy harvesting devices. Direct experimental characterization of lattice thermal conductivity in these ultra-thin systems is challenging and the impact of dopant atoms and hetero-phase interfaces, introduced unintentionally during synthesis or as part of deliberate material design, on thermal transport properties is not understood. Here, we use non-equilibrium molecular dynamics simulations to calculate lattice thermal conductivity of $${\mathrm {(Mo|W)Se_2}}$$ ( Mo | W ) Se 2 monolayer crystals including $${\mathrm {Mo}}_{1-x}{\mathrm {W}}_x{\mathrm {Se_2}}$$ Mo 1 - x W x Se 2 alloys with substitutional point defects, periodic $${\mathrm {MoSe_2}|\mathrm {WSe_2}}$$ MoSe 2 | WSe 2 heterostructures with characteristic length scales and scale-free fractal $${\mathrm {MoSe_2}}|{\mathrm {WSe_2}}$$ MoSe 2 | WSe 2 heterostructures. Each of these features has a distinct effect on phonon propagation in the crystal, which can be used to design fractal and periodic alloy structures with highly tunable thermal conductivities. This control over lattice thermal conductivity will enable applications ranging from thermal barriers to thermoelectrics.

42 ENGINEERING↗

DOE OSTI2025

Investigating Kinetic Mechanisms of Soot Formation in Plasma Pyrolysis of Methane via Active Learning (Final Technical Report)

Plasma pyrolysis of methane is an effective route for zero-carbon hydrogen production. Yet, soot generated from pyrolysis of hydrocarbons is detrimental to the climate and human health. There is ample experimental and theoretical evidence that suggests polycyclic aromatic hydrocarbons (PAHs) are the molecular precursors to soot particles. The reaction pathways of PAH formation are intricately dependent on a multitude of process parameters, whose kinetic mechanisms are not well-understood in plasma pyrolysis. This project aims to leverage advances in the kinetic modeling of soot formation in combustion, as well as in surrogate modeling and active learning, to systematically investigate the effects of process parameter on the kinetics of PAH formation in plasma pyrolysis of methane. To this end, we propose to use the PAH formation kinetics model developed by the PPPL/PU group based on the well-established ABF and HACA mechanisms, coupled with low-temperature plasma models. We will develop an active learning (AL) framework based on Bayesian optimization to systematically and data-efficiently explore the complex and multivariable parameter space of plasma pyrolysis in order to quantify the effects of plasma and feed parameters on the ABF and HACA kinetic pathways. AL is the branch of machine learning concerned with systematically querying samples from a system (experimental or computational) to train a data-driven model that maps design parameters to a performance criterion. We will use the data generated via AL to perform global sensitivity analysis, combined with uncertainty quantification, to elucidate the impact of different reaction pathways on minimizing formation of soot precursors. This study will result in an improved understanding of kinetics of PAH formation in plasma pyrolysis and can pave the way for more advanced mechanistic studies (e.g., soot nucleation mechanisms). Additionally, the findings will be useful for establishing practical strategies for increasing the pyrolysis efficiency and producing high-grade carbon for synthesis of nanomaterials.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Ultra-large-scale ab initio quantum chemical computation of bio-molecular systems: The case of spike protein of SARS-CoV-2 virus

The COVID-19 pandemic poses a severe threat to human health with an unprecedented social and economic disruption. Spike (S) glycoprotein of the SARS-CoV-2 virus is pivotal in understanding the virus anatomy, since it initiates the first contact with the ACE2 receptor in the human cell. We report results of ab initio computation of the spike protein, the largest ab initio quantum chemical computation to date on any bio-molecular system, using a divide and conquer strategy by focusing on individual structural domains. In this approach we divided the S-protein into seven structural domains: N-terminal domain (NTD), receptor binding domain (RBD), subdomain 1 (SD1), subdomain 2 (SD2), fusion peptide (FP), heptad repeat 1 with central helix (HR1-CH) and connector domain (CD). The entire Chain A has 14,488 atoms including the hydrogen atoms but excluding the amino acids with missing coordinates based on the PDB data (ID: 6VSB). The results include structural refinement, ab initio calculation of intramolecular bonding mechanism, 3- dimensional non-local inter-amino acid interaction with implications for the inter-domain interaction. Details of the electronic structure, interatomic bonding, partial charge distribution and the role played by hydrogen bond network are discussed. In the interaction among structural domains, we present new insights for crucial hinge-like movement and fusion process. Extension of such calculation to the interface between the S-protein binding domain and ACE2 receptor can provide a pathway for computational understanding of mutations and the design of therapeutic drugs to combat the COVID-19 pandemic.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2025

Ca X ML: Chemistry‐informed machine learning explains mutual changes between protein conformations and calcium ions in calcium‐binding proteins using structural and topological features

Proteins' flexibility is a feature in communicating changes in cell signaling instigated by binding with secondary messengers, such as calcium ions, associated with the coordination of muscle contraction, neurotransmitter release, and gene expression. When binding with the disordered parts of a protein, calcium ions must balance their charge states with the shape of calcium-binding proteins and their versatile pool of partners depending on the circumstances they transmit. Accurately determining the ionic charges of those ions is essential for understanding their role in such processes. However, it is unclear whether the limited experimental data available can be effectively used to train models to accurately predict the charges of calcium-binding protein variants. Here, we developed a chemistry-informed, machine-learning algorithm that implements a game theoretic approach to explain the output of a machine-learning model without the prerequisite of an excessively large database for high-performance prediction of atomic charges. We used the ab initio electronic structure data representing calcium ions and the structures of the disordered segments of calcium-binding peptides with surrounding water molecules to train several explainable models. Network theory was used to extract the topological features of atomic interactions in the structurally complex data dictated by the coordination chemistry of a calcium ion, a potent indicator of its charge state in protein. Our design created a computational tool of Ca X ML, which provided a framework of explainable machine learning model to annotate ionic charges of calcium ions in calcium-binding proteins in response to the chemical changes in an environment. Our framework will provide new insights into protein design for engineering functionality based on the limited size of scientific data in a genome space.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Knowledge-driven design of solid-electrolyte interphases on lithium metal via multiscale modelling

Abstract Due to its high energy density, lithium metal is a promising electrode for future energy storage. However, its practical capacity, cyclability and safety heavily depend on controlling its reactivity in contact with liquid electrolytes, which leads to the formation of a solid electrolyte interphase (SEI). In particular, there is a lack of fundamental mechanistic understanding of how the electrolyte composition impacts the SEI formation and its governing processes. Here, we present an in-depth model-based analysis of the initial SEI formation on lithium metal in a carbonate-based electrolyte. Thereby we reach for significantly larger length and time scales than comparable molecular dynamic studies. Our multiscale kinetic Monte Carlo/continuum model shows a layered, mostly inorganic SEI consisting of LiF on top of Li 2 CO 3 and Li after 1 µs. Its formation is traced back to a complex interplay of various electrolyte and salt decomposition processes. We further reveal that low local Li + concentrations result in a more mosaic-like, partly organic SEI and that a faster passivation of the lithium metal surface can be achieved by increasing the salt concentration. Based on this we suggest design strategies for SEI on lithium metal and make an important step towards knowledge-driven SEI engineering.

25 ENERGY STORAGE↗

DOE OSTI2025

Experimental and Computational Insights into the Structural Dynamics of the Fc Fragment of IgG1 Subtype from Biosimilar VEGF‐Trap

The constant fragment (Fc) of the immunoglobulin G1 (IgG1) subtype is increasingly recognized as a crucial scaffold in the development of advanced therapeutics due to its enhanced specificity, efficacy, and extended half‐life. A prime example is VEGF‐Trap (Aflibercept), a recombinant fusion protein that merges the Fc region of the IgG1 subtype with the binding domains of vascular endothelial growth factor receptors (VEGFR)‐1 and VEGFR‐2. The Fc region's role in N‐glycosylation is particularly important, as it significantly influences protein stability. Herein, the first near‐physiological temperature structures of the N‐glycan‐bound Fc fragment of IgG1 subtype from a biosimilar VEGF‐Trap are presented, determined using the SPring‐8 Angstrom Compact free electron LAser (SACLA) and the Turkish Light Source (Turkish DeLight). Comparative analysis with cryogenic structures, including existing data, reveals alternate conformations within the glycan‐binding pocket. Furthermore, molecular dynamics simulations indicate the presence of a high degree of structural plasticity, explaining how the protein adapts its structure through conformational changes. The observed structural fluctuations/conformational changes demonstrate the effect of N‐glycans on protein stability. These findings offer new insights into the molecular basis of Fc‐mediated functions and provide valuable information for the design of next‐generation therapeutics.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2020

Water Structure and Properties at Hydrophilic and Hydrophobic Surfaces

The properties of water on both molecular and macroscopic surfaces critically influence a wide range of physical behaviors, with applications spanning from membrane science to catalysis to protein engineering. Yet, our current understanding of water interfacing molecular and material surfaces is incomplete, in part because measurement of water structure and molecular-scale properties challenges even the most advanced experimental characterization techniques and computational approaches. This review highlights progress in the ongoing development of tools working to answer fundamental questions on the principles that govern the interactions between water and surfaces. One outstanding and critical question is what universal molecular signatures capture the hydrophobicity of different surfaces in an operationally meaningful way, since traditional macroscopic hydrophobicity measures like contact angles fail to capture even basic properties of molecular or extended surfaces with any heterogeneity at the nanometer length scale. Resolving this grand challenge will require close interactions between state-of-the-art experiments, simulations, and theory, spanning research groups and using agreed-upon model systems, to synthesize an integrated knowledge of solvation water structure, dynamics, and thermodynamics.

catalysis (heterogeneous)↗

DOE OSTI2021

Solvent and Temperature Effects on Photoinduced Proton-Coupled Electron Transfer in the Marcus Inverted Region

Concerted proton-coupled electron transfer (PCET) in the Marcus inverted region was recently demonstrated. Understanding the requirements for such reactivity is fundamentally important and holds promise as a design principle for solar energy conversion systems. Herein, we investigate the solvent polarity and temperature dependence of photoinduced proton-coupled charge separation (CS) and charge recombination (CR) in anthracene–phenol–pyridine triads: 1 (10-(4-hydroxy-3-(4-methylpyridin-2-yl)benzyl)anthracene-9-carbonitrile) and 2 (10-(4-hydroxy-3-(4-methoxypyridin-2-yl)benzyl)anthracene-9-carbonitrile). Both the CS and CR rate constants increased with increasing polarity in acetonitrile:n-butyronitrile mixtures. The kinetics were semi-quantitatively analyzed where changes in dielectric and refractive index, and thus consequently changes in driving force (–ΔG°) and reorganization energy (λ), were accounted for. The results were further validated by fitting the temperature dependence, from 180 to 298 K, in n-butyronitrile. The analyses support previous computational work where transitions to proton vibrational excited states dominate the CR reaction with a distinct activation free energy (ΔG* CR ~ 140 meV). However, the solvent continuum model fails to accurately describe the changes in ΔG° and λ with temperature via changes in dielectric constant and refractive index. Satisfactory modeling was obtained using the results of a molecular solvent model, which predicts that λ decreases with temperature, opposite to that of the continuum model. To further assess the solvent polarity control in the inverted region, the reactions were studied in toluene. Nonpolar solvents decrease both ΔG° CR and λ, slowing CR into the nanosecond time regime for 2 in toluene at 298 K. This demonstrates how PCET in the inverted region may be controlled to potentially use proton-coupled CS states for efficient solar fuel production and photoredox catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Charge carrier dynamics in two-dimensional hybrid perovskites: Dion–Jacobson vs. Ruddlesden–Popper phases

Two-dimensional (2D) halide perovskites are promising materials for environmentally stable next-generation optoelectronic device applications. Besides the widely investigated Ruddlesden–Popper (RP) phases, Dion–Jacobson (DJ) phases are attracting considerable attention due to their rapid emergence as efficient solar cell materials. However, there is very little atomistic understanding of the charge carrier dynamics under ambient conditions for these DJ-phases, limiting the possibilities to tune their optoelectronic performances through compositional engineering routes. In this study, by combining nonadiabatic molecular dynamics with time-domain density functional theory methods at room temperature, we compare the dominant non-radiative carrier recombination and dephasing processes in RP and DJ monolayered lead halide perovskites. Our systematic study demonstrates that performance-limiting nonradiative carrier recombination processes greatly depends on the electron–phonon interactions induced by structural fluctuations and instantaneous charge localization in these materials. The stiffer interlayer packing due to the presence of single spacer dications, which separates the lead iodide slabs, reduces the thermal fluctuations in the DJ-phase to a greater extent than that in the RP-phase 2D-perovskites. Specific electronic coupling between the closely spaced lead iodide layers enhances the delocalization of band-edge charge densities in DJ-phase systems. Compared to the RP-phase, reduced inelastic electron–phonon scattering in DJ-phase perovskites significantly limits intrinsic non-radiative recombination processes. The consequent enhancement in the photogenerated charge carrier lifetime makes DJ-phase perovskites potentially suitable for various optoelectronic devices. The computational insights gained from this study allow us to outline a set of robust design principles for DJ-phase perovskites to strategically tune their optoelectronic properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2005

Applications of Quantum Chemistry to the Study of Carbon Nanotubes

For several years, scientists at NASA Ames have been studying the properties of carbon nanotubes using various experimental and computational methods. In this talk, I will compare different strategies for using quantum chemistry calculations to describe the electronic structure, deformation and chemical functionalization of single wall carbon nanotubes (SWNT) and the physisorption of small molecules on nanotube surfaces. The SWNT can be treated as an infinite (periodic) or finite length carbon cylinder or as a polycyclic aromatic hydrocarbon (PAH) molecule with an imposed curvature maintained by external constraints (as if it were cut out of the SWNT surface). Calculations are carried out using DFT and MP2 methods and a variety of atomic orbital basis sets from minimal (STO-3G) to valence triple zeta. The optimal approach is based on the particular SWNT property of interest. Examples to be discussed include: nanotube fluorination and other functionalization reactions; coating of nanotubes by water vapor and low-molecular weight organic molecules; and the nature of the interface between SWNT and liquids such as water and amines. In many cases, the quantum chemistry calculations are used to parameterize or validate force fields for molecular dynamics simulations. The results of these calculations have helped explain experimental data and contributed to the design of novel materials and sensors based on carbon nanotubes. Some of this research is described in the following papers:

Jaffe, Richard L.↗

DOE OSTI2026

Finite-element-based simulations of electrodes for CO 2 cascade reduction reactions

The multielectron reduction of CO 2 to liquid fuels could be a path to scalable energy storage, but reaching this goal requires major advances in catalysis and systems engineering. Cascade catalysis, which couples sequential reactions without isolating intermediates, has emerged as a promising route to enhance selectivity and efficiency in CO 2 reduction (CO 2 R). In this review, we examine how finite-element-based simulations of continuum model [finite element method (FEM)] approaches are being used to analyze and guide CO 2 R cascade systems. We first outline the fundamentals of cascade catalysis and recent advances in catalytic materials (metallic, molecular, and hybrid architectures). We then focus on FEM developments at the electrode and device scales, emphasizing how these models capture transport phenomena, local microenvironments, and geometry-dependent effects. To clarify design principles, we present case studies of cascade electrodes organized in systems without and with integrated semiconductors. We further emphasize the integration of FEM with multiscale frameworks (density functional theory, molecular dynamics, kinetic Monte Carlo) and its role in bridging atomic-level insights with device-level performance. Finally, we identify current limitations and future prospects, including improved boundary conditions, coupling with operando experiments, and machine learning-accelerated model development. Together, these insights provide design principles for next-generation CO 2 R cascade systems for efficient solar fuel production.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Nonlinear thermodynamic computing out of equilibrium

We present the design for a thermodynamic computer that can perform arbitrary nonlinear calculations in or out of equilibrium. Simple thermodynamic circuits, fluctuating degrees of freedom in contact with a thermal bath and confined by a quartic potential, display an activity that is a nonlinear function of their input. Such circuits can therefore be regarded as thermodynamic neurons, and can serve as the building blocks of networked structures that act as thermodynamic neural networks, universal function approximators whose operation is powered by thermal fluctuations. We simulate a digital model of a thermodynamic neural network, and show that its parameters can be adjusted by genetic algorithm to perform nonlinear calculations at specified observation times, regardless of whether the system has attained thermal equilibrium. This work expands the field of thermodynamic computing beyond the regime of thermal equilibrium, enabling fully nonlinear computations, analogous to those performed by classical neural networks, at specified observation times.

Whitelam, Stephen [Lawrence Berkeley National Labo↗

DOE OSTI2021

Mechanical Behavior of Hybrids with Hyper-Connected Molecular Networks

Our research is focused on fundamental studies related to the molecular design, synthesis, characterization and modeling of molecular-reinforced hybrid glass films for superior mechanical and fracture properties. Molecular-reinforced hybrid glass films exhibit unique electro-optical properties while maintaining excellent thermal stability. They have important technological application for emerging nanoscience and energy technologies including anti-reflective and ultra-barrier layers in photovoltaics, display technologies, as well as quantum technologies. However, they are often inherently brittle in nature, do not adhere well to adjacent substrates, and exhibit poor mechanical properties. Furthermore, a robust understanding of the underlying mechanisms of deformation, fracture and fatigue has been lacking particularly given their mechanically fragile nature. The related lack of structure-property relationships between molecular structure, mechanical properties and processing represents a fundamental challenge to their integration into device technologies, compromises their thermomechanical reliability and has thwarted efforts to improve mechanical and, particularly, fracture properties.

36 MATERIALS SCIENCE↗

DOE OSTI2024

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↗

DOE OSTI2023

Automated reaction kinetics and network exploration (Arkane): A statistical mechanics, thermodynamics, transition state theory, and master equation software

The open-source statistical mechanics software described here, Arkane–Automated Reaction Kinetics and Network Exploration–facilitates computations of thermodynamic properties of chemical species, high-pressure limit reaction rate coefficients, and pressure-dependent rate coefficient over multi-well molecular potential energy surfaces (PES) including the effects of collisional energy transfer on phenomenological kinetics. Arkane can use estimates to fill in information for molecules or reactions where quantum chemistry information is missing. The software solves the internal energy master equation for complex unimolecular reaction systems. Inputs to the software include converged electronic structure computations performed by the user using a variety of supported software packages (Gaussian, Molpro, Orca, TeraChem, Q-Chem, Psi4). The software outputs high-pressure limit rate coefficients and pressure-dependent phenomenological rate coefficients, as well as computed thermodynamic properties (enthalpy, entropy, and constant pressure heat capacity) with added energy corrections. Some of the key features of Arkane include treatment of 1D, 2D or ND hindered internal rotation modes, treatment of free internal rotation modes, quantum tunneling effect consideration, transition state theory (TST) and Rice-Ramsperger-Kassel-Marcus (RRKM) rate coefficient computations, master equation solution with four implemented methods, inverse-Laplace transform of high-pressure limit rate coefficients into the energy domain, energy corrections based on bond-additivity or isodesmic reactions, automated and efficient PES exploration, and PES sensitivity analysis. The present work describes the design of Arkane, how it should be used, and refers to the theory that it employs. Arkane is distributed via the RMG-Py software suite (https://github.com/ReactionMechanismGenerator/RMG-Py).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Revealing EDL-driven reduction mechanisms in binary, ternary, and quaternary fluorinated electrolytes via an integrated MD–DFT–ML framework

Accurately predicting solid electrolyte interphase (SEI) formation requires explicitly resolving the electric double layer (EDL) structure, which deviates significantly from that of the bulk electrolyte. Although an established molecular dynamics (MD) and Density Functional Theory (DFT) framework can model SEI formation by evaluating reduction reactions of local clusters in the EDL, it suffers from a combinatorial computational bottleneck. To overcome this limitation, we introduce a machine-learning-accelerated simulation workflow (MD–DFT–ML), integrating a gradient-boosted regression model trained on EDL composition data to efficiently predict reduction potentials. We apply this framework to seven fluorinated electrolytes comprising fluorinated anions, a fluorinated ester solvent, two types of diluent (ion-solvating ester vs. non-solvating ether), and an FEC additive. The analysis shows that the EDL selectively accumulates cation-binding species; consequently, the non–cation-binding ether diluent rarely enters the EDL and makes minimal contributions to SEI formation. DFT calculations on statistically representative EDL clusters provide reduction potentials and fluorine-release pathways, while the ML model, which substantially reduces the DFT workload, predicts cluster reduction energies with a mean absolute error of 0.1 eV. The combined MD–DFT–ML approach also quantifies contributions from different sources to LiF formation in the SEI. This methodology establishes a generalizable route for multiscale modeling electrolyte and interphase design for next-generation electrochemical energy-storage systems.

DFT-MD-ML workflow↗

DOE OSTI2022

Analyzing the effect of pressure on the properties of point defects in $\gamma\mathrm{U–Mo}$ through atomistic simulations

Uranium–molybdenum (U–Mo) alloys in monolithic fuel foil are the primary candidate for the conversion of high-performance research reactors in the USA. Monolithic fuel is utilized in a plate-type design with a zirconium diffusion barrier and aluminum cladding. These fuel types are unique in that they contain no plenum for the release of fission gases, which, in conjunction with the aluminum cladding, can lead to large stress states within the fuel. The nature of how fundamental processes of radiation damage, including the evolution of point defects, under such stresses occur is unknown. In this work, we present molecular dynamics simulations of the formation energy of point defects under applied stress. We report this work will allow for the implementation of stress-dependent microstructural evolution models of nuclear fuels, including those for both fission gas bubble growth and creep, which are critical to ensure the stable and predictable behavior of research reactor fuels.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗