Interfacial solvation-structure regulation for stable Li metal anode by a desolvation coating technique.

Rechargeable lithium (Li) metal batteries face challenges in achieving stable cycling due to the instability of the solid electrolyte interphase (SEI). The Li-ion solvation structure and its desolvation process are crucial for the formation of a stable SEI on Li metal anodes and improving Li plating/stripping kinetics. This research introduces an interfacial desolvation coating technique to actively modulate the Li-ion solvation structure at the Li metal interface and regulate the participation of the electrolyte solvent in SEI formation. Through experimental investigations conducted using a carbonate electrolyte with limited compatibility to Li metal, the optimized desolvation coating layer, composed of 12-crown-4 ether-modified silica materials, selectively displaces strongly coordinating solvents while simultaneously enriching weakly coordinating fluorinated solvents at the Li metal/electrolyte interface. This selective desolvation and enrichment effect reduce solvent participation to SEI and thus facilitate the formation of a LiF-dominant SEI with greatly reduced organic species on the Li metal surface, as conclusively verified through various characterization techniques including XPS, quantitative NMR, operando NMR, cryo-TEM, EELS, and EDS. The interfacial desolvation coating technique enables excellent rate cycling stability (i.e., 1C) of the Li metal anode and prolonged cycling life of the Li||LiCoO2 pouch cell in the conventional carbonate electrolyte (E/C 2.6 g/Ah), with 80% capacity retention after 333 cycles.

The local structure in the BmimPF6/acetonitrile mixture: the charge distribution effect.

The changes of the local structure in the binary mixture of 1-butyl-3-methylimidazolium hexafluorophosphate (BmimPF6) ionic liquid and acetonitrile are investigated over the entire composition range. Two charge distribution models of the ions are considered: in the first one, the atomic fractional charges of the cations and anions are kept equal with those in the neat ionic liquid, and hence they are independent from the mole fraction of the ionic liquid, while in the second one the charge distribution is scaled up by a mole fraction dependent factor. The sum of these charges converge to +1e and -1e on the cation and anion, respectively, at infinite dilution. All the other interactions and geometry parameters of the ions (i.e., Lennard-Jones, bond stretching, angle bending and dihedral parameters) are identical in the two cases. The effect of the fractional charge distribution on the hydrogen bonding between the counterions themselves and between the ions and solvent molecules, as well as on the stacking interactions between the cations, is analyzed. To this end, two distances, characteristic of the hydrogen bond formed by the donor moiety and its nearest neighbor acceptor, as well as a coordinate system that defines unambiguously the orientation between a reference cation and its nearest neighbor, are introduced. It is shown that, with the variable charge model, the neighboring cation-anion pairs maintain their relative arrangement similar to the neat ionic liquid down to an ionic liquid mole fraction of xIL = 0.10, whereas in the case of the constant charge model such changes occur already at xIL = 0.20. Furthermore, the analysis of the first and the second nearest neighbor distance distributions of an anion around a reference cation indicates that, at this mole fraction range, there are two different preferred arrangements of the anions around the cations. In the first one, similarly to the local structure around a reference cation in the neat ionic liquid, the anion forms a distorted hydrogen bond with the cation, while in the second one the anion is located farther from the cation, forming no hydrogen bond with it. The relative population of these two types of preferred nearest neighbor cation-anion arrangements is found to be sensitive to further decrease of the ionic liquid mole fraction. These findings correlate with experimental results concerning the behavior of many physical chemical properties (e.g., excess volume, excess viscosity, chemical shift, infrared and Raman vibrational mode shifts, diffusion, etc.) that were found to undergo a drastic change in this mole fraction range. Our results show that in this composition range a transition occurs from the situation where the macroscopic physical chemical properties of the system are determined primarily by the cation-anion hydrogen bonding interactions to that where they are determined by the solvation of the cation and the anion by the molecular solvent.

Revealing the Anion-Solvent Interaction for Ultralow Temperature Lithium Metal Batteries.

Anion solvation in electrolytes can largely change the electrochemical performance of the electrolytes, yet has been rarely investigated. Herein, three anions of bis(trifluoromethanesulfonyl)imide (TFSI), bis(fluorosulfonyl)imide (FSI), and derived asymmetric (fluorosulfonyl)(trifluoro-methanesulfonyl)imide (FTFSI) are systematically examined in a weakly Li+ cation solvating solvent of bis(3-fluoropropyl)ether (BFPE). In-situ liquid secondary ion mass spectrometry demonstrates that FTFSI- and FSI- anions are associated with BFPE solvent, while weak TFSI- /BFPE cluster signals are detected. Molecular modeling further reveals that the anion-solvent interaction is accompanied by the formation of H-bonding-like interactions. Anion solvation enhances the Li+ cation transfer number and reduces the organic component in solid electrolyte interphase, which enhances the Li plating/stripping Coulombic efficiency at a low temperature of -30 °C from 42.4% in TFSI-based electrolytes to 98.7% in 1.5 m LiFTFSI and 97.9% in LiFSI-BFPE electrolytes. The anion-solvent interactions, especially asymmetric anion solvation also accelerate the Li+ desolvation kinetics. The 1.5 m LiFTFSI-BFPE electrolyte with strong anion-solvent interaction enables LiNi0.8 Mn0.1 Co0.1 O2 (NMC811)||Li (20 µm) full cell with stable cyclability even under -40 °C, retaining over 92% of initial capacity (115 mAh g-1 , after 100 cycles). The anion-solvent interactions insights allow to rational design the electrolyte for lithium metal batteries and beyond to achieve high performance.

Insight to the Local Structure of Mixtures of Imidazolium-Based Ionic Liquids and Molecular Solvents from Molecular Dynamics Simulations and Voronoi Analysis.

While the physicochemical properties as well as the NMR and vibration spectroscopic data of the mixtures of ionic liquids (ILs) with molecular solvents undergo a drastic change around the IL mole fraction of 0.2, the local structure of the mixtures pertaining to this behavior remains unclear. In this work, the local structure of 12 mixtures of 1-butyl-3-methylimidazolium cation (C4mim+) combined with perfluorinated anions, such as tetrafluoroborate (BF4-), hexafluorophosphate (PF6-), trifluoromethylsulfonate (TFO-), and bis(trifluoromethanesulfonyl)imide, (TFSI-), and aprotic dipolar solvents, such as acetonitrile (AN), propylene carbonate (PC), and gamma butyrolactone (γ-BL) is studied by molecular dynamics simulations in the entire composition range, with an emphasis on the IL mole fractions around 0.2. Distributions of metric properties corresponding to the Voronoi polyhedra of the particles (volume assigned to the particles, local density, radius of spherical voids) are determined, using representative sites of the cations, anions, and the solvent molecules, to characterize the changes in the local structure of these mixtures. By analyzing the mole fraction dependence of the average value, fluctuation, and skewness parameter of these distributions, the present study reveals that, around the IL mole fraction of 0.2, the local structure of the mixture undergoes a transition between that determined by the interionic interactions and that determined by the interactions between the ions and solvent molecules. It should be noted that the strength of the interactions between the ions and the solvent molecules, modulated by the change in the composition of the mixture, plays an important role in the occurrence of this transition. The signature of the change in the local structure is traced back to the nonlinear change of the mean values, fluctuations, and skewness values of the metric Voronoi polyhedra distributions.

Local Structure of DMF-Water Mixtures, as Seen from Computer Simulations and Voronoi Analysis.

Molecular dynamics simulations of mixtures of N,N-dimethylformamide (DMF) with water of various compositions, covering the entire composition range, are performed on the canonical (N,V,T) ensemble. The local structure of the mixtures is analyzed in terms of radial distribution functions and the contributions of the first five neighbors to them, various order parameters of the water molecules around each other, and properties of the Voronoi polyhedra of the molecules. The analyses lead to the following main conclusions. The two molecules are mixing with each other even on the molecular scale; however, small self-aggregates of both components persist even at their small mole fraction values. In particular, water-water H-bonds exist in the entire composition range, while water clusters larger than 3 and 2 molecules disappear above the DMF mole fraction values of about 0.7 and 0.9, respectively. The O atoms of the DMF molecules can well replace water O atoms in the hydrogen-bonding network. Further, the H-bonding structure is enhanced by the presence of the hydrophobic CH3 groups of the DMF molecules. On the other hand, the H-bonding network of the molecules gradually breaks down upon the addition of DMF to the system due to the lack of H-donating groups of the DMF molecules. Finally, in neat DMF, the molecules form weak, CH-donated H-bonds with each other; however, these H-bonds disappear upon the addition of water due to the increasing competition with the considerably stronger OH-donated H-bonds DMF can form with the water molecules.

Voronoi Polyhedra as a Tool for the Characterization of Inhomogeneous Distribution in 1-Butyl-3-methylimidazolium Cation-Based Ionic Liquids.

The inhomogeneity distribution in four imidazolium-based ionic liquids (ILs) containing the 1-butyl-3-methylimidazolium (C4mim) cation, coupled with tetrafluoroborate (BF4), hexafluorophosphate (PF6), bis(trifluoromethanesulfonyl)amide (TFSA), and trifluoromethanesulfonate (TfO) anions, was characterized using Voronoi polyhedra. For this purpose, molecular dynamic simulations have been performed on the isothermal-isobaric (NpT) ensemble. We checked the ability of the potential models to reproduce the experimental density, heat of vaporization, and transport properties (diffusion and viscosity) of these ionic liquids. The inhomogeneity distribution of ions around the ring, methyl, and butyl chain terminal hydrogen atoms of the C4mim cation was investigated by means of Voronoi polyhedra analysis. For this purpose, the position of the C4mim cation was described successively by the ring, methyl, and butyl chain terminal hydrogen atoms, while that of the anions was described by their F or O atom. We calculated the Voronoi polyhedra distributions of the volume, the density, and the asphericity parameters as well as that of the radius of the spherical intermolecular voids. We carried out the analysis in two steps. In the first step, both ions were taken into account. The calculated distributions gave information on the neighboring ions around a reference one. In the second step, to distinguish between like and oppositely charged ions and then to get information on the inhomogeneity distribution of the like ions, we repeated the same calculations on the same sample configurations and removed one of the ions and considered only the other one. Detailed analysis of these distributions has revealed that the ring hydrogen atoms are mainly solvated by the anions, while the methyl and butyl terminal H atoms are surrounded by like atoms. The extent of this inhomogeneity was assessed by calculating the cluster size distribution that shows that the dimers are the most abundant ones.

Distance Angle Descriptors of the Interionic and Ion-Solvent Interactions in Imidazolium-Based Ionic Liquid Mixtures with Aprotic Solvents: A Molecular Dynamics Simulation Study.

The aim of this paper is to quantify the changes of the interionic and ion-solvent interactions in mixtures of imidazolium-based ionic liquids, having tetrafluoroborate (BmimBF4), hexafluorophosphate (BmimPF6), trifluoromethylsulfonate (BmimTFO), or bis(trifluoromethanesulfonyl)imide (BmimTFSI), anions, and polar aprotic molecular solvents, such as acetonitrile (AN), γ-butyrolactone (GBL), and propylene carbonate (PC). For this purpose, we calculate, using the nearest-neighbor approach, the average distance between the imidazolium ring H atom in positions 2, 4, and 5 (H2,4,5) and the nearest high-electronegativity atom of the solvent or anion (X) as distance descriptors, and the mean angle formed by the C2,4,5-H2,4,5 bond and the H2,4,5···X axis around the H2,4,5 atom as angular descriptors of the cation-anion and cation-solvent interactions around the ring C-H groups. The behavior of these descriptors as a function of the ionic liquid mole fraction is analyzed in detail. The obtained results show that the extent of the change of these descriptors with respect to their values in the neat ionic liquid depends both on the nature of the anion and on the mixture composition. Thus, in the case of the mixtures of the molecular solvents with BmimBF4 and BmimTFO, a small change of the distance and a drastic increase of the angular descriptor corresponding to the cation-anion interactions are observed with decreasing mole fraction of the ionic liquid, indicating that the anion moves from the above/below position (with respect to the imidazolium ring plane) to a position that is nearly linearly aligned with the C2-H2 bond and hinders the possible interaction between the C2-H2 group and the solvent molecules. On the other hand, in the case of mixtures of BmimTFSI and BmimPF6 with the molecular solvents, both the observed increase of the distance descriptor and the slight change of the angular descriptor with decreasing ionic liquid mole fraction are compatible with the direct interactions of the solvent with the C2-H2 group. The behavior of these descriptors is correlated with the experimentally observed 1H chemical shift of the C2-H2 group and the red shift of the C2-H2 vibrational mode, particularly at low ionic liquid mole fractions. The present results are thus of great help in interpreting these experimental observations.

New Force Field Model for Propylene Glycol: Insight to Local Structure and Dynamics.

In this work we developed a new force field model (FFM) for propylene glycol (PG) based on the OPLS all-atom potential. The OPLS potential was refined using quantum chemical calculations, taking into account the densities and self-diffusion coefficients. The validation of this new FFM was carried out based on a wide range of physicochemical properties, such as density, enthalpy of vaporization, self-diffusion coefficients, isothermal compressibility, surface tension, and shear viscosity. The molecular dynamics (MD) simulations were performed over a large range of temperatures (293.15-373.15 K). The comparison with other force field models, such as OPLS, CHARMM27, and GAFF, revealed a large improvement of the results, allowing a better agreement with experimental data. Specific structural properties (radial distribution functions, hydrogen bonding and spatial distribution functions) were then analyzed in order to support the adequacy of the proposed FFM. Pure propylene glycol forms a continuous phase, displaying no microstructures. It is shown that the developed FFM gives rise to suitable results not only for pure propylene glycol but also for mixtures by testing its behavior for a 50 mol % aqueous propylene glycol solution. Furthermore, it is demonstrated that the addition of water to the PG phase produces a homogeneous solution and that the hydration interactions prevail over the propylene glycol self-association interactions.

Local Structure in Terms of Nearest-Neighbor Approach in 1-Butyl-3-methylimidazolium-Based Ionic Liquids: MD Simulations.

Description of the local microscopic structure in ionic liquids (ILs) is a prerequisite to obtain a comprehensive understanding of the influence of the nature of ions on the properties of ILs. The local structure is mainly determined by the spatial arrangement of the nearest neighboring ions. Therefore, the main interaction patterns in ILs, such as cation-anion H-bond-like motifs, cation-cation alkyl tail aggregation, and ring stacking, were considered within the framework of the nearest-neighbor approach with respect to each particular interaction site. We employed classical molecular dynamics (MD) simulations to study in detail the spatial, radial, and orientational relative distribution of ions in a set of imidazolium-based ILs, in which the 1-butyl-3-methylimidazolium (C4mim(+)) cation is coupled with the acetate (OAc(-)), chloride (Cl(-)), tetrafluoroborate (BF4(-)), hexafluorophosphate (PF6(-)), trifluoromethanesulfonate (TfO(-)), or bis(trifluoromethanesulfonyl)amide (TFSA(-)) anion. It was established that several structural properties are strongly anion-specific, while some can be treated as universally applicable to ILs, regardless of the nature of the anion. Namely, strongly basic anions, such as OAc(-) and Cl(-), prefer to be located in the imidazolium ring plane next to the C-H(2/4-5) sites. By contrast, the other four bulky and weakly coordinating anions tend to occupy positions above/below the plane. Similarly, the H-bond-like interactions involving the H(2) site are found to be particularly enhanced in comparison with the ones at H(4-5) in the case of asymmetric and/or more basic anions (C4mimOAc, C4mimCl, C4mimTfO, and C4mimTFSA), in accordance with recent spectroscopic and theoretical findings. Other IL-specific details related to the multiple H-bond-like binding and cation stacking issues are also discussed in this paper. The secondary H-bonding of anions with the alkyl hydrogen atoms of cations as well as the cation-cation alkyl chain aggregation turned out to be poorly sensitive to the nature of the anion.