Direct simulation of proton-coupled electron transfer reaction dynamics and mechanisms

Proton-coupled electron transfer (PCET) reactions are ubiquitous throughout chemistry and biology. However, challenges arise in both the the experimental and theoretical investigation of PCET reactions; the rare-event nature of the reactions and the coupling between quantum mechanical electron- and proton-transfer with the slower classical dynamics of the surrounding environment necessitates the development of robust simulation methodology. In the following dissertation, novel path-integral based methods are developed and employed for the direct simulation of the reaction dynamics and mechanisms of condensed-phase PCET.

Proton uptake controls electron transfer in cytochrome c oxidase

In cytochrome c oxidase, a requirement for proton pumping is a tight coupling between electron and proton transfer, which could be accomplished if internal electron-transfer rates were controlled by uptake of protons. During reaction of the fully reduced enzyme with oxygen, concomitant with the “peroxy” to “oxoferryl” transition, internal transfer of the fourth electron from CuA to heme a has the same rate as proton uptake from the bulk solution (8,000 s−1). The question was therefore raised whether the proton uptake controls electron transfer or vice versa. To resolve this question, we have studied a site-specific mutant of the Rhodobacter sphaeroides enzyme in which methionine 263 (SU II), a CuA ligand, was replaced by leucine, which resulted in an increased redox potential of CuA. During reaction of the reduced mutant enzyme with O2, a proton was taken up at the same rate as in the wild-type enzyme (8,000 s−1), whereas electron transfer from CuA to heme a was impaired. Together with results from studies of the EQ(I-286) mutant enzyme, in which both proton uptake and electron transfer from CuA to heme a were blocked, the results from this study show that the CuA → heme a electron transfer is controlled by the proton uptake and not vice versa. This mechanism prevents further electron transfer to heme a3–CuB before a proton is taken up, which assures a tight coupling of electron transfer to proton pumping.

Oxygen reduction and proton translocation by cytochrome c oxidase

Energy conversion by living organisms is central dogma of bioenergetics. The effectiveness of the energy extraction by aerobic organisms is much greater than by anaerobic ones. In aerobic organisms the final stage of energy conversion occurs in respiratory chain that is located in the inner membrane of mitochondria or cell membrane of some aerobic bacteria. The terminal complex of the respiratory chain is cytochrome c oxidase (CcO) - the subject of this study. The primary function of CcO is to reduce oxygen to water. For this, CcO accepts electrons from a small soluble enzyme cytochrome c from one side of the membrane and protons from another side. Moreover, CcO translocates protons across the membrane. Both oxygen reduction and proton translocation contributes to generation of transmembrane electrochemical gradient that is used for ATP synthesis and different types of work in the cell. Although the structure of CcO is defined with a relatively high atomic resolution (1.8 Å), its function can hardly be elucidated from the structure. The electron transfer route within CcO and its steps are very well defined. Meanwhile, the proton transfer roots were predicted from the site-specific mutagenesis and later proved by X-ray crystallography, however, the more strong proof of the players of the proton translocation machine is still required. In this work we developed new methods to study CcO function based on FTIR (Fourier Transform Infrared) spectroscopy. Mainly with use of these methods we answered several questions that were controversial for many years: [i] the donor of H+ for dioxygen bond splitting was identified and [ii] the protolytic transitions of Glu-278 one of the key amino acid in proton translocation mechanism was shown for the first time.

An ab-initio study of the proton affinity of conjugated schiff-base and related nitrogen compounds: an analysis of the triggering site of bacteriorhodopsin

The chemical groups which take part in the proton transfer reaction in bacteriorhodopsin have been studied by ab initio quantum chemical methods. The various factors such as conjugation with a linear system, electron delocalization of the guanidine type, cis-trans isomerism, geometry distortion and hydrogen bonding with charged groups can influence the properties of a given chemical group. Several systems are studied at 4-31G and STO-3G levels. Some of the Schiff-base analogues and guanidine type molecules are characterized by their molecular orbital diagrams, energy levels and the nature of charge distribution. Also, the effects of the above-mentioned factors on proton affinity are studied. It is hoped that the values thus obtained can be helpful in evaluating various structural models for proton transfer.

Drug-drug co-crystals: Temperature-dependent proton mobility in the molecular complex of isoniazid with 4-aminosalicylic acid

Two drug-drug co-crystals of the anti-tuberculosis drugs isoniazid (INH), pyrazinamide (PYR) and 4-aminosalicylic acid (PAS) are reported. The first is the 1 : 1 molecular complex of INH and PAS. The second is the monohydrate of the 1 : 1 complex of PYR and PAS. The crystal structures of both co-crystals are characterized by a number of hydrogen bonded synthons. Hydrogen bonding of the COOH center dot center dot center dot N-pyridine type is found in both cases. In the INH : PAS co-crystal, there are two symmetry independent COOH center dot center dot center dot center dot N-pyridine hydrogen bonds. In one of these, the H-atom is located on the carboxylic group and is indicative of a co-crystal. In the second case, partial proton transfer occurs across the hydrogen bond, and the extent of proton transfer depends on the temperature. This is more indicative of a salt. Drug-drug co-crystals may have some bearing in the treatment of tuberculosis.

Proton transfers at the air-water interface

Proton transfer reactions at the interface of water with hydrophobic media, such as air or lipids, are ubiquitous on our planet. These reactions orchestrate a host of vital phenomena in the environment including, for example, acidification of clouds, enzymatic catalysis, chemistries of aerosol and atmospheric gases, and bioenergetic transduction. Despite their importance, however, quantitative details underlying these interactions have remained unclear. Deeper insight into these interfacial reactions is also required in addressing challenges in green chemistry, improved water quality, self-assembly of materials, the next generation of micro-nanofluidics, adhesives, coatings, catalysts, and electrodes. This thesis describes experimental and theoretical investigation of proton transfer reactions at the air-water interface as a function of hydration gradients, electrochemical potential, and electrostatics. Since emerging insights hold at the lipid-water interface as well, this work is also expected to aid understanding of complex biological phenomena associated with proton migration across membranes.

Based on our current understanding, it is known that the physicochemical properties of the gas-phase water are drastically different from those of bulk water. For example, the gas-phase hydronium ion, H₃O⁺(g), can protonate most (non-alkane) organic species, whereas H₃O⁺(aq) can neutralize only relatively strong bases. Thus, to be able to understand and engineer water-hydrophobe interfaces, it is imperative to investigate this fluctuating region of molecular thickness wherein the ‘function’ of chemical species transitions from one phase to another via steep gradients in hydration, dielectric constant, and density. Aqueous interfaces are difficult to approach by current experimental techniques because designing experiments to specifically sample interfacial layers (< 1 nm thick) is an arduous task. While recent advances in surface-specific spectroscopies have provided valuable information regarding the structure of aqueous interfaces, but structure alone is inadequate to decipher the function. By similar analogy, theoretical predictions based on classical molecular dynamics have remained limited in their scope.

Recently, we have adapted an analytical electrospray ionization mass spectrometer (ESIMS) for probing reactions at the gas-liquid interface in real time. This technique is direct, surface-specific,and provides unambiguous mass-to-charge ratios of interfacial species. With this innovation, we have been able to investigate the following:

1. How do anions mediate proton transfers at the air-water interface?

2. What is the basis for the negative surface potential at the air-water interface?

3. What is the mechanism for catalysis ‘on-water’?

In addition to our experiments with the ESIMS, we applied quantum mechanics and molecular dynamics to simulate our experiments toward gaining insight at the molecular scale. Our results unambiguously demonstrated the role of electrostatic-reorganization of interfacial water during proton transfer events. With our experimental and theoretical results on the ‘superacidity’ of the surface of mildly acidic water, we also explored implications on atmospheric chemistry and green chemistry. Our most recent results explained the basis for the negative charge of the air-water interface and showed that the water-hydrophobe interface could serve as a site for enhanced autodissociation of water compared to the condensed phase.

Quantum Mechanics Studies of Fuel Cell Catalysts and Proton Conducting Ceramics with Validation by Experiment

We carried out quantum mechanics (QM) studies aimed at improving the performance of hydrogen fuel cells. This led to predictions of improved materials, some of which were subsequently validated with experiments by our collaborators.

In part I, the challenge was to find a replacement for the Pt cathode that would lead to improved performance for the Oxygen Reduction Reaction (ORR) while remaining stable under operational conditions and decreasing cost. Our design strategy was to find an alloy with composition Pt3M that would lead to surface segregation such that the top layer would be pure Pt, with the second and subsequent layers richer in M. Under operating conditions we expect the surface to have significant O and/or OH chemisorbed on the surface, and hence we searched for M that would remain segregated under these conditions. Using QM we examined surface segregation for 28 Pt₃M alloys, where M is a transition metal. We found that only Pt₃Os and Pt₃Ir showed significant surface segregation when O and OH are chemisorbed on the catalyst surfaces. This result indicates that Pt₃Os and Pt3Ir favor formation of a Pt-skin surface layer structure that would resist the acidic electrolyte corrosion during fuel cell operation environments. We chose to focus on Os because the phase diagram for Pt-Ir indicated that Pt-Ir could not form a homogeneous alloy at lower temperature. To determine the performance for ORR, we used QM to examine all intermediates, reaction pathways, and reaction barriers involved in the processes for which protons from the anode reactions react with O₂ to form H₂O. These QM calculations used our Poisson-Boltzmann implicit solvation model include the effects of the solvent (water with dielectric constant 78 with pH 7 at 298K). We found that the rate determination step (RDS) was the O_ad hydration reaction (O_ad + H₂O_ad -> OH_ad + OH_ad) in both cases, but that the barrier for pure Pt of 0.50 eV is reduced to 0.48 eV for Pt₃Os, which at 80 degrees C would increase the rate by 218%. We collaborated with the Pu-Wei Wu’s group to carry out experiments, where we found that the dealloying process-treated Pt2Os catalyst showed two-fold higher activity at 25 degrees C than pure Pt and that the alloy had 272% improved stability, validating our theoretical predictions.

We also carried out similar QM studies followed by experimental validation for the Os/Pt core-shell catalyst fabricated by the underpotential deposition (UPD) method. The QM results indicated that the RDS for ORR is a compromise between the OOH formation step (0.37 eV for Pt, 0.23 eV for Pt_2ML/Os core-shell) and H₂O formation steps (0.32 eV for Pt, 0.22 eV for Pt_2ML/Os core-shell). We found that Pt_2ML/Os has the highest activity (compared to pure Pt and to the Pt₃Os alloy) because the 0.37 eV barrier decreases to 0.23 eV. To understand what aspects of the core shell structure lead to this improved performance, we considered the effect on ORR of compressing the alloy slab to the dimensions of pure Pt. However this had little effect, with the same RDS barrier 0.37 eV. This shows that the ligand effect (the electronic structure modification resulting from the Os substrate) plays a more important role than the strain effect, and is responsible for the improved activity of the core- shell catalyst. Experimental materials characterization proves the core-shell feature of our catalyst. The electrochemical experiment for Pt_2ML/Os/C showed 3.5 to 5 times better ORR activity at 0.9V (vs. NHE) in 0.1M HClO₄ solution at 25 degrees C as compared to those of commercially available Pt/C. The excellent correlation between experimental half potential and the OH binding energies and RDS barriers validate the feasibility of predicting catalyst activity using QM calculation and a simple Langmuir–Hinshelwood model.

In part II, we used QM calculations to study methane stream reforming on a Ni-alloy catalyst surfaces for solid oxide fuel cell (SOFC) application. SOFC has wide fuel adaptability but the coking and sulfur poisoning will reduce its stability. Experimental results suggested that the Ni4Fe alloy improves both its activity and stability compared to pure Ni. To understand the atomistic origin of this, we carried out QM calculations on surface segregation and found that the most stable configuration for Ni₄Fe has a Fe atom distribution of (0%, 50%, 25%, 25%, 0%) starting at the bottom layer. We calculated that the binding of C atoms on the Ni4Fe surface is 142.9 Kcal/mol, which is about 10 Kcal/mol weaker compared to the pure Ni surface. This weaker C binding energy is expected to make coke formation less favorable, explaining why Ni₄Fe has better coking resistance. This result confirms the experimental observation. The reaction energy barriers for CHx decomposition and C binding on various alloy surface, Ni₄X (X=Fe, Co, Mn, and Mo), showed Ni₄Fe, Ni₄Co, and Fe₄Mn all have better coking resistance than pure Ni, but that only Ni₄Fe and Fe₄Mn have (slightly) improved activity compared to pure Ni.

In part III, we used QM to examine the proton transport in doped perovskite-ceramics. Here we used a 2x2x2 supercell of perovskite with composition Ba₈X₇M₁(OH)₁O₂₃ where X=Ce or Zr and M=Y, Gd, or Dy. Thus in each case a 4⁺ X is replace by a 3⁺ M plus a proton on one O. Here we predicted the barriers for proton diffusion allowing both includes intra-octahedron and inter-octahedra proton transfer. Without any restriction, we only observed the inter-octahedra proton transfer with similar energy barrier as previous computational work but 0.2 eV higher than experimental result for Y doped zirconate. For one restriction in our calculations is that the O_donor-O_acceptor atoms were kept at fixed distances, we found that the barrier difference between cerates/zirconates with various dopants are only 0.02~0.03 eV. To fully address performance one would need to examine proton transfer at grain boundaries, which will require larger scale ReaxFF reactive dynamics for systems with millions of atoms. The QM calculations used here will be used to train the ReaxFF force field.

Ultrasonic relaxation kinetics on fast deuteron transfer reaction

Propylamine has been selected to investigate the isotope effect of a fast deuteron transfer reaction by ultrasonic relaxation method. Ultrasonic absorption coefficients of propylamine in heavy water (D2O) at 25 degrees C in the concentration range from 0.0107 to 0.6300 mol dm(-3) have been measured by pulse and resonance methods over the frequency range from 0.8 to 220 MHz. A Debye-type single relaxation absorption has been observed in the solution. From the dependence of the ultrasonic relaxation parameters on the concentration and solution pH, the source of the observed relaxation has been attributed to a perturbation of the chemical equilibrium associated with the deuteron transfer reaction. The rate and equilibrium constants have been determined by the measurement of the deuteroxyl ion concentration dependence of the relaxation frequency. Also the standard volume change of the reaction has been determined from the concentration dependence of the maximum absorption per wavelength and the adiabatic compressibility has been calculated from the density and the sound velocity in the solution. These results have then been compared with those obtained for propylamine in light water (H2O). The forward rate constant is greater and the reverse rate constant is smaller in DO than in H2O. The standard volume change for deuteron transfer is greater than that for proton transfer reaction, and the adiabatic compressibility shows a similar trend. These data support an argument that there exists a stronger hydrogen bond in D2O than in H2O. The difference of the stability in the intermediate states, R-ND3+... OD- and R-NH3+... OH-, has also been considered from the results of the isotope effects.

Air-sea exchange of methanol and acetone during HiWinGS: Estimation of air phase, water phase gas transfer velocities

The air-sea fluxes of methanol and acetone were measured concurrently using a proton-transfer-reaction mass spectrometer (PTR-MS) with the eddy covariance (EC) technique during the High Wind Gas Exchange Study (HiWinGS) in 2013. The seawater concentrations of these compounds were also measured twice daily with the same PTR-MS coupled to a membrane inlet. Dissolved concentrations near the surface ranged from 7 to 28 nM for methanol and from 3 to 9 nM for acetone. Both gases were consistently transported from the atmosphere to the ocean as a result of their low sea surface saturations. The largest influxes were observed in regions of high atmospheric concentrations and strong winds (up to 25 m s(-1)). Comparison of the total air-sea transfer velocity of these two gases (K-a), along with the in situ sensible heat transfer rate, allows us to constrain the individual gas transfer velocity in the air phase (k(a)) and water phase (k(w)). Among existing parameterizations, the scaling of k(a) from the COARE model is the most consistent with our observations. The k(w) we estimated is comparable to the tangential (shear driven) transfer velocity previously determined from measurements of dimethyl sulfide. Lastly, we estimate the wet deposition of methanol and acetone in our study region and evaluate the lifetimes of these compounds in the surface ocean and lower atmosphere with respect to total (dry plus wet) atmospheric deposition.

Hydrogen bonding and collective proton modes in clusters and periodic layers of squaric acid: A density functional study

Hydrogen bonding in clusters and extended layers of squaric acid molecules has been investigated by density functional computations. Equilibrium geometries, harmonic vibrational frequencies, and energy barriers for proton transfer along hydrogen bonds have been determined using the Car-Parrinello method. The results provide crucial parameters for a first principles modeling of the potential energy surface, and highlight the role of collective modes in the low-energy proton dynamics. The importance of quantum effects in condensed squaric acid systems has been investigated, and shown to be negligible for the lowest-energy collective proton modes. This information provides a quantitative basis for improved atomistic models of the order-disorder and displacive transitions undergone by squaric acid crystals as a function of temperature and pressure. (C) 2001 American Institute of Physics.

Association and liquid structure of pyridine-acetic acid mixtures determined from neutron scattering using a ‘free proton’ EPSR simulation model

The liquid structure of pyridine-acetic acid mixtures have been investigated using neutron scattering at various mole fractions of acetic acid, χHOAc = 0.33, 0.50, and 0.67, and compared to the structures of neat pyridine and acetic acid. Data has been modelled using Empirical Potential Structure Refinement (EPSR) with a ‘free proton’ reference model, which has no prejudicial weighting towards either the existence of molecular or ionised species. Analysis of the neutron scattering results shows the existence of hydrogen-bonded acetic acid chains with pyridine inclusions, rather than the formation of an ionic liquid by proton transfer.

Electrochemical promotion of catalysis using solid-state proton-conducting membranes

A solid-state electrochemical reactor with ceramic proton-conducting membrane has been used to study the effect of electrochemically induced hydrogen spillover on the catalytic activity of platinum during ethylene oxidation. Suitable proton-conducting electrolyte membranes (Gd-doped BaPrO ₃ (BPG) and Y-doped BaZrO₃ (BZY)) were fabricated. These materials were chosen because of their protonic conductivity in the operational temperature region of the reaction (400-700 °C). The BZY-based electrochemical cell was used to investigate the open-circuit voltage (OCV) dependence on H₂ partial pressure with comparison being made to the theoretical OCV as predicted by the Nernst equation. Furthermore, the BZY pellets were used to study the effect of proton transfer of the catalytic activity of platinum during ethylene oxidation. The reaction was found to exhibit electrochemical promotion at 400 °C and to be electrophilic in nature, i.e. proton addition to the platinum surface resulted in an increase in reaction rate. At higher temperatures, the rate was not affected, within experimental error, by proton addition or removal. Under similar conditions, AC impedance showed that there was a large overall cell resistance at 400 °C with significantly decreased resistance at higher temperatures. It is possible that there could be a relationship between large cell resistances and the onset of electrochemical promotion in this system but there is, as yet, no conclusive evidence for this. © 2003 Elsevier B.V. All rights reserved.

Transfer reactions in the investigation of light-nuclei nucleosynthesis

Cross sections for the (6)Li(p,gamma)(7)Be, (7)Li(n,gamma)(8)Li (8)Li(n,gamma)(9)Li and (8)Li(p,gamma)(9)Be capture reactions have been investigated in the framework of the potential model. The main ingredients of the potential model are the potentials used to generate the continuum and bound-state wave functions and spectroscopic factors of the corresponding bound systems. The spectroscopic factors for the (7)Li circle times n=(8)Li(gs), (8)Li circle times n=(9)Li(gs) bound systems were obtained from a FR-DWBA analysis of neutron transfer reactions induced by (8)Li radioactive beam on a (9)Be target, while spetroscopic factor for the (8)Li circle times n=(9)Be(gs) bound system were obained from a proton transfer reaction. From the obtained capture reaction cross section, reaction rate for the (8)Li(n,gamma)(9)Li and (8)Li(p,gamma)(9)Be direct neutron and proton capture were determined and compared with other experimental and calculated values.