Role of Solvent Rearrangement on Mg2+ Solvation Structures in Dimethoxyethane Solutions using Multimodal NMR Analysis

From Bulk to Interface: Solvent Exchange Dynamics and Their Role in Ion Transport and the Interfacial Model of Rechargeable Magnesium Batteries

Multivalent battery chemistries have been explored in response to the increasing demand for high-energy rechargeable batteries utilizing sustainable resources. Solvation structures of working cations have been recognized as a key component in the design of electrolytes; however, most structure-property correlations of metal ions in organic electrolytes usually build upon favorable static solvation structures, often overlooking solvent exchange dynamics. We here report the ion solvation structures and solvent exchange rates of magnesium electrolytes in various solvents by using multimodal nuclear magnetic resonance (NMR) analysis and molecular dynamics/density functional theory (MD/DFT) calculations. These magnesium solvation structures and solvent exchange dynamics are correlated to the combined effects of several physicochemical properties of the solvents. Moreover, Mg2+ transport and interfacial charge transfer efficiency are found to be closely correlated to the solvent exchange rate in the binary electrolytes where the solvent exchange is tunable by the fraction of diluent solvents. Our primary findings are (1) most battery-related solvents undergo ultraslow solvent exchange coordinating to Mg2+ (with time scales ranging from 0.5 μs to 5 ms), (2) the cation transport mechanism is a mixture of vehicular and structural diffusion even at the ultraslow exchange limit (with faster solvent exchange leading to faster cation transport), and (3) an interfacial model wherein organic-rich regions facilitate desolvation and inorganic regions promote Mg2+ transport is consistent with our NMR, electrochemistry, and cryogenic X-ray photoelectron spectroscopy (cryo-XPS) results. This observed ultraslow solvent exchange and its importance for ion transport and interfacial properties necessitate the judicious selection of solvents and informed design of electrolyte blends for multivalent electrolytes.

Nuclear Magnetic Resonance Relaxation Pathways in Electrolytes for Energy Storage

Nuclear Magnetic Resonance (NMR) spin relaxation times have been an instrumental tool in deciphering the local environment of ionic species, the various interactions they engender and the effect of these interactions on their dynamics in conducting media. Of particular importance has been their application in studying the wide range of electrolytes for energy storage, on which this review is based. Here we highlight some of the research carried out on electrolytes in recent years using NMR relaxometry techniques. Specifically, we highlight studies on liquid electrolytes, such as ionic liquids and organic solvents; on semi-solid-state electrolytes, such as ionogels and polymer gels; and on solid electrolytes such as glasses, glass ceramics and polymers. Although this review focuses on a small selection of materials, we believe they demonstrate the breadth of application and the invaluable nature of NMR relaxometry.

Coordination-Dependent Chemical Reactivity of TFSI Anions at a Mg Metal Interface

Charge transfer across the electrode-electrolyte interface is a highly complex and convoluted process involving diverse solvated species with varying structures and compositions. Despite recent advances in in situ and operando interfacial analysis, molecular specific reactivity of solvated species is inaccessible due to a lack of precise control over the interfacial constituents and/or an unclear understanding of their spectroscopic fingerprints. However, such molecular-specific understanding is critical to the rational design of energy-efficient solid-electrolyte interphase layers. We have employed ion soft landing, a versatile and highly controlled method, to prepare well-defined interfaces assembled with selected ions, either as solvated species or as bare ions, with distinguishing molecular precision. Equipped with precise control over interfacial composition, we employed in situ multimodal spectroscopic characterization to unravel the molecular specific reactivity of Mg solvated species comprising (i.e., bis(trifluoromethanesulfonyl)imide, TFSI^-) anions and solvent molecules (i.e., dimethoxyethane, DME/G1) on a Mg metal surface relevant to multivalent Mg batteries. In situ multimodal spectroscopic characterization revealed higher reactivity of the undercoordinated solvated species [Mg-TFSI-G1]⁺ compared to the fully coordinated [Mg-TFSI-(G1)₂]⁺ species or even the bare TFSI^-. These results were corroborated by the computed reaction pathways and energy barriers for decomposition of the TFSI^- within Mg solvated species relative to bare TFSI^-. Finally, we evaluated the TFSI reactivity under electrochemical conditions using Mg(TFSI)₂-DME-based phase-separated electrolytes representing different solvated constituents. Based on our multimodal study, we report a detailed understanding of TFSI^- decomposition processes as part of coordinated solvated species at a Mg-metal anode that will aid the rational design of improved sustainable electrochemical energy technologies.

Tailoring Coordination in Conventional Ether-Based Electrolytes for Reversible Magnesium-Metal Anodes

Rechargeable magnesium (Mg) batteries based on conventional electrolytes are seriously plagued by the formation of the ion-blocking passivation layer on the Mg metal anode. By tracking the Mg²⁺ solvation sheath, this work links the passivation components to the Mg²⁺ -solvents (1,2-dimethoxyethane, DME) coordination and the consequent thermodynamically unstable DME molecules. On this basis, we propose a methodology to tailor solvation coordination by introducing the additive solvent with extreme electron richness. Oxygen atoms in phosphorus-oxygen groups compete with that in carbon-oxygen groups of DME for the coordination with Mg²⁺ , thus softening the solvation sheath deformation. Meanwhile, the organophosphorus molecules in the rearranged solvation sheath decompose on the Mg surface, increasing the Mg²⁺ transport and electrical resistance by three and one orders of magnitude, respectively. Consequently, the symmetric cells exhibit superior cycling performance of over 600 cycles with low polarization.

Understanding the Solvation-Dependent Properties of Cyclic Ether Multivalent Electrolytes Using High-Field NMR and Quantum Chemistry

Efforts to expand the technological capability of batteries have generated increased interest in divalent cationic systems. Electrolytes used for these electrochemical applications often incorporate cyclic ethers as electrolyte solvents; however, the detailed solvation environments within such systems are not well-understood. To foster insights into the solvation structures of such electrolytes, Ca(TFSI)2 and Zn(TFSI)2 dissolved in tetrahydrofuran (THF) and 2-methyl-tetrahydrofuran were investigated through multi-nuclear magnetic resonance spectroscopy (17O, 43Ca, and 67Zn NMR) combined with quantum chemistry modeling of NMR chemical shifts. NMR provides spectroscopic fingerprints that readily couple with quantum chemistry to identify a set of most probable solvation structures based on the best agreement between the theoretically predicted and experimentally measured values of chemical shifts. The multi-nuclear approach significantly enhances confidence that the correct solvation structures are identified due to the required simultaneous agreement between theory and experiment for multiple nuclear spins. Furthermore, quantum chemistry modeling provides a comparison of the solvation cluster formation energetics, allowing further refinement of the preferred solvation structures. It is shown that a range of solvation structures coexist in most of these electrolytes, with significant molecular motion and dynamic exchange among the structures. This level of solvation diversity correlates with the solubility of the electrolyte, with Zn(TFSI)2/THF exhibiting the lowest degree of each. Comparisons of analogous Ca2+ and Zn2+ solvation structures reveal a significant cation size effect that is manifested in significantly reduced cation-solvent bond lengths and thus stronger solvent bonding for Zn2+ relative to Ca2+. The strength of this bonding is further reduced by methylation of the cyclic ether ring. Solvation shells containing anions are energetically preferred in all the studied electrolytes, leading to significant quantities of contact ion pairs and consequently neutrally charged clusters. It is likely that the transport and interfacial de-solvation/re-solvation properties of these electrolytes are directed by these anion interactions. These insights into the detailed solvation structures, cation size, and solvent effects, including the molecular dynamics, are fundamentally important for the rational design of electrolytes in multivalent battery electrolyte systems.

An automated framework for high-throughput predictions of NMR chemical shifts within liquid solutions

Identifying stable speciation in multi-component liquid solutions is fundamentally important to areas from electrochemistry to organic chemistry and biomolecular systems. Here we introduce a fully automated, high-throughput computational framework for the accurate prediction of stable species in liquid solutions by computing the nuclear magnetic resonance (NMR) chemical shifts. The framework automatically extracts and categorizes hundreds of thousands of atomic clusters from classical molecular dynamics simulations, identifies the most stable species in solution and calculates their NMR chemical shifts via density functional theory calculations. Additionally, the framework creates a database of computed chemical shifts for liquid solutions across a wide chemical and parameter space. We compare our computational results to experimental measurements for magnesium bis(trifluoromethanesulfonyl)imide Mg(TFSI)2 salt in dimethoxyethane solvent. Our analysis of the Mg2+ solvation structural evolutions reveals key factors that influence the accuracy of NMR chemical shift predictions in liquid solutions. Furthermore, we show how the framework reduces the performance of over 300 13C and 600 1H density functional theory chemical shift predictions to a single submission procedure.

Amin M, Yamauchi Y, Xu M, Pan L. MoS2 Nanosheets with Expanded Interlayer Spacing for Ultrastable Aqueous Mg-Ion Hybrid Supercapacitor

Aqueous magnesium ion supercapacitors (MISs) have attracted attention due to their safety, low cost and environmental friendliness. However, the cycling stability of MISs is usually not ideal due to magnesium...

Role of a Multivalent Ion-Solvent Interaction on Restricted Mg2+ Diffusion in Dimethoxyethane Electrolytes

The diffusion behavior of Mg²⁺ in electrolytes is not as readily accessible as that from Li⁺ or Na⁺ utilizing PFG NMR, due to the low sensitivity, poor resolution, and rapid relaxation encountered when attempting ²⁵Mg NMR. In MgTFSI₂/DME solutions, "bound" DME (coordinating to Mg²⁺) and "free" DME (bulk) are distinguishable from ¹H NMR. With the exchange rates between them obtained from 2D ¹H EXSY NMR, we can extract the self-diffusivities of free DME and bound DME (which are equal to that of Mg²⁺) before the exchange occurs using PFG diffusion NMR measurements coupled with analytical formulas describing diffusion under two-site exchange. The high activation enthalpy for exhange (65-70 kJ/mol) can be explained by the structural change of bound DME as evidenced by its reduced C-H bond length. Comparison of the diffusion behaviors of Mg²⁺, TFSI^-, DME, and Li⁺ reveals a relative restriction to Mg²⁺ diffusion that is caused by the long-range interaction between Mg²⁺ and solvent molecules, especially those with suppressed motions at high concentrations and low temperatures.

Anion-Assisted Delivery of Multivalent Cations to Inert Electrodes

To understand and control key electrochemical processes-metal plating, corrosion, intercalation, etc.-requires molecular-scale details of the active species at electrochemical interfaces and their mechanisms for desolvation from the electrolyte. Using free energy sampling techniques we reveal the interfacial speciation of divalent cations in ether-based electrolytes and mechanisms for their delivery to an inert graphene electrode interface. Surprisingly, we find that anion solvophobicity drives a high population of anion-containing species to the interface that facilitate the delivery of divalent cations, even to negatively charged electrodes. Our simulations indicate that cation desolvation is greatly facilitated by cation-anion coupling. We propose anion solvophobicity as a molecular-level descriptor for rational design of electrolytes with increased efficiency for electrochemical processes limited by multivalent cation desolvation.

Ion Solvation Engineering: How to Manipulate the Multiplicity of the Coordination Environment of Multivalent Ions

Free energy analysis of solvation structures of free divalent cations, their ion pairs, and neutral aggregates in low dielectric solvents reveals the multiplicity of thermodynamically stable cation solvation configurations and identifies the micro- and macroscopic factors responsible for this phenomenon. Specifically, we show the role of ion-solvent interactions and solvent mixtures in determining the cation solvation free energy landscapes. We show that it is the entropic contribution of solvent degrees of freedom that is responsible for the solvation multiplicity, and the mutual balance between enthalpic and entropic forces or their concerted contributions is what ultimately defines the most stable ion solvation configuration and creates new ones. We show general consequences of ion solvation multiplicity on thermodynamics of complex electrolytes, specifically in the context of homogeneous or interfacial charge transfer. Identified factors and their interplay provide a pathway to formulation of solvation design rules that can be used to control bulk solvation, interfacial chemistry, and charge transfer. Our findings also suggest experimentally testable predictions.