# Parametrizing reaction probabilities for proton transfers in protic ionic liquids

## TL;DR

This paper studies proton transfer in protic ionic liquids using quantum mechanics to model reaction probabilities for use in simulations.

## Contribution

The study introduces a classical kinetic model for proton transfer probabilities in protic ionic liquids based on Morse and hyperbolic tangent functions.

## Key findings

- Proton transfer reactions in 1-methylimidazolium carboxylates have small or negligible energy barriers.
- Reaction probabilities are primarily governed by thermally activated hopping rather than quantum tunneling.
- Morse potentials and hyperbolic tangent functions accurately model the energy profiles and distance-dependent probabilities.

## Abstract

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