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

Exploiting Cation Structure and Water Content in Modulating the Acidity of Ammonium Hydrogen Sulfate Protic Ionic Liquids

In this paper, we investigated the effect of cation structure and water content on proton dissociation in alkylammonium [HSO4]- protic ionic liquids (ILs) doped with 20 wt % water and correlated this with experimental Hammett acidities. For pure systems, increased cation substitution resulted in a reduction in the number of direct anion-anion neighbors leading to larger numbers of small aggregates, which is further enhanced with addition of water. We also observed spontaneous proton dissociation from [HSO4]- to water only for primary amine-based protic ILs, preceded by the formation of an anion trimer motif. Investigation using DFT calculations revealed spontaneous proton dissociation from [HSO4]- to water can occur for each of the protic ILs investigated; however, this is dependent on the size of the anion aggregates. These findings are important in the fields of catalysis and lignocellulosic biomass, where solvent acidity is a crucial parameter in biomass fractionation and lignin chemistry.

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.

On the Rich Chemistry of Pseudo-Protic Ionic Liquid Electrolytes

Mixing weak acids and bases can produce highly complicated binary mixtures, called pseudo-protic ionic liquids, in which a complex network of effects determines the physicochemical properties that are currently impossible to predict. In this joint computational-experimental study, we investigated 1-methylimidazole-acetic acid mixtures through the whole concentration range. Effects of the varying ionization and excess of either components on the properties, such as density, diffusion coefficients, and overall hydrogen bonding structure were uncovered. A special emphasis was put on understanding the multiple factors that govern the conductivity of the system. In the presence of an excess of acetic acid, the 1-methylimidazolium acetate ion pairs dissociate more efficiently, resulting in a higher concentration of independently moving, conducting ions. However, the conductivity measurements showed that higher concentrations of acetic acid improve the conductivity beyond this effect, suggesting in addition to standard dilution effects the occurrence of Grotthuss diffusion in high acid-to-base ratios. The results here will potentially help designing novel electrolytes and proton conducting systems, which can be exploited in a variety of applications.

Ultrafast Proton Transfer Pathways Mediated by Amphoteric Imidazole

Imidazole, being an amphoteric molecule, can act both as an acid and as a base. This property enables imidazole, as an essential building block, to effectively facilitate proton transport in high-temperature proton exchange membrane fuel cells and in proton channel transmembrane proteins, enabling those systems to exhibit high energy conversion yields and optimal biological function. We explore the amphoteric properties of imidazole by following the proton transfer exchange reaction dynamics with the bifunctional photoacid 7-hydroxyquinoline (7HQ). We show with ultrafast ultraviolet-mid-infrared pump-probe spectroscopy how for imidazole, in contrast to expectations based on textbook knowledge of acid-base reactivity, the preferential reaction pathway is that of an initial proton transfer from 7HQ to imidazole, and only at a later stage a transfer from imidazole to 7HQ, completing the 7HQ tautomerization reaction. An assessment of the molecular distribution functions and first-principles calculations of proton transfer reaction barriers reveal the underlying reasons for our observations.

Search for a Grotthuss mechanism through the observation of proton transfer

The transport of protons is critical in a variety of bio- and electro-chemical processes and technologies. The Grotthuss mechanism is considered to be the most efficient proton transport mechanism, generally implying a transfer of protons between 'chains' of host molecules via elementary reactions within the hydrogen bonds. Although Grotthuss proposed this concept more than 200 years ago, only indirect experimental evidence of the mechanism has been observed. Here we report the first experimental observation of proton transfer between the molecules in pure and 85% aqueous phosphoric acid. Employing dielectric spectroscopy, quasielastic neutron, and light scattering, and ab initio molecular dynamic simulations we determined that protons move by surprisingly short jumps of only ~0.5-0.7 Å, much smaller than the typical ion jump length in ionic liquids. Our analysis confirms the existence of correlations in these proton jumps. However, these correlations actually reduce the conductivity, in contrast to a desirable enhancement, as is usually assumed by a Grotthuss mechanism. Furthermore, our analysis suggests that the expected Grotthuss-like enhancement of conductivity cannot be realized in bulk liquids where ionic correlations always decrease conductivity.

Dynamically Interfacial pH-Buffering Effect Enabled by N-Methylimidazole Molecules as Spontaneous Proton Pumps toward Highly Reversible Zinc-Metal Anodes

Aqueous zinc-metal batteries have attracted extensive attention due to their outstanding merits of high safety and low cost. However, the intrinsic thermodynamic instability of zinc in aqueous electrolyte inevitably results in hydrogen evolution, and the consequent generation of OH^- at the interface will dramatically exacerbate the formation of dead zinc and dendrites. Herein, a dynamically interfacial pH-buffering strategy implemented by N-methylimidazole (NMI) additive is proposed to remove the detrimental OH^- at zinc/electrolyte interface in real-time, thus eliminating the accumulation of by-products fundamentally. Electrochemical quartz crystal microbalance and molecular dynamics simulation results reveal the existence of an interfacial absorption layer assembled by NMI and protonated NMI (NMIH⁺ ), which acts as an ion pump for replenishing the interface with protons constantly. Moreover, an in situ interfacial pH detection method with micro-sized spatial resolution based on the ultra-microelectrode technology is developed to probe the pH evolution in diffusion layer, confirming the stabilized interfacial chemical environment in NMI-containing electrolyte. Accordingly, with the existence of NMI, an excellent cumulative plating capacity of 4.2 Ah cm^-2 and ultrahigh Coulombic efficiency of 99.74% are realized for zinc electrodes. Meanwhile, the NMI/NMIH⁺ buffer additive can accelerate the dissolution/deposition process of MnO₂ /Mn²⁺ on the cathode, leading to enhanced cycling capacity.

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.

Understanding Proton Transfer in Non-aqueous Biopolymers based on Helical Peptides: A Quantum Mechanical Study

Histidine (an imidazole-based amino acid) is a promising building block for short aromatic peptides containing a proton donor/acceptor moiety. Previous studies have shown that polyalanine helical peptides substituted at regular intervals with histidine residues exhibit both structural stability as well as high proton affinity and high conductivity. Here, we present first-principle calculations of non-aqueous histidine-containing 310-, α- and π-helices and show that they are able to form hydrogen-bonded networks mimicking proton wires that have the ability to shuttle protons via the Grotthuss shuttling mechanism. The formation of these wires enhances the stability of the helices, and our structural characterizations confirm that the secondary structures are conserved despite distortions of the backbones. In all cases, the helices exhibit high proton affinity and proton transfer barriers on the order of 1~4 kcal/mol. Zero-point energy calculations suggest that for these systems, ground state vibrational energy can provide enough energy to cross the proton transport energy barrier. Additionally, ab initio molecular dynamics results suggests that the protons are transported unidirectionally through the wire at a rate of approximately 2 Å every 20 fs. These results demonstrate that efficient deprotonation-controlled proton wires can be formed using non-aqueous histidine-containing helical peptides.

Monitoring the photochemistry of a formazan over 15 orders of magnitude in time

2,3,5-triphenyltetrazolium chloride (TTC) may convert into phenyl-benzo[c]tetrazolocinnolium chloride (PTC) and 1,3,5-triphenylformazan (TPF) under irradiation with light. The latter reaction, albeit enzymatically rather than photochemically, is used in so-called TTC assays indicating cellular respiration and cell growth. In this paper, we address the photochemistry of TPF with time-resolved spectroscopy on various time scales. TPF is stabilized by an intramolecular hydrogen bond and switches photochemically via an E-Z isomerization around an N=N double bond into another TPF-stereoisomer, from which further isomerizations around the C=N double bond of the phenylhydrazone group are possible. We investigate the underlying processes by time-resolved spectroscopy in dependence on excitation wavelength and solvent environment, resolving several intermediates over a temporal range spanning 15 orders of magnitude (hundreds of femtoseconds to hundreds of seconds) along the reaction path. In a quantum-chemical analysis, we identify 16 stable ground-state isomers and discuss which ones are identified in the experimental data. We derive a detailed scheme how these species are thermally and photochemically interconnected and conclude that proton transfer processes are involved.

Evolution of microscopic heterogeneity and dynamics in choline chloride-based deep eutectic solvents

Deep eutectic solvents (DESs) are an emerging class of non-aqueous solvents that are potentially scalable, easy to prepare and functionalize for many applications ranging from biomass processing to energy storage technologies. Predictive understanding of the fundamental correlations between local structure and macroscopic properties is needed to exploit the large design space and tunability of DESs for specific applications. Here, we employ a range of computational and experimental techniques that span length-scales from molecular to macroscopic and timescales from picoseconds to seconds to study the evolution of structure and dynamics in model DESs, namely Glyceline and Ethaline, starting from the parent compounds. We show that systematic addition of choline chloride leads to microscopic heterogeneities that alter the primary structural relaxation in glycerol and ethylene glycol and result in new dynamic modes that are strongly correlated to the macroscopic properties of the DES formed.

Tailoring the macroscopic properties of deep eutectic solvents requires knowing how these depend on the local structure and microscopic dynamics. The authors, with computational and experimental tools spanning a wide range of space- and timescales, shed light into the relationship between micro and macroscopic properties in glyceline and ethaline.

Quantum mechanical effects in acid–base chemistry

Acid–base chemistry has immense importance for explaining and predicting the chemical products formed by an acid and a base when mixed together. However, the traditional chemistry theories used to describe acid–base reactions do not take into account the effect arising from the quantum mechanical nature of the acidic hydrogen shuttling potential and its dependence on the acid base distance. Here, infrared and NMR spectroscopies, in combination with first principles simulations, are performed to demonstrate that quantum mechanical effects, including electronic and nuclear quantum effects, play an essential role in defining the acid–base chemistry when 1-methylimidazole and acetic acid are mixed together. In particular, it is observed that the acid and the base interact to form a complex containing a strong hydrogen bond, in which the acidic hydrogen atom is neither close to the acid nor to the base, but delocalized between them. In addition, the delocalization of the acidic hydrogen atom in the complex leads to characteristic IR and NMR signatures. The presence of a hydrogen delocalized state in this simple system challenges the conventional knowledge of acid–base chemistry and opens up new avenues for designing materials in which specific properties produced by the hydrogen delocalized state can be harvested.

Acid-based theories do not consider the quantum mechanical nature of the acidic hydrogen shuttling potential. Here, it is demonstrated that this particularity is needed to explain the formation acid-base complex with a delocalized acidic hydrogen.

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.

p-Aminobenzoic acid protonation dynamics in an evaporating droplet by ab initio molecular dynamics

Protonation equilibria are known to vary from the bulk to microdroplet conditions, which could induce many chemical and physical phenomena. Protonated p-aminobenzoic acid (PABA + H⁺) can be considered a model for probing the protonation dynamics in an evaporating droplet, as its protonation equilibrium is highly dependent on the formation conditions from solution via atmospheric pressure ionization sources. Experiments using diverse experimental techniques have shown that protic solvents allow formation of the O-protomer (PABA protonated in the carboxylic acid group) stable in the gas phase, while aprotic solvents yield the N-protomer (protonated in the amino group) that is the most stable protomer in solution. In this work, we explore the protonation equilibrium of PABA solvated by different numbers of water molecules (n = 0 to 32) using ab initio molecular dynamics. For n = 8-32, the protonation is either at the NH₂ group or in the solvent network. The solvent network interacts with the carboxylic acid group, but there is no complete proton transfer to form the O-protomer. For smaller clusters, however, solvent-mediated proton transfers to the carboxylic acid were observed, both via the Grotthuss mechanism and the vehicle or shuttle mechanism (for n = 1 and 2). Thermodynamic considerations allowed a description of the origins of the kinetic trapping effect, which explains the observation of the solution structure in the gas phase. This effect likely occurs in the final evaporation steps, which are outside the droplet size range covered by previous classical molecular dynamics simulations of charged droplets. These results may be considered relevant in determining the nature of the species observed in the ubiquitous ESI based mass spectrometry analysis, and in general for droplet chemistry, explaining how protonation equilibria are drastically changed from bulk to microdroplet conditions.

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.

Hopping in High Concentration Electrolytes - Long Time Bulk and Single-Particle Signatures, Free Energy Barriers, and Structural Insights

Although ion-hopping is believed to be a significant mode of transport for small ions in liquid high concentration electrolytes (HCE), its bulk signatures over sufficiently long time intervals are yet to be shown. We computationally establish the long and short time imprints of hopping in HCEs using LiBF₄-in-sulfolane mixtures as models. The high viscosity of this electrolyte leads to significant dynamic heterogeneity in Li-ion transport. Li-ions exhibit a preference to transit to previously occupied Li-ion-sites, bridged through anion or solvent molecules. Hopping in the liquid matrix was found to be an activated process, whose free energy barrier and transition state structure have been determined. Evidence for nanoscale compositional heterogeneity at high salt concentrations is also presented. The simulations shed light on the composition, stiffness, and lifetime of the solvation shell of Li ions. The understanding of HCEs gleaned from this study will spearhead the choice, engineering and applicability of this class of electrolytes.