Simultaneous proton transfer along a linear water polymer

SCF CNDO calculations were performed for species H₃O⁺·(H₂O)_n·OH^? where n was varied from one to three. The position of the intervening protons was changed simultaneously while the oxygens and remaining hydrogens were kept fixed. It was found that only one minimum occurs when n is one or two while an asymmetric double minimum potential is found when n is equal to three. A barrier of 10.4 kcal/mole was found.     

Excited-state hydrogen-atom transfer along solvent wires: water molecules stop the transfer

Excited-state hydrogen-atom transfer (ESHAT) along a hydrogen-bonded solvent wire occurs for the supersonically cooled n = 3 ammonia-wire cluster attached to the scaffold molecule 7-hydroxyquinoline (7HQ) [Tanner, C.; et al. Science 2003, 302, 1736]. Here, we study the analogous three-membered solvent-wire clusters 7HQ.(NH3)n.(H2O)m, n + m = 3, using resonant two-photon ionization (R2PI) and UV-UV hole-burning spectroscopies. Substitution of H2O for NH3 has a dramatic effect on the excited-state H-atom transfer: The threshold for the ESHAT reaction is approximately 200 cm(-1) for 7HQ.(NH3)3, approximately 350 cm(-1) for both isomers of the 7HQ.(NH3)2.H2O cluster, and approximately 600 cm(-1) for 7HQ.NH3.(H2O)2 but increases to approximately 2000 cm(-1) for the pure 7HQ.(H2O)3 water-wire cluster. To understand the effect of the chemical composition of the solvent wire on the H-atom transfer, the reaction profiles of the low-lying electronic excited states of the n = 3 pure and mixed solvent-wire clusters are calculated with the configuration interaction singles (CIS) method. For those solvent wires with an NH3 molecule at the first position, injection of the H atom into the wire can occur by tunneling. However, further H-atom transfer is blocked by a high barrier at the first (and second) H2O molecule along the solvent wire. H-atom transfer along the entire length of the solvent wire, leading to formation of the 7-ketoquinoline (7KQ) tautomer, cannot occur for any of the H2O-containing clusters, in agreement with experimentally observed absence of 7KQ fluorescence.     

Bimodal proton transfer in acid-base reactions in water

We investigate one of the fundamental reactions in solutions, the neutralization of an acid by a base. We use a photoacid, 8-hydroxy-1,3,6-trisulfonate-pyrene (HPTS; pyranine), which upon photoexcitation reacts with acetate under transfer of a deuteron (solvent: deuterated water). We analyze in detail the resulting bimodal reaction dynamics between the photoacid and the base, the first report on which was recently published. We have ascribed the bimodal proton-transfer dynamics to contributions from preformed hydrogen bonding complexes and from initially uncomplexed acid and base. We report on the observation of an additional (6 ps)(-1) contribution to the reaction rate constant. As before, we analyze the slower part of the reaction within the framework of the diffusion model and the fastest part by a static, sub-150 fs reaction rate. Adding the second static term considerably improves the overall modeling of the experimental results. It also allows to connect experimentally the diffusion controlled bimolecular reaction models as defined by Eigen-Weller and by Collins-Kimball. Our findings are in agreement with a three-stage mechanism for liquid phase intermolecular proton transfer: mutual diffusion of acid and base to form a "loose" encounter complex, followed by reorganization of the solvent shells and by "tightening" of the acid-base encounter complex. These rearrangements last a few picoseconds and enable a prompt proton transfer along the reaction coordinate, which occurs faster than our time resolution of 150 fs. Alternative models for the explanation of the slower "on-contact" reaction time of the loose encounter complex in terms of proton transmission through a von Grotthuss mechanism are also discussed.     

Energetics of electron transfer reactions in polar solvents

The free energy change of an electron transfer reaction in a polar solvent is rigorously analyzed within the framework of the dielectric continuum model. An appropriate expression for the electrostatic energy between the two product ions separated by R is derived. The present result does not support a recent claim by Suppan that, if R is close to the contact distance, the electrostatic energy should be much larger in magnitude than estimated from the usual expression −e²/ε_sR.     

Time-resolved studies of water dynamics and proton transfer at the alumina-air interface

The present study aims to understand the dynamical properties of water and OH groups layered on an alumina surface mainly by means of femtosecond IR-pump IR-probe transient absorption spectroscopy. The experimental results obtained demonstrate the existence of several kinds of O-H vibrators on the surface of alumina membranes, distinguishing them by their behavior on the femtosecond time scale and by the anisotropy of their spectral response. In the high-frequency region (>3400 cm-1), the absorption is due to well-packed aluminol groups and to physisorbed water patches on the surface. When pumping at 3200 cm-1, physisorbed water hydrogen-bonded to AlOH2+ groups is observed. The anisotropy measurements demonstrate the existence of an efficient energy-transfer mechanism among the water molecules characterized by a time constant of 400 +/- 100 fs. The persisting anisotropy at long times, especially in the case of AlOH groups and of the structured physisorbed water layer on top of them, proves the anisotropic structuring induced by the surface. The excitation at 3000 cm-1 enables the detection of a photon-induced proton-transfer reaction. The proton back-transfer reaction time constant is 350 +/- 50 fs. From anisotropy measurements, we estimate the proton hopping time to be 900 +/- 100 fs in a locally extended water network lying on the surface.     

Exploration on regulating factors for proton transfer along hydrogen-bonded water chains.

Proton transfer along a single-file hydrogen-bonded water chain is elucidated with a special emphasis on the investigation of chain length, side water, and solvent effects, as well as the temperature and pressure dependences. The number of water molecules in the chain varies from one to nine. The proton can be transported to the acceptor fragment through the single-file hydrogen-bonded water wire which contains at most five water molecules. If the number of water molecule is more than five, the proton is trapped by the chain in the hydroxyl-centered H(7)O(3) (+) state. The farthest water molecule involved in the formation of H(7)O(3) (+) is the fifth one away from the donor fragment. These phenomena reappear in the molecular dynamics simulations. The energy of the system is reduced along with the proton conduction. The proton transfer mechanism can be altered by excess proton. The augmentation of the solvent dielectric constant weakens the stability of the system, but favors the proton transfer. NMR spin-spin coupling constants can be used as a criterion in judging whether the proton is transferred or not. The enhancement of temperature increases the thermal motion of the molecule, augments the internal energy of the system, and favors the proton transfer. The lengthening of the water wire increases the entropy of the system, concomitantly, the temperature dependence of the Gibbs free energy increases. The most favorable condition for the proton transfer along the H-bonded water wire is the four-water contained chain with side water attached near to the acceptor fragment in polar solvent under higher temperature.     

Ab initio molecular dynamics study of glycine intramolecular proton transfer in water

We use ab initio molecular-dynamics simulations to quantify structural and thermodynamic properties of a model proton transfer reaction that converts a neutral glycine molecule, stable in the gas phase, to the zwitterion that predominates in aqueous solution. We compute the potential of mean force associated with the direct intramolecular proton transfer event in glycine. Structural analyses show that the average hydration number (N(w)) of glycine is not constant along the reaction coordinate, but rather progresses from N(w) = 5 in the neutral molecule to N(w) = 8 for the zwitterion. We report the free-energy difference between the neutral and charged glycine molecules, and the free-energy barrier to proton transfer. Finally, we identify the approximations inherent in our method and estimate the corresponding corrections to our reported thermodynamic predictions.     

Excited-state proton transfer in 7-hydroxy-4-methylcoumarin along a hydrogen-bonded water wire

TDDFT, RI-CC2, and CIS calculations have been performed for the nondissociative excited-state proton transfer (ESPT) in the S1 state of 7-hydroxy-4-methylcoumarin (7H4MC) along a H-bonded water wire of three water molecules bridging the proton donor (OH) and the proton acceptor (C[double bond]O) groups (7H4MC.(H2O)3). The observed structural reorganization in the water-wire cluster is interpreted as a proton-transfer (PT) reaction along the H2O solvent wire. The shift of electron density within the organic chromophore 7H4MC due to the optical excitation appears to be the driving force for ESPT. All the methods used show that the reaction path occurs in the 1pipi* state, and no crossing with a Rydberg-type 1pisigma* state is found. TDDFT and RI-CC2 calculations predict an exoergic reaction of the excited-state enol-to-keto transformation. The S1 potential energy curve reveals well-defined Cs minima of enol- and keto-clusters, separated by a single barrier with a height of 17-20 kcal/mol. After surmounting this barrier, spontaneous PT along the water wire is observed, leading without any further barrier to the keto structure. The TDDFT and RI-CC2 methods appear to be reliable approaches to describe the energy surfaces of ESPT. The CIS method predicts an endoergic ESPT reaction and an energy barrier, which is too high.     

Energetics and dynamics of electron transfer and proton transfer in dissociation of metal(III)(salen)-peptide complexes in the gas phase

Time- and collision energy-resolved surface-induced dissociation (SID) of ternary complexes of Co(III)(salen)+, Fe(III)(salen)+, and Mn(III)(salen)+ with several angiotensin peptide analogues was studied using a Fourier transform ion cyclotron resonance mass spectrometer (FT-ICR MS) specially equipped to perform SID experiments. Time-resolved fragmentation efficiency curves (TFECs) were modeled using an RRKM-based approach developed in our laboratory. The approach utilizes a very flexible analytical expression for the internal energy deposition function that is capable of reproducing both single-collision and multiple-collision activation in the gas phase and excitation by collisions with a surface. The energetics and dynamics of competing dissociation pathways obtained from the modeling provides important insight on the competition between proton transfer, electron transfer, loss of neutral peptide ligand, and other processes that determine gas-phase fragmentation of these model systems. Similar fragmentation behavior was obtained for various Co(III)(salen)-peptide systems of different angiotensin analogues. In contrast, dissociation pathways and relative stabilities of the complexes changed dramatically when cobalt was replaced with trivalent iron or manganese. We demonstrate that the electron-transfer efficiency is correlated with redox properties of the metal(III)(salen) complexes (Co > Fe > Mn), while differences in the types of fragments formed from the complexes reflect differences in the modes of binding between the metal-salen complex and the peptide ligand. RRKM modeling of time- and collision-energy-resolved SID data suggests that the competition between proton transfer and electron transfer during dissociation of Co(III)(salen)-peptide complexes is mainly determined by differences in entropy effects while the energetics of these two pathways are very similar.     

Excited-state proton transfer reactions of 10-hydroxycamptothecin

Time-resolved and steady-state emission characterization of 10-hydroxycamptothecin reveals a rich but less complex proton-transfer behavior than its parent hydroxyquinoline. The electronic effect of the additional electron-withdrawing ring makes the excited-state both less basic and more acidic than the parent and adds to the class of high-acidity excited-state proton donors in photochemistry and photobiology.     

Reversible proton transfer dynamics in bacteriorhodopsin

In a previous study of ab initio dynamics, the proton transfer in bacteriorhodopsin from protonated asp96 in the cytoplasmic region toward the deprotonated Schiff base was investigated. A quantum mechanics/molecular mechanics model was constructed from the X-ray structure of bacteriorhodopsin E204Q mutant. In this model, asp96, asp85, and thr89 as well as most of the retinal chromophore and the Schiff base link of lys216 were treated quantum mechanically while the rest of the atoms were treated molecular mechanically. A channel was found in the X-ray structure allowing a water chain to form between the asp96 and Schiff base. In the present study, a chain of four waters from asp96 to the Schiff base N coupled with one branching water supports proton transfer as a concerted event in about 3.5 ps. With both a neutral asp85 and a branched water, the dynamics is now found to be more complicated than observed in the initial study for the transition from the photocycle late M state to the N state. Proton transfer is also observed from the Schiff base back to asp96 demonstrating that there is no effective barrier to proton transfer larger than kT in a strong H-bonded network. The binding of the branched water to the four water chains can dynamically hinder the proton transfer.     

Proton-coupled electron transfer across benzimidazole bridges in bioinspired proton wires

Designing molecular platforms for controlling proton and electron movement in artificial photosynthetic systems is crucial to efficient catalysis and solar energy conversion. The transfer of both protons and electrons during a reaction is known as proton-coupled electron transfer (PCET) and is used by nature in myriad ways to provide low overpotential pathways for redox reactions and redox leveling, as well as to generate bioenergetic proton currents. Herein, we describe theoretical and electrochemical studies of a series of bioinspired benzimidazole-phenol (BIP) derivatives and a series of dibenzimidazole-phenol (BI₂P) analogs with each series bearing the same set of terminal proton-accepting (TPA) groups. The set of TPAs spans more than 6 pK_a units. These compounds have been designed to explore the role of the bridging benzimidazole(s) in a one-electron oxidation process coupled to intramolecular proton translocation across either two (the BIP series) or three (the BI₂P series) acid/base sites. These molecular constructs feature an electrochemically active phenol connected to the TPA group through a benzimidazole-based bridge, which together with the phenol and TPA group form a covalent framework supporting a Grotthuss-type hydrogen-bonded network. Infrared spectroelectrochemistry demonstrates that upon oxidation of the phenol, protons translocate across this well-defined hydrogen-bonded network to a TPA group. The experimental data show the benzimidazole bridges are non-innocent participants in the PCET process in that the addition of each benzimidazole unit lowers the redox potential of the phenoxyl radical/phenol couple by 60 mV, regardless of the nature of the TPA group. Using a series of hypothetical thermodynamic steps, density functional theory calculations correctly predicted the dependence of the redox potential of the phenoxyl radical/phenol couple on the nature of the final protonated species and provided insight into the thermodynamic role of dibenzimidazole units in the PCET process. This information is crucial for developing molecular “dry proton wires” with these moieties, which can transfer protons via a Grotthuss-type mechanism over long distances without the intervention of water molecules.

Experimental and theoretical methods characterize the thermodynamics of electrochemically driven proton-coupled electron transfer processes in bioinspired constructs involving multiple proton translocations over Grotthus-type proton wires.     

Faster proton transfer dynamics of water on SnO2 compared to TiO2

Proton jump processes in the hydration layer on the iso-structural TiO(2) rutile (110) and SnO(2) cassiterite (110) surfaces were studied with density functional theory molecular dynamics. We find that the proton jump rate is more than three times faster on cassiterite compared with rutile. A local analysis based on the correlation between the stretching band of the O-H vibrations and the strength of H-bonds indicates that the faster proton jump activity on cassiterite is produced by a stronger H-bond formation between the surface and the hydration layer above the surface. The origin of the increased H-bond strength on cassiterite is a combined effect of stronger covalent bonding and stronger electrostatic interactions due to differences of its electronic structure. The bridging oxygens form the strongest H-bonds between the surface and the hydration layer. This higher proton jump rate is likely to affect reactivity and catalytic activity on the surface. A better understanding of its origins will enable methods to control these rates.     

Quantum-classical Liouville dynamics of nonadiabatic proton transfer

A proton transfer reaction in a linear hydrogen-bonded complex dissolved in a polar solvent is studied using mixed quantum-classical Liouville dynamics. In this system, the proton is treated quantum mechanically and the remainder of the degrees of freedom is treated classically. The rates and mechanisms of the reaction are investigated using both adiabatic and nonadiabatic molecular dynamics. We use a nonadiabatic dynamics algorithm which allows the system to evolve on single adiabatic surfaces and on coherently coupled pairs of adiabatic surfaces. Reactive-flux correlation function expressions are used to compute the rate coefficients and the role of the dynamics on the coherently coupled surfaces is elucidated.     

Long-range proton transfer in aqueous Acid-base reactions

We study the mechanism of proton transfer (PT) in the aqueous acid-base reaction between the photoacid 8-hydroxy-1,3,6-pyrenetrisulfonic acid (HPTS) and acetate by probing the vibrational resonances of HPTS, acetate, and the hydrated proton with femtosecond mid-infrared laser pulses. We find that PT takes place in a distribution of hydrogen-bound reaction complexes that differ in the number of water molecules separating the acid and the base. The number of intervening water molecules ranges from 0 to 5, which, together with a strongly distance-dependent PT rate, explains the observed highly nonexponential reaction kinetics. The kinetic isotope effect for the reaction is determined to be 1.5, indicating that tunneling does not play a significant role in the transfer of the proton. Rather, the transfer mechanism is best described in terms of the adiabatic PT picture as it has been formulated by Hynes and co-workers [Staib, A.; Borgis, D.; Hynes, J. T. J. Chem. Phys. 1995, 102, 2487. Ando, K.; Hynes, J. T. J. Phys. Chem. B 1997, 101, 10464.], where solvent fluctuations play an essential role in forming the correct hydrogen-bond configuration and solvent polarization to facilitate PT.     

Ultrafast proton transfer dynamics of hydroxystilbene photoacids

The effects of 4-cyano and 3-cyano substituents on the spectroscopic properties and photoacidity of 3- and 4-hydroxystilbene have been investigated. In nonpolar solvents, the 3-hydroxycyanostilbenes have much longer singlet lifetimes and larger fluorescence quantum yields than do the 4-hydroxycyanostilbenes. The longer lifetimes of 3-hydroxystilbene and its cyano derivatives are attributed to a "meta effect" on the stilbene torsional barrier, similar to that previously observed for the aminostilbenes. The cyano substituent causes a marked increase in both ground state and excited-state acidity of the hydroxystilbenes in aqueous solution. The dynamics of excited-state proton transfer in methanol-water solution have been investigated by means of femtosecond time-resolved transient absorption spectroscopy. Assignment of the transient absorption spectra is facilitated by comparison to the spectra of the corresponding potassium salts of the conjugate bases and the methyl ethers, which do not undergo excited-state proton transfer. The 4-cyanohydroxystilbenes undergo excited-state proton transfer with rate constants of 5 x 10(11) s(-1). These rate constants are comparable to the fastest that have been reported to date for a hydroxyaromatic photoacid and approach the theoretical limit for water-mediated proton transfer. The isotope effect for proton transfer in deuterated methanol-water is 1.3 +/- 0.2, similar to the isotope effect for the dielectric response of water. The barrier for excited state double bond torsion of the conjugate bases is small for 4-cyano-4-hydroxystilbene but large for 4-cyano-3-hydroxystilbene. Thus the "meta effect" is observed for the singlet states of both the neutral and conjugate base.     

Gas-phase proton transfer reactions involving multiply charged cytochrome c ions and water under thermal conditions

Investigations of gas-phase proton transfer reactions have been performed on protein molecular ions generated by electrospray ionization (ESI). Their reactions were studied in a heated capillary inlet/reactor prior to expansion into a quadrupole mass spectrometer. Results from investigations involving protonated horse heart cytochrome c and H, O suggest that Coulombit effects can lower reaction barriers as well as aid in entropically driven reactions. For example, the charge state distribution observed by a quadrupole mass spectrometer for multiply protonated cytochrome c without the addition of any reactive gas ranges from 9+ to 19+ , with the [M + 15H]¹⁵⁺ ion being the most intense peak. With the addition of H₂O (proton affinity approximately 170.3±2 kcal/mol) to the capillary reactor at 120°C, the charge state distribution shifts to a lower charge, ranging from 13+ to less than 9+. Under the same conditions with argon (proton affinity approximately 100 kcal/mol) as the reactive gas, no shift in the charge state distribution is observed. The results demonstrate that proton transfer to water can occur for highly protonated molecular ions, a process that would be expected to be highly endothermic for singly protonated molecules (for which Coulombic destabilization is not significant). The results imply that the charge state distribution from ESI is somewhat dependent upon the mechanism and speed of the droplet evaporation/ion desolvation process, which may vary substantially with the ESI/mass spectrometry interface design.     

Energetics during proton transfer process in hydrated zündel ion complex

Zündel ion (H₅O) is one of the two important structures formed during the proton transfer process in aqueous system. This work reports microsolvation of Zündel ion using density functional theory based B3LYP method with aug‐cc‐pVTZ basis set. Interaction of Zündel ion with four water molecules in its first solvation shell is studied using many‐body analysis approach. A change in many‐body energies and their contribution to the binding energy of a complex during the proton transfer process from donor to acceptor water molecule in Zündel ion‐4H₂O complex is obtained. For the hydrated Zündel ion complex, the contribution from total two‐body, three‐body, four‐body, five‐body, and relaxation energy to the binding energy is 84.7, 14, 6.87, 1.6, and 4%, respectively, at B3LYP/aug‐cc‐pVTZ level. Relaxation energy and total five‐body energy have repulsive contribution to the binding energy of a hydrated Zündel ion complex. It is found that the relaxation energy and binding energy of a Zündel‐4H₂O complex is the maximum and minimum, respectively, when a shared proton is at equal distance from oxygen atom of donor and acceptor water molecules. A significant change in two‐body, three‐body, and four‐body energies for which Zündel ion is one of the many‐body terms is observed during the proton transfer process. A change in total two‐body, total three‐body, total four‐body, and relaxation energy is about 2.6, 1.8, 0.4, and 1.1%, respectively, during the proton transfer process. A change in two‐body, three‐body, and four‐body interaction energies between water molecules is very small during the proton transfer process. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

Energetics and electronic rearrangements of proton transfer in (H3NHOH2)+

Proton transfer between N and O in the hydrogen-bonded system (H₃NHOH₂)⁺ is studied by ab initio molecular orbital methods. Potential energy curves are calculated at the hartree–Fock level using the 4–31G basis set for hydrogen bond lengths R(NO) varying from the equilibrium value of 2.664 to 3.10 Å. Short hydrogen bonds are associated with asymmetric single-well potentials in which the minimum corresponds to the NH? O configuration. For longer R(NO) separations, the potential is of double-well form, including both N? HO and NH? O as minima. It is found that the height of the energy barrier to proton transfer is sensitive to both stretches and bends of the hydrogen bond. Continuous changes in the electron density are monitored at various stages of proton transfer via density difference maps and Mulliken population analyses. The initial loss of density from the proton-accepting molecule during the first half of the transfer is accelerated during the second half. A correlation is drawn between the energetics of transfer in a number of systems and the net charge lost by the proton-acceptor group.