Proton transfer in imidazole-based molecular crystals

Proton Conduction in Chiral Molecular Assemblies of Azolium-Camphorsulfonate Salts

Chiral molecular assemblies have attracted considerable attention because of their interesting physical properties, such as spin-selective electron transport. Cation-anion salts of three azolium cations, imidazolium (HIm+), triazolium (HTrz+), and thiazolium (HThz+), in combination with a chiral camphorsulfonate (1S-CS-) and their racemic compounds (rac-CS-) were prepared and compared in terms of phase transitions, crystal structures, dynamics of constituent molecules, dielectric responses, and proton conductivities. The cation-anion crystals containing HIm+ showed no significant difference in proton conductivity between the homochiral and racemic crystals, whereas the HTrz+-containing crystals showed higher proton conductivity and lower activation energy in the homochiral form than in the racemic form. A two-dimensional hydrogen-bonding network consisting of HTrz+ and -SO3- groups and similar in-plane rotational motion was observed in both crystals; however, the HTrz+ cation in the homochiral crystal exhibited the rotational motion modulated with translational motion, whereas the HTrz+ cation in the racemic crystal exhibited almost steady in-plane rotational motion. The different motional degrees of freedom were confirmed by crystal structure analyses and temperature- and frequency-dependent dielectric constants. In contrast, steady in-plane rotational motion with the thermally activated fluctuating motion of CS- was observed both in homochiral and racemic crystals containing HIm+, which averaged the motional space of protons resulting in similar dielectric responses and proton conductivities. The control of motional degrees of freedom in homochiral crystals affects the proton conductivity and is useful for the design of molecular proton conductors.

Mutual influence of non-covalent interactions formed by imidazole: A systematic quantum-chemical study

Imidazole is a five-membered heterocycle that is part of a number of biologically important molecules such as the amino acid histidine and the hormone histamine. Imidazole has a unique ability to participate in a variety of non-covalent interactions involving the NH group, the pyridine-like nitrogen atom or the π-system. For many biologically active compounds containing the imidazole moiety, its participation in formation of hydrogen bond NH⋯O/N and following proton transfer is the key step of mechanism of their action. In this work a systematic study of the mutual influence of various paired combinations of non-covalent interactions (e.g., hydrogen bonds and π-interactions) involving the imidazole moiety was performed by means of quantum chemistry (PW6B95-GD3/def2-QZVPD) for a series of model systems constructed based on analysis of available x-ray data. It is shown that for considered complexes formation of additional non-covalent interactions can only enhance the proton-donating ability of imidazole. At the same time, its proton-accepting ability can be both enhanced and weakened, depending on what additional interactions are added to a given system. The mutual influence of non-covalent interactions involving imidazole can be classified as weak geometric and strong energetic cooperativity-a small change in the length of non-covalent interaction formed by imidazole can strongly influence its strength. The latter can be used to develop methods for controlling the rate and selectivity of chemical reactions involving the imidazole fragment in larger systems. It is shown that the strong mutual influence of non-covalent interactions involving imidazole is due to the unique ability of the imidazole ring to effectively redistribute electron density in non-covalently bound systems with its participation.

Structures of the (Imidazole)nH+ ... Ar (n=1,2,3) complexes determined from IR spectroscopy and quantum chemical calculations

Here, we present new cryogenic infrared spectra of the (Imidazole)$$_{n}\mathrm{H}^{+}$$ n H + (n=1,2,3) ions. The data was obtained using helium tagging infrared predissociation spectroscopy. The new results were compared with the data obtained by Gerardi et al. (Chem. Phys. Lett. 501:172–178, 2011) using the same technique but with argon as a tag. Comparison of the two experiments, assisted by theoretical calculations, allowed us to evaluate the preferable attachment positions of argon to the (Imidazole)$$_{n}\mathrm{H}^{+}$$ n H + frame. Argon attaches to nitrogen-bonded hydrogen in the case of the (Imidazole)H$$^+$$ + ion, while in (Imidazole)$$_{2}\mathrm{H}^{+}$$ 2 H + and (Imidazole)$$_{3}\mathrm{H}^{+}$$ 3 H + the preferred docking sites for the argon are in the center of the complex. This conclusion is supported by analyzing the spectral features attributed to the N–H stretching vibrations. Symmetry adapted perturbation theory (SAPT) analysis of the non-covalent forces between argon and the (Imidazole)$$_{n}\mathrm{H}^{+}$$ n H + (n=1,2,3) frame revealed that this switch of docking preference with increasing complex size is caused by an interplay between induction and dispersion interactions.

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.

A dynamic and electrostatic potential prediction of the prototropic tautomerism between imidazole 3-oxide and 1-hydroxyimidazole in external electric field

In order to obtain an optimum scheme for separating the proton-transfer tautomer, a dynamic investigation into the effect of the external electric field on the proton-transfer tautomeric conversion in imidazole 3-oxide and 1-hydroxyimidazole was carried out at the M06-2X/6-311++G** and CCSD(T)/6-311++G(2d,p) level, accompanied by the analysis of the surface electrostatic potentials. The results show that, for both the forward reaction "imidazole 3-oxide → N-hydroxyimidazole free radical → 1-hydroxyimidazole" and its reverse reaction processes, the fields parallel to the N→O or N-OH bond axis affect the barrier heights and rate constants considerably more than those parallel to the other orientations. As the field strength is increased along the orientation from the O to N atom, the chemical equilibrium moves toward the direction for the formation of 1-hydroxyimidazole, while the amount of imidazole 3-oxide is increased with the increased field strength along the opposite orientation. In the fields along the orientation consistent with the dipole moment, the electrostatic potentials and their variances "abnormally" increase for the transition states with the N→O bond in comparison with those in no field (they decrease generally), which enhances the nucleophilicity of the coordination O atom and the electrophilicity of the activated H atom. The analyses of the AIM (atoms in molecules) and NICS (nucleus-independent chemical shift) were used to explain the above anomaly. Graphical Abstract Electrostatic potentials and their variances "abnormally" increase in the external electric field, which greatly affects tautomeric conversion.

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.

Liquid Water Confined in Cellulose with Variable Interfacial Hydrophilicity

We investigate liquid water confined within nanoscale cellulose slabs by means of molecular dynamics simulations. Depending on the construction of the cellulose–water interface, two different surface structures with distinct levels of hydrophilicity are exposed to the water. The different philicities are reflected in the response of the water phase to this geometric confinement, both in terms of the density profile and in the strength of the aqueous hydrogen bonding network. At the smooth surface cut along the (010) axis of the cellulose crystal, water shows typical properties of a hydrophilic confinement: the density shows fluctuations that disappear further away from the wall, the water molecules orient themselves and the coordination numbers increases at the interface. As a consequence, the water becomes “harder” at the interface, with a considerably increased local ordering. At the zigzag-shaped surface along the (111) axis, the degree of hydrophilicity is reduced, and only small effects can be seen: the density shows weak fluctuations, and the orientation of the water molecules is closer to that of bulk water than to the smooth surface. The local coordination numbers remains constant over the whole confinement. Our work shows that the nature of the exposed cellulose interface has a strong influence on how the structure of adjacent water is modified. The different ways of surface construction yield distinct degrees of hydrophilicity and spatial accessibility regarding the hydrogen bond network, resulting in a notably different interfacial water structure.

Ab initio molecular dynamics study of H-bonding and proton transfer in the phosphoric acid-N,N-Dimethylformamide system

Car-Parrinello molecular dynamics simulations of phosphoric acid (H₃PO₄)-N,N-dimethylformamide (DMF) mixtures over the whole composition range have been carried out. It has been found that the neutral molecules are the dominant species in this system. The concentration dependences of the average number of H-bonds per proton acceptor atom in P=O and C=O groups as well as per proton donor atom in DMFH⁺ ions towards phosphate species have been discussed. The H-bonding between components in all investigated mixtures of H₃PO₄ and DMF is possible. A significant fraction of the protonated DMF forms appears at phosphoric acid mole fraction higher than 0.37, indicating a high probability of proton transfer from phosphate species to oxygen atoms in C=O groups. The intermolecular proton transfer between phosphate species themselves is mainly observed when x_H3PO4 > 0.19. Satisfactory agreement with available experimental data for structural characteristics of the investigated system was obtained.

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.

Maleimide: a potential building block for the design of proton exchange membranes studied by ab initio molecular dynamics simulations

Ab initio molecular dynamics (AIMD) simulations are applied to the study of proton transport in solid state maleimide.

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.

Proton conduction in discotic mesogens

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