Search

Neat liquid phosphoric acid (H(3)PO(4)) has the highest intrinsic proton conductivity of any known substance and is a useful model for understanding proton transport in other phosphate-based systems in biology and clean energy technologies. Here, we present an ab initio molecular dynamics study that reveals, for the first time, the microscopic mechanism of this high proton conductivity. Anomalously fast proton transport in hydrogen-bonded systems involves a structural diffusion mechanism in which intramolecular proton transfer is driven by specific hydrogen bond rearrangements in the surrounding environment. Aqueous media transport excess charge defects through local hydrogen bond rearrangements that drive individual proton transfer reactions. In contrast, strong, polarizable hydrogen bonds in phosphoric acid produce coupled proton motion and a pronounced protic dielectric response of the medium, leading to the formation of extended, polarized hydrogen-bonded chains. The interplay between these chains and a frustrated hydrogen-bond network gives rise to the high proton conductivity.

Trifluoromethanesulfonic acid hydrates provide a well-defined system to study proton dissociation and transport in perfluorosulfonic acid membranes, typically used as the electrolyte in hydrogen fuel cells, in the limit of minimal water. The triflic acid pentahydrate crystal (CF(3)SO(3)H·5H(2)O) is sufficiently aqueous that it contains an extended three-dimensional water network. Despite it being extended, however, long-range proton transport along the network is structurally unfavorable and would require considerable rearrangement. Nevertheless, the triflic acid pentahydrate crystal system can provide a clear picture of the preferred locations of local protonic defects in the water network, which provides insights about related structures in the disordered, low-hydration environment of perfluorosulfonic acid membranes. Ab initio molecular dynamics simulations reveal that the proton defect is most likely to transfer to the closest water that has the expected presolvation and only contains water in its first solvation shell. Unlike the tetrahydrate of triflic acid (CF(3)SO(3)H·4H(2)O), there is no evidence of the proton preferentially transferring to a water molecule bridging two of the sulfonate groups. However, this could be an artifact of the crystal structure since the only such water molecule is separated from the proton by long O-O distances. Hydrogen bonding criteria, using the two-dimensional potential of mean force, are extracted. Radial distribution functions, free energy profiles, radii of gyration, and the root-mean-square displacement computed from ab initio path integral molecular dynamics simulations reveal that quantum effects do significantly extend the size of the protonic defect and increase the frequency of proton transfer events by nearly 15%. The calculated IR spectra confirm that the dominant protonic defect mostly exists as an Eigen cation but contains some Zundel ion characteristics. Chain lengths and ring sizes determined from the hydrogen bond network, counted using graph theory techniques, are only moderately sensitive to quantum effects. Deliberately introducing a structural defect into the native crystal yields a protonic defect with one hydrogen bond to a sulfonate group that was found to be metastable for at least 10 ps., (© 2011 American Chemical Society)

Extraction of relaxation times, lifetimes, and rates associated with the transport of topological charge defects in hydrogen-bonded networks from molecular dynamics simulations is a challenge because proton transfer reactions continually change the identity of the defect core. In this paper, we present a statistical mechanical theory that allows these quantities to be computed in an unbiased manner. The theory employs a set of suitably defined indicator or population functions for locating a defect structure and their associated correlation functions. These functions are then used to develop a chemical master equation framework from which the rates and lifetimes can be determined. Furthermore, we develop an integral equation formalism for connecting various types of population correlation functions and derive an iterative solution to the equation, which is given a graphical interpretation. The chemical master equation framework is applied to the problems of both hydronium and hydroxide transport in bulk water. For each case it is shown that the theory establishes direct links between the defect's dominant solvation structures, the kinetics of charge transfer, and the mechanism of structural diffusion. A detailed analysis is presented for aqueous hydroxide, examining both reorientational time scales and relaxation of the rotational anisotropy, which is correlated with recent experimental results for these quantities. Finally, for OH(-)(aq) it is demonstrated that the "dynamical hypercoordination mechanism" is consistent with available experimental data while other mechanistic proposals are shown to fail. As a means of going beyond the linear rate theory valid from short up to intermediate time scales, a fractional kinetic model is introduced in the Appendix in order to describe the nonexponential long-time behavior of time-correlation functions. Within the mathematical framework of fractional calculus the power law decay ∼t(-σ), where σ is a parameter of the model and depends on the dimensionality of the system, is obtained from Mittag-Leffler functions due to their long-time asymptotics, whereas (stretched) exponential behavior is found for short times.

Concerted proton-electron transfer (CPET) is documented for the homogeneous reduction of O(2) to HO(2)(*) in water by the one-electron-reduced heteropolytungstate anion, alpha-PW(12)O(40)(4-) (1(1e)). At 0.01-0.3 M H(+), O(2) reduction occurs via outer-sphere electron transfer followed by proton transfer (ETPT, with rate constant k(ET)). Between 0.30 and 1.9 M H(+), rates increase linearly with [H(+)] due to a parallel CPET pathway in which H(+) is now a reactant: (1/2)k(obs) = k(ET) + k(CPET)[H(+)] (k(ET) = 1.2 M(-1) s(-1); k(CPET) = 0.8 M(-2) s(-1)). Control experiments rule out preassociation between H(+) and 1(1e). Analysis of plausible rate expressions shows that the first-order dependence on [H(+)] is uniquely consistent with multisite CPET, and a deuterium kinetic isotope effect of 1.7 is observed. Reductions of O(2) by alpha-SiW(12)O(40)(5-) confirm theoretical predictions that CPET decreases in significance as ET becomes less endergonic. Marcus analysis, including the temperature dependence of DeltaG(o), gives reorganization energies, lambda(ET) = 41.5 kcal mol(-1) and lambda(CPET) = 52.4 kcal mol(-1). At 1.5 M H(+), approximately 75% of the (1(1e),O(2)) encounter pairs form within 6 A of H(+) ions. This value (6 +/- 1 A) is the "reaction distance" for proton diffusion and probably close to that for CPET. Even so, the 70-200 ps lifetimes of the (1(1e),O(2)) pairs provide additional time for H(+) to diffuse closer to O(2). CPET is first-order in [H(+)] because k(e) for "cage escape" from (1(1e),O(2)) pairs is much larger than k(CPET), such that the rate expression for CPET becomes -(1/2)d[1(1e)]/dt = (k(d)/k(e))k(CPET)[1(1e)][O(2)][H(+)], where k(d) is the rate constant for (1(1e),O(2)) pair formation. Overall, the findings suggest that the emergence of CPET, with hydronium ion as the proton donor, may prove a general feature of sufficiently endergonic reductions of dioxygen by otherwise "outer-sphere" complexes (or electrode reactions) at sufficiently low pH values in water.

Proton dissociation and transfer were examined with ab initio molecular dynamics (AIMD) simulations of carbon nanotubes (CNT) functionalized with perfluorosulfonic acid (-CF(2)SO(3)H) groups with 1-3 H(2)O/SO(3)H. The CNT systems were constructed both with and without fluorine atoms covalently bound to the walls to elucidate the effects of the presence of a strongly hydrophobic environment, the fluorine, on proton dissociation, hydration, and stabilization. The simulations revealed that the dissociated proton was preferentially stabilized as a hydrated hydronium cation (i.e., Eigen like) in the fluorinated CNTs but as a Zundel (H(5)O(2)(+)) cation in the nonfluorinated CNTs. This feature is attributed to the fluorine atoms forming hydrogen bonds with the water molecules coordinated to the central hydronium ion.

First-principles molecular dynamics calculations performed in a fully converged basis set are used to reveal new details about the mechanism of the anomalous proton-transport process in water, a fundamental question dating back over 200 years. By separating actual structural diffusion from simple rattling events, wherein a proton shuttles forth and back in a hydrogen bond, it is found that the former are driven by a concerted mechanism in which hydronium begins to accept a hydrogen bond from a donor water molecule while the proton-receiving water molecule simultaneously loses one of its acceptor hydrogen bonds. The kinetics of the process are found to be in good agreement with recent experiments.

A multistate empirical valence bond model for proton transport in water, which explicitly includes solvent polarization, is presented. Polarization is included for each valence-bond state via induced point dipoles, and the model is parametrized to be used with an effective path integral derived potential surface, so as to include quantum effects of the transferring proton. The new model is shown to reproduce ab initio geometries and energetics for small protonated clusters. It is also shown that the new model gives a diffusion constant for the excess proton in water, which is in good agreement with experiment, and that the qualitative features of ab initio path integral simulations [D. Marx, M. E. Tuckerman, J. Hutter, and M. Parrinello, Nature (London) 397, 601 (1999)] are well reproduced. [ABSTRACT FROM AUTHOR]

The transport of protons through aqueous, partially aqueous, or nonaqueous hydrogen-bonded media is a fundamental process in many biologically and technologically important systems. Liquid methanol is an example of a hydrogen-bonded system that, like water, supports anomalously fast proton transport. Using the methodology of ab initio molecular dynamics, in which internuclear forces are computed directly from electronic structure calculations as the simulation proceeds, we have investigated the microscopic mechanism of the proton transport process in liquid methanol at 300 K. It is found that the defect structure associated with an excess proton in liquid methanol is a hydrogen-bonded cationic chain whose length generally exceeds the average chain length in pure liquid methanol. Hydrogen bonds in the first and second solvation shells of the excess proton are considerably shorter and stronger than ordinary methanol-methanol hydrogen bonds. Along this chain, proton transfer reactions occur in an essentially random manner described by Poisson statistics. Structural diffusion of the defect structure is possible if the proton migrates toward an end of the defect chain, which causes a weakening of the hydrogen bonds at the opposite end. The latter can, therefore, be easily ruptured by ordinary thermal fluctuations. At the end of the chain where the proton resides, new hydrogen bonds are likely to form due to the strong associative nature of the excess proton. It is through this "snake-like" mechanism that the defect structure is able to diffuse through the hydrogen-bond network of the liquid. The estimated activation enthalpy of this proposed mechanism is found to be in reasonable agreement with the experimentally determined activation enthalpy. [ABSTRACT FROM AUTHOR]

Neat liquid phosphoric acid (H3PO4) has the highest intrinsic proton conductivity of any known substance and is a useful model for understanding proton transport in other phosphate-based systems in biology and clean energy technologies. Here, we present an ab initio molecular dynamics study that reveals, for the first time, the microscopic mechanism of this high proton conductivity. Anomalously fast proton transport in hydrogen-bonded systems involves a structural diffusion mechanism in which intramolecular proton transfer is driven by specific hydrogen bond rearrangements in the surrounding environment. Aqueous media transport excess charge defects through local hydrogen bond rearrangements that drive individual proton transfer reactions. In contrast, strong, polarizable hydrogen bonds in phosphoric acid produce coupled proton motion and a pronounced protic dielectric response of the medium, leading to the formation of extended, polarized hydrogen-bonded chains. The interplay between these chains and a frustrated hydrogen-bond network gives rise to the high proton conductivity. [ABSTRACT FROM AUTHOR]

Examines the solution structures and transport mechanisms of hydrated hydroxide. Variation in the transport mechanism for protons and hydroxide ions; Mechanism of transport of hydrated proton in water; Assessment of the role of solvation complexes of hydroxide ions.

Addresses questions concerning the character and anomalously high mobility of the hydrated excess proton in liquid water using computer simulation techniques. Widely held views of the proton's structure and transport; Conclusion that the proton's quantum-mechanical character plays a more important role than previously thought; Implications of results in areas related to aqueous acid-base chemistry, enzymatic proton transfers, proton transport in biological channels and atmospheric proton transfers on ice surfaces.