A Computer Simulation Model for Proton Transport in Liquid Imidazole

Electronic Fingerprint of the Protonated Imidazole Dimer Probed by X-ray Absorption Spectroscopy

Protons in low-barrier superstrong hydrogen bonds are typically delocalized between two electronegative atoms. Conventional methods to characterize such superstrong hydrogen bonds are vibrational spectroscopy and diffraction techniques. We introduce soft X-ray spectroscopy to uncover the electronic fingerprints for proton sharing in the protonated imidazole dimer, a prototypical building block enabling effective proton transport in biology and high-temperature fuel cells. Using nitrogen core excitations as a sensitive probe for the protonation status, we identify the X-ray signature of a shared proton in the solvated imidazole dimer in a combined experimental and theoretical approach. The degree of proton sharing is examined as a function of structural variations that modify the shape of the low-barrier potential in the superstrong hydrogen bond. We conclude by showing how the sensitivity to the quantum distribution of proton motion in the double-well potential is reflected in the spectral signature of the shared proton.

Ab initio molecular dynamics study of proton transport in imidazolium-based ionic liquids with added imidazole

Development of efficient anhydrous proton-conducting materials would expand the operational temperature ranges of hydrogen fuels cells (HFCs) and eliminate their dependence on maintaining sufficient hydration levels to function efficiently. Protic ionic liquids (PILs), which have high ionic densities and low vapor pressures, have emerged as a potential material for proton conducting layers in HFCs. In this work, we investigate proton transport via the Grotthuss mechanism in 1-ethylimidazolium bis-(trifluoromethanesulfonyl)imide ([C₂HIm][TFSI]) protic ionic liquids with added imidazole (Im⁰) using ab initio molecular dynamics. In particular, we vary the composition of the systems studied from pure [C₂HIm][TFSI] to those where the mole fraction of Im⁰ is 0.67. Given the large difference in pK_a between C₂HIm⁺ and HTFSI, TFSI^- does not accept acidic protons from C₂HIm⁺; conversely, imidazolium (HIm⁺) and C₂HIm⁺ have very similar pK_a values, and thus Im⁰ can readily accept protons. We find that the unprotonated nitrogen on Im⁰ dominates solvation of the labile protons on C₂HIm⁺ and other Im⁰ species, resulting in formation of robust imidazole wires. Given the amphoteric nature of Im⁰, i.e. its ability to accept and donate protons, these wires provide conduits along which protons can rapidly traverse via the Grotthuss mechanism, thereby greatly increasing the proton coefficient of self-diffusion. We find that the average length of the wires increases with added Im⁰, and thus as the mole fraction of Im⁰ increases so too does the proton diffusion constant. Lastly, we analyze our trajectories to determine the energy and time scales associated with proton transfer.

Origin of low melting point of ionic liquids: dominant role of entropy

Large structural entropy makes salts liquid at room temperature.

Characterizing and Contrasting Structural Proton Transport Mechanisms in Azole Hydrogen Bond Networks Using Ab Initio Molecular Dynamics

Imidazole and 1,2,3-triazole are promising hydrogen-bonded heterocycles that conduct protons via a structural mechanism and whose derivatives are present in systems ranging from biological proton channels to proton exchange membrane fuel cells. Here, we leverage multiple time-stepping to perform ab initio molecular dynamics of imidazole and 1,2,3-triazole at the nanosecond time scale. We show that despite the close structural similarities of these compounds, their proton diffusion constants vary by over an order of magnitude. Our simulations reveal the reasons for these differences in diffusion constants, which range from the degree of hydrogen-bonded chain linearity to the effect of the central nitrogen atom in 1,2,3-triazole on proton transport. In particular, we uncover evidence of two "blocking" mechanisms in 1,2,3-triazole, where covalent and hydrogen bonds formed by the central nitrogen atom limit the mobility of protons. Our simulations thus provide insights into the origins of the experimentally observed 10-fold difference in proton conductivity.

Physics-based, neural network force fields for reactive molecular dynamics: Investigation of carbene formation from [EMIM+][OAc-]

Reactive molecular dynamics simulations enable a detailed understanding of solvent effects on chemical reaction mechanisms and reaction rates. While classical molecular dynamics using reactive force fields allows significantly longer simulation time scales and larger system sizes compared with ab initio molecular dynamics, constructing reactive force fields is a difficult and complex task. In this work, we describe a general approach following the empirical valence bond framework for constructing ab initio reactive force fields for condensed phase simulations by combining physics-based methods with neural networks (PB/NNs). The physics-based terms ensure the correct asymptotic behavior of electrostatic, polarization, and dispersion interactions and are compatible with existing solvent force fields. NNs are utilized for a versatile description of short-range orbital interactions within the transition state region and accurate rendering of vibrational motion of the reacting complex. We demonstrate our methodology for a simple deprotonation reaction of the 1-ethyl-3-methylimidazolium cation with acetate to form 1-ethyl-3-methylimidazol-2-ylidene and acetic acid. Our PB/NN force field exhibits ∼1 kJ mol^-1 mean absolute error accuracy within the transition state region for the gas-phase complex. To characterize the solvent modulation of the reaction profile, we compute potentials of mean force for the gas-phase reaction as well as the reaction within a four-ion cluster and benchmark against ab initio molecular dynamics simulations. We find that the surrounding ionic environment significantly destabilizes the formation of the carbene product, and we show that this effect is accurately captured by the reactive force field. By construction, the PB/NN potential may be directly employed for simulations of other solvents/chemical environments without additional parameterization.

Proton Conduction Mechanism for Anhydrous Imidazolium Hydrogen Succinate Based on Local Structures and Molecular Dynamics

Anhydrous organic crystalline materials incorporating imidazolium hydrogen succinate (Im-Suc), which exhibit high proton conduction even at temperatures above 100 °C, are attractive for elucidating proton conduction mechanisms toward the development of solid electrolytes for fuel cells. Herein, quantum chemical calculations were used to investigate the proton conduction mechanism in terms of hydrogen-bonding (H-bonding) changes and restricted molecular rotation in Im-Suc. The local H-bond structures for proton conduction were characterized by vibrational frequency analysis and compared with corresponding experimental data. The calculated potential energy surface involving proton transfer (PT) and imidazole (Im) rotational motion showed that PT between Im and succinic acid was a rate-limiting step for proton transport in Im-Suc and that proton conduction proceeded via the successive coupling of PT and Im rotational motion based on a Grotthuss-type mechanism. These findings provide molecular-level insights into proton conduction mechanisms for Im-based (or -incorporated) H-bonding organic proton conductors.

Elucidating the Proton Transport Pathways in Liquid Imidazole with First-Principles Molecular Dynamics

Imidazole is a promising anhydrous proton conductor with a high conductivity comparable to that of water at a similar temperature relative to its melting point. Previous theoretical studies of the mechanism of proton transport in imidazole have relied either on empirical models or on ab initio trajectories that have been too short to draw significant conclusions. Here, we present the results of multiple time-step ab initio molecular dynamics simulations of an excess proton in liquid imidazole reaching 1 ns in total simulation time. We find that the proton transport is dominated by structural diffusion, with the diffusion constant of the proton defect being ∼8 times higher than that of self-diffusion of the imidazole molecules. By using correlation function analysis, we decompose the mechanism for proton transport into a series of first-order processes and show that the proton transport mechanism occurs over three distinct time and length scales. Although the mechanism at intermediate times is dominated by hopping along pseudo-one-dimensional chains, at longer times the overall rate of diffusion is limited by the re-formation of these chains. These results provide a more complete picture of the traditional idealized Grotthuss structural diffusion mechanism.

Proton Propensity and Orientation of Imidazolium Cation at Liquid Imidazole-Vacuum Interface: A Molecular Dynamics Simulation

Imidazole has gained attention as an alternative to anhydrous proton conductor in high-temperature proton exchange membrane fuel cells. A detailed investigation of proton propensity and the orientation of the imidazolium cation at the liquid-vacuum interface is important for understanding the interfacial properties of imidazole-based proton-conductive materials. Here, we perform all-atom molecular dynamics simulation on a slab model of the liquid imidazole-vacuum interface. Proton transportation process between the imidazolium cation and neutral imidazole molecules is described by the multistate empirical valence bond model of imidazole developed previously. The imidazolium cation shows a tendency to stay in the bulk region rather than at the outermost surface, and the NN vectors and norm vectors of both the imidazolium cation and imidazole molecules are more probable to be perpendicular to the surface normal vector at the interface than in the bulk. The orientation of the hydrogen bond cluster shows the same tendency as the NN vectors, which indicates that proton transportation along the direction of the surface normal vector is hindered. The instantaneous surface analyses show that the fluctuation is depressed when the imidazolium cation is near the outermost surface, which makes it less favorable for the cation appearing at the interface.

Equilibrium Structure, Hydrogen Bonding, and Proton Conductivity in Half-Neutralized Diamine Ionic Liquids

Recent experiments on proton conducting ionic liquids point to half-neutralized diamine-triflate salts as promising candidates for applications in power generation and energy conversion electrochemical devices. Structural and dynamical properties of the simplest among these compounds are investigated by a combination of density functional theory (DFT) and molecular dynamics (MD) simulations based on an empirical force field. Three different cations have been considered, consisting of a pair of amine-ammonium terminations joined by a short aliphatic segment -(CH₂) _n- with n = 2, 3, and 4. First, the ground state structure, vibrational eigenstates, and hydrogen-bonding properties of single ions, neutral ion pairs, small neutral aggregates of up to eight ions, and molecularly thin hydrogen bonded wires have been investigated by DFT computations. Second, structural and dynamical properties of homogeneous liquid and amorphous phases are investigated by MD simulations over the temperature range of 200 ≤ T ≤ 440 K. Structure factors, radial distribution functions, diffusion coefficient, and electrical conductivity are computed and discussed, highlighting the inherent structural heterogeneity of these compounds. The core investigation, however, is the characterization of connected paths consisting of cation chains that could support proton transport via a Grotthuss-type mechanism. Since simulations are carried out using a force field of fixed bonding topology, this analysis is based on the equilibrium structure only, using geometrical criteria to identify potential paths for proton conduction. Paths of connected cations can reach a length of 80 cations and 30 Å, provided that bridging oxygen atoms from triflate anions are taken into account. The effects of water contamination at 1% weight concentration on the structure, dynamics, and paths for proton transport are discussed.

Local structure and hydrogen bond characteristics of imidazole molecules for proton conduction in acid and base proton-conducting composite materials

Protons in composite materials of acidic polymers and imidazole molecules transport with rotational motion of imidazole in hydrogen bonds.

Potential energy construction in the diabatic picture for quantum mechanical rate constants of intermolecular proton transfer

We propose a simple method for potential construction in the diabatic picture and the estimation of thermal rate constants for intermolecular proton transfer reactions using quantum dynamics simulations carried out on the constructed potentials.

Ab initio and density functional theory (DFT) studies on triflic acid with water and protonated water clusters

The structure, stability and infrared spectral signatures of triflic acid (TA) with water clusters (W_n) and protonated water clusters (TAH⁺W_n, n = 1 - 6) were computed using DFT and MP2 methods. Our calculations show that a minimum of three water molecules are necessary to stabilize the dissociated zwitterionic form of TA. It can be seen from the results that there is no significant movement of protons in smaller (n = 1 and 2) and linear (n = 1 - 6) types of water clusters. Further, the geometries of TAW_n clusters first form a neutral pair (NP) to contact ion pair (CIP), then form a solvent separated ion pair (SSIP) in a water hexamer. These findings reveal that proton transfer may take place through NP to CIP and then CIP to SSIP. The calculated binding energies (BEs) of ion pair clusters is always higher than that of NP clusters (i.e., more stable than the NP). Existing excess proton linear chain clusters transfer a proton to adjacent water molecules via a Grotthuss mechanism, whereas the same isomers in the branched motifs do not conduct protons. Examination of geometrical parameters and infrared frequencies reveals hydronium ion (H₃O⁺ also called Eigen cation) formation in both TAW_n and protonated TAW_n clusters. The stability of Eigen water clusters is three times higher than that of other non-Eigen water clusters. Our study shows clearly that formation of ion pairs in TAW_n and TAH⁺W_n clusters greatly favors proton transfer to neighboring water molecules and also enhances the stability of these complexes.

Dynamics of proton transfer in imidazole hydrogen-bond chains

Potential energy curve of H-bond (2) in G4-[1]^(ω1=−92) is nearly identical to H-bond (1) in G2-[1]^(ω1=−94), the protonated H-bond (2) becomes a new precursor for the next transfer and confirms the Eigen–Zundel–Eigen mechanism.

Efficient Multistate Reactive Molecular Dynamics Approach Based on Short-Range Effective Potentials

Nonbonded interactions between molecules usually include the van der Waals force and computationally expensive long-range electrostatic interactions. This article develops a more efficient approach: the effective-interaction multistate empirical-valence-bond (EI-MS-EVB) model. The EI-MS-EVB method relies on a mapping of all interactions onto a short-range and thus, computationally efficient effective potential. The effective potential is tabulated by matching its force to known trajectories obtained from the full-potential empirical multistate empirical-valence-bond (MS-EVB) model. The effective pairwise interaction depends on and is uniquely determined by the atomic configuration of the system, varying only with respect to the hydrogen-bonding topology. By comparing the EI-MS-EVB and full MS-EVB calculations of several equilibrium and dynamic properties important to hydrated excess proton solvation and transport, we show that the EI-MS-EVB model produces very accurate results for the specific system in which the tabulated potentials were generated. The EI-MS-EVB potential also transfers reasonably well to similar systems with different temperatures and box sizes. The EI-MS-EVB method also reduces the computational cost of the nonbonded interactions by about 1 order of magnitude in comparison with the full algorithm.

Proton Transport in Imidazoles: Unraveling the Role of Supramolecular Structure

The impact of supramolecular hydrogen bonded networks on dynamics and charge transport in 2-ethyl-4-methylimidazole (2E4MIm), a model proton-conducting system, is investigated by broadband dielectric spectroscopy, depolarized dynamic light scattering, viscometry, and calorimetry. It is observed that the slow, Debye-like relaxation reflecting the supramolecular structure in neat 2E4MIm is eliminated upon the addition of minute amounts of levulinic acid. This is attributed to the dissociation of imidazole molecules and the breaking down of hydrogen-bonded chains, which leads to a 10-fold enhancement of ionic conductivity.

Proton transport under external applied voltage

Proton transport through an electrolyte layer between platinum electrodes under a range of applied voltages is explored using reactive molecular dynamics simulation. The proton transport process is decomposed into vehicular and Grotthuss hopping components, and the two mechanisms and their correlation are investigated as a function of applied voltage. At higher applied voltages, the effect of the hopping mechanism is much larger as compared with the vehicular mechanism. As the voltage is increased, the net correlation between the two mechanisms goes from negative to positive, and both the hopping frequency as well as the number of consecutive forward hops increases. This behavior results in a larger total diffusion constant at higher values of the voltage. The behavior of the hydrated excess proton is therefore substantially different under an applied external voltage than in the normal bulk water environment.

Proton dissociation and transfer in a phosphoric acid doped imidazole system

Fluctuations of local-dielectric environment and H-bond chain lengths lead to intermediate complexes and proton transfer along the Im H-bond chains.

Ion transport in a class of imidazole-based liquid/solid protic ionics

A class of protic ionic-compounds were prepared by Brønsted acid-base reaction of imidazole or benzimidazole with one of the following acids: trifluoromethanesulfonic, nonafluorobutanesulfonic, para-toluenesulfonic and trifluoroacetic. Except those based on trifluoroacetic acid, all prepared compounds are thermally stable up to at least 270 °C. They are solid up to temperatures between 134 and 220 °C, depending on their constituent acid and base. A simple physico-mathematical model of ion motion in the lattice was developed and implemented to correctly interpret frequency-dependent electrical response of these materials, particularly in the solid state, and investigate their ion-conducting behavior as a function of temperature. These ionic compounds display sensible ionic conductivity up to ca. 5 × 10(-4) and 5 × 10(-2) S cm(-1) in the solid and molten state, respectively, under fully anhydrous conditions. The presence of absorbed water, after brief exposure to an ambient atmosphere, enhances conduction properties remarkably. Conductivity values up to 10(-3) and 10(-1) S cm(-1) were registered, respectively, in the solid and molten state, after short exposure to (humid) ambient air. It is argued how absorbed water molecules may remove protons from (ImH)(+) or (BImH)(+) groups, thereby enabling a chain mechanism of proton-hopping through non-protonated Im or BIm sites. It is discussed how these results and methods may inspire designing protic ionic-materials at the solid-state, with enhanced proton conduction even under fully-anhydrous conditions.

Proton conduction in discotic mesogens

We show that discotic mesogens can be used to lower the activation energy barrier for proton transport.