An official website of the United States government
Proton transport in aqueous systems occurs by making and breaking covalent bonds, a process that classical force fields cannot reproduce. Various attempts have been made to remedy this deficiency, by valence bond theory or instantaneous proton transfers, but the ability of such methods to provide a realistic picture of this fundamental process has not been fully evaluated. Here we compare an ab initio molecular dynamics (AIMD) simulation of an excess proton in water to a simulation of a classical H3O+ in TIP3P water. The energy gap upon instantaneous proton transfer from H3O+ to an acceptor water molecule is much higher in the classical simulation than in the AIMD configurations evaluated with the same classical potential. The origins of this discrepancy are identified by comparing the solvent structures around the excess proton in the two systems. One major structural difference is in the tilt angle of the water molecules that accept an hydrogen bond from H3O+. The lack of lone pairs in TIP3P produces a tilt angle that is too large and generates an unfavorable geometry after instantaneous proton transfer. This problem can be alleviated by the use of TIP5P, which gives a tilt angle much closer to the AIMD result. Another important factor that raises the energy gap is the different optimal distance in water-water vs H3O+-water H-bonds. In AIMD the acceptor is gradually polarized and takes a hydronium-like configuration even before proton transfer actually happens. Ways to remedy some of these problems in classical simulations are discussed.
more »
« less
Structural aspects of the topological model of the hydrogen bond in water on auto-dissociation via proton transfer
Molecular dynamics (MD) simulations were used to investigate the structure and lifetimes of hydrogen bonds and auto dissociation via proton transfer in bulk water using a reactive and dissociative all-atom potential that has previously been shown to match a variety of water properties and proton transfer. Using the topological model, each molecule's donated and accepted hydrogen bonds were labeled relative to the other hydrogen bonds on neighboring waters, providing a description of the effect of these details on the structure, dynamics and autoionization of water molecules. In agreement with prior data, asymmetric bonding at the sub-100 femtosecond timescale is observed, as well as the existence of linear, bifurcated, and dangling hydrogen bonds. The lifetime of the H-bond, 2.1 ps, is consistent with experimental data, with short time librations on the order of femtoseconds. The angular correlation functions, the presence of a second shell water entering the first shell, and OH vibrational stretch frequencies were all consistent with experiment or ab initio calculations. The simulations show short-lived (femtoseconds) dissociation of a small fraction of water molecules followed by rapid recombination. The role of the other H-bonds to the acceptor and on the donor plays an important part in proton transfer between the molecules in auto dissociation and is consistent with the role of a strong electric field caused by local (first and second shell) waters on initiating dissociation. The number of H-bonds to the donor water is 4.3 per molecule in the simulations, consistent with previous data regarding the number of hydrogen bonds required to generate this strong local electric field that enhances dissociation. The continuous lifetime autocorrelation function of the H-bond for those molecules that experience dissociation is considerably longer than that for all molecules that show no proton transfer.
more »
« less
- Award ID(s):
- 1609044
- PAR ID:
- 10090791
- Date Published:
- Journal Name:
- Physical Chemistry Chemical Physics
- Volume:
- 20
- Issue:
- 24
- ISSN:
- 1463-9076
- Page Range / eLocation ID:
- 16414 to 16427
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Short hydrogen bonds (SHBs), which have donor and acceptor separations below 2.7 Å, occur extensively in small molecules and proteins. Due to their compact structures, SHBs exhibit prominent covalent characters with elongated Donor-H bonds and highly downfield (>14 ppm) 1H NMR chemical shifts. In this work, we carry out first principles simulations on a set of model molecules to assess how quantum effects determine the symmetry and chemical shift of their SHBs. From simulations that incorporate the quantum mechanical nature of both the electrons and nuclei, we reveal a universal relation between the chemical shift and the position of the proton in a SHB, and unravel the origin of the observed downfield spectral signatures. We further develop a metric that allows one to accurately and efficiently determine the proton position directly from its 1H chemical shift, which will facilitate the experimental examination of SHBs in both small molecules and biological macromolecules.more » « less
-
null (Ed.)Hydrogen bonding plays a critical role in maintaining order and structure in complex biological and synthetic systems. N -heterocyclic carbenes (NHCs) represent one of the most versatile tools in the synthetic chemistry toolbox, yet their potential as neutral carbon hydrogen bond acceptors remains underexplored. This report investigates this capability in a strategic manner, wherein carbene-based hydrogen bonding can be assessed by use of ditopic NH -containing molecules. N–H bonds are unique as there are three established reaction modes with carbenes: non-traditional hydrogen bonding adducts (X–H⋯:C), salts arising from proton transfer ([H–C] + [X] − ), or amines from insertion of the carbene into the N–H bond. Yet, there are no established rules to predict product distributions or the strength of these associations. Here we seek to correlate the hydrogen bond strength of symmetric and asymmetric ditopic secondary amines with 1,3-bis(2,6-diisopropylphenyl)imidazol-2-ylidene ( IPr , a representative NHC). In symmetric and asymmetric ditopic amine adducts both the solid-state (hydrogen bond lengths, NHC interior angles) and solution-state ( 1 H Δ δ of NH signals, 13 C signals of carbenic carbon) can be related to the p K a of the parent amine.more » « less
-
The three-dimensional architecture of biomolecules often creates specialized structural elements, notably short hydrogen bonds that have donor–acceptor separations below 2.7 Å. In this work, we statistically analyze 1663 high-resolution biomolecular structures from the Protein Data Bank and demonstrate that short hydrogen bonds are prevalent in proteins, protein–ligand complexes and nucleic acids. From these biological macromolecules, we characterize the preferred location, connectivity and amino acid composition in short hydrogen bonds and hydrogen bond networks, and assess their possible functional importance. Using electronic structure calculations, we further uncover how the interplay of the structural and chemical features determines the proton potential energy surfaces and proton sharing conditions in biological short hydrogen bonds.more » « less
-
Abstract Several renewable energy schemes aim to use the chemical bonds in abundant molecules like water and ammonia as energy reservoirs. Because the O-H and N-H bonds are quite strong (>100 kcal/mol), it is necessary to identify substances that dramatically weaken these bonds to facilitate proton-coupled electron transfer processes required for energy conversion. Usually this is accomplished through coordination-induced bond weakening by redox-active metals. However, coordination-induced bond weakening is difficult with earth’s most abundant metal, aluminum, because of its redox inertness under mild conditions. Here, we report a system that uses aluminum with a redox non-innocent ligand to achieve significant levels of coordination-induced bond weakening of O-H and N-H bonds. The multisite proton-coupled electron transfer manifold described here points to redox non-innocent ligands as a design element to open coordination-induced bond weakening chemistry to more elements in the periodic table.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1039/c8cp02592d
