Efficient Generation of Torsional Energy Profiles by Multifidelity Gaussian Processes for Hindered Rotor Corrections.

Accurate thermochemistry computations often require proper treatment of torsional modes. The one-dimensional hindered rotor model has proven to be a computationally efficient solution, given a sufficiently accurate potential energy surface. Methods that provide potential energies at various compromises of uncertainty and computational time demand can be optimally combined within a multifidelity treatment.

Automatic Potential Energy Surface Exploration by Accelerated Reactive Molecular Dynamics Simulations: From Pyrolysis to Oxidation Chemistry.

Automatic potential energy surface (PES) exploration is important to a better understanding of reaction mechanisms. Existing automatic PES mapping tools usually rely on predefined knowledge or computationally expensive on-the-fly quantum-chemical calculations. In this work, we have developed the PESmapping algorithm for discovering novel reaction pathways and automatically mapping out the PES using merely one starting species is present.

The first HyDRA challenge for computational vibrational spectroscopy.

Vibrational spectroscopy in supersonic jet expansions is a powerful tool to assess molecular aggregates in close to ideal conditions for the benchmarking of quantum chemical approaches. The low temperatures achieved as well as the absence of environment effects allow for a direct comparison between computed and experimental spectra. This provides potential benchmarking data which can be revisited to hone different computational techniques, and it allows for the critical analysis of procedures under the setting of a blind challenge.

Improved modeling of anharmonicity for furan microsolvation.

Computational benchmark data for complexes requires accurate models of anharmonic torsional motion. State-of-the-art hindered rotor treatments come with a number of difficulties, regarding discontinuities from badly converged points or coupling, oscillations, or the consideration and correction of stationary points. Their manual handling introduces a level of arbitrariness not suitable for benchmark procedures.

A Reactive Molecular Dynamics Study of Chlorinated Organic Compounds. Part II: A ChemTraYzer Study of Chlorinated Dibenzofuran Formation and Decomposition Processes.

The front cover artwork is provided by Prof. K. Leonhard's group at RWTH Aachen University.

Influence of Oxygen Atoms and Ring Strain on the Low-Temperature Oxidation Pathways of 1,3-Dioxolane.

Biohybrid fuels are a promising solution for making the transportation sector more environmentally friendly. One such interesting fuel candidate is 1,3-dioxolane, which can be produced from inedible biomass. However, very little kinetics data are available for the low-temperature oxidation of this fuel molecule.

A Reactive Molecular Dynamics Study of Chlorinated Organic Compounds. Part I: Force Field Development.

This work presents a novel parametrization for the ReaxFF formalism as a means to investigate reaction processes of chlorinated organic compounds. Force field parameters cover the chemical elements C, H, O, Cl and were obtained using a novel optimization approach involving relaxed potential energy surface scans as training targets. The resulting ReaxFF parametrization shows good transferability, as demonstrated on two independent ab initio validation sets.

A Reactive Molecular Dynamics Study of Chlorinated Organic Compounds. Part II: A ChemTraYzer Study of Chlorinated Dibenzofuran Formation and Decomposition Processes.

In our two-paper series, we first present the development of ReaxFF CHOCl parameters using the recently published ParAMS parametrization tool. In this second part, we update the reactive Molecular Dynamics - Quantum Mechanics coupling scheme ChemTraYzer and combine it with our new ReaxFF parameters from Part I to study formation and decomposition processes of chlorinated dibenzofurans. We introduce a self-learning method for recovering failed transition-state searches that improves the overall ChemTraYzer transition-state search success rate by 10 percentage points to a total of 48 %.

Completely Computational Model Setup for Spectroscopic Techniques: The Molecular Dynamics Indirect Hard Modeling Approach.

The spectroscopic quantification of mixture compositions usually requires pure compounds and mixtures of known compositions for calibration. Since they are not always available, methods to fill such gaps have evolved, which are, however, not generally applicable. Therefore, calibration can be extremely challenging, especially when multiple unstable species, for example, intermediates, exist in a system.

Efficient Reaction Space Exploration with ChemTraYzer-TAD.

The development of a reaction model is often a time-consuming process, especially if unknown reactions have to be found and quantified. To alleviate the reaction modeling process, automated procedures for reaction space exploration are highly desired. We present ChemTraYzer-TAD, a new reactive molecular dynamics acceleration technique aimed at efficient reaction space exploration.

Coupled Anharmonic Thermochemistry from Stratified Monte Carlo Integration.

This study presents configuration integral Monte Carlo integration (CIMCI), a new semiclassical method for handling fully coupled anharmonicity in gas-phase thermodynamics that promises to be black boxable, to be applicable to all kinds of anharmonicity, and to scale better at higher dimensionality than other methods for handling gas-phase molecular anharmonicity. The method does so using automatically and recursively stratified, simultaneous Monte Carlo (MC) integration of multiple functions, following a modified version of the standard MISER scheme that converges at a rate of about the square of naïve MC integration. For the small systems analyzed by this study where proper reference data is available (HO and HO), the method's anharmonic entropy corrections match reference data better than those of other black box anharmonic methods, e.

The first microsolvation step for furans: New experiments and benchmarking strategies.

The site-specific first microsolvation step of furan and some of its derivatives with methanol is explored to benchmark the ability of quantum-chemical methods to describe the structure, energetics, and vibrational spectrum at low temperature. Infrared and microwave spectra in supersonic jet expansions are used to quantify the docking preference and some relevant quantum states of the model complexes. Microwave spectroscopy strictly rules out in-plane docking of methanol as opposed to the top coordination of the aromatic ring.

Correcting Rate Constants from Anharmonic Molecular Dynamics for Quantum Effects.

Anharmonicity can greatly affect rate constants. One or even several orders of magnitude of deviation are found for obtaining rate constants using the standard rigid-rotor harmonic-oscillator model. In turn, reactive molecular dynamics (MD) simulations are a powerful way to explore chemical reaction networks and calculate rate constants from the fully anharmonic potential energy surface.

Automated Chemical Kinetic Modeling via Hybrid Reactive Molecular Dynamics and Quantum Chemistry Simulations.

An automated scheme for obtaining chemical kinetic models from scratch using reactive molecular dynamics and quantum chemistry simulations is presented. This methodology combines the phase space sampling of reactive molecular dynamics with the thermochemistry and kinetics prediction capabilities of quantum mechanics. This scheme provides the NASA polynomial and modified Arrhenius equation parameters for all species and reactions that are observed during the simulation and supplies them in the ChemKin format.

The furan microsolvation blind challenge for quantum chemical methods: First steps.

Herein we present the results of a blind challenge to quantum chemical methods in the calculation of dimerization preferences in the low temperature gas phase. The target of study was the first step of the microsolvation of furan, 2-methylfuran and 2,5-dimethylfuran with methanol. The dimers were investigated through IR spectroscopy of a supersonic jet expansion.

Assessing Statistical Uncertainties of Rare Events in Reactive Molecular Dynamics Simulations.

Reactive molecular dynamics (MD) simulations are a versatile tool which allow for studying reaction pathways and rates simultaneously. However, most reactions will be observed only a few times in such a simulation due to computational limitations or slow kinetics, and it is unclear how this will influence the obtained rate constants. Therefore, we propose a method based on the Poisson distribution to assess the statistical uncertainty of reaction rate constants obtained from reactive MD simulations.

Prediction of Chain Propagation Rate Constants of Polymerization Reactions in Aqueous NIPAM/BIS and VCL/BIS Systems.

Microgels have a wide range of possible applications and are therefore studied with increasing interest. Nonetheless, the microgel synthesis process and some of the resulting properties of the microgels, such as the cross-linker distribution within the microgels, are not yet fully understood. An in-depth understanding of the synthesis process is crucial for designing tailored microgels with desired properties.

General formulation of rovibrational kinetic energy operators and matrix elements in internal bond-angle coordinates using factorized Jacobians.

We show how inverse metric tensors and rovibrational kinetic energy operators in terms of internal bond-angle coordinates can be obtained analytically following a factorization of the Jacobian worked out by Frederick and Woywod. The structure of these Jacobians is exploited in two ways: On one hand, the elements of the metric tensor as well as its determinant all have the form ∑rsin(α)cos(β). This form can be preserved by working with the adjugate metric tensor that can be obtained without divisions.

Efficient yet accurate approximations for ab initio calculations of alcohol cluster thermochemistry.

We have calculated the binding enthalpies and entropies of gas phase alcohol clusters from ethanol to 1-decanol. In addition to the monomers, we have investigated dimers, tetramers, and pentamers. Geometries have been obtained at the B3LYP/TZVP level and single point energy calculations have been performed with the Resolution of the Identity-MP2 (RIMP2) method and basis set limit extrapolation using aug-cc-pVTZ and aug-cc-pVQZ basis sets.

Automated discovery of reaction pathways, rate constants, and transition states using reactive molecular dynamics simulations.

We provide a methodology for deducing quantitative reaction models from reactive molecular dynamics simulations by identifying, quantifying, and evaluating elementary reactions of classical trajectories. Simulations of the inception stage of methane oxidation are used to demonstrate our methodology. The agreement of pathways and rates with available literature data reveals the potential of reactive molecular dynamics studies for developing quantitative reaction models.

Hydrogen abstraction from n-butyl formate by H˙ and HO₂˙.

The combustion chemistry of esters has been elucidated in the past through the study of smaller formates and acetates. Hydrogen abstraction from the fuel as an initiation step is mostly modeled based on estimations for similar abstractions from nonoxygenated hydrocarbons. This study reports computed ab initio rates for abstractions by H˙ and HO₂˙ radicals from the recently proposed biofuel candidate n-butyl formate.