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

Markus Meuwly

Markus Meuwly contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
physics.chem-phBiological PhysicsBiomoleculesphysics.ao-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
20works
0followers
4topics
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

20 published item(s)

preprint2026arXiv

Full Reaction Pathway Dynamics for Atmospheric Decomposition Reactions: The Photodissociation of H$_2$COO

Branching ratios for fragmentation channels of important meta- and unstable species are essential for a molecular-level characterization of atmospheric chemistry. Here, the molecular product channels for the decomposition dynamics of the smallest Criegee intermediate, H$_2$COO, are quantitatively investigated. Using a high-quality, full-dimensional machine learned potential energy surface (CASPT2/aug-cc-pVTZ), the translational, rotational, and vibrational energy distributions of the CO$_2$+H$_2$, H$_2$O+CO, and HCO+OH fragmentation channels were analyzed to elucidate partitioning of the available energy. The CO$_2$ + H$_2$ product forms through two different pathways that bifurcate after formation of the OCH$_2$O intermediate. Along the direct pathway, CO$_2$ is preferentially vibrationally excited with H$_2$in its vibrational ground state, whereas for the indirect pathway going through formic acid, H$_2$ can populate levels with $v > 0$. For all product channels passing through energized formic acid, the lifetime distributions are described by stretched exponentials with $β$ ranging from 1.1 to 1.7. This is a clear signature of non-RRKM effects and suggests that the explicit molecular dynamics needs to be followed for a quantitative and realistic description of the photodissociation dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Computational Vibrational Spectroscopy

Vibrational spectroscopy is a powerful technique to characterize the near-equilibrium dynamics of molecules in the gas- and the condensed-phase. This contribution summarizes efforts from computer-based methods to gain insight into the relationship between structure and spectroscopic response. Methods for this purpose include physics-based empirical energy functions, machine-learned force fields, and methods that separate sampling conformational space and determining the data for spectral analysis such as map-based approaches.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Double Proton Transfer in Hydrated Formic Acid Dimer: Interplay of Spatial Symmetry and Solvent-Generated Force on Reactivity

The double proton transfer (DPT) reaction in hydrated formic acid dimer (FAD) is investigated at molecular-level detail. For this, a global and reactive machine learned (ML) potential energy surface (PES) is developed to run extensive (more than 100 ns) mixed ML/MM molecular dynamics (MD) simulations in explicit molecular mechanics (MM) solvent at MP2-quality for the solute. Simulations with fixed - as in a conventional empirical force field - and conformationally fluctuating - as available from the ML-based PES - charge models for FAD shows significant impact on the competition between DPT and dissociation of FAD into two formic acid monomers. With increasing temperature the barrier height for DPT in solution changes by about 10% ($\sim 1$ kcal/mol) between 300 K and 600 K. The rate for DPT is largest, $\sim 1$ ns$^{-1}$, at 350 K and decreases for higher temperatures due to destabilisation and increased probability for dissociation of FAD. The water solvent is found to promote the first proton transfer by exerting a favourable solvent-induced Coulomb force along the O-H$\cdots$O hydrogen bond whereas the second proton transfer is significantly controlled by the O-O separation and other conformational degrees of freedom. Double proton transfer in hydrated FAD is found to involve a subtle interplay and balance between structural and electrostatic factors.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Hydration Dynamics and IR Spectroscopy of 4-Fluorophenol

Halogenated groups are relevant in pharmaceutical applications and potentially useful spectroscopic probes for infrared spectroscopy. In this work, the structural dynamics and infrared spectroscopy of $para$-fluorophenol (F-PhOH) and phenol (PhOH) is investigated in the gas phase and in water using a combination of experiment and molecular dynamics (MD) simulations. The gas phase and solvent dynamics around F-PhOH and PhOH is characterized from atomistic simulations using empirical energy functions with point charges or multipoles for the electrostatics, Machine-Learning (ML) based parametrization and with full $\textit{ab initio}$ (QM) and mixed Quantum Mechanical/Molecular Mechanics (QM/MM) simulations with a particular focus on the CF- and OH-stretch region. The CF-stretch band is heavily mixed with other modes whereas the OH-stretch in solution displays a characteristic high-frequency peak around 3600 cm$^{-1}$ most likely associated with the -OH group of PhOH and F-PhOH together with a characteristic progression below 3000 cm$^{-1}$ due to coupling with water modes which is also reproduced by several of the simulations. Solvent and radial distribution functions indicate that the CF-site is largely hydrophobic except for simulations using point charges which renders them unsuited for correctly describing hydration and dynamics around fluorinated sites.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Local Hydration Control and Functional Implications Through S-Nitrosylation of Proteins: Kirsten rat sarcoma virus (KRAS) and Hemoglobin (Hb)

S-nitrosylation, the covalent addition of NO to the thiol side chain of cysteine, is an important post-transitional modification (PTM) that can affect the function of proteins. As such, PTMs extend and diversify protein functions and thus characterizing consequences of PTM at a molecular level is of great interest. Although PTMs can be detected through various direct/indirect methods, they lack the capabilities to investigate the modifications at the molecular level. In the present work local and global structural dynamics, their correlation, the hydration structure, and the infrared spectroscopy for WT and S-nitrosylated Kirsten rat sarcoma virus (KRAS) and Hemoglobin (Hb) are characterized from molecular dynamics simulations. It is found that for KRAS attaching NO to Cys118 rigidifies the protein in the Switch-I region which has functional implications, whereas for Hb nitrosylation at Cys93 at the $β_1$ chain increases the flexibility of secondary structural motives for Hb in its T$_{0}$ and R$_{4}$ conformational substates. Solvent water access decreased by 40\% after nitrosylation in KRAS, similar to Hb for which, however, local hydration of the R$_4$NO state is yet lower than for T$_0$NO. Finally, S-nitrosylation leads to detectable peaks for NO stretch, however, congested IR region will make experimental detection of these bands difficult.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Molecular Dynamics with Conformationally Dependent, Distributed Charges

Accounting for geometry-induced changes in the electronic distribution in molecular simulation is important for capturing effects such as charge flow, charge anisotropy and polarization. Multipolar force fields have demonstrated their ability to qualitatively and correctly represent chemically significant features such as sigma holes. It has also been shown that off-center point charges offer a compact alternative with similar accuracy. Here it is demonstrated that allowing relocation of charges within a minimally distributed charge model (MDCM) with respect to their reference atoms is a viable route to capture changes in the molecular charge distribution depending on geometry. The approach, referred to as ``flexible MDCM'' (fMDCM) is validated on a number of small molecules and provides accuracies in the electrostatic potential (ESP) of 0.5 kcal/mol on average compared with reference data from electronic structure calculations whereas MDCM and point charges have root mean squared errors of a factor of 2 to 5 higher. In addition, MD simulations in the $NVE$ ensemble using fMDCM for a box of flexible water molecules with periodic boundary conditions show a width of 0.1 kcal/mol for the fluctuation around the mean at 300 K on the 10 ns time scale. The accuracy in capturing the geometry dependence of the ESP together with the long-time stability in energy conserving simulations makes fMDCM a promising tool to introduce advanced electrostatics into atomistic simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Photodissociation dynamics of N$_{3}^{+}$

The photodissociation dynamics of N$_3^+$ excited from its $^3 Σ_{\rm g}^{-}$ ground to the first excited singlet and triplet states is investigated. Three dimensional potential energy surfaces for the $^1$A$'$, $^1$A$''$, and $^3$A$'$ electronic states, correlating with the $^1 Δ_{\rm g}$ and $^3 Π_{\rm u}$ states in linear geometry, for N$_3^+$ are constructed using high level electronic structure calculations and represented as reproducing kernels. The reference {\it ab initio} energies are calculated at the MRCI+Q/aug-cc-pVTZ level of theory. For following the photodissociation dynamics in the excited states, rotational and vibrational distributions $P(v')$ and $P(j')$ for the N$_2$ product are determined from vertically excited ground state distributions. Due to the different shapes of the ground state $^3$A$^{''}$ PES and the excited states, appreciable angular momentum $j' \sim 60$ is generated in the diatomic fragments. The lifetimes in the excited states extend to at least 50 ps. Notably, results from sampling initial conditions from a thermal ensemble and from the Wigner distribution of the ground state wavefunction are comparable.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Transfer learning for affordable and high quality tunneling splittings from instanton calculations

The combination of transfer learning (TL) a low level potential energy surface (PES) to a higher level of electronic structure theory together with ring-polymer instanton (RPI) theory is explored and applied to malonaldehyde. The RPI approach provides a semiclassical approximation of the tunneling splitting and depends sensitively on the accuracy of the PES. With second order Møller-Plesset perturbation theory (MP2) as the low-level (LL) model and energies and forces from coupled cluster singles, doubles and perturbative triples (CCSD(T)) as the high-level (HL) model, it is demonstrated that CCSD(T) information from only 25 to 50 judiciously selected structures along and around the instanton path suffice to reach HL-accuracy for the tunneling splitting. In addition, the global quality of the HL-PES is demonstrated through a mean average error of 0.3 kcal/mol for energies up to 40 kcal/mol above the minimum energy structure (a factor of 2 higher than the energies employed during TL) and $< 2 $ cm$^{-1}$ for harmonic frequencies compared with computationally challenging normal mode calculations at the CCSD(T) level.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Transfer Learned Potential Energy Surfaces: Accurate Anharmonic Vibrational Dynamics and Dissociation Energies for the Formic Acid Monomer and Dimer

The vibrational dynamics of formic acid monomer (FAM) and dimer (FAD) is investigated from machine-learned potential energy surfaces at the MP2 (PES$_{\rm MP2}$) and transfer-learned (PES$_{\rm TL}$) to the CCSD(T) levels of theory. The normal modes and anharmonic frequencies of all modes below 2000 cm$^{-1}$ agree favourably with experiment whereas the OH-stretch mode is challenging for FAM and FAD from normal mode analyses and finite-temperature MD simulations. VPT2 calculations on PES$_{\rm TL}$ for FAM reproduce the experimental OH frequency to within 22 cm$^{-1}$. For FAD the VPT2 calculations find the high-frequency OH stretch at 3011cm$^{-1}$, compared with an experimentally reported, broad ($\sim 100$ cm$^{-1}$) absorption band with center frequency estimated at $\sim 3050$ cm$^{-1}$. In agreement with earlier reports, MD simulations at higher temperature shift the position of the OH-stretch in FAM to the red, consistent with improved sampling of the anharmonic regions of the PES. However, for FAD the OH-stretch shifts to the blue and for temperatures higher than 1000 K the dimer partly or fully dissociates using PES$_{\rm TL}$. Including zero-point energy corrections from diffusion Monte Carlo simulations for FAM and FAD and corrections due to basis set superposition and completeness errors yield a dissociation energy of $D_0 = -14.23 \pm 0.08$ kcal/mol compared with an experimentally determined value of $-14.22 \pm 0.12$ kcal/mol.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Isomerization and Decomposition Reactions of Acetaldehyde Relevant to Atmospheric Processes from Dynamics Simulations on Neural Network-Based Potential Energy Surfaces

Acetaldehyde (AA) isomerization (to vinylalcohol, VA) and decomposition (into either CO+CH$_4$ and H$_2$+H$_2$CCO) is studied using a fully dimensional, reactive potential energy surface represented as a neural network (NN). The NN, trained on 432&#39;399 reference structures from MP2/aug-cc-pVTZ calculations has a MAE of 0.0453 kcal/mol and an RMSE of 1.186 kcal/mol for a test set of 27&#39;399 structures. For the isomerization process AA $\rightarrow$ VA the minimum dynamical path implies that the C-H vibration, and the C-C-H (with H being the transferring H-atom) and the C-C-O angles are involved to surmount the 68.2 kcal/mol barrier. Using an excess energy of 93.6 kcal/mol - the energy available in the solar spectrum and sufficient to excite to the first electronically excited state - to initialize the molecular dynamics, no isomerization to VA is observed on the 500 ns time scale. Only with excess energies of $\sim$ 127.6 kcal/mol (including the zero point energy of the AA molecule), isomerization occurs on the nanosecond time scale. Given that collisional de-excitation at atmospheric conditions in the stratosphere occurs on the 100 ns time scale, it is concluded that formation of VA following photoexcitation of AA from actinic photons is unlikely. This also limits the relevance of this reaction pathway to be a source for formic acid.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

N$_3^+$: Full-Dimensional Potential Energy Surface, Vibrational Energy Levels and Ground State Dynamics

The fundamentals and higher vibrationally excited states for the N$_3^+$ ion in its electronic ground state have been determined from quantum bound state calculations on 3-dimensional potential energy surfaces (PESs) at the CCSD(T)-F12 and MRCI+Q levels of theory. The vibrational fundamentals are at 1130 cm$^{-1}$ ($ν_1$, symmetric stretch), 807 cm$^{-1}$ ($ν_3$, asymmetric stretch), and 406 cm$^{-1}$ ($ν_2$, bend) on the higher-quality CCSD(T)-F12 surface. For $ν_1$, the calculations are close to the estimated frequency from experiment (1170 cm$^{-1}$) and previous calculations\cite{rosmus.n3:1994} which find it at 1190 cm$^{-1}$. Calculations of the vibrational states on the MRCI+Q PES are in qualitative agreement with those using the CCSD(T)-F12 PES. Analysis of the reference CASSCF wave function for the MRCI+Q calculations provides further insight into the shape of the PES and lends support for the reliability of Hartree-Fock as the reference wave function for the coupled cluster calculations. According to this, N$_3^+$ has mainly single reference character in all low-energy regions of its electronic ground state ($^3$A$&#39;&#39;$) 3d PES.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Reactive atomistic simulations of Diels-Alder-type reactions: Conformational and dynamic effects in the polar cycloaddition of 2,3-dibromobutadiene radical ions with maleic anhydride

The kinetics, dynamics and conformational specificities for the ionic Diels-Alder reaction (polar cycloaddition) of maleic anhydride with 2,3-dibromobutadiene radical ions have been studied theoretically using multisurface adiabatic reactive molecular dynamics. A competition of concerted and stepwise reaction pathways was found and both the s-cis and s-trans conformers of of the diene are reactive. The analysis of the minimum dynamic path of the reaction indicates that both, rotations and vibrations of the reactant molecules are important for driving the system towards the transition state. The rates were computed as $k = 5.1 \times 10^{-14}$ s$^{-1}$ for the s-cis and $k = 3.8 \times 10^{-14}$ s$^{-1}$ for the s-trans conformer of 2,3-dibromobutadiene at an internal temperature of 300 K. The present results are to be contrasted with the neutral variant of the title system in which only the gauche conformer of the diene was found to undergo a considerably slower, concerted and mostly synchronous reaction driven by the excitation of rotations. The results presented here inform detailed experimental studies of the dynamics of polar cycloadditions under single-collision conditions in the gas phase.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Dynamics and Infrared Spectrocopy of Monomeric and Dimeric Wild Type and Mutant Insulin

The infrared spectroscopy and dynamics of -CO labels in wild type and mutant insulin monomer and dimer are characterized from molecular dynamics simulations using validated force fields. It is found that the spectroscopy of monomeric and dimeric forms in the region of the amide-I vibration differs for residues B24-B26 and D24-D26, which are involved in dimerization of the hormone. Also, the spectroscopic signatures change for mutations at position B24 from phenylalanine - which is conserved in many organisms and known to play a central role in insulin aggregation - to alanine or glycine. Using three different methods to determine the frequency trajectories - solving the nuclear Schrödinger equation on an effective 1-dimensional potential energy curve, instantaneous normal modes, and using parametrized frequency maps - lead to the same overall conclusions. The spectroscopic response of monomeric WT and mutant insulin differs from that of their respective dimers and the spectroscopy of the two monomers in the dimer is also not identical. For the WT and F24A and F24G monomers spectroscopic shifts are found to be $\sim 20$ cm$^{-1}$ for residues (B24 to B26) located at the dimerization interface. Although the crystal structure of the dimer is that of a symmetric homodimer, dynamically the two monomers are not equivalent on the nanosecond time scale. Together with earlier work on the thermodynamic stability of the WT and the same mutants it is concluded that combining computational and experimental infrared spectroscopy provides a potentially powerful way to characterize the aggregation state and dimerization energy of modified insulins.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Thermal Activation of Methane by MgO$^+$: Temperature Dependent Kinetics, Reactive Molecular Dynamics Simulations and Statistical Modeling

The kinetics of MgO$^+$ + CH$_4$ was studied experimentally using the variable ion source, temperature adjustable selected ion flow tube (VISTA-SIFT) apparatus from 300 $-$ 600 K and computationally by running and analyzing reactive atomistic simulations. Rates and product branching fractions were determined as a function of temperature. The reaction proceeded with a rate of $k = 5.9 \pm 1.5 10^{-10}(T/300 $ K$)^{-0.5 \pm 0.2}$ cm$^3$ s$^{-1}$. MgOH$^+$ was the dominant product at all temperatures, but Mg$^+$, the co-product of oxygen-atom transfer to form methanol, was observed with a product branching fraction of $0.08 \pm 0.03 (T / 300 $ K$)^{-0.8 \pm 0.7}$. Reactive molecular dynamics simulations using a reactive force field, as well as a neural network yield rate coefficients about one order of magnitude lower. This underestimation of the rates is traced back to the multireference character of the transition state [MgOCH$_4$]$^+$. Statistical modeling of the temperature-dependent kinetics provides further insight into the reactive potential surface. The rate limiting step was found to be consistent with a four-centered activation of the C-H bond, consistent with previous calculations. The product branching was modeled as a competition between dissociation of an insertion intermediate directly after the rate-limiting transition state, and traversing a transition state corresponding to a methyl migration leading to a Mg-CH$_3$OH$^+$ complex, though only if this transition state is stabilized significantly relative to the dissociated MgOH$^+$ + CH$_3$ product channel. An alternative non-statistical mechanism is discussed, whereby a post-transition state bifurcation in the potential surface could allow the reaction to proceed directly from the four-centered TS to the Mg-CH$_3$OH$^+$ complex thereby allowing a more robust competition between the product channels.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Water Dynamics Around Proteins: T- and R-States of Hemoglobin and Melittin

The water dynamics, as characterized by the local hydrophobicity (LH), is investigated for tetrameric hemoglobin and dimeric melittin. For the T0 to R0 transition in Hb it is found that LH provides additional molecular-level insight into the Perutz mechanism, i.e., the breaking and formation of salt bridges at the alpha1 / beta2 and alpha2 / beta1 interface is accompanied by changes in LH. For Hb in cubic water boxes with 90 Aengstroem and 120 Aengstroem edge length it is observed that following a decrease in LH as a consequence of reduced water density or change of water orientation at the protein/water interface the alpha / beta interfaces are destabilized; this is a hallmark of the Perutz stereochemical model for the T to R transition in Hb. The present work thus provides a dynamical view of the classical structural model relevant to the molecular foundations of Hb function. For dimeric melittin, earlier results by Cheng and Rossky (Nature, 1998, 392, 696-699) are confirmed and interpreted on the basis of LH from simulations in which the protein structure is frozen. For the flexible melittin dimer the changes in the local hydration can be as much as 30 % than for the rigid dimer, reflecting the fact that protein and water dynamics are coupled.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Formation and Stabilization of Ground and Excited State Singlet O$_2$ upon Recombination of $^3$P Oxygen on Amorphous Solid Water

The recombination dynamics of $^3$P oxygen atoms on cold amorphous solid water to form triplet and singlet molecular oxygen (O$_2$) is followed under conditions representative for cold clouds. It is found that both, formation of ground state ($X ^3 Σ_{g}^{-}$) O$_2$ and molecular oxygen in the two lowest singlet states ($a ^1Δ_g$ and $b ^1Σ_g^+$) is possible and that the species can stabilize. The relative proportions of the species is approximately 1:1:1. These results also agree qualitatively with a kinetic model based on simplified wavepacket simulations. As the chemical reactivity of triplet and singlet O$_2$ is different it is likely that substantial amounts of $a ^1Δ_g$ and $b ^1Σ_g^+$ oxygen influences the chemical evolution of cold clouds.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

High-Dimensional Potential Energy Surfaces for Molecular Simulations

An overview of computational methods to describe high-dimensional potential energy surfaces suitable for atomistic simulations is given. Particular emphasis is put on accuracy, computability, transferability and extensibility of the methods discussed. They include empirical force fields, representations based on reproducing kernels, using permutationally invariant polynomials, and neural network-learned representations and combinations thereof. Future directions and potential improvements are discussed primarily from a practical, application-oriented perspective.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Near dissociation states for H$_2^+$-He on MRCI and FCI potential energy surfaces

A new analytical potential energy surface (PES) has been constructed for H$_2^+$-He using a reproducing kernel Hilbert space (RKHS) representation from an extensive number of $ab initio$ energies computed at the multi-reference and full configuration interaction level of theory. For the MRCI PES the long-range interaction region of the PES is described by analytical functions and is connected smoothly to the short-range interaction region, represented as a RKHS. All ro-vibrational states for the ground electronic state of H$_2^+$-He are calculated using two different methods to determine quantum bound states. Comparing transition frequencies for the near-dissociation states for $ortho$- and $para$-H$_2^+$-He allows assignment of the 15.2 GHz line to a $J=2$ $e/f$ parity doublet of $ortho$-H$_2^+$-He whereas the experimentally determined 21.8 GHz line is only consistent with a $(J=0)$ $\rightarrow$ $(J=1)$ $e/e$ transition in $para$-H$_2^+$-He.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Reactive Dynamics and Spectroscopy of Hydrogen Transfer from Neural Network-Based Reactive Potential Energy Surfaces

The in silico exploration of chemical, physical and biological systems requires accurate and efficient energy functions to follow their nuclear dynamics at a molecular and atomistic level. Recently, machine learning tools gained a lot of attention in the field of molecular sciences and simulations and are increasingly used to investigate the dynamics of such systems. Among the various approaches, artificial neural networks (NNs) are one promising tool to learn a representation of potential energy surfaces. This is done by formulating the problem as a mapping from a set of atomic positions $\mathbf{x}$ and nuclear charges $Z_i$ to a potential energy $V(\mathbf{x})$. Here, a fully-dimensional, reactive neural network representation for malonaldehyde (MA), acetoacetaldehyde (AAA) and acetylacetone (AcAc) is learned. It is used to run finite-temperature molecular dynamics simulations, and to determine the infrared spectra and the hydrogen transfer rates for the three molecules. The finite-temperature infrared spectrum for MA based on the NN learned on MP2 reference data provides a realistic representation of the low-frequency modes and the H-transfer band whereas the CH vibrations are somewhat too high in frequency. For AAA it is demonstrated that the IR spectroscopy is sensitive to the position of the transferring hydrogen at either the OCH- or OCCH$_3$ end of the molecule. For the hydrogen transfer rates it is demonstrated that the O-O vibration is a gating mode and largely determines the rate at which the hydrogen is transferred between the donor and acceptor. Finally, possibilities to further improve such NN-based potential energy surfaces are explored. They include the transferability of an NN-learned energy function across chemical species (here methylation) and transfer learning from a lower level of reference data (MP2) to a higher level of theory (pair natural orbital-LCCSD(T)).

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The ${\rm N}(^4S) +{\rm O}_2(X^3Σ^-_g) \leftrightarrow {\rm O}(^3P) + {\rm NO}(X^2Π)$ Reaction: Thermal and Vibrational Relaxation Rates for the $^{2}$A$&#39;$, $^{4}$A$&#39;$ and $^{2}$A$&#39;&#39;$ States

The kinetics and vibrational relaxation of the ${\rm N}(^4S) +{\rm O}_2(X^3Σ^-_g) \leftrightarrow {\rm O}(^3P) + {\rm NO}(X^2Π)$ reaction is investigated over a wide temperature range based on quasiclassical trajectory simulations on 3-dimensional potential energy surfaces (PESs) for the lowest three electronic states. Reference energies at the multi reference configuration interaction level are represented as a reproducing kernel and the topology of the PESs is rationalized by analyzing the CASSCF wavefunction of the relevant states. The forward rate matches one measurement at 1575 K and is somewhat lower than the high-temperature measurement at 2880 K whereas for the reverse rate the computations are in good agreement for temperatures between 3000 and 4100 K. The temperature-dependent equilibrium rates are consistent with results from JANAF and CEA results. Vibrational relaxation rates for O + NO($ν=1$) $\rightarrow$ O + NO($ν=0$) are consistent with a wide range of experiments. This process is dominated by the dynamics on the $^2$A$&#39;$ and $^4$A$&#39;$ surfaces which both contribute similarly up to temperatures $T \sim 3000$ K, and it is found that vibrationally relaxing and non-relaxing trajectories probe different parts of the potential energy surface. The total cross section depending on the final vibrational state monotonically decreases which is consistent with early experiments and previous simulations but at variance with other recent experiments which reported an oscillatory cross section.

0 citations0 reviewsNo code linkedOpen access