Halide, Ammonium, and Alkali Metal Ion Parameters for Modeling Aqueous Solutions

The influence of P-toluenesulfonyl isocyanate on the solvation structure and solid electrolyte interphase formation on graphite anode under high temperatures

The formation and stability of the solid electrolyte interphase (SEI) play a crucial role in determining the performance and lifespan of lithium-ion batteries (LIBs) at evaluated temperatures. Electrolyte additives have emerged as promising candidates for modulating the formation kinetics and stability of the SEI layer. In this study, molecular dynamics (MD) and ab initio molecular dynamics (AIMD) simulations were employed to investigate the influences of P-toluenesulfonyl isocyanate (PTSI) as an electrolyte additive in the LiPF6/EC/DMC electrolyte on the SEI formation process on the graphite anode under high temperatures. MD simulations revealed that PTSI induces modifications in the solvation structure, resulting in two ethylene carbonate (EC) and two dimethyl carbonate (DMC) molecules in the first solvation shell of Li+. Furthermore, the PTSI additive suppressed the decomposition of solvent molecules and PF6 anions in the LiPF6/EC/DMC/PTSI electrolyte under high temperatures according to the AIMD simulations, contributing to a more stable SEI layer. Moreover, this suppression can be attributed to the reduced reactivity of solvent molecules and salt anions, as evidenced by the enhanced bond strength and decreased bond length. These microscopic insights provide a multi-faceted understanding of the functionality of PTSI additives, facilitating the rational design of novel additives.

Molecular dynamics insights into the dynamical behavior of structurally modified water in aqueous deep eutectic solvents (ADES)

Recent studies have demonstrated that the presence of water in deep eutectic solvents (DESs) significantly affects their dynamics, structure, and physical properties. Although the structural changes due to the addition of water are well understood, the microscopic dynamics of these changes have been rarely studied. Here, we performed molecular dynamics simulation of 30% (v/v) (∼0.57 molar fraction) water mixture of DES containing CH3CONH2 and NaSCN/KSCN at various salt fractions to understand the microscopic structure and dynamics of water. The simulated results reveal a heterogeneous environment for water molecules in aqueous DES (ADES), which is influenced by the nature of the cation. The diffusion coefficients of water in ADESs are significantly lower than that in neat water and concentrated aqueous NaSCN/KSCN solution. When Na+ ions are replaced by K+ ions in the ADES system, the diffusion coefficient increases, which is consistent with the measured nuclear magnetic resonance data. Self-dynamic structure factor for water and other simulated dynamic quantities, such as reorientation, hydrogen-bond, and residence time correlation functions, show markedly slower dynamics inside ADES than in the neat water and aqueous salt solution. Moreover, these dynamics become faster when Na+ ions in ADES are replaced by K+ ions. The results suggest that the structural environment of water in Na+-rich ADES is rigid due to the presence of cation-bound water and geometrically constrained water. The medium becomes less rigid as the KSCN fraction increases due to the relatively weaker interaction of K+ ions with water than Na+ ions, which accelerates the dynamical processes.

Molecule Design for Non-Aqueous Wide-Temperature Electrolytes via the Intelligentized Screening Method

Operating a lithium-ion battery (LIB) in a wide temperature range is essential for ensuring a stable electricity supply amidst fluctuating temperatures caused by climate or terrain changes. Electrolyte plays a pivotal role in determining the temperature durability of batteries. However, specialized electrolytes designed for either low or high temperatures typically possess distinct features. Therefore, wide-temperature electrolytes (WTEs) are necessary as they encompass a combination of diverse properties, which complicates the clear instruction of WTE design. Here we represent an artificial intelligence (Al)-assisted workflow of WTE design through stepwise parameterizations and calculations. Linear mono-nitriles are identified as ideal wide-liquidus-range solvents that can "softly" solvate lithium ions by weak interactions. In addition, the explainable modules revealed the halogenoid similarity of cyanide as fluorine on the electrolyte properties (e.g. boiling point and dielectric constant). With the further introduction of an ether bond, 3-methoxypropionitrile (MPN) has been eventually determined as a main electrolyte solvent, enabling the battery operation from -60 to 120 °C. Particularly, a LiCoO2/Li cell using the proposed WTE can realize stable cycling with capacity retention reaching 72.3 % after 50 cycles under a high temperature of 100 °C.

Simulating the Competitive Ion Pairing of Hydrated Electrons with Chaotropic Cations

Experiments show that the absorption spectrum of the hydrated electron (ehyd-) blue-shifts in electrolyte solutions compared with what is seen in pure water. This shift has been assigned to the ehyd-'s competitive ion-pairing interactions with the salt cation relative to the salt anion based on the ions' positions on the Hofmeister series. Remarkably, little work has been done investigating the ehyd-'s behavior when the salts have chaotropic cations, which should greatly change the ion-pairing interactions given that the ehyd- is a champion chaotrope. In this work, we remedy this by using mixed quantum/classical simulations to analyze the behavior of two different models of the ehyd- in aqueous RbF and RbI electrolyte solutions as a function of salt concentration. We find that the magnitude of the salt-induced spectral blue-shift is determined by a combination of the number of chaotropic Rb+ cations near the ehyd- and the number of salt anions near those cations so that the spectrum of the ehyd- directly reflects its local environment. We also find that the use of a soft-cavity ehyd- model predicts stronger competitive interactions with Rb+ relative to I- than a more traditional hard cavity model, leading to different predicted spectral shifts that should provide a way to distinguish between the two models experimentally. Our simulations predict that at the same concentration, salts with chaotropic cations should produce larger spectral blue-shifts than salts with kosmotropic cations. We also found that at high salt concentrations with chaotropic cations, the predicted blue-shift is greater when the salt anion is kosmotropic instead of chaotropic. Our goal is for this work to inspire experimentalists to make such measurements, which will help provide a spectroscopic means to distinguish between simulations models that predict different hydration structures for the ehyd-.

Real-time structural dynamics of the ultrafast solvation process around photo-excited aqueous halides

This work investigates and describes the structural dynamics taking place following charge-transfer-to-solvent photo-abstraction of electrons from I- and Br- ions in aqueous solution following single- and 2-photon excitation at 202 nm and 400 nm, respectively. A Time-Resolved X-ray Solution Scattering (TR-XSS) approach with direct sensitivity to the structure of the surrounding solvent as the water molecules adopt a new equilibrium configuration following the electron-abstraction process is utilized to investigate the structural dynamics of solvent shell expansion and restructuring in real-time. The structural sensitivity of the scattering data enables a quantitative evaluation of competing models for the interaction between the nascent neutral species and surrounding water molecules. Taking the I0-O distance as the reaction coordinate, we find that the structural reorganization is delayed by 0.1 ps with respect to the photoexcitation and completes on a time scale of 0.5-1 ps. On longer time scales we determine from the evolution of the TR-XSS difference signal that I0: e- recombination takes place on two distinct time scales of ∼20 ps and 100 s of picoseconds. These dynamics are well captured by a simple model of diffusive evolution of the initial photo-abstracted electron population where the charge-transfer-to-solvent process gives rise to a broad distribution of electron ejection distances, a significant fraction of which are in the close vicinity of the nascent halogen atoms and recombine on short time scales.

Are NaTFSI and NaFSI Salt-Based Water-in-Salt Electrolytes Structurally Similar or Different?

Water-in-salt electrolytes (WiSEs) are a promising class of electrolytes due to their wide electrochemical stability window and nonflammability. In this study, we explore the structural organization of sodium bis(trifluoromethylsulfonyl)imide (NaTFSI) and sodium bis(fluorosulfonyl)imide (NaFSI) salt-based aqueous electrolytes, covering dilute to highly concentrated regions, by employing an all-atom molecular dynamics simulation. For the NaTFSI-based electrolyte, we observe that Na+ ions are mostly surrounded by water molecules at all the salt concentrations due to the very strong interaction between them. While TFSI anions weakly coordinate with Na+ ions and other TFSI anions, they also mostly prefer to be surrounded by water molecules. These interactions were found to have moderate dependence on the concentration of the NaTFSI salt. For the NaFSI-based electrolyte, while the Na+-water interaction is stronger at lower salt concentrations, the number of nearest neighbor FSI anions is found to be more than that of water at higher concentrations (≥20 m). This is because the increase in the salt concentration leads to expulsion of water molecules from the solvation shell of Na+ ions and enhances the interaction between Na+ ions and oxygen atoms of FSI. At the highest salt concentration (solubility limit), the bulk-like water structure is completely disrupted and dominated by an anionic network in the FSI-based electrolyte. In contrast, water-water hydrogen bonding network is still present even in the highly concentrated TFSI-based electrolyte. The simulated X-ray scattering pattern displays a low-q peak, revealing the presence of an intermediate range ordering due to alternating anion-rich and water/Na+-rich regions in both the electrolytes. However, the characteristic length scale corresponding to the low-q peak decreases with increasing the salt content in both the electrolytes.

Solvation Structures and Transport Mechanisms of Cations in Water-in-Bisalt Electrolytes for High-Concentration Lithium-Ion Batteries Revealed by Molecular Dynamics Simulations

The development of safe and cost-effective electrolytes for rechargeable batteries is currently underway. While water-based electrolytes hold promise, their restricted electrochemical stability window poses a challenge. Combining multiple ionic species emerges as a promising strategy to broaden this stability window and optimize Li-ion battery performance. This study focuses on dual-cation electrolytes, which blend lithium and potassium acetates to enhance the electrochemical characteristics of the solution at high concentrations. We investigated the solvation structure of each ion and its interactions on a molecular level. Our analysis reveals that ion clusters and aggregates are formed through shared acetate and water molecules at high salt concentrations. Furthermore, the residence time analyses of atom pairs indicate that cations diffuse in vehicular mode at low concentrations. In contrast, they switch to a structural mode at high concentrations due to diminishing water content. This study offers a comprehensive model for exploring diverse solvation structures of cations and gaining insights into their diffusion mechanisms within water-in-bisalt electrolytes for aqueous Li-ion batteries.

All-atom molecular dynamics simulations of polymer and polyelectrolyte brushes

Densely grafted polymer and polyelectrolyte (PE) brushes, owing to their significant abilities to functionalize surfaces for a plethora of applications in sensing, diagnostics, current rectification, surface wettability modification, drug delivery, and oil recovery, have attracted significant attention over the past several decades. Unfortunately, most of the attention has primarily focused on understanding the properties of the grafted polymer and the PE chains with little attention devoted to studying the behavior of the brush-supported ions (counterions needed to screen the PE chains) and water molecules. Over the past few years, our group has been at the forefront of addressing this gap: we have employed all-atom molecular dynamics (MD) simulations for studying a wide variety of polymer and PE brush systems with specific attention to unraveling the properties and behavior of the brush-supported water molecules and ions. Our findings have revealed some of the most fascinating properties of such brush-supported ions and water molecules, including the most remarkable control of nanofluidic transport afforded by the specific ion and water responses induced by the PE brushes grafted on the inner walls of the nanochannel. This feature article aims to summarize some of our key contributions associated with such atomistic simulations of polymer and PE brushes and brush-supported water molecules and counterions.

A locally solvent-tethered polymer electrolyte for long-life lithium metal batteries

Solid polymer electrolytes exhibit enhanced Li+ conductivity when plasticized with highly dielectric solvents such as N,N-dimethylformamide (DMF). However, the application of DMF-containing electrolytes in solid-state batteries is hindered by poor cycle life caused by continuous DMF degradation at the anode surface and the resulting unstable solid-electrolyte interphase. Here we report a composite polymer electrolyte with a rationally designed Hofmann-DMF coordination complex to address this issue. DMF is engineered on Hofmann frameworks as tethered ligands to construct a locally DMF-rich interface which promotes Li+ conduction through a ligand-assisted transport mechanism. A high ionic conductivity of 6.5 × 10-4 S cm-1 is achieved at room temperature. We demonstrate that the composite electrolyte effectively reduces the free shuttling and subsequent decomposition of DMF. The locally solvent-tethered electrolyte cycles stably for over 6000 h at 0.1 mA cm-2 in Li | |Li symmetric cell. When paired with sulfurized polyacrylonitrile cathodes, the full cell exhibits a prolonged cycle life of 1000 cycles at 1 C. This work will facilitate the development of practical polymer-based electrolytes with high ionic conductivity and long cycle life.

Effect of K+ Force Fields on Ionic Conductivity and Charge Dynamics of KOH in Ethylene Glycol

Predicting ionic conductivity is crucial for developing efficient electrolytes for energy storage and conversion and other electrochemical applications. An accurate estimate of ionic conductivity requires understanding complex ion-ion and ion-solvent interactions governing the charge transport at the molecular level. Molecular simulations can provide key insights into the spatial and temporal behavior of electrolyte constituents. However, such insights depend on the ability of force fields to describe the underlying phenomena. In this work, molecular dynamics simulations were leveraged to delineate the impact of force field parameters on ionic conductivity predictions of potassium hydroxide (KOH) in ethylene glycol (EG). Four different force fields were used to represent the K+ ion. Diffusion-based Nernst-Einstein and correlation-based Einstein approaches were implemented to estimate the ionic conductivity, and the predicted values were compared with experimental measurements. The physical aspects, including ion-aggregation, charge distribution, cluster correlation, and cluster dynamics, were also examined. A force field was identified that provides reasonably accurate Einstein conductivity values and a physically coherent representation of the electrolyte at the molecular level.

Reparameterization of Polarizable Force Fields for Studying Ion Transfer across Liquid-Liquid Interfaces

We have developed a general scheme for refining classical polarizable molecular dynamics (MD) force fields that can accurately describe the molecular interactions in systems with liquid-liquid interfaces. While ab initio MD (AIMD) simulations can naturally describe molecular interactions, they are often so computationally expensive that simulating large system sizes and/or long time scales is usually infeasible. To resolve this, we parameterized efficient and accurate classical polarizable force fields that use AIMD reference data by minimizing both the relative entropy and the root mean squared deviation in atomic forces. We utilized our new multiscale models to study chloride ion transfer across the water-dichloromethane (DCM) interface with and without the tetraethylammonium cation as the phase-transfer catalyst. Our calculated free-energy barrier for the water-DCM interface is consistent with the other reported simulation results. We further analyzed the ion-transfer process by studying the hydration shell structures around the chloride ion and the ion-pair formation to better understand the mechanism. We observed that electronic polarizability is an important factor for the studied phase-transfer catalyst to lower the free-energy barrier of the ion transfer. Using the water-benzene interface system as an additional example, we show that our parameterization scheme provides a general route for modeling different liquid-liquid interface systems even when the experimental data or force field parameters are not readily available.

Temperature-dependent dielectric relaxation measurements of (acetamide + K/Na SCN) deep eutectic solvents: Decoding the impact of cation identity via computer simulations

The impact of successive replacement of K+ by Na+ on the megahertz-gigahertz polarization response of 0.25[fKSCN + (1 - f)NaSCN] + 0.75CH3CONH2 deep eutectic solvents (DESs) was explored via temperature-dependent (303 ≤ T/K ≤ 343) dielectric relaxation (DR) measurements and computer simulations. Both the DR measurements (0.2 ≤ ν/GHz ≤ 50) and the simulations revealed multi-Debye relaxations accompanied by a decrease in the solution static dielectric constant (ɛs) upon the replacement of K+ by Na+. Accurate measurements of the DR response of DESs below 100 MHz were limited by the well-known one-over-frequency divergence for conducting solutions. This problem was tackled in simulations by removing the zero frequency contributions arising from the ion current to the total simulated DR response. The temperature-dependent measurements revealed a much stronger viscosity decoupling of DR times for Na+-containing DES than for the corresponding K+ system. The differential scanning calorimetry measurements indicated a higher glass transition temperature for Na+-DES (∼220 K) than K+-DES (∼200 K), implying more fragility and cooperativity for the former (Na+-DES) than the latter. The computer simulations revealed a gradual decrease in the average number of H bonds (⟨nHB⟩) per acetamide molecule and increased frustrations in the average orientational order upon the replacement of K+ by Na+. Both the measured and simulated ɛs values were found to decrease linearly with ⟨nHB⟩. Decompositions of the simulated DR spectra revealed that the cation-dependent cross interaction (dipole-ion) term contributes negligibly to ɛs and appears in the terahertz regime. Finally, the simulated collective single-particle reorientational relaxations and the structural H-bond fluctuation dynamics revealed the microscopic origin of the cation identity dependence shown by the measured DR relaxation times.

Effect of the Ion, Solvent, and Thermal Interaction Coefficients on Battery Voltage

In order to increase the adoption of batteries for sustainable transport and energy storage, improved charging and discharging capabilities of lithium-ion batteries are necessary. To achieve this, accurate data that describe the internal state of the cells are essential. Several models have been derived, and transport coefficients have been reported for use in these models. We report for the first time a complete set of transport coefficients to model the concentration and temperature polarization in a lithium-ion battery ternary electrolyte, allowing us to test common assumptions. We include effects due to gradients in chemical potentials and in temperature. We find that the voltage contributions due to salt and solvent polarization are of the same order of magnitude as the ohmic loss and must be taken into account for more accurate modeling and understanding of battery performance. We report new Soret and Seebeck coefficients and find thermal polarization to be significant in cases relevant to battery research. The analysis is suitable for electrochemical systems, in general.

Ab Initio Molecular Dynamics Simulations of Aqueous LiTFSI Solutions─Structure, Hydrogen Bonding, and IR Spectra

Simulations of Density Functional Theory-based ab initio molecular dynamics (AIMD) have been performed for a series of aqueous lithium bis(trifluoromethylsulfonyl)imide (LiTFSI) solutions with concentrations ranging from salt-in-water to water-in-salt systems. Analysis of the structure of electrolytes has revealed a preference of Li+ cations to interact with water molecules. In concentrated LiTFSI solutions, water molecules form small associates. The total number of hydrogen bonds (HBs) in the system decreases with salt concentrations, with bonds to water acceptors being only partially replaced by interactions with TFSI anions. Infrared (IR) spectra in the region of the O-H stretching frequency calculated from AIMD trajectories are in good agreement with experimental data. Statistics of oscillations of individual O-H bonds have shown correlations between vibrational frequencies and the structure of HBs formed by water. The changes in the IR spectrum have been related to the varying contributions of different local environments of the water molecules. The abundances of the three spectral components calculated from the simulations agree well with the decomposition of the experimental IR spectra reported in the literature.

Madrid-2019 force field: An extension to divalent cations Sr2+ and Ba2

In this work, we present a parameterization of Sr2+ and Ba2+ cations, which expands the alkali earth set of cations of the Madrid-2019 force field. We have tested the model against the experimental densities of eight different salts, namely, SrCl2, SrBr2, SrI2, Sr(NO3)2, BaCl2, BaBr2, BaI2, and Ba(NO3)2. The force field is able to reproduce the experimental densities of all these salts up to their solubility limit. Furthermore, we have computed the viscosities for two selected salts, finding that the experimental values are overestimated, but the predictions are still reasonable. Finally, the structural properties for all the salts have been calculated with this model and align remarkably well with experimental observations.

Variation of Surface Propensity of Halides with Droplet Size and Temperature: The Planar Interface Limit

The radial number density profiles of halide and alkali ions in aqueous clusters with equimolar radius ≲1.4 nm, which correspond to ≲255 H2O molecules, have been extensively studied by computations. However, the surface abundance of Cl-, Br-, and I- relative to the bulk interior in these smaller clusters may not be representative of the larger systems. Indeed, here we show that the larger the cluster is, the lower the relative surface abundance of chaotropic halides is. In droplets with an equimolar radius of ≈2.45 nm, which corresponds to ≈2000 H2O molecules, the polarizable halides show a clear number density maximum in the droplet's bulk-like interior. A similar pattern is observed in simulations of the aqueous planar interface with halide salts at room temperature. At elevated temperature the surface propensity of Cl- decreases gradually, while that of I- is partially preserved. The change in the chaotropic halide location at higher temperatures than the room temperature may considerably affect photochemical reactivity in atmospheric aerosols, vapor-liquid nucleation and growth mechanisms, and salt crystallization via solvent evaporation. We argue that the commonly used approach of nullifying parameters in a force field in order to find the factors that determine the ion location does not provide transferable insight into other force fields.

Stable Harsh-Temperature Lithium Metal Batteries Enabled by Tailoring Solvation Structure in Ether Electrolytes

For lithium metal batteries (LMBs), the elevated operating temperature results in severe capacity fading and safety issues due to unstable electrode-electrolyte interphases and electrolyte solvation structures. Therefore, it is crucial to construct advanced electrolytes capable of tolerating harsh environments to ensure stable LMBs. Here, we proposed a stable localized high-concentration electrolyte (LHCE) by introducing the highly solvating power solvent diethylene glycol dimethyl ether (DGDME). Computational and experimental evidence discloses that the original DGDME-LHCE shows favorable features for high-temperature LMBs, including high Li+-binding stability, electro-oxidation resistance, thermal stability, and nonflammability. The tailored solvated sheath structure achieves the preferred decomposition of anions, inducing the stable (cathode and Li anode)/interphases simultaneously, which enables a homogeneous Li plating-stripping behavior on the anode side and a high-voltage tolerance on the cathode side. For the Li||Li cells coupled with DGDME-LHCE, they showcase outstanding reversibility (a long lifespan of exceeding 1900 h). We demonstrate exceptional cyclic stability (∼95.59%, 250 cycles), high Coulombic efficiency (>99.88%), and impressive high-voltage (4.5 V) and high-temperature (60 °C) performances in Li||NCM523 cells using DGDME-LHCE. Our advances shed light on an encouraging ether electrolyte tactic for the Li-metal batteries confronted with stringent high-temperature challenges.

Explicit solvent machine-learned coarse-grained model of sodium polystyrene sulfonate to capture polymer structure and dynamics

Strongly charged polyelectrolytes (PEs) demonstrate complex solution behavior as a function of chain length, concentrations, and ionic strength. The viscosity behavior is important to understand and is a core quantity for many applications, but aspects remain a challenge. Molecular dynamics simulations using implicit solvent coarse-grained (CG) models successfully reproduce structure, but are often inappropriate for calculating viscosities. To address the need for CG models which reproduce viscoelastic properties of one of the most studied PEs, sodium polystyrene sulfonate (NaPSS), we report our recent efforts in using Bayesian optimization to develop CG models of NaPSS which capture both polymer structure and dynamics in aqueous solutions with explicit solvent. We demonstrate that our explicit solvent CG NaPSS model with the ML-BOP water model [Chan et al. Nat Commun 10, 379 (2019)] quantitatively reproduces NaPSS chain statistics and solution structure. The new explicit solvent CG model is benchmarked against diffusivities from atomistic simulations and experimental specific viscosities for short chains. We also show that our Bayesian-optimized CG model is transferable to larger chain lengths across a range of concentrations. Overall, this work provides a machine-learned model to probe the structural, dynamic, and rheological properties of polyelectrolytes such as NaPSS and aids in the design of novel, strongly charged polymers with tunable structural and viscoelastic properties.

Accurate Calculation of Solvation Properties of Lithium Ions in Nonaqueous Solutions

We perform all-atom molecular dynamics simulations of lithium triflate in 1,2-dimethoxyethane using six different literature force fields. This system is representative of many experimental studies of lithium salts in solvents and polymers. We show that multiple historically common force fields for lithium ions give qualitatively incorrect results when compared with those from experiments and quantum chemistry calculations. We illustrate the importance of correctly selecting force field parameters and give recommendations on the force field choice for lithium electrolyte applications.

Computation of Electrical Conductivities of Aqueous Electrolyte Solutions: Two Surfaces, One Property

In this work, we computed electrical conductivities under ambient conditions of aqueous NaCl and KCl solutions by using the Einstein-Helfand equation. Common force fields (charge q = ±1 e) do not reproduce the experimental values of electrical conductivities, viscosities, and diffusion coefficients. Recently, we proposed the idea of using different charges to describe the potential energy surface (PES) and the dipole moment surface (DMS). In this work, we implement this concept. The equilibrium trajectories required to evaluate electrical conductivities (within linear response theory) were obtained by using scaled charges (with the value q = ±0.75 e) to describe the PES. The potential parameters were those of the Madrid-Transport force field, which accurately describe viscosities and diffusion coefficients of these ionic solutions. However, integer charges were used to compute the conductivities (thus describing the DMS). The basic idea is that although the scaled charge describes the ion-water interaction better, the integer charge reflects the value of the charge that is transported due to the electric field. The agreement obtained with experiments is excellent, as for the first time electrical conductivities (and the other transport properties) of NaCl and KCl electrolyte solutions are described with high accuracy for the whole concentration range up to their solubility limit. Finally, we propose an easy way to obtain a rough estimate of the actual electrical conductivity of the potential model under consideration using the approximate Nernst-Einstein equation, which neglects correlations between different ions.

Closed-loop cathode recycling in solid-state batteries enabled by supramolecular electrolytes

Deconstructing solid-state batteries (SSBs) to physically separated cathode and solid-electrolyte particles remains intensive, as does the remanufacturing of cathodes and separators from the recovered materials. To address this challenge, we designed supramolecular organo-ionic (ORION) electrolytes that are viscoelastic solids at battery operating temperatures (-40° to 45°C) yet are viscoelastic liquids above 100°C, which enables both the fabrication of high-quality SSBs and the recycling of their cathodes at end of life. SSBs implementing ORION electrolytes alongside Li metal anodes and either LFP or NMC cathodes were operated for hundreds of cycles at 45°C with less than 20% capacity fade. Using a low-temperature solvent process, we isolated the cathode from the electrolyte and demonstrated that refurbished cells recover 90% of their initial capacity and sustain it for an additional 100 cycles with 84% capacity retention in their second life.

Exploring Pyrrolo-Fused Heterocycles as Promising Anticancer Agents: An Integrated Synthetic, Biological, and Computational Approach

Five new series of pyrrolo-fused heterocycles were designed through a scaffold hybridization strategy as analogs of the well-known microtubule inhibitor phenstatin. Compounds were synthesized using the 1,3-dipolar cycloaddition of cycloimmonium N-ylides to ethyl propiolate as a key step. Selected compounds were then evaluated for anticancer activity and ability to inhibit tubulin polymerization in vitro. Notably, pyrrolo[1,2-a]quinoline 10a was active on most tested cell lines, performing better than control phenstatin in several cases, most notably on renal cancer cell line A498 (GI₅₀ 27 nM), while inhibiting tubulin polymerization in vitro. In addition, this compound was predicted to have a promising ADMET profile. The molecular details of the interaction between compound 10a and tubulin were investigated through in silico docking experiments, followed by molecular dynamics simulations and configurational entropy calculations. Of note, we found that some of the initially predicted interactions from docking experiments were not stable during molecular dynamics simulations, but that configurational entropy loss was similar in all three cases. Our results suggest that for compound 10a, docking experiments alone are not sufficient for the adequate description of interaction details in terms of target binding, which makes subsequent scaffold optimization more difficult and ultimately hinders drug design. Taken together, these results could help shape novel potent antiproliferative compounds with pyrrolo-fused heterocyclic cores, especially from an in silico methodological perspective.

Dielectric relaxation and dielectric decrement in ionic acetamide deep eutectic solvents: Spectral decomposition and comparison with experiments

Frequency-dependent dielectric relaxation in three deep eutectic solvents (DESs), (acetamide+LiClO4/NO3/Br), was investigated in the temperature range, 329 ≤ T/K ≤ 358, via molecular dynamics simulations. Subsequently, decomposition of the real and the imaginary components of the simulated dielectric spectra was carried out to separate the rotational (dipole-dipole), translational (ion-ion), and ro-translational (dipole-ion) contributions. The dipolar contribution, as expected, was found to dominate all the frequency-dependent dielectric spectra over the entire frequency regime, while the other two components together made tiny contributions only. The translational (ion-ion) and the cross ro-translational contributions appeared in the THz regime in contrast to the viscosity-dependent dipolar relaxations that dominated the MHz-GHz frequency window. Our simulations predicted, in agreement with experiments, anion-dependent decrement of the static dielectric constant (ɛs ∼ 20 to 30) for acetamide (ɛs ∼ 66) in these ionic DESs. Simulated dipole-correlations (Kirkwood g factor) indicated significant orientational frustrations. The frustrated orientational structure was found to be associated with the anion-dependent damage of the acetamide H-bond network. Single dipole reorientation time distributions suggested slowed down acetamide rotations but did not indicate presence of any "rotationally frozen" molecule. The dielectric decrement is, therefore, largely static in origin. This provides a new insight into the ion dependence of the dielectric behavior of these ionic DESs. A good agreement between the simulated and the experimental timescales was also noticed.

How NaFTA salt affects the structural landscape and transport properties of Pyrr1,3FTA ionic liquid

Recently, it has been demonstrated that ionic liquids (ILs) with an asymmetric anion render a wider operational temperature range and can be used as a solvent in sodium ion batteries. In the present study, we examine the microscopic structure and dynamics of pure 1-methyl-1-propylpyrrolidinium fluorosulfonyl(trifluoromethylsulfonyl)amide (Pyrr_1,3FTA) IL using atomistic molecular dynamics simulations. How the addition of the sodium salt (NaFTA) having the same anion changes the structural landscape and transport properties of the pure IL has also been explored. The simulated x-ray scattering structure functions reveal that the gradual addition of NaFTA salt (up to 1.2 molal) suppresses the charge alternating feature of the pure IL because of the replacement of the Pyrr⁺ cations with the Na⁺ ions. The Na⁺ ions are majorly found near the oxygen atoms of the anions, but the probability of finding the Na⁺ ions near these atoms slightly decreases with increasing salt concentration. As expected, the Na⁺ ions stay away from the Pyrr⁺ cations. However, the probability of finding the anions around anions increases with increasing salt concentration. The simulated self-diffusion coefficients of the ions in the pure IL reveal slightly faster diffusion of the Pyrr⁺ cations as compared to the FTA^- anions. Interestingly, in the salt solution, despite having smaller size, the diffusion of the Na⁺ ions is found to be lesser than the Pyrr⁺ cations and the FTA^- anions. The analysis of the ionic conductivity and transport numbers reveals that the fractional contribution of the FTA^- anion to the overall conductivity remains nearly constant with increasing salt concentration, but the contribution of Pyrr⁺ cation decreases and Na⁺ ion increases.

Scaled charges for ions: An improvement but not the final word for modeling electrolytes in water

In this work, we discuss the use of scaled charges when developing force fields for NaCl in water. We shall develop force fields for Na⁺ and Cl^- using the following values for the scaled charge (in electron units): ±0.75, ±0.80, ±0.85, and ±0.92 along with the TIP4P/2005 model of water (for which previous force fields were proposed for q = ±0.85 and q = ±1). The properties considered in this work are densities, structural properties, transport properties, surface tension, freezing point depression, and maximum in density. All the developed models were able to describe quite well the experimental values of the densities. Structural properties were well described by models with charges equal to or larger than ±0.85, surface tension by the charge ±0.92, maximum in density by the charge ±0.85, and transport properties by the charge ±0.75. The use of a scaled charge of ±0.75 is able to reproduce with high accuracy the viscosities and diffusion coefficients of NaCl solutions for the first time. We have also considered the case of KCl in water, and the results obtained were fully consistent with those of NaCl. There is no value of the scaled charge able to reproduce all the properties considered in this work. Although certainly scaled charges are not the final word in the development of force fields for electrolytes in water, its use may have some practical advantages. Certain values of the scaled charge could be the best option when the interest is to describe certain experimental properties.

Difference in "Supercooling" Affinity between (Acetamide + Na/KSCN) Deep Eutectics: Reflections in the Simulated Anomalous Motions of the Constituents and Solution Microheterogeneity Features

Deep depression of freezing points of ionic amide deep eutectic solvents (DESs) is known to exhibit a significant dependence on the identity of ions present in those systems and the nature of the functional group attached to the host amide. This deep depression of the freezing point is sometimes termed as "supercooling". For (acetamide + electrolyte) DESs, experiments have revealed signatures of ion-dependent spatiotemporal heterogeneity features. The focus of this work is to provide microscopic explanations of these experimentally observed macroscopic system properties in terms of particle jumps and insights about the origin of the cation dependence. For this purpose, extensive molecular dynamics simulations have been performed employing (acetamide + Na/KSCN) deep eutectics as representative ionic systems at 303, 318, 333, and 348 K. The individual translational motions of acetamide and the ions are followed, and their connections to solution heterogeneity are explored. The center-of-mass motion for Na⁺ has been found to be more anomalous than that for K⁺. This difference corroborates well with experimental reports on heterogeneous relaxations in these systems. Simulated viscosity coefficients and dynamic heterogeneity features also reflect this difference. Moreover, simulated reorientational relaxations of acetamide molecules in these ionic DESs suggest that a Na⁺-containing DES is more heterogeneous than the corresponding K⁺-containing system. Estimated void and neck distributions for acetamide molecules differ as the alkali metal ions differ. In brief, this study provides a detailed microscopic view of the cation dependence of the microheterogeneous relaxation dynamics of these DESs reported repeatedly by different experiments.

Hydrophilic Versus Hydrophobic Coupling in the Pressure Dependence of the Chemical Potential of Alkali Metal and Halide Ions in Water

We computed the chemical potential for some alkali metal ions (K⁺, Rb⁺, and Cs⁺) and two halide ions (Br^- and I^-) in aqueous solution at ambient T and various pressures in the range 1-8000 atm. Results were obtained from classic Monte Carlo simulations in the NPT ensemble by means of the free energy perturbation method. Here, the chemical potential is computed as the sum of a term relative to a Lennard-Jones solute and a term relative to the process in which this solute is transformed into the ion. Hydrophobic and hydrophilic features of these two components of the chemical potential show opposite behaviors under isothermal compression. The increase in pressure determines an increase in the hydrophobic component, which becomes more positive with a stronger effect for larger ions. Correspondingly, the values of the hydrophilic component become more negative for alkali ions, whereas they are only slightly affected by compression for halide ions. Hydrophobic-hydrophilic quasi-compensation in the slopes is observed for Rb⁺. For a smaller ion, such as K⁺, the dependence on pressure of the hydrophilic component is slightly dominant. For a larger ion, as observed in the cases of Cs⁺, Br^-, and I^-, the hydrophobic component assumes the determinant role. Pressure dependence of the chemical potential is little affected by corrections introduced for molecular potential truncation. This view can change for possible boundary artifacts that could have affected the static electrostatic potential. Some inference is obtained from comparison with experimental data at 1 atm on the free energy of hydration. Discrepancies show the characteristic asymmetry between cations and anions. The further addition of a correction based on the static potential significantly reduces these discrepancies with important error cancellation on the sum of chemical potentials of ions of opposite charge. The correction is applied also at higher pressures, and results are compared with those obtained by adding an alternative correction that is based on the water number density. Regardless of the ion, changes of the chemical potential induced by an increase in pressure appear to be dominated by the hydrophobic component, in particular when using the alternative correction. For bromide and iodide electrolytes, the two corrections give chemical potentials in good agreement.

Ion Solvation Free Energy Calculation Based on Ab Initio Molecular Dynamics Using a Hybrid Solvent Model

Free energy calculation of small molecules or ion species in aqueous solvent is one of the most important problems in electrochemistry study. Although there are many previous approaches to calculate such free energies, they are based on ab initio methods and suffer from various limitations and approximations. In the current work, we developed a hybrid approach based on ab initio molecular dynamics (AIMD) simulations to calculate the ion solvation energy. In this approach, a small water cluster surrounding the central ion is used, and implicit solvent model is applied outside the water cluster. A dynamic potential well is used during AIMD to keep the water cluster together. Quasi-harmonic approximation is used to calculate the entropy contribution, while the total energy average is used to calculate the enthalpy term. The obtained solvation voltages of the bulk metal agree with experiments within 0.3 eV, and the simulation results for the solvation energies of gaseous ions are close to the experimental observations. Besides the free energies, radial pair distribution functions and coordination numbers of hydrated cations are also obtained. The remaining challenges of this method are also discussed.

Heterogeneity and Nanostructure of Superconcentrated LiTFSI-EmimTFSI Hybrid Aqueous Electrolytes: Beyond the 21 m Limit of Water-in-Salt Electrolyte

Ionic liquids such as EmimTFSI (1-ethyl-3-methylimidazolium bis(trifluoromethanesulfonyl)imide) have been found to improve the solubility of LiTFSI salt in water-in-salt electrolyte (WiSE) from 21 to 60 m. However, the molecular origin of such enhancement in the solubility is still unknown. In the present work, we elucidate the microscopic structures of LiTFSI-EmimTFSI-based hybrid aqueous electrolytes and compare them with the structure of LiTFSI-based WiSE using molecular dynamics simulations. Our analysis reveals the presence of alternating water-rich clusters and TFSI-rich extended domains in the WiSE. In these clusters and domains, the Li⁺ ions reside such that the total number of oxygen atoms around them is conserved to four, where water contributes about three oxygen atoms. The addition of EmimTFSI in the WiSE results in removal of water from the nearest-neighbor solvation shell of TFSI^- ions but not from the Li⁺ ions. Significant structural changes are observed when LiTFSI salt is further added to LiTFSI-EmimTFSI aqueous solution. In both the hybrid electrolytes, water and Emim⁺ cations are found to avoid each other. The simulated X-ray scattering structure factor reveals the presence of larger length-scale heterogeneity in the most concentrated solution of the hybrid aqueous electrolyte. We observe that this nanoscale heterogeneity originates from a water-TFSI-Emim-TFSI-water-TFSI-Emim-TFSI-like arrangement in which Li⁺ ions are dispersed such that the coordination number of oxygen atoms around them is enhanced to five, wherein the major contribution comes from the TFSI^- ions. We envision that the enhanced LiTFSI solubility originates from the replacement of water molecules with TFSI^- ions in the first solvation shell of Li⁺ ions.

Molecular simulation of glycerol-derived triether podands for lithium ion solvation

Solvate ionic liquids (ILs) are promising candidates for several applications due to their stability, high coulombic efficiency, and low volatility. In this work, we investigate the solvation of lithium-bistriflimide by different glycerol-derived triether solvents, using molecular dynamics simulations. Very strong interactions between Li⁺ and the solvent oxygen sites are found, leading to significant conformational changes in the solvent. By comparing the conformation of the neat solvents with their IL mixtures at different concentrations and temperatures, we find that the presence of Li⁺ induces a distinct crown-like structure in the solvent molecules. The Li⁺ cations and the surrounding solvent form a podand complex, which is stable even at elevated temperatures. These glycerol-derived solvents exhibit distinct interactions with Li⁺ cations which may be exploited in electrolytic applications or lithium recovery processes.

Structures of cationic and anionic polyelectrolytes in aqueous solutions: the sign effect

In this study, we use molecular dynamics simulation to explore the structures of anionic and cationic polyelectrolytes in aqueous solutions. We first confirm the significantly stronger solvation effects of single anions compared to cations in water at the fixed ion radii, due to the reversal orientations of asymmetric dipolar H₂O molecules around the ions. Based on this, we demonstrate that the solvation discrepancy of cations/anions and electrostatic correlations of ionic species can synergistically cause the nontrivial structural difference between single anionic and cationic polyelectrolytes. The cationic polyelectrolyte shows an extended structure whereas the anionic polyelectrolyte exhibits a collapsed structure, and their structural differences decline with increasing the counterion size. Furthermore, we corroborate that multiple cationic polyelectrolytes or multiple anionic polyelectrolytes can exhibit largely differential molecular architectures in aqueous solutions. In the solvation dominant regime, the polyelectrolyte solutions exhibit uniform structures; whereas, in the electrostatic correlation dominant regime, the polyelectrolyte solutions exhibit heterogeneous structures, in which the likely charged chains microscopically aggregate through counterion condensations. Increasing the intrinsic chain rigidity causes polyelectrolyte extension and hence moderately weakens the inter-chain clustering. Our work highlights the various, unique structures and molecular architectures of polyelectrolytes in solutions caused by the multi-body correlations between polyelectrolytes, counterions and asymmetric dipolar solvent molecules, which provides insights into the fundamental understanding of ion-containing polymers.

Overlap Concentration of Sodium Polystyrene Sulfonate in Solution

The overlap concentration c* of sodium polystyrene sulfonate in water is calculated using multichain atomistic and coarse grained (CG) simulations for a range of chain lengths. Fully atomistic molecular dynamics simulations are carried out for N = 32-192 monomers. The CG model was parameterized to match the end-to-end distance from the atomistic simulations at small N and allows us to simulate a much larger N. Treating the hydrophobic backbone by inclusion of attraction between monomers is an essential addition to the CG model. The simulation c* are in agreement with experimental data, yet at c*, the chains are not fully stretched, even for N as large as 1200. This implies that none of the experimental systems are in the scaling regime and to reach the scaling regime for NaPSS chains much longer than N = 1200 are required.

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.

Two-dimensional quantum-sheet films with sub-1.2 nm channels for ultrahigh-rate electrochemical capacitance

Dense, thick, but fast-ion-conductive electrodes are critical yet challenging components of ultrafast electrochemical capacitors with high volumetric power/energy densities^1-4. Here we report an exfoliation-fragmentation-restacking strategy towards thickness-adjustable (1.5‒24.0 μm) dense electrode films of restacked two-dimensional 1T-MoS₂ quantum sheets. These films bear the unique architecture of an exceptionally high density of narrow (sub-1.2 nm) and ultrashort (~6.1 nm) hydrophobic nanochannels for confinement ion transport. Among them, 14-μm-thick films tested at 2,000 mV s^-1 can deliver not only a high areal capacitance of 0.63 F cm^-2 but also a volumetric capacitance of 437 F cm^-3 that is one order of magnitude higher than that of other electrodes. Density functional theory and ab initio molecular dynamics simulations suggest that both hydration and nanoscale channels play crucial roles in enabling ultrafast ion transport and enhanced charge storage. This work provides a versatile strategy for generating rapid ion transport channels in thick but dense films for energy storage and filtration applications.

Computer Simulations of Deep Eutectic Solvents: Challenges, Solutions, and Perspectives

Deep eutectic solvents (DESs) are one of the most rapidly evolving types of solvents, appearing in a broad range of applications, such as nanotechnology, electrochemistry, biomass transformation, pharmaceuticals, membrane technology, biocomposite development, modern 3D-printing, and many others. The range of their applicability continues to expand, which demands the development of new DESs with improved properties. To do so requires an understanding of the fundamental relationship between the structure and properties of DESs. Computer simulation and machine learning techniques provide a fruitful approach as they can predict and reveal physical mechanisms and readily be linked to experiments. This review is devoted to the computational research of DESs and describes technical features of DES simulations and the corresponding perspectives on various DES applications. The aim is to demonstrate the current frontiers of computational research of DESs and discuss future perspectives.

Rational Design of Nonbonded Point Charge Models for Monovalent Ions with Lennard-Jones 12-6 Potential

Ions are of central importance in nature, and a variety of potential models was proposed to model ions in different phases for an in-depth exploration of ion-related systems. Here, we developed point charge models of 14 monovalent ions with the traditional 12-6 Lennard-Jones (LJ) potential for use in conjunction with 11 water models of TIP3P, OPC3, SPC/E, SPC/E_b, TIP3P-FB, a99SB-disp, TIP4P-Ew, OPC, TIP4P/2005, TIP4P-D, and TIP4P-FB. The designed models reproduced the real hydration free energy (HFE) of ions and the ion-oxygen distance (IOD) in the first hydration shell accurately and simultaneously, a performance similar to the previously reported 12-6-4 LJ-type ion models (12-6 LJ plus an attractive C₄ term for cations or a repulsive one for anions). This work, along with our previous work on di-, tri-, and tetravalent metal cations (J. Chem. Inf. Model. 2021, 61, 4031-4044; J. Chem. Inf. Model. 2021, 61, 4613-4629), demonstrates the feasibility of the simple 12-6 LJ potential in ion modeling. In order for the 12-6 LJ potential to reproduce both the HFE and IOD, the LJ R parameters need to be close to Shannon's ionic radii for the highly charged cations and to the Stokes's van der Waals (vdW) radii for the monovalent ions. With an additional C₄ term, the R parameters of 12-6-4 LJ ion models agree well with the Stokes's vdW radii and have a more physical meaning. It appears that the C₄ term can be merged into the 12-6 LJ potential by a rational tuning of R and the LJ well depth. Simulations of the osmotic coefficients of alkali chloride solutions and the properties of gaseous and solid alkali halides indicate the necessity of further optimizing ion-ion interactions via, for instance, targeting more properties or using a more physical (polarizable) model.

Temperature-Dependent Dielectric Relaxation in Ionic Acetamide Deep Eutectics: Partial Viscosity Decoupling and Explanations from the Simulated Single-Particle Reorientation Dynamics and Hydrogen-Bond Fluctuations

We report here temperature-dependent (293 ≤ T (K) ≤ 336) dielectric relaxation (DR) measurements of (acetamide + LiBr/NO₃^-/ClO₄^-) deep eutectic solvents (DESs) in the frequency window of 0.2 ≤ ν (GHz) ≤ 50 and explore, via molecular dynamics simulations, the relative roles for the collective single-particle reorientational relaxations and the H-bond dynamics of acetamide in the measured DR response. In addition, DR measurements of neat molten acetamide were performed. Recorded DR spectra of these DESs require multi-Debye fits and produce well-separated DR time scales that are spread over several picoseconds to ∼1 ns. Simulations suggest DR time scales derive contributions from both the collective reorientational (Cl(t)) relaxation and structural H-bond (C_HB(t)) dynamics of acetamide. A good correlation between the measured and simulated activation energies further reveals a strong connection between the measured DR and the simulated Cl(t) and C_HB(t). Average DR times exhibit a strong fractional viscosity dependence, suggesting substantial microheterogeneity in these media. Simulations of Cl(t) and C_HB(t) reveal strong stretched exponential relaxations with a stretching exponent, 0.4 ≤ β ≤ 0.7. The ratio between the average reorientational correlation times of first and second ranks, ⟨τ⟩l=1/⟨τ⟩l=2, deviates appreciably from Debye's l(l+1) law for homogeneous media. Importantly, a pronounced translation-rotation decoupling between the simulated reorientation and center-of-mass diffusion times was observed.

The solvation structure, transport properties and reduction behavior of carbonate-based electrolytes of lithium-ion batteries

Despite the extensive employment of binary/ternary mixed-carbonate electrolytes (MCEs) for Li-ion batteries, the role of each ingredient with regards to the solvation structure, transport properties, and reduction behavior is not fully understood. Herein, we report the atomistic modeling and transport property measurements of the Gen2 (1.2 M LiPF₆ in ethylene carbonate (EC) and ethyl methyl carbonate (EMC)) and EC-base (1.2 M LiPF₆ in EC) electrolytes, as well as their mixtures with 10 mol% fluoroethylene carbonate (FEC). Due to the mixing of cyclic and linear carbonates, the Gen2 electrolyte is found to have a 60% lower ion dissociation rate and a 44% faster Li⁺ self-diffusion rate than the EC-base electrolyte, while the total ionic conductivities are similar. Moreover, we propose for the first time the anion–solvent exchange mechanism in MCEs with identified energetic and electrostatic origins. For electrolytes with additive, up to 25% FEC coordinates with Li⁺, which exhibits a preferential reduction that helps passivate the anode and facilitates an improved solid electrolyte interphase. The work provides a coherent computational framework for evaluating mixed electrolyte systems.

The different roles of the anion, cyclic and linear carbonates, and additive in mixed-carbonate electrolytes are revealed. The anion–solvent exchange mechanism and factors influencing the solid electrolyte interphase (SEI) formation are deciphered.

How Temperature, Pressure, and Salt Concentration Affect Correlations in LiTFSI/EMIM-TFSI Electrolytes: A Molecular Dynamics Study

Classical polarizable molecular dynamics simulations have been performed for LiTFSI solutions in the EMIM-TFSI ionic liquid. Different temperature or pressure values and salt concentrations have been examined. The structure and dynamics of the solvation shell of Li+ cations, diffusion coefficients of ions, conductivities of the electrolytes, and correlations between motions of ions have been analyzed. The results indicated that regardless of the conditions, significant correlations are present in all systems. The degree of correlations depends mainly on the salt fraction in the electrolyte and is much less affected by temperature and pressure changes. A positive correlation between motions of Li+ cations and TFSI anions, leading to the occurrence of negative Li+ transference numbers, exists for all conditions, although temperature and pressure changes affect the speed of anion exchange in Li+ solvation shells.

Effects of All-Atom Molecular Mechanics Force Fields on Amyloid Peptide Assembly: The Case of PHF6 Peptide of Tau Protein

Molecular dynamics (MD) simulations play a vital role in revealing the mechanism of amyloid aggregation that is crucial to the therapeutic agent development for Alzheimer's Disease. However, the accuracy of MD simulation results strongly depends on the force field employed. In our previous benchmark for 17 all-atom force fields on modeling of amyloid aggregation using the Aβ_16-22 dimer, we showed that AMBER14SB and CHARMM36m are suitable force fields for amyloid aggregation simulation, while GROMOS54a7 and OPLSAA are not good for the task. In this work, we continue assessing the applicability of atomistic force fields on amyloid aggregation using the VQIVYK (PHF6) peptide which is essential for tau-protein aggregation. Although, both Aβ_16-22 and PHF6 peptides formed fibrils in vitro, the PHF6 fibrils are parallel β-sheets, while the Aβ_16-22 fibrils are antiparallel β-sheets. We performed an all-atom large-scale MD simulation in explicit water on the PHF6 dimer and octa-peptides systems using five mainstream force fields, including AMBER99SB-disp, AMBER14SB, CHARMM36m, GROMOS54a7, and OPLSAA. The accumulated simulation time is 0.2 ms. Our result showed that the β-sheet structures of PHF6 peptides sampled by AMBER99SB-disp, AMBER14SB, GROMOS54a7, and OPLSAA are in favor of the antiparallel β-sheets, while the dominant type of β-sheet structures is parallel β-sheet by using CHARMM36m. Among the five force fields, CHARMM36m provides the strongest CH-π interaction that was observed in an NMR study. The comparison between our results and experimental observation indicates that CHARMM36m achieved the best performance on modeling the aggregation of PHF6 peptides. In summary, CHARMM36m is currently the most suitable force field for studying the aggregation of both amyloid-β and Tau through MD simulations.