Analytic functions fit to proton transfer potentials

Kinetic Monte Carlo simulations of proton conductivity

The kinetic Monte Carlo method is used to model the dynamic properties of proton diffusion in anhydrous proton conductors. The results have been discussed with reference to a two-step process called the Grotthuss mechanism. There is a widespread belief that this mechanism is responsible for fast proton mobility. We showed in detail that the relative frequency of reorientation and diffusion processes is crucial for the conductivity. Moreover, the current dependence on proton concentration has been analyzed. In order to test our microscopic model the proton transport in polymer electrolyte membranes based on benzimidazole C(7)H(6)N(2) molecules is studied.

Analysis of Complexes Pairing Hydroperoxyl Radical with Peroxyformic Acid

Ab initio quantum calculations are used to analyze the binding of complexes pairing OOH with HOOCHO. Six minima are located on the potential energy surface, all of cyclic geometry. Of particular interest are the OH...O and CH...O H-bonds that arise in the complexes and the manner in which these interactions influence the internal properties of the subunits. The analysis is complicated by the presence of an intramolecular H-bond in the unperturbed HOOCHO molecule, which must be broken in order to form the pair of intermolecular H-bonds that are responsible for the binding in the most stable complex. The CH bond of HOOCHO is contracted, and its stretching frequency undergoes a blue shift, when this group participates in a H-bond.

Energy barriers of proton transfer reactions between amino acid side chain analogs and water from ab initio calculations

Proton transfer reactions were studied in all titratable pairs of amino acid side chains where, under physiologically reasonable conditions, one amino acid may function as a donor and the other one as an acceptor. Energy barriers for shifting the proton from donor to acceptor atom were calculated by electronic structure methods at the MP2/6-31++G(d,p) level, and the well-known double-well potentials were characterized. The energy difference between both minima can be expressed by a parabola using as argument the donor-acceptor distance R(DA). In this work, the fit parameters of the quadratic expression are determined for each donor-acceptor pair. Moreover, it was found previously that the energy barriers of the reactions can be expressed by an analytical expression depending on the distance between donor and acceptor and the energy difference between donor and acceptor bound states. The validity of this approach is supported by the extensive new data set. This new parameterization of proton transfer barriers between titratable amino acid side chains allows us to very efficiently estimate proton transfer probabilities in molecular modelling studies or during classical molecular dynamics simulation of biomolecular systems.

Theoretical Evidence for a NH···XC Blue-Shifting Hydrogen Bond: Complexes Pairing Monohalomethanes with HNO

Correlated ab initio calculations are used to analyze the interaction between nitrosyl hydride (HNO) and CH3X (X = F, Cl, Br). Three minima are located on the potential energy surface of each complex. The more strongly bound contains a NH...X bond, along with CH...O; CH...O and CH...N bonds occur in the less stable minimum. Binding energies of the global minimum lie in the range of 11-13 kJ/mol, and there is little sensitivity to the identity of the halogen atom. Unlike most other such hydrogen bonds, the NH covalent bond in this set of complexes becomes shorter, and its stretching frequency shifts to the blue, upon forming the NH...X hydrogen bond. The amount of this blue shift varies in the order F > Cl > Br.

Discrete kink dynamics in hydrogen-bonded chains: the two-component model

We study discrete topological solitary waves (kinks and antikinks) in two nonlinear diatomic chain models that describe the collective dynamics of proton transfers in one-dimensional hydrogen-bonded networks. The essential ingredients of the models are (i) a realistic (anharmonic) ion-proton interaction in the hydrogen bond, (ii) a harmonic coupling between the protons in adjacent hydrogen bonds, and (iii) a harmonic coupling between the nearest-neighbor heavy ions (an isolated diatomic chain with the lowest acoustic band) or instead a harmonic on-site potential for the heavy ions (a diatomic chain subject to a substrate with two optical bands), both providing a bistability of the hydrogen-bonded proton. Exact two-component (kink and antikink) discrete solutions for these models are found numerically. We compare the soliton solutions and their properties in both the one- (when the heavy ions are fixed) and two-component models. The effect of stability switchings, discovered previously for a class of one-component kink-bearing models, is shown to exist in these two-component models as well. However, the presence of the second component, i.e., the softness of the heavy-ion sublattice, brings principal differences, like a significant difference in the stability switchings behavior for the kinks and the antikinks. Water-filled carbon nanotubes are briefly discussed as possible realistic systems, where topological discrete (anti)kink states might exist.

Discrete kink dynamics in hydrogen-bonded chains: the one-component model

We study topological solitary waves (kinks and antikinks) in a nonlinear one-dimensional Klein-Gordon chain with the on-site potential of a double-Morse type. This chain is used to describe the collective proton dynamics in quasi-one-dimensional networks of hydrogen bonds, where the on-site potential plays the role of the proton potential in the hydrogen bond. The system supports a rich variety of stationary kink solutions with different symmetry properties. We study the stability and bifurcation structure of all these stationary kink states. An exactly solvable model with a piecewise "parabola-constant" approximation of the double-Morse potential is suggested and studied analytically. The dependence of the Peierls-Nabarro potential on the system parameters is studied. Discrete traveling-wave solutions of a narrow permanent profile are shown to exist, depending on the anharmonicity of the Morse potential and the cooperativity of the hydrogen bond (the coupling constant of the interaction between nearest-neighbor protons).

Collective proton transport with weak proton-proton coupling

A mechanism based on the peculiar dynamical properties of the hydrogen bond is suggested to describe fast proton transfers in hydrogen-bonded chains with realistic (weak) intersite proton-proton interactions. For this purpose, a two-sublattice model is suggested and studied analytically that admits an exact soliton solution. The solitons, which are topological objects (kinks and antikinks) on the one hand, but dynamical entities (like Boussinesq solitons) on the other hand, are shown to have a sufficiently large width to allow free propagation along the chain.