Field-driven Ion Pairing Dynamics in Concentrated Electrolytes

We investigate ion pairing dynamics in electrolytes driven far from equilibrium using molecular simulations and nonequilibrium rate theory. Focusing on 0.5 M $\mathrm{LiPF_6}$ in water and acetonitrile under uniform electric fields, we compute transition path theory observables including reactive fluxes and mean first-passage times of ion pairing. Moreover, we introduce a dynamical proxy of free-ion population, where its field-induced change is strongly correlated with the nonlinear enhancement of conductivity, yielding an increase of $40 \ %$ at 50 mV/Å in acetonitrile, compared to less than $10 \ %$ in aqueous electrolytes. Further kinetic analysis elucidates that Onsager’s classical theory substantially overestimates field-induced enhancement of ion pair dissociation in molecular electrolytes. This discrepancy arises from solvent-mediated dynamical pathways and field-induced dielectric decrement that suppress ion pair dissociation within explicit solvents, highlighting that a faithful description of molecular details is essential. Our results provide a molecular interpretation of nonlinear electrolyte transport beyond continuum theories and establish a general framework for quantifying nonequilibrium reaction kinetics in condensed phase systems.

💡 Research Summary

This paper investigates how strong electric fields affect ion‑pair dynamics in concentrated electrolytes, using 0.5 M LiPF₆ dissolved in water and acetonitrile as model systems. Classical continuum theories such as Onsager’s description of the second Wien effect predict a field‑enhanced dissociation constant based solely on electrostatic length scales, but they neglect solvent structure, dielectric response, and collective dynamics that become crucial at high concentration. To address these gaps, the authors perform extensive nonequilibrium molecular dynamics simulations (LAMMPS, SPC/fw water, OPLS‑AA force fields) at 300 K with uniform fields ranging from 0 to 50 mV Å⁻¹. Each trajectory is equilibrated for 10 ns and then sampled for 100 ns, providing a statistically robust steady‑state ensemble.

A central methodological advance is the definition of a cluster‑based nearest‑counterion distance, R