[論文レビュー] Reversible Ionic Aggregation Kinetics in Concentrated Electrolytes

Zachary A. H. Goodwin|arXiv (Cornell University)|Mar 7, 2026

Electrostatics and Colloid Interactions被引用数 0

ひとこと要約

The paper develops a macroscopic rate-equation formalism for reversible ionic aggregation in concentrated electrolytes (specifically salt-in-ionic liquids) and tests it against atomistic MD simulations, revealing qualitative agreement and multi-time-scale behavior.

ABSTRACT

Here we develop and test a formalism for reversible ionic aggregation kinetics in an example concentrated electrolyte. Specifically, building on previous equilibrium work of McEldrew and co-workers in the context of concentrated electrolytes, and non-equilibrium properties of thermoreversible polymers and patchy particle systems, we develop the formalism for how ionic associations in electrolytes change subject to a step-change in conditions. This is achieved through solving a macroscopic rate equation of open/occupied association sites, which is a solution of the reversible Smoluchowski aggregation equation. We compare the derived equations against atomistic molecular dynamics simulations of a salt-in-ionic liquid. Good qualitative agreement is obtained, but quantitative differences are found, which highlights the multiple time scales of the associations that exist in concentrated electrolytes. We hope this formalism acts as the starting point for investigating these properties in other electrolytes, and developing it further to investigate the non-Newtonian behaviour of concentrated electrolytes, double layer charging, and the slow dynamics of these electrolytes in confinement.

研究の動機と目的

• Extend non-equilibrium formalism from polymer/patchy-particle systems to concentrated electrolytes.
• Link ionic association kinetics to observable properties in concentrated electrolytes such as SiILs.
• Provide a framework to study time-dependent aggregation after step-changes in conditions.
• Assess the validity and limits of a single-rate-constant description for aggregation kinetics.

提案手法

• Adopt a thermoreversible aggregation formalism for Cayley-tree ionic clusters forming in a salt-in-ionic liquid.
• Derive a macroscopic rate equation for association-site occupancy p(t) from the reversible Smoluchowski equation.
• Express cluster distributions with equilibrium constants K_lm and detailed free-energy contributions (comb, conf, bind).
• Solve for time-dependent association probabilities using mass-action-like relations and initial conditions after step-changes.
• Extend symmetric and asymmetric cases to handle different association-site pools (ψ_+, ψ_-).
• Test theoretical predictions against atomistic MD simulations of NaTFSI in EMIMTFSI to assess qualitative/quantitative agreement.

Figure 1: Schematic of non-equilibrium kinetics of ionic aggregation in the studied salt-in-ionic liquid (SIIL). The alkali metal cation and IL anion bind together to form aggregates. From step-changing a property of the system, we investigate how they evolve in time. The alkali metal cation and IL

実験結果

リサーチクエスチョン

• RQ1 How do reversible ionic aggregation kinetics evolve in concentrated electrolytes after a step-change in conditions?
• RQ2 Can a macroscopic rate equation for association-site occupancy capture the non-equilibrium dynamics of ionic aggregation?
• RQ3 What are the limits of a single-rate-constant description for aggregation across time scales in concentrated electrolytes?
• RQ4 How well does the theory describe both symmetric and asymmetric association scenarios in salt-in-ionic liquids?

主な発見

• Qualitative agreement between the reversible Smoluchowski-based kinetics and MD simulations of a salt-in-ionic liquid.
• Quantitative differences indicate multiple timescales of ionic associations, challenging single-rate descriptions.
• The framework links time-dependent association probabilities to equilibrium cluster distributions, embedding non-equilibrium evolution in a quasi-equilibrium picture.
• Analytical solutions show how relaxation times depend on association rates and current association probabilities, predicting possible very long relaxation times as p approaches unity.
• Symmetric case yields a closed-form p(t) solution; asymmetric case extends to p_+- and p_-+ with analogous relaxation forms.
• The study highlights the multi-time-scale nature of concentrated electrolytes and provides a starting point for exploring non-Newtonian behavior and interfacial dynamics in these systems.

Figure 2: Coordination numbers between Na + and TFSI - in NaTFSI 0.75 EMIMTFSI 0.25 as a function of time, for the charge rescaling case, at the indicated temperatures, from MD simulations and theory.