Search

The hydrogen abstraction reaction between H and H2S, yielding HS and H2 as products, has been studied within the framework of interstellar surface chemistry. High-temperature rate constants below 2000 K are calculated in the gas phase and are in agreement with previously reported values. Subsequently, low-temperature rate constants down to 55 K are presented for the first time that are of interest to astrochemistry, i.e., covering both bimolecular and unimolecular reaction mechanisms. For this, a so-called implicit surface model is used. In strict terms, this is a structural gas-phase model in which the restriction of the rotation in the solid state is taken into account. The calculated kinetic isotope effects are explained in terms of the difference in activation and delocalization. All rate constants are calculated at the UCCSD(T)-F12/cc-VTZ-F12 level of theory. Finally, we show that the energetics of the reaction is affected to an only small extent by the presence of H2O or H2S molecular clusters that simulate an ice surface, calculated at the MPWB1K/def2-TZVP level of theory. [ABSTRACT FROM AUTHOR]

Taurine/α-ketoglutarate dioxygenase is one of the most studied α-ketoglutarate-dependent dioxygenases (αKGDs), involved in several biotechnological applications. We investigated the key step in the catalytic cycle of the αKGDs, the hydrogen transfer process, by a quantum mechanics/molecular mechanics approach (B3LYP/CHARMM22). Analysis of the charge and spin densities during the reaction demonstrates that a concerted mechanism takes place, where the H atom transfer happens simultaneously with the electron transfer from taurine to the FeO cofactor. We found the quantum tunneling of the hydrogen atom to increase the rate constant by a factor of 40 at 5 °C. As a consequence, a quite high kinetic isotope effect close to 60 is obtained, which is consistent with the experimental value. [ABSTRACT FROM AUTHOR]

The radical mechanism of the conversion of glutamate to methylaspartate catalyzed by glutamate mutase is studied with quantum mechanical/molecular mechanical (QM/MM) simulations based on density functional theory (DFT/MM). The hydrogen transfer between the substrate and the cofactor is found to be rate limiting with a barrier of 101.1 kJ mol-1. A careful comparison to the uncatalyzed reaction in water is performed. The protein influences the reaction predominantly electrostatically and to a lesser degree sterically. Our calculations shed light on the atomistic details of the reaction mechanism. The well-known arginine claw and Glu 171 (Clostridium cochlearium notation) are found to have the strongest influence on the reaction. However, a catalytic role of Glu 214, Lys 322, Gln 147, Glu 330, Lys 326, and Met 294 is found as well. The arginine claw keeps the intermediates in place and is probably responsible for the enantioselectivity. Glu 171 temporarily accepts a proton from the glutamyl radical intermediate and donates it back at the end of the reaction. We relate our results to experimental data when available. Our simulations lead to further understanding of how glutamate mutase catalyzes the carbon skeleton rearrangement of glutamate. [ABSTRACT FROM AUTHOR]

Biological nitrogen fixation has been investigated beginning with the monoprotonated dinitrogen bound to the FeMo cofactor of nitrogenase up to the formation of the two ammonia molecules. The energy differences of the relevant intermediates, the reaction barriers, and potentially relevant side branches are presented. During the catalytic conversion, nitrogen bridges two Fe atoms of the central cage, replacing a sulfur bridge present before dinitrogen binds to the cofactor. A transformation from cis- to trans-diazene has been found. The strongly exothermic cleavage of the dinitrogen bond takes place, while the Fe atoms are bridged by a single nitrogen atom. The dissociation of the second ammonia from the cofactor is facilitated by the closing of the sulfur bridge following an intramolecular proton transfer. This closes the catalytic cycle. [ABSTRACT FROM AUTHOR]

The catalytic cycle of acetylene reduction at the FeMo cofactor of nitrogenase has been investigated on the basis of density functional theory. C2H2 binds to the same site as N2, but it binds to a less reduced state of the cofactor. In a manner similar to that of N2 binding, one of the sulfur bridges opens during acetylene binding. The model explains the strong noncompetitive inhibition of N2 reduction by C2H2 and the weak competitive inhibition of C2H2 reduction by N2. Our proposed mechanism is consistent with experimentally observed stereoselectivity and the ability of C2H2 to suppress H2 production by nitrogenase. [ABSTRACT FROM AUTHOR]

Organometallic complexes are frequently deposited on solid surfaces, but little is known about how the resulting complex-solid interactions alter their properties. Here, a series of complexes of the type Cu(dppf)(L x ) + (dppf = 1,1'-bis(diphenylphosphino)ferrocene, L x = mono- and bidentate ligands) were synthesized, physisorbed, ion-exchanged, or covalently immobilized on solid surfaces and investigated by 31 P MAS NMR spectroscopy. Complexes adsorbed on silica interacted weakly and were stable, while adsorption on acidic γ-Al 2 O 3 resulted in slow complex decomposition. Ion exchange into mesoporous Na-[Al]SBA-15 resulted in magnetic inequivalence of 31 P nuclei verified by 31 P- 31 P RFDR and 1 H- 31 P FSLG HETCOR. DFT calculations verified that a MeCN ligand dissociates upon ion exchange. Covalent immobilization via organic linkers as well as ion exchange with bidentate ligands both lead to rigidly bound complexes that cause broad 31 P CSA tensors. We thus demonstrate how the interactions between complexes and functional surfaces determine and alter the stability of complexes. The applied Cu(dppf)(L x ) + complex family members are identified as suitable solid-state NMR probes for investigating the influence of support surfaces on deposited inorganic complexes.

Recently, searches were made for H 2 CCS and HCCSH in a variety of interstellar environments─all of them resulted in nondetections of these two species. Recent findings have indicated the importance of destruction pathways, e.g., with atomic hydrogen, in explaining the consistent nondetection of other species, such as the H 2 C 3 O family of isomers. We have thus performed ab initio calculations looking at reactions of H 2 CCS, HCCSH, and related species with atomic hydrogen. Our results show that H 2 CCS and HCCSH are both destroyed barrierlessly by atomic hydrogen, thus providing a plausible explanation for the nondetections. We further find that subsequent reactions with atomic hydrogen can barrierlessly lead to CH 3 CH 2 SH, which has been detected. Astrochemical simulations including these reactions result not only in reproducing the observed abundance of H 2 CCS in TMC-1 but also show that CH 3 CH 2 SH, produced via our H-addition pathways and subsequently trapped on grains, can desorb in warmer sources up to abundances that match previous observations of CH 3 CH 2 SH in Orion KL. These results, taken together, point to the importance of grain-surface H-atom addition reactions and highlight the chemical links between cold prestellar cores and their subsequent, warmer evolutionary stages.

We report new computational and experimental evidence of an efficient and astrochemically relevant formation route to formaldehyde (H 2 CO). This simplest carbonylic compound is central to the formation of complex organics in cold interstellar clouds and is generally regarded to be formed by the hydrogenation of solid-state carbon monoxide. We demonstrate H 2 CO formation via the reaction of carbon atoms with amorphous solid water. Crucial to our proposed mechanism is a concerted proton transfer catalyzed by the water hydrogen bonding network. Consequently, the reactions 3 C + H 2 O → 3 HCOH and 1 HCOH → 1 H 2 CO can take place with low or without barriers, contrary to the high-barrier traditional internal hydrogen migration. These low barriers (or the absence thereof) explain the very small kinetic isotope effect in our experiments when comparing the formation of H 2 CO to D 2 CO. Our results reconcile the disagreement found in the literature on the reaction route C + H 2 O → H 2 CO.

Artificial neural networks (NNs) are one of the most frequently used machine learning approaches to construct interatomic potentials and enable efficient large-scale atomistic simulations with almost ab initio accuracy. However, the simultaneous training of NNs on energies and forces, which are a prerequisite for, e.g., molecular dynamics simulations, can be demanding. In this work, we present an improved NN architecture based on the previous GM-NN model [Zaverkin V.; Kästner, J. J. Chem. Theory Comput . 2020, 16, 5410-5421], which shows an improved prediction accuracy and considerably reduced training times. Moreover, we extend the applicability of Gaussian moment-based interatomic potentials to periodic systems and demonstrate the overall excellent transferability and robustness of the respective models. The fast training by the improved methodology is a prerequisite for training-heavy workflows such as active learning or learning-on-the-fly.

Covalent organic frameworks (COFs) with a pore size beyond 5 nm are still rarely seen in this emerging field. Besides obvious complications such as the elaborated synthesis of large linkers with sufficient solubility, more subtle challenges regarding large-pore COF synthesis, including pore occlusion and collapse, prevail. Here we present two isoreticular series of large-pore imine COFs with pore sizes up to 5.8 nm and correlate the interlayer interactions with the structure and thermal behavior of the COFs. By adjusting interlayer interactions through the incorporation of methoxy groups acting as pore-directing "anchors", different stacking modes can be accessed, resulting in modified stacking polytypes and, hence, effective pore sizes. A strong correlation between stacking energy toward highly ordered, nearly eclipsed structures, higher structural integrity during thermal stress, and a novel, thermally induced phase transition of stacking modes in COFs was found, which sheds light on viable design strategies for increased structural control and stability in large-pore COFs.

Geometry optimization based on Gaussian process regression (GPR) was extended to internal coordinates. We used delocalized internal coordinates composed of distances and several types of angles and compared two methods of including them. In both cases, the GPR surrogate surface is trained on geometries in internal coordinates. In one case, it predicts the gradient in Cartesian coordinates and in the other, in internal coordinates. We tested both methods on a set of 30 small molecules and one larger Rh complex taken from the study of a catalytic mechanism. The former method is slightly more efficient, while the latter method is somewhat more robust. Both methods reduce the number of required optimization steps compared to GPR in Cartesian coordinates or the standard L-BFGS optimizer. We found it advantageous to use automatically adjusted hyperparameters to optimize them.

Atom probe tomography allows us to measure the three-dimensional composition of materials with up to atomic resolution by evaporating the material using high electric fields. Initially developed for metals, it is increasingly used for covalently bound structures. To aid the interpretation of the obtained fragmentation pattern, we modeled the fragmentation and desorption of self-assembled monolayers of thiolate molecules on a gold surface in strong electrostatic fields using density functional theory. We used a cluster model and a periodic model of amino-undecanethiolate, NH 2 (CH 2 ) 11 S, and fluoro-decanethiolate, CF 3 (CF 2 ) 7 (CH 2 ) 2 S. In the former molecule, the fragment CH 2 NH 2 + was found to evaporate at fields of 5.4-7.7 V/nm. It was followed by different hydrocarbon fragments. Fluoro-decanethiolate evaporates CF 3 + at fields of 5.7-6.7 V/nm in the cluster model and at 15.4-23.1 V/nm in the periodic model, followed by CF 2 + and C 2 F 4 2+ . Detailed analysis of the electronic structure during the evaporation process revealed a stepwise accumulation of the charge in the head groups exposed to the strongest fields, followed by dissociation of covalent bonds. These observations will facilitate the analysis of atom probe experiments of covalently bound structures.

We calculated reaction rate constants including atom tunneling for the hydrogen abstraction reaction CH 3 OH + H → CH 2 OH + H 2 with the instanton method. The potential energy was fitted by a neural network that was trained to UCCSD(T)-F12/VTZ-F12 data. Bimolecular gas-phase rate constants were calculated using microcanonic instanton theory. All H/D isotope patterns on the CH 3 group and the incoming H atom are studied. Unimolecular reaction rate constants, representing the reaction on a surface, down to 30 K, are presented for all isotope patterns. At 30 K, they range from 4100 for the replacement of the abstracted H by D to ∼8 for the replacement of the abstracting H to ∼2 to 6 for secondary KIEs. The 12 C/ 13 C kinetic isotope effect is 1.08 at 30 K, while the 16 O/ 18 O kinetic isotope effect is extremely small. A simple kinetic surface model using these data predicts high abundances of the deuterated forms of methanol.

Canonical instanton theory is a widespread approach to describe the dynamics of chemical reactions in low temperature environments when tunneling effects become dominant. It is a semiclassical theory which requires locating classical periodic orbits on the upside-down potential energy surface, so-called instantons, and the computation of second order quantum corrections. The calculation of these corrections usually involves a matrix diagonalization. In this paper we present an alternative approach, which requires to solve only linear systems of equations involving sparse matrices. Furthermore, the proposed method provides a reliable and numerically stable way to obtain stability parameters in multidimensional systems, which are of particular interest in the context of microcanonical instanton theory.

We implemented a gradient-based algorithm for transition state search which uses Gaussian process regression. Besides a description of the algorithm, we provide a method to find the starting point for the optimization if only the reactant and product minima are known. We perform benchmarks on 27 test systems against the dimer method and partitioned rational function optimization as implemented in the DL-FIND library. We found the new optimizer to significantly decrease the number of required energy and gradient evaluations.

Instanton theory is an established method to calculate rate constants of chemical reactions including atom tunneling. Technical and methodological improvements increased its applicability. Still, a large number of energy and gradient calculations is necessary to optimize the instanton tunneling path, and second derivatives of the potential energy along the tunneling path have to be evaluated, restricting the range of suitable electronic structure methods. To enhance the applicability of instanton theory, we present a dual-level approach in which instanton optimizations and Hessian calculations are performed using an efficient but approximate electronic structure method, and the potential energy along the tunneling path is recalculated using a more accurate method. This procedure extends the applicability of instanton theory to high-level electronic structure methods for which analytic gradients may not be available, like local linear-scaling approaches. We demonstrate for the analytical Eckart barrier and three molecular systems how the dual-level instanton approach corrects for the largest part of the error caused by the inaccuracy of the efficient electronic structure method. This reduces the error of the calculated rate constants significantly.

The hydrogen abstraction reaction between H and H 2 S, yielding HS and H 2 as products, has been studied within the framework of interstellar surface chemistry. High-temperature rate constants below 2000 K are calculated in the gas phase and are in agreement with previously reported values. Subsequently, low-temperature rate constants down to 55 K are presented for the first time that are of interest to astrochemistry, i.e., covering both bimolecular and unimolecular reaction mechanisms. For this, a so-called implicit surface model is used. In strict terms, this is a structural gas-phase model in which the restriction of the rotation in the solid state is taken into account. The calculated kinetic isotope effects are explained in terms of the difference in activation and delocalization. All rate constants are calculated at the UCCSD(T)-F12/cc-VTZ-F12 level of theory. Finally, we show that the energetics of the reaction is affected to an only small extent by the presence of H 2 O or H 2 S molecular clusters that simulate an ice surface, calculated at the MPWB1K/def2-TZVP level of theory.

The flavin-containing enzyme monoamine oxidase (MAO) is essential for the enzymatic decomposition of amine neurotransmitters. The exact mechanism of the oxidative deamination of amines to aldehydes by the enzyme has not yet been fully understood despite extensive research on the area. The rate limiting step is the reductive half-reaction where the Hα together with two electrons of the amine substrate is transferred to the flavin cofactor. However, it is still not known whether the hydrogen is transferred as a proton or a hydride. Experimental results cannot be fully explained by either of those mechanisms. In our previous work, theoretical results based on QM/MM calculations of the full enzyme show an intermediate situation between these two cases. In this paper, we report on an in-depth computational analysis concerning the role of the enzymatic environment for the reaction mechanism of human MAO-B with different p-substituted benzylamines as substrates. Our results show that steric and electrostatic effects from the active site environment turn the mechanism closer to an asynchronous polar nucleophilic mechanism. We found indications that the protein environment of MAO-A enhances the polar nucleophilic character of the mechanism compared to that of MAO-B.

Quantum mechanical tunneling of atoms plays a significant role in many chemical reactions. The crossover temperature between classical and quantum movement is a convenient preliminary indication of the importance of tunneling for a particular reaction. Here we show, using instanton theory, that quantum tunneling is possible significantly above this crossover temperature for specific forms of the potential energy surface. We demonstrate the effect on an analytic potential as well as a chemical system. While protons move asynchronously along a Grotthuss chain in the classical high-temperature range, the onset of tunneling results in a synchronization of their movement.

The flavoenzyme monoamine oxidase (MAO) is essential for the enzymatic decomposition of neurotransmitters. While it is commonly accepted that the rate limiting step of the reaction is the stereoselective abstraction of a hydrogen from the substrate, the precise mechanism is unknown. We modeled the reaction of human MAO-B with benzylamine by means of QM/MM calculations based on density functional theory. Oxidation of the unprotonated substrate was found to proceed with rates in good agreement with experimental values, while the protonated substrate does not react at room temperature. Our results support a concerted asynchronous polar nucleophilic mechanism. The lone pair of the amine-nitrogen interacts with a carbon atom of the flavin cofactor. During the reaction, this lone pair, as well as a proton, are transferred to the cofactor. Analysis of the electronic structure during the reaction rules out a radical mechanism.

Nudged elastic band (NEB) is a well established chain-of-states method to locate the minimum energy path in configuration space. Unfortunately, existing search algorithms suffer from slow convergence. We provide an analytic derivative of the nudged elastic band force, enabling a full Newton-Raphson optimization. For molecular systems, the components of the step belonging to translations and rotations are removed with an efficient algorithm. Minimization of the NEB force is ensured by reversing components for which the Newton-Raphson step would increase the force. We achieve quadratic convergence of this optimizer when applied to simple test cases where analytic Hessians are available: one analytic two-dimensional potential and a system of Lennard-Jones particles.

The proper description of explicit water shells is of enormous importance for all-atom calculations. We propose a new approach for the setup of water shells around proteins based on Tabu-Search global optimization and compare its efficiency with standard molecular dynamics protocols using the chignolin protein as a test case. Both algorithms generate reasonable water shells, but the new approach provides solvated systems with an increased water-enzyme interaction and offers further advantages. It enables a stepwise buildup of the solvent shell, so that the more important inner part can be prepared more carefully. It also allows the generation of solute structures which can be biased either toward the (experimental) starting structure or the underlying theoretical model, i.e., the employed force field.

The role of quantum mechanical atom tunneling during the conversion of glutamate to methylaspartate catalyzed by glutamate mutase is investigated by quantum mechanical/molecular mechanical (QM/MM) simulations based on coupled cluster and density functional calculations. The use of instanton theory allows us to calculate the tunneling contributions of up to 78 atoms in the active site. We calculate kinetic isotope effects (KIEs) and compare them to experimental data. The simulations lead to deuterium KIEs of 10 for the hydrogen abstraction from glutamate substrate and 16 for the hydrogen abstraction from methylaspartate substrate, which are consistent with the experimental results. The hydrogen abstraction from methylaspartate has higher primary deuterium and tritium (46.1) KIEs than the abstraction from glutamate. The tunneling effect increases the reaction rate by a factor of 12.3 for the hydrogen abstraction from methylaspartate at 0. Tunneling is supported by the environment by preparing the enzyme through classical motions. Consideraton of the tunneling contributions of more and more atoms around the active center shows that the motions at the ribose ring play a central role during the tunneling enhancement of the hydrogen transfers. Our simulations give new insight into the catalytic process in glutamate mutase and the way enzymes use tunneling effects for a successful catalysis.

We calculate, down to low temperature and for different isotopes, the reaction rate constants for the hydrogen abstraction reaction H + H(3)COH → H(2) + CH(2)OH/CH(3)O. These explain the known abundances of deuterated forms of methanol in interstellar clouds, where CH(2)DOH can be almost as abundant as CH(3)OH. For abstraction from both the C- and the O-end of methanol, the barrier-crossing motion involves the movement of light hydrogen atoms. Consequently, tunneling plays a dominant role already at relatively high temperature. Our implementation of harmonic quantum transition state theory with on the fly calculation of forces and energies accounts for these tunneling effects. The results are in good agreement with previous semiclassical and quantum dynamics calculations (down to 200 K) and experimental studies (down to 295 K). Here we extend the rate calculations down to lower temperature: 30 K for abstraction from the C-end of methanol and 80 K for abstraction from the OH-group. At all temperatures, abstraction from the C-end is preferred over abstraction from the O-end, more strongly so at lower temperature. Furthermore, the tunneling behavior strongly affects the kinetic isotope effects (KIEs). D + H(3)COH → HD + CH(2)OH has a lower vibrationally adiabatic barrier than H + H(3)COH → H(2) + CH(2)OH, giving rise to an inverse KIE (k(H)/k(D) < 1) at high temperature, in accordance with previous experiments and calculations. However, since tunneling is more facile for the light H atom, abstraction by H is favored over abstraction by D below ~135 K, with a KIE k(H)/k(D) of 11.2 at 30 K. The H + D(3)COD → HD + CD(2)OD reaction is calculated to be much slower than the D + H(3)COH → HD + CH(2)OH, in agreement with low-temperature solid-state experiments, which suggests the preference for H (as opposed to D) abstraction from the C-end of methanol to be the mechanism by which interstellar methanol is deuterium-enriched.

We implemented and compared four algorithms to locate instantons, i.e., the most likely tunneling paths at a given temperature. These allow to calculate reaction rates, including atom tunneling, down to very low temperature. An instanton is a first-order saddle point of the Euclidean action in the space of closed Feynman paths. We compared the Newton-Raphson method to the partitioned rational function optimization (P-RFO) algorithm, the dimer method, and a newly proposed mode-following algorithm, where the unstable mode is directly estimated from the instanton path. We tested the algorithms on three chemical systems, each including a hydrogen transfer, at different temperatures. Overall, the Newton-Raphson turned out to be the most promising method, with our newly proposed mode following, being the fall-back option.