Parametrizing reaction probabilities for proton transfers in protic ionic liquids

Abstract

Protic ionic liquids are promising electrolytes for electrochemical applications owing to their 3 intrinsic proton conductivity, but quantitative understanding of the underlying proton transfer 4 processes remains limited. Here, we present a systematic quantum-mechanical investigation 5 of proton transfer in a series of 1-methylimidazolium carboxylates, with the specific goal of 6 parametrizing reaction probabilities for use in reactive molecular dynamics simulations. Density 7 functional theory scans were performed to map relaxed potential energy surfaces along two 8 collective variables, the donor–acceptor distance and the proton transfer coordinate. The resulting 9 energy profiles were accurately represented by Morse potentials. 10 From the donor–acceptor distance scans, the distance-dependent reaction probabilities were 11 fitted using a hyperbolic tangent function. Analysis of the proton transfer coordinate revealed 12 small or even negligible energy barriers for the proton transfer reactions, which in turn resulted in 13 low empirical valence bond coupling energies between the reactant and product states. Quantum 14 tunneling effects appear to play only a minor role in these processes. Consequently, the proton 15 transfer reaction probabilities are predominantly governed by thermally activated hopping events, 16 which are captured within a classical kinetic model framework.