Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “electrolyte theory”

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 361 records · Page 20

Effect of Added Salt on Disordered Poly(ethylene oxide)-Block-Poly(methyl methacrylate) Copolymer Electrolytes

In this report we studied the effect of salt addition on a diblock copolymer system with a negative Flory-Huggins interaction parameter, χ, indicative of attractive interactions between the two blocks. The system studied is poly(ethylene oxide)-block-poly(methyl methacrylate) (PEO-PMMA) with added lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) salt. We studied two asymmetric block copolymers, PEO-PMMA(10-33) and PEO-PMMA(10-64), where the numbers refer to the molar masses of the blocks in kg mol -1 . The small-angle X-ray scattering (SAXS) profiles for PEO-PMMA(10-33) were featureless at all salt concentrations. In contrast, PEO-PMMA(10-64) exhibited SAXS peaks when the salt concentration was between 0.22 ≤ m (mol Li/kg polymer) ≤ 0.44. The appearance of SAXS peaks only in PEO-PMMA(10-64) is consistent with the predictions of ionic self-consistent field theory developed by de la Cruz and co-workers, which predicts that in systems with negative χ, ordered phases are only found when the volume fraction of the ionic block is about 10%.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

From Micro-environments to Macroscopic Effects: How the Alkaline Hydrogen Evolution Reaction Drives Cu Cathodic Corrosion

Cathodic corrosion of copper (Cu) has posed a significant challenge for over a century, impeding various technological progresses such as electrochemical conversion of CO 2 (eCO 2 RR) into fuels and other value-added carbon products. Here, in this study, employing a combined Density Functional Theory (DFT) and kinetic Monte Carlo (kMC) simulation approach, we delve into the atomistic level mechanism driving this phenomenon in Cu. Our hypothesis posits the pivotal role of alkaline hydrogen evolution reaction (HER) in facilitating cathodic corrosion in Cu. We rigorously develop a pH-dependent hydroxide (OH) adsorption mechanism and calculate the equilibrium OH coverage (𝜃 OH ) at varying pH levels, the thermodynamic stability of subsurface oxygen (O sub ), as well as the Cu-vacancy mediated diffusion of subsurface oxygens (O sub ). Through comprehensive analysis, we establish correlation among various microenvironments, including oxygen diffusion in subsurface layers, pH-dependent OH adsorption, and Cu dissolution into the electrolyte as (Cu-OH) complexes. Furthermore, our investigation explores the correlation between surface coordination environment of active sites and cathodic corrosion of Cu. Finally, by integrating DFT-derived thermodynamic data into a kMC model, we successfully predict the formation of experimentally observed corrosion pits on Cu-surfaces. This combined approach not only advances our fundamental understanding of Cu cathodic corrosion but also offers insights crucial for developing effective corrosion mitigation strategies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Accurate simulation of surfaces and interfaces of ten FCC metals and steel using Lennard–Jones potentials

Abstract The earlier integration of validated Lennard–Jones (LJ) potentials for 8 fcc metals into materials and biomolecular force fields has advanced multiple research fields, for example, metal–electrolyte interfaces, recognition of biomolecules, colloidal assembly of metal nanostructures, alloys, and catalysis. Here we introduce 12-6 and 9-6 LJ parameters for classical all-atom simulations of 10 further fcc metals (Ac, Ca (α), Ce (γ), Es (β), Fe (γ), Ir, Rh, Sr (α), Th (α), Yb (β)) and stainless steel. The parameters reproduce lattice constants, surface energies, water interfacial energies, and interactions with (bio)organic molecules in 0.1 to 5% agreement with experiment, as well as qualitative mechanical properties under standard conditions. Deviations are reduced up to a factor of one hundred in comparison to earlier Lennard–Jones parameters, embedded atom models, and density functional theory. We also explain a quantitative correlation between atomization energies from experiments and surface energies that supports parameter development. The models are computationally very efficient and applicable to an exponential space of alloys. Compatibility with a wide range of force fields such as the Interface force field (IFF), AMBER, CHARMM, COMPASS, CVFF, DREIDING, OPLS-AA, and PCFF enables reliable simulations of nanostructures up to millions of atoms and microsecond time scales. User-friendly model building and input generation are available in the CHARMM-GUI Nanomaterial Modeler. As a limitation, deviations in mechanical properties vary and are comparable to DFT methods. We discuss the incorporation of reactivity and features of the electronic structure to expand the range of applications and further increase the accuracy.

Kanhaiya, Krishan (ORCID:0000000236222655)↗

Tracking ion intercalation into layered Ti 3 C 2 MXene films across length scales

Enhancing the energy stored and power delivered by layered materials relies strongly on improved understanding of the intricate interplay of electrolyte ions, solvents, and electrode interactions as well as the role of confinement. Here we report a highly integrated study with multiscale theory/modelling and experiments to track the intercalation of aqueous Li + , Na + , K + , Cs + , and Mg 2+ ions into Ti 3 C 2 MXene. The integrated analysis of experiments assisted by theory/modelling allows for a deep understanding of energy storage processes highlighting the importance of the dynamics of cations, their positionings between MXene sheets, their effects on mechanical properties and capacitive energy storage. Computational simulations and operando calorimetry measurements prove the processes involving cation dehydration and H+ rehydration, showing a good correlation for heat variations between experiments and theory. Operando liquid AFM mapped energy dissipation of ions appears non-uniformly across the MXene surface, indicating heterogeneities of ions inside the MXene and confirming partially the ion behaviour obtained in theory. We directly demonstrate that the average distance between the cation and MXene surface follows a modified two-sided Helmholtz model when plotted versus the open circuit potential capacitance, revealing a different electrical double layer mechanism in confinement. This new fundamental understanding lays the foundation for improved functional devices utilizing electrodes and membranes made of two-dimensional materials.

36 MATERIALS SCIENCE↗

Review—Development of Highly Active and Durable Hybrid Compressive Platinum Lattice Catalysts for Polymer Electrolyte Membrane Fuel Cells: Mathematical Modeling and Experimental Work

This review provides a comprehensive overview on the development of highly active and durable platinum catalysts with ultra-low Pt loadings for polymer electrolyte membrane fuel cells (PEMFCs) through a combined mathematical modeling and experimental work. First, simulation techniques were applied to evaluate the validity of the Tafel approximation for the calculation of the mass activity (MA) and specific activity (SA). A one-dimensional agglomeration model was developed and solved to understand the effects of exchange current density, porosity, agglomerate size, Nafion ® film thickness, and Pt loading on the MA and SA. High porosity (> 60%) and agglomerations at high Pt loadings cause the loss of the Tafel approximation and consequently the decrease in MA and SA. A new structure parameter was introduced to estimate the real porous structure using the fractal theory. The volumetric catalyst density was corrected by the fractal dimension (measured by Hg porosimetry), which gave a good agreement with the experimental values. The loading-dependent Tafel equation was then derived, which contains both the utilization and the non-linear scaling factor. Second, activated carbon composite support (ACCS) with optimized surface area, porosity, pore size, and pore size distribution was developed. The hydrophilic/hydrophobic ratio, structural properties (amorphous/crystalline ratio), and the number of active sites were optimized through metal-catalyzed pyrolysis. Stability of ACCS and Pt/ACCS were evaluated using an accelerated stress test (AST). The results indicated that Pt/ACCS showed no significant loss of MA and power density after 5,000 cycles at 1.0–1.5 V, while the commercial Pt/C catalysts showed drastic losses of MA and power density. Finally, monolayers of compressed Pt (core–shell-type Pt 3 Co 1 ) catalysts were structured by diffusing Co atoms (previously embedded in ACCS) into Pt. Compressive Pt lattice (Pt * ) catalysts were synthesized through an annealing procedure developed at the University of South Carolina (USC). The Pt * /ACCS catalyst showed high initial power density (rated) of 0.174 g Pt kW −1 and high stability (24 mV loss) at 0.8 A cm −2 after 30,000 cycles (0.6–1.0 V). The outstanding performance of Pt * /ACCS is due to the synergistic effect of ACCS and compressive Pt * lattice.

25 ENERGY STORAGE↗

A metal-supported single-atom catalytic site enables carbon dioxide hydrogenation

Nitrogen-doped graphene-supported single atoms convert CO 2 to CO, but fail to provide further hydrogenation to methane – a finding attributable to the weak adsorption of CO intermediates. To regulate the adsorption energy, here we investigate the metal-supported single atoms to enable CO 2 hydrogenation. We find a copper-supported iron-single-atom catalyst producing a high-rate methane. Density functional theory calculations and in-situ Raman spectroscopy show that the iron atoms attract surrounding intermediates and carry out hydrogenation to generate methane. The catalyst is realized by assembling iron phthalocyanine on the copper surface, followed by in-situ formation of single iron atoms during electrocatalysis, identified using operando X-ray absorption spectroscopy. The copper-supported iron-single-atom catalyst exhibits a CO 2 -to-methane Faradaic efficiency of 64% and a partial current density of 128 mA cm –2 , while the nitrogen-doped graphene-supported one produces only CO. The activity is 32 times higher than a pristine copper under the same conditions of electrolyte and bias.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Molybdenum Carbide Electrocatalyst In Situ Embedded in Porous Nitrogen–Rich Carbon Nanotubes Promotes Rapid Kinetics in Sodium–Metal–Sulfur Batteries

This work is the first report of a molybdenum carbide-based electrocatalyst for sulfur-based sodium metal batteries (SMBs/NMBs). MoC/Mo 2 C is in-situ grown on nitrogen-doped carbon nanotubes in parallel with formation of extensive nanoporosity. Sulfur impregnation (50 wt% S) results in unique triphasic architecture termed MoC/Mo 2 C@PCNT-S. Quasi-solid-state phase transformation to Na 2 S is promoted in carbonate electrolyte, with in-situ time-resolved Raman, XPS and optical analysis demonstrating minimal soluble polysulfides. MoC/Mo 2 C@PCNT-S cathodes delivered among the most promising rate performance characteristics in literature, achieving 987 mAh g -1 at 1 Ag -1 , 818 mAh g -1 at 3 A g -1 , and 621 mAh g -1 at 5 A g -1 . The cells deliver superior cycling stability, retaining 650 mAh g -1 after 1000 cycles at 1.5 Ag -1 , corresponding to 0.028% capacity decay per cycle. High mass loading cathodes (64 wt% S, 12.7 mg cm -2 ) also show cycling stability, with anode degradation due to deep plating/stripping driving capacity decay. Density functional theory (DFT) demonstrates that formation energy of Na 2 S x (1 ≤ x ≤ 4) on surface of MoC/Mo 2 C is significantly lowered compared to analogous redox in liquid. Strong binding of Na 2 S x (1 ≤ x ≤ 4) on MoC/Mo 2 C surfaces results from charge transfer between the sulfur and Mo sites on carbides' surface.

25 ENERGY STORAGE↗

Accelerating amorphous polymer electrolyte screening by learning to reduce errors in molecular dynamics simulated properties

Polymer electrolytes are promising candidates for the next generation lithium-ion battery technology. Large scale screening of polymer electrolytes is hindered by the significant cost of molecular dynamics (MD) simulation in amorphous systems: the amorphous structure of polymers requires multiple, repeated sampling to reduce noise and the slow relaxation requires long simulation time for convergence. Here, we accelerate the screening with a multi-task graph neural network that learns from a large amount of noisy, unconverged, short MD data and a small number of converged, long MD data. We achieve accurate predictions of 4 different converged properties and screen a space of 6247 polymers that is orders of magnitude larger than previous computational studies. Further, we extract several design principles for polymer electrolytes and provide an open dataset for the community. Our approach could be applicable to a broad class of material discovery problems that involve the simulation of complex, amorphous materials.

36 MATERIALS SCIENCE↗

Mechanism of ion pairing–induced co-ion penetration into the stern layer

In classical electric double layer (EDL) theory, co-ions—ions carrying the same charge as the electrode—are assumed to be excluded from the Stern layer because of electrostatic repulsion. However, ion pairing with adsorbed counterions may enable co-ions to access the EDL, a phenomenon that remains underexplored. Key questions therefore remain unresolved, including how deeply co-ions penetrate, how ion-specific properties govern them, and how ion pairing within the EDL differs from bulk solution. Here, in this study, we employ molecular dynamics simulations combined with well-tempered metadynamics to quantify co-ion penetration and ion-pairing thermodynamics at positively charged Mg-Al layered double hydroxide interfaces (σ = +20.6 μC/cm 2 ) in aqueous chloride electrolytes across 0.005–3 M range. We reveal a concentration-driven crossover: at dilute conditions, counterion–co-ion contact ion pairing (CIP) is thermodynamically unfavorable or only weakly favorable (ΔG CIP = +0.12 kcal/mol for Li + to −0.42 kcal/mol for K + ), leading to co-ion exclusion. As concentration increases (≥0.5 M), enhanced screening reduces surface repulsion and stabilizes CIP thermodynamics (ΔG CIP = −0.89 kcal/mol for Li + to −1.05 kcal/mol for K + at 3 M), enabling monovalent co-ions to penetrate the Stern layer following the hydration hierarchy K + > Na + > Li + , while divalent co-ions (Mg 2+ , Ca 2+ ) remain in the diffuse layer, forming only solvent-separated pairs. Comparison with bulk solution shows interfacial ion pairing is suppressed at low concentrations but converges to bulk-like thermodynamics at high concentrations, indicating surface effects vanish under strong screening. These results establish a quantitative framework linking ion-pairing thermodynamics to co-ion penetration and providing strategies to design electrochemical interfaces.

Co-ion penetration↗

Ir(hkl) Surface Electrochemistry in a Nonadsorbing Acidic Medium

The fundamental properties of electrochemical materials depend on the multiple and often complex interactions between electrode surface sites and electrolyte species at the electrochemical interface. Despite Iridium use in electrolyzer systems, much of its surface electrochemistry remains underexplored. This study investigates the surface electrochemistry of Ir(111), Ir(100), and Ir(110) surfaces in acidic media. Using cyclic voltammetry and CO charge displacement experiments, we establish the charge states and adsorbate coverages as a function of the electrode potential, revealing the presence of hydrogen and hydroxyl co-adsorption at low potentials on (111), and almost no coverage of H ad on (110) facet. In situ Shell Isolated Nanoparticle Enhanced Raman Spectroscopy experiments provide direct evidence of the formation of key adsorbate species, such as hydrogen, hydroxyl, and oxygen, but most importantly, their interactions with interfacial water, confirmed by Density Functional Theory calculations. Our findings highlight the role of co-adsorption and interspecies interactions, with microkinetic adsorption voltammetry simulations corroborating the influence of lateral interactions on adsorption dynamics, particularly for Ir(100) where the OHad formation occurs as a sharp adsorption/desorption current. Our results underscores the importance of interfacial water and hydrogen bonding networks in shaping the electrochemical behavior on Ir surfaces, refining our baseline understanding of the Ir surface electrochemistry necessary for the development of advanced Ir-based electrochemical materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Theory-Enabled High-Throughput Screening of Ion Dissociation Explains Conductivity Enhancements in Diluted Ionic Liquid Mixtures

The growing demand for room-temperature ionic liquids (RTILs) for energy applications necessitates the development of an efficient screening platform. Here, in this study, we successfully developed a fully automated high-throughput RTIL screening platform specifically designed for assessing ionic conductivity. By utilizing the 96 wells of a microtiter plate as individual electrolysis cells, we measured the ionic conductivity of 22 different RTILs, encompassing various combinations of cations and anions, and benchmarked the values with existing literature. We also employed the screening platform to investigate the conductivities of RTIL mixtures with a nonaqueous solvent, ethylene glycol (EG). Specific combinations of RTILs with EG result in approximately 200% enhancement in the conductivity values compared to the pure RTILs. To understand the underlying mechanisms responsible for this enhancement, we developed a theoretical framework for ionic conductivity that considers factors such as the degree of dissociation, viscous forces, and molal volume of the RTIL-EG mixtures. The optimized electrolyte mixture was then employed in the migration-assisted moisture gradient (MAMG) CO 2 capture process to study the effects of improved ionic conductivity on the energy efficiency of the process. Notably, the enhanced conductivity of the RTIL-EG mixture led to a nearly 50% reduction in energy consumption for capturing CO 2 . These outcomes highlight the effectiveness of our strategy in screening RTILs and improving existing processes. Moreover, this fully automated high-throughput setup, combined with the developed theoretical framework, provides a comprehensive platform for screening and studying RTIL mixtures with different solvents, enabling their application in various fields.

25 ENERGY STORAGE↗

On the Statistical Mechanics of Mass Accommodation at Liquid–Vapor Interfaces

Here we propose a framework for describing the dynamics associated with the adsorption of small molecules to liquid-vapor interfaces using an intermediate resolution between traditional continuum theories that are bereft of molecular detail and molecular dynamics simulations that are replete with them. In particular, we develop an effective single particle equation of motion capable of describing the physical processes that determine thermal and mass accommodation probabilities. The effective equation is parametrized with quantities that vary through space away from the liquid-vapor interface. Of particular importance in describing the early time dynamics is the spatially dependent friction, for which we propose a numerical scheme to evaluate from molecular simulation. Taken together with potentials of mean force computable with importance sampling methods, we illustrate how to compute the mass accommodation coefficient and residence time distribution. Throughout, we highlight the case of ozone adsorption in aqueous solutions and its dependence on electrolyte composition.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Real-Space Local Dynamics in 1,2,3-Triazole Using Inelastic Neutron Scattering

Enhancing proton transport in polymer electrolytes is crucial for advancing next-generation solid-state batteries, yet our understanding of proton conductivity in nonaqueous environments remains limited due to a lack of atomic-scale insights. Here, in this study, we investigated the atomic-scale dynamics of 1,2,3-triazole, a small molecule capable of dynamic hydrogen bonding, as a model system for proton hopping in nonaqueous environments. Using the real-space correlation function determined by the double Fourier transformation of inelastic neutron scattering spectra, we identified that the self-motion of protons and intermolecular dynamics occur on comparable time scales. Furthermore, we observed that the activation energy associated with the intermolecular dynamics matches the energy barrier for molecular rotations determined through Density Functional Theory calculations. These findings underscore the importance of controlling molecular dynamics at the atomic scale to control proton transport. Additionally, we demonstrated that intermolecular dynamics in systems involving protons can be studied using inelastic neutron scattering even without deuteration, thereby providing a broader avenue for studying atomic-scale dynamics in soft matter systems.

Shinohara, Yuya [Oak Ridge National Laboratory (OR↗

Understanding of three different polyvinylpyrrolidone (PVP) based battery binders blends on graphene surfaces from first principles via DFT simulations

Binders play a crucial role in binding electrodes, current collectors, and conductive agents during battery manufacturing. The selection of the binders and their blends considerably affect various properties such as mechanical durability, adhesion performance, ionic/electronic conductivities, and solid electrolyte interface stability for lithium-ion batteries. In this work, polyvinylidene di-fluoride (PVDF), polyacrylic acid (PAA), polyvinylpyrrolidone (PVP), and lithiated polyacrylic acid (Li-PAA) were selected as the model binders, and their blends were studied on graphene surfaces. Further, among them, PVDF is hardly soluble in water unlike the other binders, so in practice, the blend of PVP:PVDF is not thought to be a proper combination. However, the PVP:PVDF binder has also been included because the purpose of this work is to calculate several electronic properties such as binding energies, intermolecular interaction energies, bond critical points, electron density, and the Laplacian of electron density at critical bonding regions of all these blends at the molecular level via density functional theory (DFT) simulations in order to evaluate and compare how is the interaction strength and bonding type of three different lithium-ion batteries (LIBs) binder blends on graphene surfaces. Most stable binder pairs and their binding mechanisms on graphene surfaces were studied. The delocalization of lithium in Li-PAA was studied with the presence of water.

25 ENERGY STORAGE↗

Resolving the Solvation Structure and Transport Properties of Aqueous Zinc Electrolytes from Salt-in-Water to Water-in-Salt Using Neural Network Potential

Zn Cl 2 solutions are promising electrolytes for aqueous zinc-ion batteries. Here, we report a joint computational and experimental study of the structural and dynamic properties of aqueous Zn Cl 2 electrolytes with concentrations ranging from salt-in-water to water-in-salt (WIS). By developing a neural network potential (NNP) model, we perform molecular dynamics (MD) simulations with accuracy but at much larger lengths and longer timescales. The NNP predicted structures are validated by the structure factors measured by X-ray total scattering experiments. The MD trajectories provide a comprehensive and quantitative picture of the Zn 2 + solvation shell structures. Additionally, we find that the O − H covalent bonds in water are strengthened with increasing salt concentration, thus expanding the electrochemical stability window of aqueous electrolytes. In terms of dynamic properties, the calculated and experimentally measured conductivities are in good agreement. Through the analysis of the calculated cation transference number, we propose a three-stage charge carrier transport mechanism with increasing concentration: independent ion transport, strongly correlated ion transport, and small positive charge carrier diffusion through negatively charged polymeric clusters. Our study provides fundamental atomic scale insights into the structure and transport properties of the Zn Cl 2 electrolyte that can aid the optimization and development of WIS electrolytes. Published by the American Physical Society 2025

25 ENERGY STORAGE↗

Excitonic Effects in X-ray Absorption Spectra of Fluoride Salts and Their Surfaces

Given their natural abundance and thermodynamic stability, fluoride salts may appear as evolving components of electrochemical interfaces in Li-ion batteries and emergent multivalent ion cells. This is due to the practice of employing electrolytes with fluorine-containing species (salt, solvent, or additives) that electrochemically decompose and deposit on the electrodes. Operando X-ray absorption spectroscopy (XAS) can probe the electrode–electrolyte interface with a single-digit nanometer depth resolution and offers a wealth of insights into the evolution and Coulombic efficiency or degradation of prototype cells, provided that the spectra can be reliably interpreted in terms of local oxidation state, atomic coordination, and electronic structure about the excited atoms. Here we explore fluorine K-edge XAS of mono- (Li, Na, and K) and di-valent (Mg, Ca, and Zn) fluoride salts from a theoretical standpoint and discover a surprising level of detailed electronic structure information about these materials despite the relatively predictable oxidation state and ionicity of the fluoride anion and the metal cation. Utilizing a recently developed many-body approach based on the ΔSCF method, we calculate the XAS using density functional theory and experimental spectral profiles are well reproduced despite some experimental discrepancies in energy alignment within the literature, which we can correct for in our simulations. We outline a general methodology to explain shifts in the main XAS peak energies in terms of a simple exciton model and explain line-shape differences resulting from the mixing of core-excited states with metal d character (for K and Ca specifically). Given ultimate applications to evolving interfaces, some understanding of the role of surfaces and their terminations in defining new spectral features is provided to indicate the sensitivity of such measurements to changes in interfacial chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A nitrogen- and carbonyl-rich conjugated small-molecule organic cathode for high-performance sodium-ion batteries

Organic-based cathode materials have attracted considerable attention for sustainable Na-ion batteries due to their great promise to overcome the issues arising from the insertion of large Na + into the rigid structures of conventional transition-metal-containing inorganic electrodes. The structural flexibility as a result of weak intermolecular interactions and simple electron- and ion-storage mechanisms found in organic compounds ensure facile and reversible transport of Na ions. In order to address the general drawbacks of small-molecule organic electrodes such as material dissolution in commonly used organic electrolytes and poor electronic conductivity, we report the use of a nitrogen- and carbonyl-rich highly extended π-conjugated small molecule, hexaazatrianthranylene (HATA) embedded quinone (HATAQ), as a cathode for sodium-ion batteries. The unique hydrogen bonds between highly functionalized conjugated HATAQ molecules allow supramolecular graphite-like 2D layered arrangements in the solid state which help facilitate the structural stability during long-term cycling and promote charge transfer. The compound delivers a capacity as high as 460 mA h g -1 at 500 mA g -1 and an excellent capacity retention of 99% (~138 mA h g -1 ) after 5000 cycles at an extremely high rate of 60 A g -1 . The reaction kinetics and redox mechanism of the material have been elucidated by several characterization techniques together with density functional theory (DFT) studies. The insights gained in this work could pave the way for ultra-high-performance small-molecule organic cathodes for sustainable energy storage.

25 ENERGY STORAGE↗

Visualizing the SEI formation between lithium metal and solid-state electrolyte

The solid electrolyte interphase (SEI) is regarded as the most important factor affecting the durability of lithium-metal anode in all-solid-state batteries (ASSBs). Despite its significance, the nucleation and growth mechanism of SEI is not yet well understood. Here, we elucidate the thermodynamics and kinetics governing SEI formation at the Li|β-Li 3 PS 4 interface at the atomic scale via thermodynamic phase equilibrium analysis and machine-learning-potential-assisted molecular dynamics (MD) simulations. An accurate moment tensor potential using the machine-learning method is developed for a reactive model of Li|β-Li 3 PS 4 . This potential enabled us to perform large-scale MD simulations with the model size expanded to the experimental dimensions (~40 nm) while maintaining the same level of accuracy as density functional theory calculations. The results reveal a four-stage evolution process at the Li|β-Li 3 PS 4 interface, namely (i) fast ion diffusion, (ii) nucleation, (iii) Li 2 S growth, and (iv) stabilization. Notably, we demonstrate that the SEI can be categorized into crystalline and amorphous regions. The simulated SEI thickness, structure, and composition closely match experimental findings, validating the accuracy of the MD simulations. We further disclose the significant impact of ion diffusion kinetic limitations on the phase formation and crystallization of interfacial products. Furthermore, we shed light on the detailed potential energy (PE) distribution of lithium along the direction perpendicular to the interface. Finally, this information is crucial for better understanding interfacial ion mobility.

25 ENERGY STORAGE↗