Biomedical subjects

Christofer S Tautermann

Publications and source records attributed to Christofer S Tautermann.

9 recordsLinked to original sources

PubMed2006

Comparative study of cluster- and supercell-approaches for investigating heterogeneous catalysis by electronic structure methods: tunneling in the reaction N + H --> NH on Ru(0001).

Different ruthenium clusters of various sizes are constructed with the aim to model the Ru(0001) surface with a sufficient accuracy for predicting catalysis by hybrid density functional methods (B3LYP). As an example reaction the hydrogenation step N(ads) + H(ads) --> NH(ads) from the catalytic production cycle of ammonia is chosen. A cluster of 12 ruthenium atoms is found to reproduce experimental geometries and frequencies of the various reactants on the surface satisfyingly. To get the geometries of adsorbed hydrogen qualitatively correct it is shown that second layer atoms have to be included in the model cluster. Boundary effects are believed to have minor effects on optimized geometries, whereas the effects on reaction barriers are significant. A comparison of model cluster calculations to a periodic supercell approach employing plane waves and density functional methods (RPBE) reveals similar barriers for reaction. The influence of tunneling in this reaction is determined by the small curvature tunneling approach on the electronic surfaces.

Journal Article↗

PubMed2005

The importance of tunneling in the first hydrogenation step in ammonia synthesis over a Ru(0001) surface.

The hydrogenation of nitrogen (N(ads)+H(ads)-->NH(ads)) on metal surfaces is an important step in ammonia catalysis. We investigate the reaction dynamics of this hydrogenation step by time independent scattering theory and variational transition state theory (VTST) including tunneling corrections. The potential energy surface is derived by hybrid density functional theory on a model cluster composed of 12 ruthenium atoms resembling a Ru(0001) surface. The scattering calculations are performed on a reduced dimensionality potential energy hypersurface, where two dimensions are treated explicitly and all others are included implicitly by the zero-point correction. The VTST calculations include quantum effects along the reaction coordinate by applying the small curvature tunneling scheme. Even at room temperature (where ruthenium already shows catalytic activity) we find rate enhancement by tunneling by a factor of approximately 70. Inspection of the reaction probabilities shows that the major contribution to reactivity comes from the vibrational ground state of the reactants into vibrationally excited product states. The reaction rates are higher than determined in previous studies, and are compatible with experimental overall rates for ammonia synthesis.

Journal Article↗

PubMed2004

Double hydrogen tunneling revisited: the breakdown of experimental tunneling criteria.

Formic acid dimer was chosen as a model system to investigate synchronous double proton transfer by means of variational transition state theory (VTST) for various isotopically modified hydrogen species. The electronic barrier for the double proton transfer was evaluated to be 7.9 kcal/mol, thus being significantly lower than it was determined in previous studies. The tunneling probabilities were evaluated at temperatures from 100 up to 400 K and typical Arrhenius behavior with enhancement by tunneling is observed. When comparing the transmission factors kappa in dependence of the mass of the tunneling hydrogen, it was found that there are two maxima, one at very low masses (e.g., 0.114 amu, corresponding to the muonium entity) and one maximum at around 2 amu (corresponding to deuterium). With the knowledge of the VTST-hydrogen transfer rates and the corresponding tunneling corrections, various tunneling criteria were tested (e.g., Swain-Schaad exponents) and were shown to fail in this reaction in predicting the extent of tunneling. This finding adds another aspect in the ongoing "Tunneling-Enhancement by Enzymes" discussion, as the used tunneling criteria based on experimental reaction rates may fail to predict tunneling behavior correctly.

Journal Article↗

PubMed2004

The ground-state tunneling splitting of various carboxylic acid dimers.

Carboxylic acid dimers in gas phase reveal ground-state tunneling splittings due to a double proton transfer between the two subunits. In this study we apply a recently developed accurate semiclassical method to determine the ground-state tunneling splittings of eight different carboxylic acid derivative dimers (formic acid, benzoic acid, carbamic acid, fluoro formic acid, carbonic acid, glyoxylic acid, acrylic acid, and N,N-dimethyl carbamic acid) and their fully deuterated analogs. The calculated splittings range from 5.3e-4 to 0.13 cm(-1) (for the deuterated species from 2.8e-7 to 3.3e-4 cm(-1)), thus indicating a strong substituent dependence of the splitting, which varies by more than two orders of magnitude. One reason for differences in the splittings could be addressed to different barriers heights, which vary from 6.3 to 8.8 kcal/mol, due to different mesomeric stabilization of the various transition states. The calculated splittings were compared to available experimental data and good agreement was found. A correlation could be found between the tunneling splitting and the energy barrier of the double proton transfer, as the splitting increases with increased strength of the hydrogen bonds. From this correlation an empirical formula was derived, which allows the prediction of the ground-state tunneling splitting of carboxylic acid dimers at a very low cost and the tunneling splittings for parahalogen substituted benzoic acid dimers is predicted.

Journal Article↗

PubMed2003

Influence of backbone conformations of human carbonic anhydrase II on carbon dioxide hydration: hydration pathways and binding of bicarbonate.

In this study, the hydration of carbon dioxide and the formation of bicarbonate in human carbonic anhydrase II have been examined. From semiempirical QM/MM molecular dynamics studies, dominant conformations of the protein backbone, possibly contributing to the catalytic activity, have been isolated and further examined by means of density functional QM/MM methods. In agreement with experimental observations, a binding site for cyanate, which acts as an inhibitor, has been located, whereas for carbon dioxide, depending on the conformation of the protein environment, either a different binding site or no binding site has been found. In the latter case, carbon dioxide diffuses barrierless to the zinc-bound oxygen, and then a weakly bound bicarbonate complex is formed. The formed complex is characterized by a long C-O bond to the zinc-bound hydroxide. The nature of the calculated stationary points was verified by determination of vibrational frequencies. Finally, the dissociation of the formed bicarbonate from zinc has been considered. Therefore, a water molecule was included in the QM zone of the QM/MM hybrid potential, and minimization yielded a pentacoordinated intermediate. From a potential energy scan, an activation energy of 6.2 kcal/mol for dissociation of bicarbonate from Zn has been found.

Bicarbonates↗

PubMed2003

Modeling anhydrous and aqua copper(II) amino acid complexes: a new molecular mechanics force field parametrization based on quantum chemical studies and experimental crystal data.

This paper presents the vacuum structures of aquacopper(II) bis(amino acid) complexes with glycine, sarcosine, N,N-dimethylglycine, and N-tert-butyl-N-methylglycine estimated using the B3LYP method. The differences between the B3LYP vacuum structures and experimental crystal structures suggested considerable influence of crystal lattice packing effects on the changes in the complexes' geometries. A previously developed molecular mechanics force field for modeling anhydrous copper(II) amino acidates was reoptimized to simulate these changes and predict the properties of both trans and cis anhydrous and aqua copper(II) amino acid complexes. The modeling included experimental molecular and crystal structures of 13 anhydrous and 10 aqua copper(II) amino acidates with the same atom types (Cu(II), C, H, N, and O) but various copper(II) coordination polyhedron geometries, crystal symmetries, and intermolecular interactions. The empirical parameters of the selected potential energy functions were optimized on the B3LYP vacuum copper(II) coordination geometries of three anhydrous copper(II) amino acidates and on experimental crystalline internal coordinates and unit cell dimensions of six anhydrous and six aqua copper(II) amino acid complexes. The respective equilibrium structures were calculated in vacuo and in simulated crystalline environment. The efficacy of the final force field, FFW, was examined. The total root-mean-square deviations between the experimental and theoretical crystal values were 0.018 A in the bond lengths, 2.2 degrees in the valence angles, 5.5 degrees in the torsion angles, and 0.395 A in the unit cell lengths. FFW reproduced the unit cell volumes in the range from -8.1 to 9.6%. The means of Cu to axial water oxygen distances were 2.4 +/- 0.1 A (experiment) and 2.6 +/- 0.1 A (FFW). This paper describes the ability of the molecular mechanics model and FFW force field to simulate the flexibility of the metal coordination polyhedron. The new force field proved effective in predicting the most stable molecular conformation of copper(II) amino acidato systems in vacuo.

Journal Article↗

PubMed2003

Extended method for adiabatic mode reordering.

The task of vibrational mode reordering is very important for reaction valley studies and for the determination of small curvature tunneling effects. An extended algorithm for adiabatic mode reordering is presented. It is based on the method introduced by Konkoli et al. [J Comput Chem 1997, 18, 1282], which is shown to suffer from numerical problems in the region of frequency-crossings and avoided crossings. One improvement is the use of cubic splines for the interpolation of the projected matrix of force constants, which allows larger step sizes between the discrete points along the reaction path, where vibrational analysis is performed. The main improvement is the use of perturbation theory to resolve crossings and avoided crossings. Within this theoretical framework it becomes clear why the method of the maximal overlap between the normal modes cannot work properly, as eigenvectors associated with nearby eigenvalues tend to become "wobbly". Thus a perturbative procedure is designed that is used for all cases where two harmonic frequencies approach each other and the overlap of the associated normal modes is of no practical use. Advantages of the new procedure are the use of larger step sizes along the minimum energy path and the much more reliable resolution of mode-crossings and avoided crossings independent of the systems symmetry. In addition to that it is shown that one should be very cautious in all computational situations when working with eigenvectors associated with nearby eigenvalues.

Journal Article↗

PubMed2002

About the stability of sulfurous acid (H2SO3) and its dimer.

The characterization and isolation of sulfurous acid (H2SO3) have never been accomplished and thus still remain one of the greatest open challenges of inorganic chemistry. It is known that H2SO3 is thermodynamically unstable. In this study, however, we show that a Ci-symmetric dimer of sulfurous acid (H2SO3)2 is 3.5 kcal mol-1 more stable than its dissociation products SO2 and H2O at 77 K. Additionally, we have investigated the kinetic stability of the sulfurous acid monomer with respect to dissociation into SO2 and H2O and the kinetic isotope effect (KIE) on this reaction by transition-state theory. At 77 K, the half-life of H2SO3 is 15 x 10(9) years, but for the deuterated molecule (D2SO3) it increases to 7.9 x 10(26) years. At room temperature, the half-life of sulfurous acid is only 24 hours; however, a KIE of 3.2 x 10(4) increases it to a remarkable 90 years. Water is an efficient catalyst for the dissociation reaction since it reduces the reaction barrier tremendously. With the aid of two water molecules, one can observe a change in the reaction mechanism for sulfurous acid decomposition with increasing temperature. The most likely mechanism below 170 K is via an eight-membered transition-state ring; yet, above 170 K, a mechanism with a six-membered transition state ring becomes the predominant one. For deuterated sulfurous acid, this change in reaction mechanism can be observed at 120 K. Consequently, between 120 and 170 K, different predominant reaction mechanisms occur for the decomposition of normal and deuterated sulfurous acid when assisted by two water molecules. However, the much longer half-life of deuterated sulfurous acid and the stability of the sulfurous acid dimer at 77 K are encouraging for future synthesis and characterization under laboratory conditions.

Journal Article↗

PubMed2002

Towards the experimental decomposition rate of carbonic acid (H2CO3) in aqueous solution.

Dry carbonic acid has recently been shown to be kinetically stable even at room temperature. Addition of water molecules reduces this stability significantly, and the decomposition (H2CO3 + nH2O --> (n+1)H2O + CO2) is extremely accelerated for n = 1, 2, 3. By including two water molecules, a reaction rate that is a factor of 3000 below the experimental one (10 s(-1)) at room temperature was found. In order to further remove the gap between experiment and theory, we increased the number of water molecules involved to 3 and took into consideration different mechanisms for thorough elucidation of the reaction. A mechanism whereby the reaction proceedes via a six-membered transition state turns out to be the most efficient one over the whole examined temperature range. The determined reaction rates approach experimental values in aqueous solution reasonably well; most especially, a significant increase in the rates in comparison to the decomposition reaction with fewer water molecules is found. Further agreement with experiment is found in the kinetic isotope effects (KIE) for the deuterated species. For water-free carbonic acid, the KIE (i.e., kH2CO3/kD2CO3) for the decomposition reaction is predicted to be 220 at 300 K, whereas it amounts to 2.2-3.0 for the investigated mechanisms including three water molecules. This result is therefore reasonably close to the experimental value of 2 (at 300 K). These KIEs are in much better accordance with the experiment than the KIE for decomposition with fewer water entities.

Carbonic Acid↗