Local Structure of Li+ in Concentrated Ethylene Carbonate Solutions Studied by Low-Frequency Raman Scattering and Neutron Diffraction with 6Li/7Li Isotopic Substitution Methods

Effects of Me-Solvent Interactions on the Structure and Infrared Spectra of MeTFSI (Me = Li, Na) Solutions in Carbonate Solvents-A Test of the GFN2-xTB Approach in Molecular Dynamics Simulations

We investigated the performance of the computationally effective GFN2-xTB approach in molecular dynamics (MD) simulations of liquid electrolytes for lithium/sodium batteries. The studied systems were LiTFSI and NaTFSI solutions in ethylene carbonate or fluoroethylene carbonate and the neat solvents. We focused on the structure of the electrolytes and on the manifestations of ion-solvent interactions in the vibrational spectra. The IR spectra were calculated from MD trajectories as Fourier transforms of the dipole moment. The results were compared to the data obtained from ab initio MD. The spectral shifts of the carbonyl stretching mode calculated from the GFN2-xTB simulations were in satisfactory agreement with the ab initio MD data and the experimental results for similar systems. The performance in the region of molecular ring vibrations was significantly worse. We also found some differences in structural data, suggesting that the GFN2-xTB overestimates interactions of Me ions with TFSI anions and Na+ binding to solvent molecules. We conclude that the GFN2-xTB method is an alternative worth considering for MD simulations of liquids, but it requires testing of its applicability for new systems.

Extension of the TraPPE Force Field for Battery Electrolyte Solvents

Optimizing electrolyte formulations is key to improving performance of Li-/Na-ion batteries, where transport properties (diffusion coefficient, viscosity) and permittivity need to be predicted as functions of temperature, salt concentration and solvent composition. More efficient and reliable simulation models are urgently needed, owing to the high cost of experimental methods and the lack of united-atom molecular dynamics force fields validated for electrolyte solvents. Here the computationally efficient TraPPE united-atom force field is extended to be compatible with carbonate solvents, optimizing the charges and dihedral potential. Computing the properties of electrolyte solvents, ethylene carbonate (EC), propylene carbonate (PC), dimethyl carbonate (DMC), diethyl carbonate (DEC), and dimethoxyethane (DME), we observe that the average absolute errors in the density, self-diffusion coefficient, permittivity, viscosity, and surface tension are approximately 15% of the corresponding experimental values. Results compare favorably to all-atom CHARMM and OPLS-AA force fields, offering computational performance improvement of at least 80%. We further use TraPPE to predict the structure and properties of LiPF₆ salt in these solvents and their mixtures. EC and PC form complete solvation shells around Li⁺ ions, while the salt in DMC forms chain-like structures. In the poorest solvent, DME, LiPF₆ forms globular clusters despite DME's higher permittivity than DMC.

Local Structure of Li+ in Superconcentrated Aqueous LiTFSA Solutions

It has been reported that aqueous lithium ion batteries (ALIBs) can operate beyond the electrochemical window of water by using a superconcentrated electrolyte aqueous solution. The liquid structure, particularly the local structure of the Li⁺, which is rather different from conventional dilute solution, plays a crucial role in realizing the ALIB. To reveal the local structure around Li⁺, the superconcentrated LiTFSA (TFSA: bis(trifluoromethylsulfonil)amide) aqueous solutions were investigated by means of Raman spectroscopic experiments, high-energy X-ray total scattering measurements, and the neutron diffraction technique with different isotopic composition ratios of ⁶Li/⁷Li and H/D. The Li⁺ local structure changes with the increase of the LiTFSA concentration; the oligomer ([Li_p(TFSA)_q]^(p-q)+ (q > 2) forms at the molar fraction of LiTFSA (x_LiTFSA) > 0.25. The average structure can be determined in which two water molecules and two oxygen atoms of TFSA anion(s) coordinate to the Li⁺ in the superconcentrated LiTFSA aqueous solution (LiTFSA)_0.25(H₂O)_0.75. In addition, the intermolecular interaction between the neighboring water molecules was not found, and the hydrogen-bonded interaction in the solution should be significantly weak. According to the coordination number of the oxygen atom (TFSA or H₂O), a variety of TFSA^- and H₂O coordination manners would exist in this solution; in particular, the oligomer is formed in which the monodentate TFSA cross-links Li⁺.

Elucidating the mechanism behind the infrared spectral features and dynamics observed in the carbonyl stretch region of organic carbonates interacting with lithium ions

Ultrafast infrared spectroscopy has become a very important tool for studying the structure and ultrafast dynamics in solution. In particular, it has been recently applied to investigate the molecular interactions and motions of lithium salts in organic carbonates. However, there has been a discrepancy in the molecular interpretation of the spectral features and dynamics derived from these spectroscopies. Hence, the mechanism behind spectral features appearing in the carbonyl stretching region was further investigated using linear and nonlinear spectroscopic tools and the co-solvent dilution strategy. Lithium perchlorate in a binary mixture of dimethyl carbonate (DMC) and tetrahydrofuran was used as part of the dilution strategy to identify the changes of the spectral features with the number of carbonates in the first solvation shell since both solvents have similar interaction energetics with the lithium ion. Experiments showed that more than one carbonate is always participating in the lithium ion solvation structures, even at the low concentration of DMC. Moreover, temperature-dependent study revealed that the exchange of the solvent molecules coordinating the lithium ion is not thermally accessible at room temperature. Furthermore, time-resolved IR experiments confirmed the presence of vibrationally coupled carbonyl stretches among coordinated DMC molecules and demonstrated that this process is significantly altered by limiting the number of carbonate molecules in the lithium ion solvation shell. Overall, the presented experimental findings strongly support the vibrational energy transfer as the mechanism behind the off-diagonal features appearing on the 2DIR spectra of solutions of lithium salt in organic carbonates.

Computing the frequency fluctuation dynamics of highly coupled vibrational transitions using neural networks

The description of frequency fluctuations for highly coupled vibrational transitions has been a challenging problem in physical chemistry. In particular, the complexity of their vibrational Hamiltonian does not allow us to directly derive the time evolution of vibrational frequencies for these systems. In this paper, we present a new approach to this problem by exploiting the artificial neural network to describe the vibrational frequencies without relying on the deconstruction of the vibrational Hamiltonian. To this end, we first explored the use of the methodology to predict the frequency fluctuations of the amide I mode of N-methylacetamide in water. The results show good performance compared with the previous experimental and theoretical results. In the second part, the neural network approach is used to investigate the frequency fluctuations of the highly coupled carbonyl stretch modes for the organic carbonates in the solvation shell of the lithium ion. In this case, the frequency fluctuation predicted by the neural networks shows a good agreement with the experimental results, which suggests that this model can be used to describe the dynamics of the frequency in highly coupled transitions.

MeTFSI (Me = Li, Na) Solvation in Ethylene Carbonate and Fluorinated Ethylene Carbonate: A Molecular Dynamics Study

Classical and ab initio molecular dynamics (MD) simulations have been performed for electrolytes based on LiTFSI and NaTFSI solutions in ethylene carbonate and its mono- and difluoro derivatives. Differences between electrolytes with Li⁺ or Na⁺ ions and the effect of fluorination on the structure and transport properties have been analyzed. The observed differences are related to the strength of Me⁺-carbonate binding, which is weaker for the Na⁺ cation and/or fluorinated solvents. Infrared spectra have been computed from ab initio MD and density functional tight binding (DFTB) MD trajectories. The changes of vibrational frequencies have been related to the local structure of the electrolyte and to interactions between salt cations and solvent molecules. The frequency shifts obtained from the AIMD simulations agree with experimental data, whereas DFTB underestimates Na⁺-carbonate interactions.

Solvation Structure of Li+ in Concentrated Acetonitrile and N,N-Dimethylformamide Solutions Studied by Neutron Diffraction with 6Li/7Li Isotopic Substitution Methods

Neutron diffraction measurements on ⁶Li/⁷Li isotopically substituted 10 and 33 mol % *LiTFSA (lithium bis(trifluoromethylsulfonyl)amide)-AN-d₃ (acetonitrile-d₃) and 10 and 33 mol % *LiTFSA-DMF-d₇(N,N-dimethylformamide-d₇) solutions have been carried out in order to obtain structural insights on the first solvation shell of Li⁺ in highly concentrated organic solutions. Structural parameters concerning the local structure around Li⁺ have been determined from the least squares fitting analysis of the first-order difference function derived from the difference between carefully normalized scattering cross sections observed for ⁶Li-enriched and natural abundance solutions. In 10 mol % LiTFSA-AN-d₃ solution, 3.25 ± 0.04 AN molecules are coordinated to Li⁺ with a intermolecular Li⁺···N(AN) distance of 2.051 ± 0.007 Å. It has been revealed that 1.67 ± 0.07 AN molecules and 2.00 ± 0.01 TFSA^- are involved in the first solvation shell of Li⁺ in the 33 mol % LiTFSA-AN solution. The nearest neighbor Li⁺···N_AN and Li⁺···O_TFSA^- distances are obtained to be r(Li⁺···N) = 2.09 ± 0.01 Å and r(Li⁺···O) = 1.88 ± 0.01 Å, respectively. The first solvation shell of Li⁺ in the 10 mol % LiTFSA-DMF-d₇ solutions contains 3.4 ± 0.1 DMF molecules with an intermolecular Li⁺···O_DMF distance of 1.95 ± 0.02 Å. In highly concentrated 33 mol % LiTFSA-DMF-d₇ solutions, there are 1.3 ± 0.2 DMF molecules and 3.2 ± 0.2 TFSA^- in the first solvation shell of Li⁺ with intermolecular distances of r(Li⁺···O_DMF) = 1.90 ± 0.02 Å and r(Li⁺···O_TFSA^-) = 2.01 ± 0.01 Å, respectively. The Li⁺···TFSA^- contact ion pair stably exists in highly concentrated 33 mol % LiTFSA-AN and -DMF solutions.

Physicochemical and Structural Properties of a Hydrophobicity/Hydrophilicity Switchable Ionic Liquid

Herein, we report the physicochemical and structural properties of a new solubility-switchable ionic liquid (IL) comprising the glycerammonium (GA) cation with a hydrophilic group, the GA cation attached to an acetal-based protective group [protected GA (PGA)], and bis(trifluoromethanesulfonyl)amide (TFSA). The interionic volumes (V_inter) of the hydrophobic [PGA][TFSA] and hydrophilic [GA][TFSA] ILs were evaluated based on solution density, revealing weaker ion-ion interactions in these relative to conventional ILs. The [PGA][TFSA] and [GA][TFSA] also exhibit poor ion-conducting properties, with up to an order of magnitude lower ionic conductivity (σ) and self-diffusion coefficient (D), as compared with conventional ILs. Radial distribution functions derived from high-energy X-ray total scattering experiments [G^exp(r)] and molecular dynamics (MD) simulations [G^MD(r)] indicate that nearest-neighbor ion-ion interactions in the [PGA][TFSA] and [GA][TFSA] are comparable to those in imidazolium-based IL. Conversely, these are appreciably weakened at the second- and third-neighbors and thus less structured in the long range (r > 12 Å) and very different from the highly ordered imidazolium IL. The atom-atom pair correlation function derived from the MD simulations disclose that at a local scale, specific interactions are absent, with only an electrostatic interaction in the [PGA][TFSA], whereas the GA cations interact with TFSA anions via hydrogen bonding of diol groups in the GA and O atoms in the TFSA. No hydrogen bonding group within the PGA cation leads to weak ion-hydration resulting in a phase separation of [PGA][TFSA] and water; in contrast, the GA cations are easily hydrogen-bonded with water molecules to be miscible in aqueous solutions.

Molecular Structure, Chemical Exchange, and Conductivity Mechanism of High Concentration LiTFSI Electrolytes

High concentration lithium electrolytes have been found to be good candidates for high energy density and high voltage lithium batteries. Recent studies have shown that limiting the free solvent molecules in the electrolytes prevents the degradation of the battery electrodes. However, the molecular level knowledge of the structure and dynamics of such an electrolyte system is limited, especially for electrolytes based on typical organic carbonates. In this article, the interactions and motions involved in lithium bis(trifluoromethanesulfonyl)imide in carbonyl-containing solvents are investigated using linear and time-resolved vibrational spectroscopies and computational methods. Our results suggest that the overall structure and the speciation of the three high concentration electrolytes are similar. However, the cyclic carbonate-based electrolyte presents an additional interaction as a result of dimer formation. Time-resolved studies reveal similar and fast dynamics for the structural motions of solvent molecules in electrolytes composed of linear molecules, while the electrolyte made of cyclic solvent molecules shows slower structural changes as a result of the dimer formation. Additionally, a picosecond time scale process is observed and assigned to the coordination and decoordination of solvent molecules from a lithium-ion solvation shell. This process of solvent exchange is found to be directly correlated to the making and breaking of structures between the lithium-ion and the anion and, consequently, to the conduction mechanism. Overall, our data show that the molecular structure of the solvent does not significantly affect the speciation and distribution of the lithium-ion solvation shells. However, the presence of dimerization between solvent molecules of two neighboring lithium-ions appears to produce a microscopic ordering that it is manifested macroscopically in properties of the electrolyte, such as its viscosity.

Solvation Structure of Li+ in Methanol and 2-Propanol Solutions Studied by ATR-IR and Neutron Diffraction with 6Li/7Li Isotopic Substitution Methods

Neutron diffraction measurements have been carried out on 10 mol % LiTFSA (TFSA: bis(trifluoromethylsulfonil)amide) solutions in methanol- d₄ and 2-propanol- d₈ to obtain information on the solvation structure of Li⁺. The detailed coordination structure of solvent molecules within the first solvation shell of Li⁺ was determined through the least-squares fitting analysis of the difference function between normalized scattering cross sections observed for ⁶Li/⁷Li isotopically substituted sample solutions. The nearest-neighbor Li⁺···O distance and coordination number determined for the 10 mol % LiTFSA-methanol- d₄ solution are r_LiO = 1.98 ± 0.02 Å and n_LiO = 3.8 ± 0.6, respectively. In the 2-propanol- d₈ solution, it has been revealed that 2-propanol- d₈ molecules within the first solvation shell of Li⁺ take at least two different coordination geometries with the intermolecular nearest-neighbor Li⁺···O distance of r_LiO = 1.93 ± 0.04 Å. The Li⁺···O coordination number, n_LiO = 3.3 ± 0.3, is determined. Ion-pair formation in the LiTFSA-methanol and LiTFSA-2-propanol solutions has been investigated by the attenuated total reflection infrared spectroscopic method. Mole fractions of free, Li⁺-bound, and aggregated TFSA^- are derived from the peak deconvolution analysis of vibrational bands observed for TFSA^-.