BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Johannes Kästner

Johannes Kästner contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
physics.chem-phastro-ph.GAcond-mat.mtrl-sciphysics.comp-phphysics.app-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
19works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

19 published item(s)

preprint2022arXiv

Hydrogen abstraction reactions in formic and thioformic acid isomers by hydrogen and deuterium atoms. Insights on isomerism and deuteration

The isomerism of molecules in the interstellar medium and the mechanisms behind it are essential questions in the chemistry of organic molecules in space. In particular, for the simple formic and thioformic acids, the low temperatures found in molecular clouds indicate that cis-trans isomerization in the gas-phase must be impeded. Reactions happening on top of interstellar dust grains may explain the isomer interconversion at low temperatures. We studied the isomerization processes of formic and thioformic acid susceptible to happen on the surface of interstellar dust grains and initiated by H abstraction reactions. Similarly, deuterium enrichment of the acids can occur by the same mechanism. Our objective is to shed light on both topics to increase our understanding of key precursors of organic molecules in space.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Status and direction of atom probe analysis of frozen liquids

Imaging of liquids and cryogenic biological materials by electron microscopy has been recently enabled by innovative approaches for specimen preparation and the fast development of optimised instruments for cryo-enabled electron microscopy (cryo-EM). Yet, Cryo-EM typically lacks advanced analytical capabilities, in particular for light elements. With the development of protocols for frozen wet specimen preparation, atom probe tomography (APT) could advantageously complement insights gained by cryo-EM. Here, we report on different approaches that have been recently proposed to enable the analysis of relatively large volumes of frozen liquids from either a flat substrate or the fractured surface of a wire. Both allowed for analysing water ice layers which are several microns thick consisting of pure water, pure heavy-water and aqueous solutions. We discuss the merits of both approaches, and prospects for further developments in this area. Preliminary results raise numerous questions, in part concerning the physics underpinning field evaporation. We discuss these aspects and lay out some of the challenges regarding the APT analysis of frozen liquids.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Atom Tunneling in the Hydroxylation Process of Taurine/$α$-Ketoglutarate Dioxygenase Identified by Quantum Mechanics/Molecular Mechanics Simulations

TauD 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 QM/MM 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 quantum tunneling of the hydrogen atom to increase the rate constant by a factor of 40 at 5\textcelsius{}. As a consequence a quite high KIE value of close to 60 is obtained, which is consistent with the experimental value.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Atom-Tunneling in Chemistry

Quantum mechanical tunneling of atoms is increasingly found to play an important role in many chemical transformations. Experimentally, atom-tunneling can be indirectly detected by temperature-independent rate constants at low temperature or by enhanced kinetic isotope effects. On the contrary, using computational investigations the influence of tunneling on the reaction rates can directly be monitored. The tunnel effect, for example, changes reaction paths and branching ratios, enables chemical reactions in an astrochemical environment that would be impossible by thermal transition, and influences biochemical processes.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Calculation of Reaction Rate Constants in the Canonical and Microcanonical Ensemble

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Comparison of classical reaction paths and tunneling paths studied with the semiclassical instanton theory

Atom tunneling in the hydrogen atom transfer reaction of the 2,4,6-tri-tert-butylphenyl radical to 3,5-di-tert-butylneophyl, which has a short but strongly curved reaction path, was investigated using instanton theory. We found the tunneling path to deviate qualitatively from the classical intrinsic reaction coordinate, the steepest-descent path in mass-weighted Cartesian coordinates. To perform that comparison, we implemented a new variant of the predictor-corrector algorithm for the calculation of the intrinsic reaction coordinate. We used the reaction force analysis method as a means to decompose the reaction barrier into structural and electronic components. Due to the narrow energy barrier, atom tunneling is important in the abovementioned reaction, even above room temperature. Our calculated rate constants between 350 K and 100 K agree well with experimental values. We found a H/D kinetic isotope effect of almost 106 at 100 K. Tunneling dominates the protium transfer below 400 K and the deuterium transfer below 300 K. We compared the lengths of the tunneling path and the classical path for the hydrogen atom transfer in the reaction HCl + Cl and quantified the corner cutting in this reaction. At low temperature, the tunneling path is about 40% shorter than the classical path.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Dual-Level Approach to Instanton Theory

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 2$^{\text{nd}}$ 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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Efficient Training of ANN Potentials by Including Atomic Forces via Taylor Expansion and Application to Water and a Transition-Metal Oxide

Artificial neural network (ANN) potentials enable the efficient large-scale atomistic modeling of complex materials with near first-principles accuracy. For molecular dynamics simulations, accurate energies and interatomic forces are a prerequisite, but training ANN potentials simultaneously on energies and forces from electronic structure calculations is computationally demanding. Here, we introduce an efficient alternative method for the training of ANN potentials on energy and force information based on an extrapolation of the total energy via a Taylor expansion. By translating the force information to approximate energies, the quadratic scaling with the number of atoms exhibited by conventional force-training methods can be avoided, which enables the training on reference data sets containing complex atomic structures. We demonstrate for different materials systems, clusters of water molecules, bulk liquid water, and a lithium transition-metal oxide, that the proposed force-training approach provides substantial improvements over schemes that train on energies only. Including force information for training reduces the size of the reference data sets required for ANN potential construction, increases the transferability of the potential, and generally improves the force prediction accuracy. For a set of water clusters, the Taylor-expansion approach achieves around 50% of the force error improvement compared to the explicit training on all force components, at a much smaller computational cost. The alternative force training approach thus simplifies the construction of general ANN potentials for the prediction of accurate energies and interatomic forces for diverse types of materials, as demonstrated here for water and a transition-metal oxide.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Gaussian Process Regression for Geometry Optimization

We implemented a geometry optimizer based on Gaussian process regression (GPR) to find minimum structures on potential energy surfaces. We tested both a two times differentiable form of the Matérn kernel and the squared exponential kernel. The Matérn kernel performs much better. We give a detailed description of the optimization procedures. These include overshooting the step resulting from GPR in order to obtain a higher degree of interpolation vs. extrapolation. In a benchmark against the L-BFGS optimizer of the DL-FIND library on 26 test systems, we found the new optimizer to generally reduce the number of required optimization steps.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Gaussian Process Regression for Transition State Search

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Instanton Rate Constant Calculations Close to and Above the Crossover Temperature

Canonical instanton theory is known to overestimate the rate constant close to a system-dependent crossover temperature and is inapplicable above that temperature. We compare the accuracy of the reaction rate constants calculated using recent semi-classical rate expressions to those from canonical instanton theory. We show that rate constants calculated purely from solving the stability matrix for the action in degrees of freedom orthogonal to the instanton path is not applicable at arbitrarily low temperatures and employ two methods to overcome this. Furthermore, as a by-product of the developed methods, we derive a simple correction to canonical instanton theory that can alleviate this known overestimation of rate constants close to the crossover temperature. The combined methods accurately reproduce the rate constants of the canonical theory along the whole temperature range without the spurious overestimation near the crossover temperature. We calculate and compare rate constants on three different reactions: H in the Müller--Brown potential, methylhydroxycarbene $\to$ acetaldehyde and H$_2$ + OH $\to$ H + H$_2$O.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Instanton rate constant calculations using interpolated potential energy surfaces in non-redundant, rotationally and translationally invariant coordinates

A trivial flaw in the utilization of artificial neural networks in interpolating chemical potential energy surfaces (PES) whose descriptors are Cartesian coordinates is their dependence on simple translations and rotations of the molecule under consideration. A different set of descriptors can be chosen to circumvent this problem, internuclear distances, inverse internuclear distances or z-matrix coordinates are three such descriptors. The objective is to use an interpolated PES in instanton rate constant calculations, hence information on the energy, gradient and Hessian is required at coordinates in the vicinity of the tunneling path. Instanton theory relies on smoothly fitted Hessians, therefore we use energy, gradients and Hessians in the training procedure. A major challenge is presented in the proper back-transformation of the output gradients and Hessians from internal coordinates to Cartesian coordinates. We perform comparisons between our method, a previous approach and on-the-fly rate constant calcuations on the hydrogen abstraction from methanol and on the hydrogen addition to isocyanic acid.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Isomers in Interstellar Environments (I): The Case of Z- and E-Cyanomethanimine

In this work, we present the results of our investigation into the chemistry of Z- and E-cyanomethanimine (HNCHCN), both of which are possible precursors to the nucleobase adenine. Ab initio quantum chemical calculations for a number of reactions with atomic hydrogen were carried out. We find that the reaction H + Z/E-HNCHCN leading both to H-addition as well as H$_2$-abstraction proceed via similar short-range barriers with bimolecular rate coefficients on the order of $\sim10^{-17}$ cm$^{3}$ s$^{-1}$. These results were then incorporated into astrochemical models and used in simulations of the giant molecular cloud G+0.693. The calculated abundances obtained from these models were compared with previous observational data and found to be in good agreement, with a predicted [Z/E] ratio of $\sim3$ - somewhat smaller than the previously derived value of $6.1\pm2.4$. We find that the [Z/E] ratio in our simulations is due mostly to ion-molecule destruction rates driven by the different permanent dipoles of the two conformers. Based on these results, we propose a general rule-of-thumb for estimating the abundances of isomers in interstellar environments.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Low-Temperature Kinetic Isotope Effects in CH3OH+H -> CH2OH+H2 Shed Light on the Deuteration of Methanol in Space

We calculated reaction rate constants including atom tunneling for the hydrogen abstraction reaction CH3OH+H -> CH2OH+H2 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 CH3 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 about 2--6 for secondary KIEs. The $^\text{12}$C/$^\text{13}$C kinetic isotope effect is 1.08 at 30 K, while the $^\text{16}$O/$^\text{18}$O kinetic isotope effect is vanishingly small. A simple kinetic surface model using these data predicts high abundances of the deuterated forms of methanol.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Neural-Network Assisted Study of Nitrogen Atom Dynamics on Amorphous Solid Water. I. Adsorption & Desorption

Dynamics of adsorption and desorption of (4S)-N on amorphous solid water are analyzed using molecular dynamics simulations. The underlying potential energy surface was provided by machine-learned interatomic potentials. Binding energies confirm the latest available theoretical and experimental results. The nitrogen sticking coefficient is close to unity at dust temperatures of 10 K but decreases at higher temperatures. We estimate a desorption time scale of 1 μs at 28 K. The estimated time scale allows chemical processes mediated by diffusion to happen before desorption, even at higher temperatures. We found that the energy dissipation process after a sticking event happens on the picosecond timescale at dust temperatures of 10 K, even for high energies of the incoming adsorbate. Our approach allows the simulation of large systems for reasonable time scales at an affordable computational cost and ab-initio accuracy. Moreover, it is generally applicable for the study of adsorption dynamics of interstellar radicals on dust surfaces.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Potential energy surface interpolation with neural networks for instanton rate calculations

Artificial neural networks are used to fit a potential energy surface. We demonstrate the benefits of using not only energies, but also their first and second derivatives as training data for the neural network. This ensures smooth and accurate Hessian surfaces, which are required for rate constant calculations using instanton theory. Our aim was a local, accurate fit rather than a global PES, because instanton theory requires information on the potential only in the close vicinity of the main tunneling path. Elongations along vibrational normal modes at the transition state are used as coordinates for the neural network. The method is applied to the hydrogen abstraction reaction from methanol, calculated on a coupled-cluster level of theory. The reaction is essential in astrochemistry to explain the deuteration of methanol in the interstellar medium.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Rate constants from instanton theory via a microcanonical approach

Microcanonical instanton theory offers the promise of providing rate constants for chemical reactions including quantum tunneling of atoms over the whole temperature range. We discuss different rate expressions, which require the calculation of stability parameters of the instantons. The traditional way of obtaining these stability parameters is shown to be numerically unstable in practical applications. We provide three alternative algorithms to obtain such stability parameters for non-separable systems, i.e., systems in which the vibrational modes perpendicular to the instanton path couple to movement along the path. We show the applicability of our algorithms on two molecular systems: H$_2$ + OH $\rightarrow$ H$_2$O + H using a fitted potential energy surface and HNCO + H $\rightarrow$ NH$_2$CO using a potential obtained on-the-fly from density functional calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Tunneling Rate Constants for H2CO+H on Amorphous Solid Water Surfaces

Formaldehyde (H2CO) is one of the most abundant molecules observed in the icy mantle covering interstellar grains. Studying its evolution can contribute to our understanding of the formation of complex organic molecules in various interstellar environments. In this work, we investigated the hydrogenation reactions of H2CO yielding CH3O, CH2OH, and the hydrogen abstraction resulting in H2+HCO on an amorphous solid water (ASW) surface using a quantum mechanics/molecular mechanics (QM/MM) model. The binding energies of H2CO on the ASW surface vary broadly, from 1000 to 9370 K. No correlation was found between binding energies and activation energies of hydrogenation reactions. Combining instanton theory with QM/MM modeling, we calculated rate constants for the Langmuir Hinshelwood and the Eley Rideal mechanisms for the three product channels of H+H2CO surface reactions down to 59 K. We found that the channel producing CH2OH can be ignored, owing to its high activation barrier leading to significantly lower rates than the other two channels. The ASW surface influences the reactivity in favor of formation of CH3O (branching ratio 80%) and hinders the H2CO dissociation into H2+HCO. In addition, kinetic isotope effects are strong in all reaction channels and vary strongly between the channels. Finally, we provide fits of the rate constants to be used in astrochemical models.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Vibrational analysis of methyl cation - rare gas atom complexes: CH$_3^+$-Rg (Rg=He, Ne, Ar, Kr)

The vibrational spectra of simple CH$_3^+$-Rg (Rg=He, Ne, Ar, Kr) complexes have been studied by vibrational configuration interaction (VCI) theory relying on multidimensional potential energy surfaces (PES) obtained from explicitly correlated coupled cluster calculations, CCSD(T)-F12a. In agreement with experimental results, the series of rare gas atoms leads to rather unsystematic results and indicates huge zero point vibrational energy effects for the helium complex. In order to study these sensitive complexes more consistently, we also introduce configuration averaged vibrational self-consistent field theory (CAVSCF), which is a generalization of standard VSCF theory to several configurations. The vibrational spectra of the complexes are compared to that of the methyl cation, for which corrections due to scalar-relativistic effects, high-level coupled-cluster terms, i.e. CCSDTQ, and core-valence correlation have explicitly been accounted for. The occurrence of tunneling splittings for the vibrational ground-state of CH$_3^+$-He has been investigated on the basis of semiclassical instanton theory. These calculations and a direct comparison of the energy profiles along the intrinsic reaction coordinates (IRC) with that of the hydronium cation, H$_3$O$^+$, suggest that tunneling effects for vibrationally excited states should be very small.

0 citations0 reviewsNo code linkedOpen access