Abstract
High concentration water-in-salt electrolyte (WiSE) systems can expand the electrochemical stability window of water, thereby enabling the application of water-based electrolytes in Li-ion batteries. However, the solvation structure and the dynamics of the ions are not yet fully resolved, and prior molecular-mechanics-based molecular dynamics (MMMD) simulation studies present contrasting viewpoints. In the present work, we utilize first-principles molecular dynamics (FPMD) simulations to study the structure and dynamics of high-concentration (10 and 20 m) LiTFSI electrolyte solutions at 298 and 373 K. Although computationally more expensive than the MMMD simulations, the FPMD simulations, in which the forces on the nuclei are obtained from Kohn-Sham density functional theory reflecting the instantaneous arrangement of the electron density, may offer a more accurate representation of WiSE systems where polarization and charge transfer are important. The FPMD simulations demonstrate disruption of the water hydrogen bonding environment and concurrent formation of an anionic network upon increasing the LiTFSI concentration from 10 to 20 m. However, nanoscale spatial heterogeneity is not observed. Analysis of the Li+ cation dynamics obtained from both FPMD and MMMD simulations indicates that ion transport proceeds predominantly via a mixed-mode mechanism, with contributions from both vehicular motion and hopping depending on concentration and temperature.
| Original language | English |
|---|---|
| Article number | 144504 |
| Journal | Journal of Chemical Physics |
| Volume | 163 |
| Issue number | 14 |
| DOIs | |
| State | Published - Oct 14 2025 |
Funding
This research was supported by a collaborative National Science Foundation Grant No. OAC-1835067 to the University of Minnesota, Grant Nos. OAC-1835874 and DMR-2119575 to Vanderbilt University, and Grant No. OAC-1835630 to the University of Notre Dame. Computational resources were provided by the Minnesota Supercomputing Institute at the University of Minnesota, Twin Cities. Y.Z. and E.M. thank the Center for Research Computing (CRC) at the University of Notre Dame for providing computational resources. The authors thank Co D. Quach and Matthew W. Thompson for help with the GROMACS simulations and Saumil Chheda for the helpful discussions.