Solubilities of Ethylene and Carbon Dioxide Gases in Lithium-Ion Battery Electrolyte.

During Li-ion battery operation, (electro)chemical side reactions occur within the cell that can promote or degrade performance. These complex reactions produce byproducts in the solid, liquid, and gas phases. Studying byproducts in these three phases can help optimize battery lifetimes. To relate the measured gas-phase byproducts to species dissolved in the liquid-phase, equilibrium proprieties such as the Henry's law constants are required. The present work implements a pressure decay experiment to determine the thermodynamic equilibrium concentrations between the gas and liquid phases for ethylene (C2H4) and carbon dioxide (CO2), which are two gases commonly produced in Li-ion batteries, with an electrolyte of 1.2 M LiPF6 in 3:7 wt/wt ethylene carbonate/ethyl methyl carbonate and 3 wt % fluoroethylene carbonate (15:25:57:3 wt % total composition). The experimentally measured pressure decay curve is fit to an analytical dissolution model and extrapolated to predict the final pressure at equilibrium. The relationship between the partial pressures and concentration of dissolved gas in electrolyte at equilibrium is then used to determine Henry's law constants of 2.0 × 104 kPa for C2H4 and k CO2 = 1.1 × 104 kPa for CO2. These values are compared to Henry's law constants predicted from density functional theory and show good agreement within a factor of 3.

Design Principles and Routes for Calcium Alkoxyaluminate Electrolytes.

Multivalent-ion battery technologies are increasingly attractive options for meeting diverse energy storage needs. Calcium ion batteries (CIB) are particularly appealing candidates for their earthly abundance, high theoretical volumetric energy density, and relative safety advantages. At present, only a few Ca-ion electrolyte systems are reported to reversibly plate at room temperature: for example, aluminates and borates, including Ca[TPFA]2, where [TPFA]- = [Al(OC(CF3)3)4]- and Ca[B(hfip)4]2, [B(hfip)4]2- = [B(OCH(CF3)2)4]-. Analyzing the structure of these salts reveals a common theme: the prevalent use of a weakly coordinating anion (WCA) consisting of a tetracoordinate aluminum/boron (Al/B) center with fluorinated alkoxides. Leveraging the concept of theory-aided design, we report an innovative, one-pot synthesis of two new calcium-ion electrolyte salts (Ca[Al(tftb)4]2, Ca[Al(hftb)4]2) and two reported salts (Ca[Al(hfip)4]2 and Ca[TPFA]2) where hfip = (-OCH(CF3)2), tftb = (-OC(CF3)(Me)2), hftb = (-OC(CF3)2(Me)), [TPFA]- = [Al(OC(CF3)3)4]-. We also reveal the dependence of Coulombic efficiency on their inherent propensity for cation-anion coordination.

CoeffNet: predicting activation barriers through a chemically-interpretable, equivariant and physically constrained graph neural network.

Activation barriers of elementary reactions are essential to predict molecular reaction mechanisms and kinetics. However, computing these energy barriers by identifying transition states with electronic structure methods (e.g., density functional theory) can be time-consuming and computationally expensive. In this work, we introduce CoeffNet, an equivariant graph neural network that predicts activation barriers using coefficients of any frontier molecular orbital (such as the highest occupied molecular orbital) of reactant and product complexes as graph node features. We show that using coefficients as features offer several advantages, such as chemical interpretability and physical constraints on the network's behaviour and numerical range. Model outputs are either activation barriers or coefficients of the chosen molecular orbital of the transition state; the latter quantity allows us to interpret the results of the neural network through chemical intuition. We test CoeffNet on a dataset of SN2 reactions as a proof-of-concept and show that the activation barriers are predicted with a mean absolute error of less than 0.025 eV. The highest occupied molecular orbital of the transition state is visualized and the distribution of the orbital densities of the transition states is described for a few prototype SN2 reactions.

A Critical Analysis of Chemical and Electrochemical Oxidation Mechanisms in Li-Ion Batteries.

Electrolyte decomposition limits the lifetime of commercial lithium-ion batteries (LIBs) and slows the adoption of next-generation energy storage technologies. A fundamental understanding of electrolyte degradation is critical to rationally design stable and energy-dense LIBs. To date, most explanations for electrolyte decomposition at LIB positive electrodes have relied on ethylene carbonate (EC) being chemically oxidized by evolved singlet oxygen (1O2) or electrochemically oxidized. In this work, we apply density functional theory to assess the feasibility of these mechanisms. We find that electrochemical oxidation is unfavorable at any potential reached during normal LIB operation, and we predict that previously reported reactions between the EC and 1O2 are kinetically limited at room temperature. Our calculations suggest an alternative mechanism in which EC reacts with superoxide (O2-) and/or peroxide (O22-) anions. This work provides a new perspective on LIB electrolyte decomposition and motivates further studies to understand the reactivity at positive electrodes.

Chemical Reaction Networks Explain Gas Evolution Mechanisms in Mg-Ion Batteries.

Out-of-equilibrium electrochemical reaction mechanisms are notoriously difficult to characterize. However, such reactions are critical for a range of technological applications. For instance, in metal-ion batteries, spontaneous electrolyte degradation controls electrode passivation and battery cycle life. Here, to improve our ability to elucidate electrochemical reactivity, we for the first time combine computational chemical reaction network (CRN) analysis based on density functional theory (DFT) and differential electrochemical mass spectroscopy (DEMS) to study gas evolution from a model Mg-ion battery electrolyte─magnesium bistriflimide (Mg(TFSI)2) dissolved in diglyme (G2). Automated CRN analysis allows for the facile interpretation of DEMS data, revealing H2O, C2H4, and CH3OH as major products of G2 decomposition. These findings are further explained by identifying elementary mechanisms using DFT. While TFSI- is reactive at Mg electrodes, we find that it does not meaningfully contribute to gas evolution. The combined theoretical-experimental approach developed here provides a means to effectively predict electrolyte decomposition products and pathways when initially unknown.

Assessing the Accuracy of Density Functional Approximations for Predicting Hydrolysis Reaction Kinetics.

Hydrolysis reactions are ubiquitous in biological, environmental, and industrial chemistry. Density functional theory (DFT) is commonly employed to study the kinetics and reaction mechanisms of hydrolysis processes. Here, we present a new data set, Barrier Heights for HydrOlysis - 36 (BH2O-36), to enable the design of density functional approximations (DFAs) and the rational selection of DFAs for applications in aqueous chemistry. BH2O-36 consists of 36 diverse organic and inorganic forward and reverse hydrolysis reactions with reference energy barriers ΔE‡ calculated at the CCSD(T)/CBS level. Using BH2O-36, we evaluate 63 DFAs. In terms of mean absolute error (MAE) and mean relative absolute error (MRAE), ωB97M-V is the best-performing DFA tested, while MN12-L-D3(BJ) is the best-performing pure (nonhybrid) DFA. Broadly, we find that range-separated hybrid DFAs are necessary to approach chemical accuracy (0.043 eV). Although the best-performing DFAs include a dispersion correction to account for long-range interactions, we find that dispersion corrections do not generally improve MAE or MRAE for this data set.

Chemical reaction networks and opportunities for machine learning.

Chemical reaction networks (CRNs), defined by sets of species and possible reactions between them, are widely used to interrogate chemical systems. To capture increasingly complex phenomena, CRNs can be leveraged alongside data-driven methods and machine learning (ML). In this Perspective, we assess the diverse strategies available for CRN construction and analysis in pursuit of a wide range of scientific goals, discuss ML techniques currently being applied to CRNs and outline future CRN-ML approaches, presenting scientific and technical challenges to overcome.

Quantum chemical calculations of lithium-ion battery electrolyte and interphase species.

Lithium-ion batteries (LIBs) represent the state of the art in high-density energy storage. To further advance LIB technology, a fundamental understanding of the underlying chemical processes is required. In particular, the decomposition of electrolyte species and associated formation of the solid electrolyte interphase (SEI) is critical for LIB performance. However, SEI formation is poorly understood, in part due to insufficient exploration of the vast reactive space. The Lithium-Ion Battery Electrolyte (LIBE) dataset reported here aims to provide accurate first-principles data to improve the understanding of SEI species and associated reactions. The dataset was generated by fragmenting a set of principal molecules, including solvents, salts, and SEI products, and then selectively recombining a subset of the fragments. All candidate molecules were analyzed at the ωB97X-V/def2-TZVPPD/SMD level of theory at various charges and spin multiplicities. In total, LIBE contains structural, thermodynamic, and vibrational information on over 17,000 unique species. In addition to studies of reactivity in LIBs, this dataset may prove useful for machine learning of molecular and reaction properties.

A chemically consistent graph architecture for massive reaction networks applied to solid-electrolyte interphase formation.

Modeling reactivity with chemical reaction networks could yield fundamental mechanistic understanding that would expedite the development of processes and technologies for energy storage, medicine, catalysis, and more. Thus far, reaction networks have been limited in size by chemically inconsistent graph representations of multi-reactant reactions (e.g. A + B → C) that cannot enforce stoichiometric constraints, precluding the use of optimized shortest-path algorithms. Here, we report a chemically consistent graph architecture that overcomes these limitations using a novel multi-reactant representation and iterative cost-solving procedure. Our approach enables the identification of all low-cost pathways to desired products in massive reaction networks containing reactions of any stoichiometry, allowing for the investigation of vastly more complex systems than previously possible. Leveraging our architecture, we construct the first ever electrochemical reaction network from first-principles thermodynamic calculations to describe the formation of the Li-ion solid electrolyte interphase (SEI), which is critical for passivation of the negative electrode. Using this network comprised of nearly 6000 species and 4.5 million reactions, we interrogate the formation of a key SEI component, lithium ethylene dicarbonate. We automatically identify previously proposed mechanisms as well as multiple novel pathways containing counter-intuitive reactions that have not, to our knowledge, been reported in the literature. We envision that our framework and data-driven methodology will facilitate efforts to engineer the composition-related properties of the SEI - or of any complex chemical process - through selective control of reactivity.

Aqueous Diels-Alder reactions for thermochemical storage and heat transfer fluids identified using density functional theory.

Thermal storage and transfer fluids have important applications in industrial, transportation, and domestic settings. Current thermal fluids have relatively low specific heats, often significantly below that of water. However, by introducing a thermochemical reaction to a base fluid, it is possible to enhance the fluid's thermal properties. In this work, density functional theory (DFT) is used to screen Diels-Alder reactions for use in aqueous thermal fluids. From an initial set of 52 reactions, four are identified with moderate aqueous solubility and predicted turning temperature near the liquid region of water. These reactions are selectively modified through 60 total functional group substitutions to produce novel reactions with improved solubility and thermal properties. Among the reactions generated by functional group substitution, seven have promising predicted thermal properties, significantly improving specific heat (by as much as 30.5%) and energy storage density (by as much as 4.9%) compared to pure water.

BonDNet: a graph neural network for the prediction of bond dissociation energies for charged molecules.

A broad collection of technologies, including e.g. drug metabolism, biofuel combustion, photochemical decontamination of water, and interfacial passivation in energy production/storage systems rely on chemical processes that involve bond-breaking molecular reactions. In this context, a fundamental thermodynamic property of interest is the bond dissociation energy (BDE) which measures the strength of a chemical bond. Fast and accurate prediction of BDEs for arbitrary molecules would lay the groundwork for data-driven projections of complex reaction cascades and hence a deeper understanding of these critical chemical processes and, ultimately, how to reverse design them. In this paper, we propose a chemically inspired graph neural network machine learning model, BonDNet, for the rapid and accurate prediction of BDEs. BonDNet maps the difference between the molecular representations of the reactants and products to the reaction BDE. Because of the use of this difference representation and the introduction of global features, including molecular charge, it is the first machine learning model capable of predicting both homolytic and heterolytic BDEs for molecules of any charge. To test the model, we have constructed a dataset of both homolytic and heterolytic BDEs for neutral and charged (-1 and +1) molecules. BonDNet achieves a mean absolute error (MAE) of 0.022 eV for unseen test data, significantly below chemical accuracy (0.043 eV). Besides the ability to handle complex bond dissociation reactions that no previous model could consider, BonDNet distinguishes itself even in only predicting homolytic BDEs for neutral molecules; it achieves an MAE of 0.020 eV on the PubChem BDE dataset, a 20% improvement over the previous best performing model. We gain additional insight into the model's predictions by analyzing the patterns in the features representing the molecules and the bond dissociation reactions, which are qualitatively consistent with chemical rules and intuition. BonDNet is just one application of our general approach to representing and learning chemical reactivity, and it could be easily extended to the prediction of other reaction properties in the future.