Reversible Ionic Aggregation Kinetics in Concentrated Electrolytes
Pith reviewed 2026-05-15 14:48 UTC · model grok-4.3
The pith
A macroscopic rate equation for open and occupied association sites solves the reversible Smoluchowski aggregation equation and tracks ionic cluster changes after step changes in concentrated electrolytes.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The derived macroscopic rate equation of open/occupied association sites is a solution of the reversible Smoluchowski aggregation equation and predicts how ionic associations change subject to a step-change in conditions, with good qualitative agreement to atomistic MD simulations but quantitative differences highlighting multiple time scales.
What carries the argument
Macroscopic rate equation tracking the fractions of open and occupied association sites, which solves the reversible Smoluchowski aggregation equation for ionic clusters.
If this is right
- The formalism can be applied to investigate non-Newtonian flow behaviour in concentrated electrolytes.
- It provides a route to model double-layer charging dynamics under changing conditions.
- It describes the slow relaxation of electrolyte structure in confined geometries.
- It enables study of non-equilibrium behaviour without full atomistic simulation at each instant.
Where Pith is reading between the lines
- The same rate-equation structure may apply directly to related systems such as thermoreversible polymers and patchy particles.
- The observed quantitative gaps indicate that future refinements could incorporate at least two distinct relaxation times to improve accuracy for specific electrolyte compositions.
- The model offers a practical starting point for simulating electrolyte response in devices where sudden perturbations occur, such as during rapid charging or shear flow.
Load-bearing premise
The assumption that a single macroscopic rate equation for association sites captures the essential reversible aggregation kinetics without additional microscopic details or separate fitted time scales.
What would settle it
Time-resolved measurement of the fraction of occupied association sites after an abrupt change in ion concentration or temperature that deviates from the functional form predicted by the single rate equation beyond the reported quantitative mismatch.
Figures
read the original abstract
Here we develop a formalism for reversible ionic aggregation kinetics in an example concentrated electrolyte, building on previous equilibrium work of McEldrew and co-workers, and thermoreversible polymers and patchy particle systems. This is achieved through solving a macroscopic rate equation of open/occupied association sites, shown to be a solution of the reversible Smoluchowski aggregation equation, which predicts how ionic associations in electrolytes change subject to a step-change in conditions. We test the derived equations against atomistic molecular dynamics simulations of a salt-in-ionic liquid, where good qualitative agreement is obtained, but quantitative differences are found. This highlights the multiple time scales that exist in concentrated electrolytes, with a fast timescale preceding a longer timescale. We hope this formalism opens new avenues in understanding the dynamics and non-equilibrium behaviour in electrolytes. For example, the formalism can be developed further to investigate the non-Newtonian behaviour of concentrated electrolytes, double layer charging, and the slow dynamics of these electrolytes in confinement.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper develops a formalism for reversible ionic aggregation kinetics in concentrated electrolytes by deriving a macroscopic rate equation tracking open and occupied association sites. This equation is presented as a solution to the reversible Smoluchowski aggregation equation and is used to predict the response of ionic associations to step changes in conditions. The model is tested against atomistic MD simulations of a salt-in-ionic liquid, yielding qualitative agreement but quantitative differences that highlight multiple timescales in the system.
Significance. If the central derivation is shown to hold under the relevant approximations, the work supplies a useful reduced-order model for non-equilibrium ionic aggregation, extending equilibrium theories to dynamics. The parameter-free character of the derivation and the direct, falsifiable predictions for step-change responses are strengths. The qualitative match to MD supports potential applications to double-layer charging and non-Newtonian flow, though the quantitative discrepancies indicate the model captures only part of the multi-scale physics.
major comments (2)
- [Theory/derivation section (around the rate-equation introduction)] The claim that the macroscopic rate equation is a solution of the reversible Smoluchowski aggregation equation is central but lacks explicit derivation steps showing the required closure (e.g., size-independent kernel or neglect of higher-order correlations). In concentrated electrolytes, long-range Coulomb forces and chain/network aggregates likely violate standard closures, as implied by the reported quantitative mismatch with MD; this needs to be demonstrated or the claim qualified as approximate.
- [Results and MD comparison] The MD comparison (results section) reports good qualitative agreement but unresolved quantitative differences in relaxation timescales. This directly tests the model's predictive power; without addressing why a single ODE fails to capture the observed multi-timescale behavior (e.g., via additional terms or explicit justification of the closure), the validation remains incomplete for the claimed regime.
minor comments (3)
- [Abstract] The abstract states 'good qualitative agreement' without specifying the comparison metric or the magnitude of deviations; adding this would improve clarity.
- [Figures] Figure captions for the step-change predictions should explicitly list the initial and final conditions (concentration, temperature) and the observable plotted.
- [References] Citations to the Smoluchowski framework and McEldrew equilibrium work should include the precise references used for the reversible kernel.
Simulated Author's Rebuttal
We thank the referee for their careful reading of our manuscript and for providing constructive feedback. We address each of the major comments point by point below, and have made revisions to the manuscript where necessary to strengthen the presentation of the derivation and the discussion of the model's limitations.
read point-by-point responses
-
Referee: [Theory/derivation section (around the rate-equation introduction)] The claim that the macroscopic rate equation is a solution of the reversible Smoluchowski aggregation equation is central but lacks explicit derivation steps showing the required closure (e.g., size-independent kernel or neglect of higher-order correlations). In concentrated electrolytes, long-range Coulomb forces and chain/network aggregates likely violate standard closures, as implied by the reported quantitative mismatch with MD; this needs to be demonstrated or the claim qualified as approximate.
Authors: We appreciate the referee pointing out the need for more explicit steps in the derivation. In the revised manuscript, we have added a new appendix that provides a detailed derivation from the reversible Smoluchowski equation to the macroscopic rate equation, explicitly stating the closure approximations, including the use of a size-independent kernel and mean-field neglect of higher-order correlations. We agree that in concentrated electrolytes, long-range Coulomb forces and the possibility of chain or network aggregates can violate these standard closures, which likely contributes to the quantitative differences with MD. We have therefore qualified the claim in the main text as an approximate solution under the mean-field assumption, valid for capturing the primary aggregation kinetics. revision: yes
-
Referee: [Results and MD comparison] The MD comparison (results section) reports good qualitative agreement but unresolved quantitative differences in relaxation timescales. This directly tests the model's predictive power; without addressing why a single ODE fails to capture the observed multi-timescale behavior (e.g., via additional terms or explicit justification of the closure), the validation remains incomplete for the claimed regime.
Authors: We concur that the single ODE does not fully capture the multi-timescale relaxation seen in the MD simulations. The model is intended to describe the slower, reversible aggregation process, while the faster timescales arise from local ion dynamics and solvent effects not included in this reduced description. In the revised results section, we have added explicit discussion justifying the closure and explaining the origin of the timescale separation. We also outline how the formalism could be extended with additional terms or coupled equations to account for multiple modes, thereby improving the validation for the claimed regime of concentrated electrolytes. revision: yes
Circularity Check
Derivation chain is self-contained with no circular reductions
full rationale
The paper derives the macroscopic rate equation for open/occupied association sites and states that it is shown to be a solution of the reversible Smoluchowski aggregation equation, building on prior equilibrium work of McEldrew et al. and thermoreversible polymer frameworks. Validation occurs against independent atomistic MD simulations of a salt-in-ionic liquid, yielding qualitative agreement with noted quantitative differences attributed to multiple timescales. No load-bearing step reduces by construction to parameters fitted from the same validation data, no self-citation chain supplies an unverified uniqueness theorem, and no ansatz is smuggled via citation. The derivation remains independent of the target simulation results, satisfying the criteria for a non-circular finding.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption The macroscopic rate equation of open/occupied association sites is a solution of the reversible Smoluchowski aggregation equation
- domain assumption Qualitative features of ionic aggregation kinetics can be captured by tracking only the fraction of associated sites without resolving full microscopic configurations
Reference graph
Works this paper leans on
- [1]
-
[2]
L. Suo, O. Borodin, T. Gao, M. Olguin, J. Ho, X. Fan, C. Luo, C. Wang, and K. Xu, Science350, 938 (2015). 14
work page 2015
-
[3]
J. Wang, Y. Yamada, K. Sodeyama, C. H. Chiang, Y. Tateyama, and A. Yamada, Nat. Commun.7, 12032 (2016)
work page 2016
-
[4]
F. Wang, O. Borodin, T. Gao, X. Fan, W. Sun, F. Han, A. Faraone, J. A. Dura, K. Xu, and C. Wang, Nature Materials17, 543 (2018)
work page 2018
- [5]
-
[6]
L. Suo, O. Borodin, Y. Wang, X. Rong, W. Sun, X. Fan, S. Xu, M. A. Schroeder, A. V. Cresce, F. Wang,et al., Advanced Energy Materials7, 1701189 (2017)
work page 2017
-
[7]
Q. Dou, S. Lei, D.-W. Wang, Q. Zhang, D. Xiao, H. Guo, A. Wang, H. Yang, Y. Li, S. Shi, and X. Yan, Energy Environ. Sci.11, 3212 (2018)
work page 2018
-
[8]
A. A. Kornyshev, J. Phys. Chem. B111, 5545 (2007)
work page 2007
-
[9]
M. V. Fedorov and A. A. Kornyshev, Chem. Rev.114, 2978 (2014)
work page 2014
- [10]
-
[11]
J. P. Hallett and T. Welton, Chem. Rev.111, 3508–3576 (2011)
work page 2011
-
[12]
N. Molinari, J. P. Mailoa, and B. Kozinsky, J. Phys. Chem. Lett.10, 2313 (2019)
work page 2019
-
[13]
N. Molinari, J. P. Mailoa, N. Craig, J. Christensen, and B. Kozinsky, J. Power Sources428, 27 (2019)
work page 2019
- [14]
- [15]
-
[16]
C. M. Efaw, Q. Wu, N. Gao, Y. Zhang, H. Zhu, K. Ger- ing, M. F. Hurley, H. Xiong, E. Hu, X. Cao,et al., Nature Materials22, 1531 (2023)
work page 2023
- [17]
-
[18]
S. C. Kim, J. Wang, R. Xu, P. Zhang, Y. Chen, Z. Huang, Y. Yang, Z. Yu, S. T. Oyakhire, W. Zhang, L. C. Greenburg, M. S. Kim, D. T. Boyle, P. Sayavong, Y. Ye, J. Qin, Z. Bao, and Y. Cui, Nature Energy8, 814 (2023)
work page 2023
-
[19]
Y. Xie, J. Wang, B. H. Savitzky, Z. Chen, Y. Wang, S. Betzler, K. Bustillo, K. Persson, Y. Cui, L.-W. Wang, C. Ophus, P. Ercius, and H. Zheng, Sci. Adv.9, eadc9721 (2023)
work page 2023
- [20]
-
[21]
The Li-ion recharge- able battery: A perspective,
J. B. Goodenough and K. S. Park, “The Li-ion recharge- able battery: A perspective,” (2013)
work page 2013
-
[22]
M. Li, C. Wang, Z. Chen, K. Xu, and J. Lu, Chem. Rev. (2020)
work page 2020
-
[23]
M. A. Gebbie, B. Liu, W. Guo, S. R. Anderson, and S. G. Johnstone, ACS Catal.13, 16222 (2023)
work page 2023
- [24]
-
[25]
C. M. Phelan, E. Bj¨ orklund, J. Singh, M. Fraser, P. N. Didwal, G. J. Rees, Z. Ruff, P. Ferrer, D. C. Grin- ter, C. P. Grey,et al., Chemistry of Materials36, 3334 (2024)
work page 2024
-
[26]
C. Phelan, A. Bhandari, J. Singh, M. F. nad Erik Bj¨ orklund, J. Swallow, P. Bencok, C. Grey, Z. Goodwin, C. Skylaris, and R. Weatherup, https://chemrxiv.org/doi/pdf/10.26434/chemrxiv- 2025-k903g (2025)
-
[27]
Z. Yu, L. A. Curtiss, R. E. Winans, Y. Zhang, T. Li, and L. Cheng, The Journal of Physical Chemistry Letters 11, 1276 (2020)
work page 2020
-
[28]
O. Borodin, L. Suo, M. Gobet, X. Ren, F. Wang, A. Faraone, J. Peng, M. Olguin, M. Schroeder, M. S. Ding,et al., ACS nano11, 10462 (2017)
work page 2017
-
[29]
P. Lannelongue, R. Bouchal, E. Mourad, C. Bodin, M. Olarte, S. le Vot, F. Favier, and O. Fontaine, Jour- nal of The Electrochemical Society165, A657 (2018)
work page 2018
- [30]
-
[31]
L. G. Chagas, S. Jeong, I. Hasa, and S. Passerini, ACS Appl. Mater. Interfaces11, 22278 (2019)
work page 2019
-
[32]
T. C. Lourenco, L. G. Dias, and J. L. Da Silva, ACS Appl. Energy Mater.4, 4444 (2021)
work page 2021
-
[33]
R.-S. K¨ uhnel, D. Reber, and C. Battaglia, ACS Energy Lett.2, 2005 (2017)
work page 2005
- [34]
-
[35]
Z. Yu, N. P. Balsara, O. Borodin, A. A. Gewirth, N. T. Hahn, E. J. Maginn, K. A. Persson, V. Srinivasan, M. F. Toney, K. Xu, K. R. Zavadil, L. A. Curtiss, and L. Cheng, ACS Energy Lett.7, 461 (2022)
work page 2022
-
[36]
M. A. Gebbie, M. Valtiner, X. Banquy, E. T. Fox, W. A. Henderson, and J. N. Israelachvili, Proceedings of the National Academy of Sciences110, 9674 (2013)
work page 2013
-
[37]
M. A. Gebbie, H. A. Dobes, M. Valtiner, and J. N. Israelachvili, Proceedings of the National Academy of Sciences112, 7432–7437 (2015)
work page 2015
-
[38]
M. Han, H. Kim, C. Leal, M. Negrito, J. D. Batteas, and R. M. Espinosa-Marzal, Adv Mater Interfaces7, 2001313 (2020)
work page 2020
-
[39]
M. Han, R. Zhang, A. A. Gewirth, and R. M. Espinosa- Marzal, Nano Lett.21, 2304 (2021)
work page 2021
-
[40]
T. S. Groves, C. S. Perez-Martinez, R. Lhermerout, and S. Perkin, The Journal of Physical Chemistry Letters 12, 1702 (2021)
work page 2021
-
[41]
A. Lee, C. S. Perez-Martinez, A. M. Smith, S. Perkin, et al., Faraday discussions199, 239 (2017)
work page 2017
-
[42]
A. M. Smith, A. A. Lee, and S. Perkin, J. Phys. Chem. Lett.7, 2157 (2016)
work page 2016
- [43]
-
[44]
A. H¨ artel, M. B¨ ultmann, and F. Coupette, Physical Review Letters130, 108202 (2023)
work page 2023
-
[45]
Kjellander, Physical Chemistry Chemical Physics18, 18985 (2016)
R. Kjellander, Physical Chemistry Chemical Physics18, 18985 (2016)
work page 2016
-
[46]
Y. Avni, R. M. Adar, and D. Andelman, Physical Re- view E101, 010601 (2020)
work page 2020
-
[47]
R. M. Adar, S. A. Safran, H. Diamant, and D. Andel- man, Physical Review E100, 042615 (2019)
work page 2019
-
[48]
J. P. de Souza, Z. A. Goodwin, M. McEldrew, A. A. Kornyshev, and M. Z. Bazant, Phys. Rev. Lett.125, 116001 (2020)
work page 2020
-
[49]
E. Krucker-Velasquez and J. W. Swan, The Journal of Chemical Physics155, 134903 (2021)
work page 2021
-
[50]
T. Yamaguchi, A. Dukhin, Y.-J. Ryu, D. Zhang, O. Borodin, M. A. Gonz` alez, O. Yamamuro, D. L. Price, and M.-L. Saboungi, J. Phys. Chem. Lett.15, 76 (2024)
work page 2024
-
[51]
J. Lim, K. Park, H. Lee, J. Kim, K. Kwak, and M. Cho, Journal of the American Chemical Society140, 15661 (2018)
work page 2018
-
[52]
Y. S. Meng, V. Srinivasan, and K. Xu, Science378, 1065 (2022)
work page 2022
- [53]
-
[54]
O. Borodin, X. Ren, J. Vatamanu, A. von Wald Cresce, J. Knap, and K. Xu, Acc. Chem. Res.50, 2886 (2017)
work page 2017
-
[55]
Q. Wu, M. T. McDowell, and Y. Qi, JACS145, 2473 (2023)
work page 2023
- [56]
-
[57]
M. McEldrew, Z. A. Goodwin, S. Bi, M. Z. Bazant, and A. A. Kornyshev, J. Chem. Phys.152, 234506 (2020)
work page 2020
-
[58]
M. McEldrew,Ion Aggregation, Correlated Ion Trans- port and the Double Layer in Super-Concentrated Elec- trolytes, Ph.D. thesis, MIT (2021)
work page 2021
-
[59]
M. McEldrew, Z. A. H. Goodwin, H. Zhao, M. Z. Bazant, and A. A. Kornyshev, J. Phys. Chem B125, 2677–2689 (2021)
work page 2021
-
[60]
M. McEldrew, Z. A. H. Goodwin, N. Molinari, B. Kozin- sky, A. A. Kornyshev, and M. Z. Bazant, J. Phys. Chem. B125, 13752 (2021)
work page 2021
-
[61]
M. McEldrew, Z. A. Goodwin, S. Bi, A. Kornyshev, and M. Z. Bazant, J. Electrochem. Soc.168, 050514 (2021)
work page 2021
-
[62]
Z. A. H. Goodwin, M. McEldrew, B. Kozinsky, and M. Z. Bazant, PRX Energy2, 013007 (2023)
work page 2023
-
[63]
Z. A. Goodwin, M. McEldrew, J. P. de Souza, M. Z. Bazant, and A. A. Kornyshev, J. Chem. Phys.157, 094106 (2022)
work page 2022
-
[64]
Z. A. Goodwin and A. A. Kornyshev, Electrochim. Acta 434, 141163 (2022)
work page 2022
-
[65]
D. M. Markiewitz, Z. A. Goodwin, Q. Zheng, M. McEl- drew, R. M. Espinosa-Marzal, and M. Z. Bazant, ACS Applied Materials & Interfaces17, 29515 (2025)
work page 2025
-
[66]
D. M. Markiewitz, Z. A. Goodwin, M. McEldrew, J. P. de Souza, X. Zhang, R. M. Espinosa-Marzal, and M. Z. Bazant, Faraday Discussions , 365 (2024)
work page 2024
-
[67]
P. J. Flory, The Journal of chemical physics10, 51 (1942)
work page 1942
-
[68]
W. H. Stockmayer, The Journal of chemical physics11, 45 (1943)
work page 1943
-
[69]
W. H. Stockmayer, The Journal of Chemical Physics 12, 125 (1944)
work page 1944
-
[70]
W. H. Stockmayer, J. Polym. Sci.9, 69 (1952)
work page 1952
- [71]
- [72]
- [73]
- [74]
- [75]
-
[76]
Tanaka, Physica A: Statistical Mechanics and its Ap- plications257, 245 (1998)
F. Tanaka, Physica A: Statistical Mechanics and its Ap- plications257, 245 (1998)
work page 1998
- [77]
- [78]
-
[79]
F. Sciortino and E. Zaccarelli, Current Opinion in Solid State and Materials Science15, 246 (2011)
work page 2011
-
[80]
F. Sciortino, C. D. Michele, S. Corezzi, J. Russo, E. Za- ccarelli, and P. Tartaglia, Soft Matter5, 2571 (2009)
work page 2009
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.