New Insights into the Dynamics of Concerted Proton Tunneling in Cyclic Water and Hydrogen Fluoride Clusters

Systematic analysis of electronic barrier heights and widths for concerted proton transfer in cyclic hydrogen bonded clusters: (HF)n, (HCl)n and (H2O)n where n = 3, 4, 5

The MP2 and CCSD(T) methods are paired with correlation consistent basis sets as large as aug-cc-pVQZ to optimize the structures of the cyclic minima for (HF)n, (HCl)n and (H2O)n where n = 3-5, as well as the corresponding transition states (TSs) for concerted proton transfer (CPT). MP2 and CCSD(T) harmonic vibrational frequencies confirm the nature of each minimum and TS. Both conventional and explicitly correlated CCSD(T) computations are employed to assess the electronic dissociation energies and barrier heights for CPT near the complete basis (CBS) limit for all 9 clusters. Results for (HF)n are consistent with prior studies identifying Cnh and Dnh point group symmetry for the minima and TSs, respectively. Our computations also confirm that CPT proceeds through Cs TS structures for the C1 minima of (H2O)3 and (H2O)5, whereas the process goes through a TS with D2d symmetry for the S4 global minimum of (H2O)4. This work corroborates earlier findings that the minima for (HCl)3, (HCl)4 and (HCl)5 have C3h, S4 and C1 point group symmetry, respectively, and that the Cnh structures are not minima for n = 4 and 5. Moreover, our computations show the TSs for CPT in (HCl)3, (HCl)4 and (HCl)5 have D3h, D2d, and C2 point group symmetry, respectively. At the CCSD(T) CBS limit, (HF)4 and (HF)5 have the smallest electronic barrier heights for CPT (≈15 kcal mol-1 for both), followed by the HF trimer (≈21 kcal mol-1). The barriers are appreciably higher for the other clusters (around 27 kcal mol-1 for (H2O)4 and (HCl)3; roughly 30 kcal mol-1 for (H2O)3, (H2O)5 and (HCl)4; up to 38 kcal mol-1 for (HCl)5). At the CBS limit, MP2 significantly underestimates the CCSD(T) barrier heights (e.g., by ca. 2, 4 and 7 kcal mol-1 for the pentamers of HF, H2O and HCl, respectively), whereas CCSD overestimates these barriers by roughly the same magnitude. Scaling the barrier heights and dissociation energies by the number of fragments in the cluster reveals strong linear relationships between the two quantities and with the magnitudes of the imaginary vibrational frequency for the TSs.

Unveiling the role of intra and interatomic interactions in the energetics of reaction schemes: a quantum chemical topology analysis

In this work we present a detailed analysis of selected reaction schemes in terms of the atomic components of the electronic energy defined by the quantum theory of atoms in molecules and the interacting quantum atoms method. The aim is to provide an interpretation tool for the energy change involved in a chemical reaction by means of the atomic and interaction contributions to the energies of the molecules involved. Ring strain in cyclic alkanes, the resonance energy of aromatic and antiaromatic molecules, local aromaticity in polycyclic aromatic hydrocarbons, intermolecular bonding in hydrogen fluoride clusters, and hydration of d-block metal dications were selected for the study. It was found that in addition to the changes in the strong C-C interactions in the carbon skeleton of the organic molecular rings, other contributions not usually considered to be important such as those between C and H atoms (either bonded or not) need to be considered in order to account for the net energy changes. The analysis unveils the role of the ionic and covalent contributions to the hydrogen bonding in HF clusters and the energetic origin and extent of cooperative effects involved. Moreover, the "double-hump" behavior observed for the hydration energy trend of [M(H2O)6]2+ complexes is explained in terms of the deformation energy of the metal cation and the increasingly covalent metal-water interactions. In addition, proper comparisons with the description provided by other methodologies are briefly discussed. The topological approach proposed in this contribution proves to be useful for the description of energy changes of apposite reaction schemes in chemically meaningful terms.

The mechanism of proton transfer between adjacent sites on the molecular surface

The surface of a protein, or a membrane, is spotted with a multitude of proton binding sites, some of which are only few A apart. When a proton is released from one site, it propagates through the water by a random walk under the bias of the local electrostatic potential determined by the distribution of the charges on the protein. Eventually, the released protons are dispersed in the bulk, but during the first few nanoseconds after the dissociation, the protons can be trapped by encounter with nearby acceptor sites. While the study of this reaction on the surface of a protein suffers from experimental and theoretical difficulties, it can be investigated with simple model compounds like derivatives of fluorescein. In the present study, we evaluate the mechanism of proton transfer reactions that proceed, preferentially, inside the Coulomb cage of the dye molecules. Kinetic analysis of the measured dynamics reveals the role of the dimension of the Coulomb cage on the efficiency of the reaction and how the ordering of the water molecules by the dye affects the kinetic isotope effect.

Dynamics of the Concerted Triple Proton Transfer in Cyclic Water Trimer Using the Multiconfiguration Molecular Mechanics Algorithm

Cyclic water clusters are important molecular species to understand the nature of hydrogen bonded networks. Theoretical studies for the dynamics of triple proton transfer in the cyclic water trimer were performed. The potential energy surface (PES) of triple proton transfer is generated by the multiconfiguration molecular mechanics (MCMM) algorithm. We have used the MP2/6-31G(d,p) level for high-level ab initio data (energies, gradients, and Hessians), which are used in the Shepard interpolation. Eight high-level reference points were added step by step, including two points for the critical configurations of the large curvature tunneling paths. The more high-level points are used, the better the potential energy surfaces become. The rate constant and kinetic isotope effect (KIE) for the triple proton transfer at 300 K, which have been calculated by the canonical variational transition-state theory with microcanonical optimized multidimensional semiclassical tunneling approximation, are 1.6 x 10(-3) s(-1) and 230, respectively. Tunneling is very important not only for the triple proton transfer but also for the triple deuterium transfer. The MCMM results show good agreement with those from the direct ab initio dynamics calculations.

An ab initio study of tunneling splittings in the water trimer

Tunneling splittings in the water trimer have been determined by the semiclassical WKB method, based on pathways characterized at the computational level of second-order Møller-Plesset theory with basis sets of aug-cc-pVTZ quality. This calculation takes into account the single-flip and bifurcation tunneling rearrangements of the water trimer. The predicted splittings are 37.93 cm(-1) for the flip and 6.50x10(-3) cm(-1) for bifurcation, which agree quite well with the corresponding experimental values of 43.52 cm(-1) and 9.63x10(-3) cm(-1).

The Mechanism of Proton Transfer between Adjacent Sites Exposed to Water

The surface of a protein, or a membrane, is spotted with a multitude of proton-binding sites, some of which are only a few angstroms apart. When a proton is released from one site, it propagates through the water by a random walk under the bias of the local electrostatic potential determined by the distribution of the charges on the protein. Some of the released protons disperse into the bulk, but during the first few nanoseconds, the released protons can be trapped by encounter with nearby acceptor sites. This process resembles a scenario which corresponds with the time-dependent Debye-Smoluchowski equation. In the present study, we investigated the mechanism of proton transfer between sites that are only a few angstroms apart, using as a model the proton exchange between sites on a small molecule, fluorescein, having two, spectrally distinguishable, proton-binding sites. The first site is the oxyanion on the chromophore ring structure. The second site is the carboxylate moiety on the benzene ring of the molecule. Through our experiments, we were able to reconstruct the state of protonation at each site and the velocity of proton transfer between them. The fluorescein was protonated by a few nanosecond long proton pulse under specific conditions that ensured that the dye molecules would be protonated only by a single proton. The dynamics of the protonation of the chromophore were measured under varying initial conditions (temperature, ionic strength, and different solvents (H(2)O or D(2)O)), and the velocity of the proton transfer between the two sites was extracted from the overall global analysis of the signals. The dynamics of the proton transfer between the two proton-binding sites of the fluorescein indicated that the efficiency of the site-to-site proton transfer is very sensitive to the presence of the screening electrolyte and has a very high kinetic isotope effect (KIE = 55). These two parameters clearly distinguish the mechanism from proton diffusion in bulk water. The activation energy of the reaction (E(a) = 11 kcal mol(-1)) is also significantly higher than the activation energy for proton dissociation in bulk water (E(a) approximately 2.5 kcal mol(-1)). These observations are discussed with respect to the effect of the solute on the water molecules located within the solvation layer.

An ab initio study of tunneling splittings in the water dimer

Tunneling splittings in the water dimer have been determined by the semiclassical WKB method, based on pathways characterized at the computational level of second-order Møller-Plesset (MP2) theory with basis sets of aug-cc-pVTZ quality. This calculation takes into account all three acceptor tunneling, donor-acceptor interchange, and bifurcation tunneling rearrangements of the water dimer. The tunneling splittings were evaluated as 7.73 cm(-1) (large splitting) and 0.42 cm(-1) (small splitting), which are in good agreement with the corresponding experimental values of 11.18 cm(-1) and 0.70 cm(-1), respectively.