Inferring Thermal Dissociation Kinetics of Ion Clusters from Molecular Dynamics
Abstract
Electrospray ionic liquid ion sources operating in the pure ion regime emit metastable, low solvation number ion clusters whose post-emission dissociation modifies plume composition, energy distributions, and diagnostic observables. While field-free measurements support first-order, thermally activated decay on microsecond timescales, the acceleration region is characterized by rapidly varying electric fields and nanosecond residence times such that dissociation becomes strongly field-enhanced. In this work, microcanonical ensemble molecular dynamics (MD) trajectories are used to infer dissociation kinetics of EMI-BF4 ion clusters under controlled internal energy and applied electric field. Positive and negative dimers, trimers, and tetramers are simulated at temperatures between 600 to 1000 K and uniform fields of 10^6 to 10^9 V/m. Dissociation lifetimes are extracted using a connectivity-based fragmentation criterion, and product channels are classified to obtain pathway-resolved branching probabilities. Across all solvation numbers, lifetimes are field independent in the weak field limit, but collapse by orders of magnitude above 10^8 V/m. Cluster polarity dependence is most prominent in the temperature-controlled regime and diminish as electrostatic work dominates. The MD simulations further reveal that field-driven changes in dissociation topology, including a transition in trimer breakup from neutral pair evaporation to charged core ion ejection, and three competing tetramer pathways involving single neutral loss, double neutral emission, and core ion ejection. The MD dataset is reduced to compact parameterizations of {τ_n}(E,T), yielding transferable kinetic inputs for multiscale plume transport models.
Related
- Modelling Energy Transfer Efficiency in Collision Induced Dissociation and Unfolding Using Molecular Dynamics Simulations
- Machine learning kinetics from molecular dynamics data
- Beyond CO₂ dissociation: The pivotal role of O* surface species in reverse water-gas shift kinetics over nickel
- Integer electron-count cuts in variational second-order reduced-density-matrix theory. I. Molecular dissociation with fixed real-space and orbital partitions
- Computational Insights into Mechanostability and Dissociation Dynamics of the Dengue Virus Envelope Protein Ectodomain Dimer Across pH and Temperature Gradients
BibTeX
@misc{laws2026inferring,
title = {{Inferring Thermal Dissociation Kinetics of Ion Clusters from Molecular Dynamics}},
author = {Nicholas Laws and Adler Smith and Elaine Petro},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.25070},
eprint = {2609.25070},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.25070},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - Inferring Thermal Dissociation Kinetics of Ion Clusters from Molecular Dynamics
AU - Nicholas Laws
AU - Adler Smith
AU - Elaine Petro
PY - 2026
DA - 2026/09/17
PB - arXiv
DO - 10.48550/arxiv.2609.25070
UR - https://arxiv.org/abs/2609.25070
AB - Electrospray ionic liquid ion sources operating in the pure ion regime emit metastable, low solvation number ion clusters whose post-emission dissociation modifies plume composition, energy distributions, and diagnostic observables. While field-free measurements support first-order, thermally activated decay on microsecond timescales, the acceleration region is characterized by rapidly varying electric fields and nanosecond residence times such that dissociation becomes strongly field-enhanced. In this work, microcanonical ensemble molecular dynamics (MD) trajectories are used to infer dissociation kinetics of EMI-BF4 ion clusters under controlled internal energy and applied electric field. Positive and negative dimers, trimers, and tetramers are simulated at temperatures between 600 to 1000 K and uniform fields of 10^6 to 10^9 V/m. Dissociation lifetimes are extracted using a connectivity-based fragmentation criterion, and product channels are classified to obtain pathway-resolved branching probabilities. Across all solvation numbers, lifetimes are field independent in the weak field limit, but collapse by orders of magnitude above 10^8 V/m. Cluster polarity dependence is most prominent in the temperature-controlled regime and diminish as electrostatic work dominates. The MD simulations further reveal that field-driven changes in dissociation topology, including a transition in trimer breakup from neutral pair evaporation to charged core ion ejection, and three competing tetramer pathways involving single neutral loss, double neutral emission, and core ion ejection. The MD dataset is reduced to compact parameterizations of {τ_n}(E,T), yielding transferable kinetic inputs for multiscale plume transport models.
ER -