SEARCH · Engineering Papers

Results for “electrochemical simulation”

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 127 records · Page 7

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Understanding the Surprising Ionic Conductivity Maximum in Zn(TFSI) 2 Water/Acetonitrile Mixture Electrolytes

Aqueous electrolytes composed of 0.1 M zinc bis-(trifluoromethyl-sulfonyl)-imide (Zn-(TFSI) 2 ) and acetonitrile (ACN) were studied using combined experimental and simulation techniques. The electrolyte was found to be electrochemically stable when the ACN V% is higher than 74.4. In addition, it was found that the ionic conductivity of the mixed solvent electrolytes changes as a function of ACN composition, and a maximum was observed at 91.7 V% of ACN although the salt concentration is the same. This behavior was qualitatively reproduced by molecular dynamics (MD) simulations. Detailed analyses based on experiments and MD simulations show that at high ACN composition the water network existing in the high water composition solutions breaks. As a result, the screening effect of the solvent weakens and the correlation among ions increases, which causes a decrease in ionic conductivity at high ACN V%. Furthermore, this study provides a fundamental understanding of this complex mixed solvent electrolyte system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Electrification and decarbonization of spent Li-ion batteries purification by using an electrochemical membrane reactor

The expanding electric vehicle market brings with it exponential growth in the use of lithium (Li)-ion batteries (LIB) for which a wave of spent LIB is expected to come within the next 5 to 10 years. Due to the economic and strategic value imbedded within the metals contained in LIB, different recycling technologies, including hydrometallurgy, pyrometallurgy and direct recycling, are under development. Being different from previous hydrometallurgical methods, which may have high chemical consumption and negative environmental impact, an electrochemical membrane reactor is designed and validated for the first time, to electrify and decarbonize the impurity removal process. This reactor electroplates copper (Cu) and electrochemically precipitates aluminum (Al) and iron (Fe) from simulated spent LIB leachates, by consuming only air, water, and electricity, and the impurities are reduced to <1 ppm. The purified leachate maintains 99.5 % of the nickel (Ni), 95.4 % of the cobalt (Co) and 99.14 % of manganese (Mn) from the original leachate solution, and then can be directly applied for cathode precursor synthesis. Additionally, the purification process doesn’t introduce extra impurity, and the reactor restoration process generates valuable by-product hydro sulfate (H2SO4). This electrochemical process can reduce the cost, because of the much less chemical consumption and the valuable by-product generation, and mitigates the waste emissions, because of no extra impurity introduced and no greenhouse gas (GHG) produced. In conclusion, the chemical precipitation method uses significant amount of NaOH, which induced GHG emission during the manufacturing process.

25 ENERGY STORAGE↗

DOE OSTI2022

Numerical Simulation and Modeling of Hydrogen Gas Evolution on Planar and Microwire Array Electrodes

The impact of gas evolution on the electrochemical characteristics of planar electrodes and microwire array electrodes has been analyzed using modeling and simulation. The impacts can mainly be broken into three phenomena: a) a shift in the local reversible hydrogen electrode potential; b) hyperpolarization; and c) an increase in the solution resistance of the electrolyte. The local reversible hydrogen electrode potential shift was found to play the most important role, constituting >40% of the total potential drop between the cathode and reference electrode, following correction for cell resistance. Compared to planar electrodes, a microwire array structure reduces the impact of bubbles on the solution conductance, but the shift in the local reversible hydrogen electrode potential varies with distance from the actual electrode surface.

25 ENERGY STORAGE↗

DOE OSTI2024

CO 2 Conversion to Butene via a Tandem Photovoltaic–Electrochemical/Photothermocatalytic Process: A Co-design Approach to Coupled Microenvironments

Here, we developed a tandem, unassisted, solar-driven electrochemical and photothermocatalytic process for the single-pass conversion of CO 2 to butene using only simulated solar irradiation as the energetic input. The two-step process involves electrochemical CO 2 reduction (CO 2 R) to ethylene followed by ethylene dimerization to butene. We assessed two unassisted electrochemical setups to concentrate ethylene in the CO 2 R reactor, achieving concentrations up to 5.4 vol.% with 1.8% average solar-to-ethylene conversion and 5.6% average CO 2 -to-ethylene single-pass conversion under 1-sun illumination. When passed through the photothermocatalytic ethylene oligomerization reactor, we generated 600 ppm of butene under 3-sun illumination. Through analysis of this process, we identified that the presence of H 2 , CO, and H 2 O leads to rapid deactivation of the Ni-based ethylene oligomerization catalyst.

14 SOLAR ENERGY↗

DOE OSTI2025

Bubble Transport through a Porous Lattice with an Applied Inlet Flow

Within gas-evolving electrochemical systems, bubbles negatively impact performance by covering electrode active sites for reactions, blocking electric field lines, and obstructing liquid electrolyte flow causing pressure buildup. Recent additive manufacturing advances have enabled tuned porous electrode microstructures to be created, but producing systems that maximize electrochemical throughput and minimize bubble impact remains challenging. Thus, improved physical understanding of and modeling capabilities for bubble behavior are critical to improve electrolyzer design. To address this need, this study examines rising stage bubbles within a lattice with an applied liquid flow—an underexplored regime that strongly influences an electrochemical bubble’s fate. Notably, theoretical predictions and resolved bubble simulations are complemented by experiments from a 3D-printed visualization cell that matches the simulation geometry. The minimum threshold flow rate to achieve bubble breakthrough is found to be larger for higher porosities and for smaller bubbles. Different-sized bubbles decrease expected electrochemical performance in different ways; smaller bubbles tend to stay stuck but cover less solid surface, while larger bubbles more readily break through but cover more surface while in the lattice. The bubble trajectory, deformation, and contact area provide insight into these different behaviors. These findings provide design guidelines toward creating more effective electrolyzers.

Guo, Jack [Lawrence Livermore National Laboratory ↗

DOE OSTI2020

Transport Processes in a Li-ion Cell during an Internal Short-Circuit

Internal short-circuit in a lithium-ion cell causes an abrupt increase in cell temperature and triggers subsequent thermal runaway. In this work, we present a detailed electrochemical-thermal model to investigate the physical behavior during an internal short-circuit. Simulations at wide range of heat transfer coefficients and short-circuit resistances are conducted to illustrate electrochemical and thermal behavior under a wide range of conditions. The Joule heating at the shorted region promotes electrochemical reactions nearby, causing in-plane non-uniformity of electrolyte and active material transport. Furthermore, it is found that diffusion in solid active materials plays a significant role at very high shorting currents (~20 C), as electrochemical reactions rate are being controlled progressively more by availability of Li + at the interface, due to limitations in diffusion through the active material with increasing discharge rates. This diffusion limitation causes a drop in available energy, and subsequently a decrease in cell equilibrium potential and the heat generation rate at the location of the short. On the other hand, rapid depletion of lithium concentration in the electrolyte and accumulation on the electrode surface results in highly non-uniform transport properties resulting in higher heat generation rates. Hence, the heating regime shifts from "local heating" to "global heating". Based on the findings, important design parameters for battery safety are discussed.

25 ENERGY STORAGE↗

DOE OSTI2022

Review—Development of Highly Active and Stable Catalyst Supports and Platinum–Free Catalysts for PEM Fuel Cell

Metal-free, metal-containing, and template-assisted metal-containing nitrogen-modified carbon-based catalysts carbon composite supports, and highly active and electrochemically stable hybrid cathode catalysts for oxygen reduction reaction are reviewed in this manuscript. Here, novel procedures are developed for the synthesis of (i) highly catalytically active PGM-free catalyst and stable carbon composite catalyst supports. The carbon composite catalyst supports are engineered with optimized BET surface area and pore size distribution and with a well-defined kinetic and mass transfer region during the reaction and their performance is dicussed in detail in this review. The surface activation results in increasing carbon graphitization and inclusion of non-metallic active sites on the support surface. The USC catalysts exhibit the initial performance of 0.91 V and a maximum power density of 177 mW cm –2 , with well-defined kinetic and mass-transfer regions and ~2.5% H 2 O 2 production. Activated carbon composite support (ACCS) is modified to optimize its kinetic activity and its electrochemical stability shows excellent thermal stability and support stability under simulated start-up/shut-down operating conditions. The stabilities of various supports developed in this study are compared with those of a commercial Pt/C catalyst. The active sites for the ORR identified with electrochemical and physicochemical methods are pyridinic-N and quaternary-N.

25 ENERGY STORAGE↗

DOE OSTI2020

In-situ/operando X-ray absorption spectroscopic investigation of the electrode/electrolyte interface on the molecular scale

A method for performing X-ray absorption spectroscopy (XAS) in the soft X-ray region was developed to investigate the surface-sensitive electron signal of electrode/electrolyte interfaces in common multivalent based organic solvents under in situ/operando conditions. Our approach enables us to probe the molecular-scale structure of electrode interfaces by measuring total electron yield and is suitable for redox systems exhibiting low intrinsic electrochemical current. In-situ F K-edge XAS measurements in a 0.5 M Magnesium bis(trifluoromethanesulfonimide)/2-Methyltetrahydrofuran (Mg(TFSI) 2 /2-MeTHF) electrolyte were carried out to determine the evolution of interfacial species during the electrochemical charging/discharging process. Time-dependent density-functional theory (TD-DFT) simulation indicate that the F K-edge evolution is the result of interfacial chemical environment change driven by electrochemical potential. In addition, we performed the “operando XAS” which runs the CV and collects spectra simultaneously. Furthermore, by using this method, the non-equilibrium state of condensed matter interfaces and interfacial dynamical transient process can be revealed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Removal of As(V) from wastewaters using magnetic iron oxides formed by zero-valent iron electrocoagulation

Electrocoagulation of zero-valent iron has been widely applied to the removal of dissolved arsenic, but the solid-liquid separation of arsenic-containing precipitates remains technically challenging. In this work, zero-valent iron was electrochemically oxidized to magnetic iron oxides for the removal of As(V) from simulated and actual mining wastewaters. The results indicated that lepidocrocite was formed when zero-valent iron was oxidized by dissolved oxygen, but ferrihydrite and green rust were first formed and then transformed to magnetic iron oxides (mainly magnetite and maghemite) in the electrochemical oxidation from 0 to 0.9 V (vs. SCE), which facilitates the adsorption of As(V) and subsequent solid-liquid separation under a magnetic field. In simulated As(V)-containing solution with initial pH 7.0, zero-valent iron was electrochemically oxidized to magnetite and maghemite at 0.6 V (vs. SCE) for 2 h. The As(V) concentration first decreased from 5127.5 to 26.8 μg L -1 with a removal ratio of 99.5%. In actual mining wastewaters, zero-valent iron was electrochemically oxidized to maghemite at 0.6 V (vs. SCE) for 24 h, and the As(V) concentration decreased from 5486.4 to 3.6 μg L -1 with a removal ratio of 99.9%. The removal ratio of As(V) increased slightly with increasing potential, and increased first and then decreased with increasing initial pH. Compared with that of SO 4 2- and NO 3 - , the presence of Cl - significantly enhanced the removal of As(V). This work provides a highly efficient, facile and low-cost technique for the treatment of arsenic-containing wastewaters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

A Segregated Approach for Modeling the Electrochemistry in the 3-D Microstructure of Li-Ion Batteries and Its Acceleration Using Block Preconditioners

Abstract Battery performance is strongly correlated with electrode microstructure. Electrode materials for lithium-ion batteries have complex microstructure geometries that require millions of degrees of freedom to solve the electrochemical system at the microstructure scale. A fast-iterative solver with an appropriate preconditioner is then required to simulate large representative volume in a reasonable time. In this work, a finite element electrochemical model is developed to resolve the concentration and potential within the electrode active materials and the electrolyte domains at the microstructure scale, with an emphasis on numerical stability and scaling performances. The block Gauss-Seidel (BGS) numerical method is implemented because the system of equations within the electrodes is coupled only through the nonlinear Butler–Volmer equation, which governs the electrochemical reaction at the interface between the domains. The best solution strategy found in this work consists of splitting the system into two blocks—one for the concentration and one for the potential field—and then performing block generalized minimal residual preconditioned with algebraic multigrid, using the FEniCS and the Portable, Extensible Toolkit for Scientific Computation libraries. Significant improvements in terms of time to solution (six times faster) and memory usage (halving) are achieved compared with the MUltifrontal Massively Parallel sparse direct Solver. Additionally, BGS experiences decent strong parallel scaling within the electrode domains. Last, the system of equations is modified to specifically address numerical instability induced by electrolyte depletion, which is particularly valuable for simulating fast-charge scenarios relevant for automotive application.

25 ENERGY STORAGE↗

DOE OSTI2025

Cooperative Effects Associated with High Electrolyte Concentrations in Driving the Conversion of CO2 to C2H4 on Copper

Compared to a conventional electrolyte concentration of 1 M HCOOK, the use of a highly concentrated 7.1 M HCOOK electrolyte increases the Faradaic efficiency (FE) ratio of C2H4/CO from 2.2 +- 0.3 to 18.3 +- 4.8 at -1.08 V vs. reversible hydrogen electrode (RHE) on a Cu gas-diffusion electrode. Based on electrochemical analysis and ab initio molecular dynamics (AIMD) simulation, the identity and concentration of the cation and anion play more important roles in controlling the CO2R reaction pathway than the bulk CO2 solubility and the bulk pH of electrolytes. In situ attenuated reflectance surface enhanced infrared absorption spectroscopy (ATR-SEIRAS) suggests that, unlike 1 M HCOOK, the *CO-bridge-binding mode on Cu is dominant in 7.1 M HCOOK electrolyte, which potentially results in less CO release and higher yield of C2H4. This study demonstrates that although we can tailor the electrolyte composition to shift product selectivity, the factors that control the product selectivity are numerous and cannot be distilled down into one correlated property-reactivity relationship.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Spontaneous sodium ion storage behaviors of reduced graphene oxide anodes exceeding 100% Coulombic efficiency by modulated ion solvation

Rechargeable batteries are essential energy storage devices that power portable devices and electrical vehicles throughout the world. In general, it is thought that the electrochemical performance of rechargeable batteries is mostly determined by the electrodes within them and that the electrolyte plays a relatively passive role. However, ion transport and storage can be greatly influenced by the electrolyte solution structure, specifically, ion solvation within the bulk and ion desolvation across the electrode/electrolyte interfaces. Herein, we studied the role of the electrolyte as an active component of electrochemical energy storage devices. We found that with an appropriate electrolyte formulation, ion storage in disordered carbonaceous anode materials can occur spontaneously without externally supplied electrical energy. Reduced graphene oxide (RGO) in an ether-based electrolyte demonstrates ‘spontaneous' ion storage behaviors of adsorbing and inserting the solvated ions utilizing facilitated permeability and wettability of RGO, which results in Coulombic efficiency of ~145% due to additional charging capacity of ~180 mAh g -1 during electrochemical processes. The unexpected spontaneous ion storage behavior was extensively investigated using a combination of electrochemical analyses and diagnostics, advanced characterizations, and computational simulation. In conclusion, we believe the spontaneous ion storage behavior offers a new way to further improve the energy efficiency of practical rechargeable batteries.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

C–C Bond Formation during Electrochemical CO 2 Reduction on Pristine Cu(100) Unlikely to Involve Adsorbed CO at Any Potential

Formation of hydrocarbons containing two or more carbon atoms (C 2+ ) during heterogeneous electrochemical CO and CO 2 reduction (ECOR and ECO 2 R) only occurs, among pure metals, on Cu electrodes. Moreover, the activity and selectivity is facet dependent, with Cu(100) generally preferentially forming ethylene over methane. Previously, we found via quantum-mechanics-based modeling that, unlike standard density functional theory, more accurate correlated wavefunction methods predict that non-electroactive coupling pathways involving two adsorbed COs (*CO) or a *CO and a *COH to form C–C bonds on Cu(100) are kinetically inhibited, with the former also thermodynamically unfavorable. Here, we extend that embedded complete active space second order perturbation theory (ECASPT2) study, further showing that electrochemical coupling of two *COs to form an anionic dimer [OC*–*CO] (1+δ)– , followed by protonation to form [OC*–*COH] δ− , is not kinetically competitive with the reduction of *CO to *COH at relevant ECO/CO 2 R potentials. Our simulations therefore suggest that the ability of Cu(100) to electrochemically synthesize C 2+ molecules from CO and CO 2 is unlikely to be via *CO, at least on pristine Cu(100). Instead, hydrogenated CO species (*COH, *CH x OH, or *CH x ) are most likely to be the key intermediates in C–C bond formation.

Martirez, John Mark P. [Princeton Plasma Physics L↗

DOE OSTI2024

Rational Design of High-Performing Electrodes in Energy Storage Devices

Our approach in this project is to control conductivities of different ions by modulating the pore sizes in MOFs as sulfur hosts in cathodes. To validate our hypothesis that controlling the pore size of MOFs can effectively prevent LiPSs shuttling while allowing sufficient Li ion transport, which will ensure electrochemical cycling performance, we performed molecular dynamics (MD) simulations using a series of UiO MOFs with different linker sizes including LiPSs and electrolyte components (1M LiTFSI in DME/DOL 1:1 volume) in pores. Specifically, UiO-66, UiO-67, and UiO-68 are Zr-based MOFs with a linker length of 5.67 Å, 10.07 Å, and 14.35 Å, respectively, leading to different pore diameters of 8.4 Å, 12.8 Å, and 17.8 Å. We calculated the mean-squared displacements of Li and LiPS from MD simulations at 300 K for 50 ns or over to assess their mobility in MOFs.

25 ENERGY STORAGE↗

DOE OSTI2026

Microstructure, Transport, and Mechanics of Compacted Clay Simulated at the 0.1 μm Scale (1400 Smectite Clay Particles) Using a Coarse-Grained Model with Explicit Counterions

Clay-rich geomaterials play a critical role in many subsurface systems. The macroscale properties of these materials (low permeability, high ionic conductivity, high swelling pressure, etc.) are sensitive to molecular-level adsorption and hydration interactions at clay−water interfaces. Efforts to develop multiscale simulation approaches to predict these properties reveal a scale gap between atomistic simulations (typically limited to systems smaller than 10 nm) and continuum-scale models (which use computational grid elements with dimensions ≳ 10 μm). In this study, we present a coarse-grained (CG) framework that partly bridges this gap by simulating compacted smectite clay assemblages with dimensions of 0.1 μm containing 1,400 clay particles across a range of dry densities (1,050 to 1,850 kg·m −3 ) and Na/Ca counterion compositions (Na fraction ranging from 0.2 to 1). The simulated systems, along with their reconstructed binary three-dimensional pore networks, are used to evaluate the microstructure, pore size distribution, tortuosity, ion diffusivity, and swelling pressure of compacted smectite clay. Results show that our approach captures important features of the mesoscale heterogeneity of compacted clays, including tactoid formation, hierarchical porosity, and anisotropic pore networks. Results also reveal how compaction and counterion composition govern emergent behaviors, including dominant pore sizes, directional transport, and electrochemical response. This work highlights the potential of CG simulations to bridge molecular and continuum scales and to advance geotechnical and environmental applications involving clay-rich materials as well as related nanoporous media such as geopolymers and calcium-silicate-hydrate. However, the results also suggest that accurate prediction of certain microstructural and mechanical properties (e.g., swelling pressure) may require even larger-scale systems on the order of 1 μm.

02 PETROLEUM↗

DOE OSTI2025

MOSCATO Development and Integration in Fiscal Year 2025: Implementation of Multiphase, Multiphysics Modeling Capabilities for Molten Salt Systems

MOSCATO (Molten Salt Chemistry and Transport) is a multiphysics code that provides high-fidelity, coupled simulations of fluid flow, heat transfer, mass transfer, chemistry, electrochemical phenomena, and alloy corrosion for molten salt systems. In FY25, significant developments were made to the code package, enhancing its capabilities for modeling all relevant phenomena within operating moltens salt reactors (MSRs). The developments and activities in FY25 included: 1. Implementation of Level-Set methods to enable modeling of single-bubble behavior in molten salts. In FY25, the Level-Set two-phase flow modeling implementation was improved to simulate single bubble behavior with molten salt media. The large density and viscosity ratios between typical gases and molten salt liquids present challenges for these types of numerical solvers. With enhancements to the pressure projection method, MOSCATO’s Level-Set solver was able to be successfully validated to experiments related to helium bubble rise in stagnant molten salt. The simulated bubble rising velocity showed reasonable good agreement with experimental measurements. The bubble shape and dynamics were also visually compared with experimental snapshots, demonstrating a good qualitative match. 2. Generation of mass transfer correlations for multiphase flow systems. To enable calculations of the tritium transport across the interface between gas bubbles and salt, we modeled high- Schmidt-number mass transfer around a sphere across a broad range of Reynolds numbers. The mesh near the sphere surface was highly refined to resolve steep concentration gradients caused by the low diffusion coefficient. Literature-based mass transfer correlations were compared with the numerical results, and modifications were proposed to improve agreement, particularly at higher Schmidt numbers. These mass transfer correlations were subsequently provided to other national laboratories to help enable high quality mass transfer simulations using lower-order solvers under development within the NEAMS program. 3. Preliminary implementation of a bubbly flow solver. To model bubbly flow in molten salt, we implemented a bubbly flow solver for void fractions less than 5%. To do so, an algebraic relative velocity model that assumes small bubbles with rapid momentum equilibration was added to MOSCATO to compute bubble velocities. Preliminary comparisons with experimental data showed reasonable agreement, and further improvements are underway. 4. Generation of mass transfer correlations for MSRE subchannel The Molten-Salt Reactor Experiment (MSRE) was a landmark historical project that demonstrated the feasibility of molten-salt reactor technology. The MSRE campaign also generated a significant body of experimental data and reports that continue to support molten-salt–related research. In this report, the MSRE core subchannel was used as the reference geometry for a mass transfer study performed with MOSCATO. The geometry and computational mesh were adapted from a previous study, providing adequate resolution for the relatively low Reynolds number in this case. Additional mesh refinement was applied to reach higher Schmidt numbers, enabling the derivation of a reliable mass-transfer correlation for the present scenario. 5. Simulations of oxygen ingressions into molten salt. In the previous fiscal year, we initiated a study to simulate oxygen ingression in stagnant salt. As oxygen enters the salt through its surface, it reacts with Ce 3+ to form solid CeO 2 and other reaction products. To more fully capture the complex diffusion-convection-reaction mechanisms, capabilities for modeling natural convection in the salt vessel were added. These were needed as the flow of the ingressed gas induced flow in the salt caused by surface shear and non-isothermal effects. With these updated physics in place, we were able to successfully reproduce the experimental results for the rate of change of CeCl 3 concentrations versus time. 6. Flow corrosion model validation. In FY24, MOSCATO’s corrosion model was validated against static corrosion experiments. In FY25, this work was extended to a flow corrosion experiment, where FLiNaK salt was driven by natural convection, with initial salt impurities to initiate corrosion. Despite uncertainties in parameters such as elemental diffusion coefficients in the alloy and unknown H + concentrations, the simulations achieved good agreement with experimental results, especially in predicting sample mass losses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2022

An Electrochemical Hydrogen-Looping System for Low-Cost CO 2 Capture from Seawater

Carbon dioxide (CO 2 ) capture from seawater is a viable, long-term solution for addressing legacy carbon pollution on the planet. However, significant challenges to widespread applicability exist in terms of efficiency, cost, and scalability. In this work, we have developed an electrochemical hydrogen-looping (EHL) system that capitalizes on the generation of proton and OH – from seawater. A CO 2 capture efficiency as high as 91% is achieved for our EHL flow cell. We have demonstrated that the electrochemical hydrogen-looping flow cell can remove CO 2 from simulated seawater at a record low electric energy consumption of 660 kWh ton –1 CO 2 , which is a significant reduction of 56% compared with the currently used bipolar membrane electrodialysis technology. In conclusion, the electrochemical hydrogen-looping flow cell with its much lower energy consumption provides an economical approach for direct removal of CO 2 from seawater at practical scale.

Electrical energy↗

DOE OSTI2022

Understanding Solid-Electrolyte Interphase Formation at the Lithium Metal Anode of Lithium-Sulfur Batteries

Lithium-sulfur (Li-S) batteries are one of the most promising chemistries for the next generation of energy storage devices. In a standard Li-S cell, the metallic lithium anode is an essential component of the battery due to its low density, extremely high theoretical specific capacity, and very low negative electrochemical potential. However, several challenges related to the use of Li metal have prevented this battery technology from becoming commercially available. For instance, Li-anodes are highly reactive, which results in the continuous decomposition of the electrolyte and the formation of the solid-electrolyte interphase (SEI) layer. This, combined with the non-uniform deposition of Liions during plating at the anode surface, can seriously affect the performance, cycling, and safety of the battery. In addition, it has been suggested that a controlled SEI formation at the metallic Li anode can yield enhanced battery performance. Hence, a comprehensive understanding of how the SEI layer is formed can help elucidate improvements in this battery technology. In this work, the formation of the SEI layer in Li-S batteries is investigated using density functional theory (DFT) and ab initio molecular dynamics (AIMD) simulations. The stability of liquid organic-based electrolyte components, salt concentration, electron-rich environments, and the use of inorganic solid-state electrolyte (SSE) materials are explored in order to provide molecular-level fundamental insights into how the nature and composition of the electrolyte can alter the initial stages of the SEI formation. The products of electrolyte decomposition and electro-/chemical reaction pathways of liquid electrolyte component such as solvent and salts are predicted from AIMD simulations of electrolyte-anode interfaces. Energetics of reactions from DFT ground-state optimizations are also presented to confirm the decomposition mechanisms. Some significant differences are then drawn regarding the use of low and high salt concentrations. Simulations under electron-rich environments show additional multielectron electrochemical reactions of solvent and salt decomposition taking place due to the excess of electrons and the presence of radical anions in the solution. Finally, the initial stages of formation, stability, and main constituents of interfaces between Li-metal and Sbased electrodes with sulfide-based SSEs are characterized in detail.

25 ENERGY STORAGE↗