🔬 condensed matter

Reversible Ionic Aggregation Kinetics in Concentrated Electrolytes

This paper develops and validates a formalism for modeling reversible ionic aggregation kinetics in concentrated electrolytes by solving a macroscopic rate equation derived from the Smoluchowski aggregation equation, which shows qualitative agreement with atomistic molecular dynamics simulations while revealing complex multi-timescale dynamics that warrant further investigation into non-Newtonian behavior and confinement effects.

Original authors: Zachary A. H. Goodwin

Published 2026-03-10

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Zachary A. H. Goodwin

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine you are at a crowded party. In a normal room, people mingle freely, chatting with whoever they like. But in a concentrated electrolyte (like the special liquids used in next-generation batteries), the room is so packed that people are constantly bumping into each other, forming tight little groups or "cliques."

This paper is about understanding how fast these cliques form and break apart when the conditions of the party suddenly change.

Here is a breakdown of the research using simple analogies:

1. The Setting: The "Salt-in-Ionic Liquid" Party

The researchers studied a specific type of liquid called a Salt-in-Ionic Liquid (SiIL).

The Guests: Imagine two types of guests:
- The "Active" Guests (Alkali Metal Cations): These are the ones who want to form groups.
- The "Background" Guests (Ionic Liquid Cations): These are the wallflowers who mostly just watch and don't join the groups directly, but they affect the atmosphere.
- The "Connectors" (Anions): These are the social butterflies that link the Active Guests together.
The Goal: In a battery, these "Active Guests" need to move around to charge and discharge the device. But if they get stuck in too many cliques (aggregates), the battery slows down.

2. The Problem: Predicting the Chaos

Scientists already had a map of what the party looks like when it's calm and steady (equilibrium). They knew how many cliques existed and how big they were.

But what happens if you suddenly turn up the music or change the temperature? How fast do the cliques break up? How fast do new ones form?

The Old Way: Scientists usually just guessed or looked at snapshots.
The New Way (This Paper): The author, Zachary Goodwin, built a mathematical recipe (a formalism) to predict exactly how these groups evolve over time, second by second.

3. The Theory: The "Velcro" Model

The theory treats these ions like pieces of Velcro.

Each ion has a certain number of "hooks" (binding sites).
When they bump into each other, they stick (aggregate).
When they get pulled apart, they unstick (dissociate).

The math predicts that if you suddenly change the environment (like turning off the electricity in the simulation to see how the groups react), the system should settle into a new pattern at a specific speed. It's like predicting how long it takes for a crowd to reorganize after a fire alarm goes off.

4. The Experiment: The Digital Dance Floor

To test their math, the researchers ran computer simulations (Molecular Dynamics).

They created a virtual room with thousands of these ions.
They suddenly "turned off" the electrical attraction between them (making them ignore each other) and then "turned it back on" instantly.
They watched how fast the ions ran back to their friends to form cliques again.

5. The Surprise: One Speed vs. Many Speeds

Here is where the story gets interesting.

The Theory's Prediction: The math assumed everyone moves at one single speed. It predicted that the cliques would form smoothly and steadily, like a line of dominoes falling.
The Reality: The computer simulation showed something more chaotic.
- Phase 1 (The Sprint): Immediately after the change, ions that were already standing next to each other snapped together super fast (in a few picoseconds). It was a frantic sprint.
- Phase 2 (The Shuffle): After that initial burst, the process slowed down dramatically. The ions had to physically swim through the crowded liquid to find new partners. This was much slower, taking nanoseconds.

The Analogy: Imagine a group of people trying to form a conga line.

The Theory thought: "Everyone will grab the hand of the person next to them at a steady, average pace."
The Reality was: "The people standing right next to each other grabbed hands instantly (Fast Phase), but the people on the other side of the room had to weave through the crowd, which took forever (Slow Phase)."

6. Why This Matters

The researchers found that their math was qualitatively correct (it got the general shape of the curve right) but quantitatively off (it couldn't predict the exact timing because it missed the "two-speed" nature of the process).

Why does this matter for you?

Better Batteries: If we understand that ions move in "sprints" and "shuffles," we can design batteries that charge faster.
Non-Newtonian Fluids: These liquids act like weird fluids that get thicker or thinner depending on how you push them (like ketchup or oobleck). Understanding the "speed" of the cliques helps explain why.
Future Tech: This math provides a starting point to understand how these liquids behave in tiny, confined spaces (like inside a battery electrode), which is crucial for making smaller, more powerful devices.

The Bottom Line

The paper is a first draft of a new rulebook for how ions dance in crowded liquids. It successfully predicted the general rhythm of the dance but discovered that the dancers have two different tempos: a fast sprint for neighbors and a slow shuffle for everyone else. This discovery opens the door to designing better energy storage systems by accounting for these different speeds.

1. Problem Statement

Concentrated electrolytes, such as salt-in-ionic liquids (SiILs), water-in-salt electrolytes (WiSEs), and pure ionic liquids (ILs), exhibit complex non-equilibrium behaviors critical for energy storage applications (e.g., batteries, supercapacitors). While equilibrium theories (specifically the thermoreversible ionic aggregation model by McEldrew et al.) successfully describe static properties like activity coefficients and transference numbers, there is a significant gap in understanding the kinetics of ionic aggregation.

Key challenges include:

Non-Newtonian behavior: These electrolytes often display viscoelastic properties and anomalous transport mechanisms.
Lack of Microscopic Kinetic Theory: Existing models link coordination environments to equilibrium properties but fail to quantitatively describe the time-dependent evolution of ionic associations following a perturbation (e.g., a step-change in temperature or field).
Multi-scale Dynamics: It is known that desolvation and ion association occur over multiple time scales, but a unified formalism linking these to macroscopic rate equations has been missing.

2. Methodology

The authors developed a theoretical formalism and validated it against atomistic Molecular Dynamics (MD) simulations.

A. Theoretical Framework

The study extends the equilibrium theory of thermoreversible aggregation (based on Flory-Stockmayer polymer theory and Cayley tree approximations) to non-equilibrium conditions.

System: Salt-in-ionic liquids (SiILs), modeled as a mixture of alkali metal cations ( $+$ ), ionic liquid cations ( $\oplus$ ), and anions ($-$). The IL cations are treated as non-associating but influencing the solvent environment.
Core Equations: The authors derived a macroscopic rate equation for the time evolution of association probabilities ( $p(t)$ $p (t)$ ) by adapting the reversible Smoluchowski aggregation equation.
- They utilized the "RAf" (Reversible Aggregation with fixed functionality) model, where the coalescence rate depends on the number of dangling association sites ( $\sigma_{lm}$ ) on clusters.
- The time evolution is governed by the balance between the formation rate ( $k$ ) and dissociation rate ( $f$ ) of single associations.
Cases Analyzed:
1. Symmetric Case: Where the concentration of association sites for cations and anions is equal ( $\psi_+ = \psi_-$ ), leading to a single association probability $p(t)$ .
2. Asymmetric Case: Where $\psi_+ \neq \psi_-$ , requiring coupled differential equations for $p_{+-}(t)$ and $p_{-+}(t)$ .
Quasi-Equilibrium Assumption: The theory assumes that the cluster distribution adjusts instantaneously to the changing association probabilities, maintaining a quasi-equilibrium state described by the equilibrium cluster distribution formula but with time-dependent parameters.

B. Computational Validation

System: NaTFSI dissolved in EMIMTFSI (an ionic liquid).
Simulation Tool: LAMMPS with the CL&P force field (charges rescaled to 0.7 to mimic polarizability).
Non-Equilibrium Protocol: To test kinetics, the authors performed a "charge rescaling" experiment:
1. Equilibrated the system at 300 K.
2. Scaled charges to near-zero (effectively breaking ionic associations) and equilibrated.
3. Instantly restored full charges (step-change) and tracked the re-formation of coordination numbers and aggregates over time.
Metrics: Tracked coordination numbers (CN), association probabilities, cluster size distributions ( $c_{lm}$ ), and gelation criteria (percolation thresholds).

3. Key Contributions

Formalism Development: Successfully adapted the reversible Smoluchowski equation from patchy particle systems and thermoreversible polymers to concentrated electrolytes. This provides the first microscopic theory linking ionic association kinetics to macroscopic rate equations in this context.
Quasi-Equilibrium Solution: Demonstrated that the time-dependent association probabilities $p(t)$ can serve as the central variable to describe the evolution of the entire cluster distribution, effectively connecting time and temperature through the probability of bond formation.
Identification of Multi-Time Scales: The comparison between theory and simulation revealed a fundamental limitation in the single-rate-constant assumption. The study identified that ionic aggregation in concentrated electrolytes is governed by multiple time scales:
- A fast time scale ( $\sim$ ps) driven by local re-association of ions already in proximity.
- A slow time scale ( $\sim$ ns) limited by diffusion and structural rearrangement of the electrolyte.

4. Results

Qualitative Agreement: The derived equations successfully captured the qualitative trends of the MD simulations. The theory correctly predicted the rapid initial increase in coordination numbers followed by a slower approach to equilibrium, as well as the formation of percolating ionic networks (gels).
Quantitative Discrepancy:
- The theory, assuming a single rate constant ( $k$ ), under-predicted coordination numbers at short times (fast regime) and over-predicted them at long times (slow regime).
- Fitting a single decay constant ( $\tau$ ) to the simulation data yielded values that did not perfectly match the entire time evolution curve.
Gelation Dynamics:
- The theory predicted the time to reach the percolation threshold (gel point) based on the product of association probabilities ( $p_{+-}p_{-+}$ ).
- MD simulations confirmed the formation of a percolating network, but the timing varied slightly due to the multi-scale nature of the kinetics.
- Cluster distribution analysis showed that at the gel point, large aggregates (up to 30+ ions) formed, consistent with percolation theory, though the theory slightly underestimated the speed of this formation.
Temperature and Composition Dependence:
- Increasing temperature decreased the association constant ( $\lambda$ ) and the relaxation time ( $\tau$ ), as expected.
- However, the rate of association formation ( $k$ ) did not follow a simple Arrhenius (activated) behavior, suggesting the process is not dominated by a single high energy barrier but rather by complex solvent reorganization.
- Higher salt concentrations increased the association constant but decreased the relaxation time (slower kinetics), indicating stronger but more sluggish network formation.

5. Significance and Future Outlook

Bridging Theory and Simulation: This work provides a crucial link between microscopic ionic associations and non-equilibrium structural changes, offering a framework to interpret residence times and transport properties.
Transport Mechanisms: The identification of multi-time scales has profound implications for understanding ion transport. The fast time scale suggests structural diffusion (hopping) may dominate over vehicular motion in the short term, while the slow time scale relates to bulk diffusion.
Applications:
- Battery Kinetics: The formalism can be extended to study desolvation dynamics at electrode interfaces, a critical factor in intercalation rates.
- Double Layer Charging: It offers a pathway to model the charging dynamics of electrical double layers (EDL) in concentrated electrolytes, accounting for the relaxation of ionic aggregates.
- Non-Newtonian Behavior: The theory lays the groundwork for predicting the frequency-dependent viscoelastic response of these electrolytes, potentially explaining the non-Newtonian behavior observed in WiSEs.
- Confinement Effects: The authors suggest extending this formalism to confined geometries (e.g., nanopores) to investigate "anomalous underscreening" and confinement-induced gelation.

In conclusion, while the single-rate-constant approximation has quantitative limitations due to the complex, multi-scale nature of concentrated electrolytes, the developed formalism serves as a robust starting point for investigating non-equilibrium phenomena, linking coordination environments to macroscopic transport and rheological properties.