Ab initio based potential energy surface and kinetics study of the OH + NH3 hydrogen abstraction reaction

Caracal: A Versatile Ring Polymer Molecular Dynamics Simulation Package

A new open-source program package named Caracal covering simulations of molecular systems with ring polymer molecular dynamics (RPMD) is presented. It combines a powerful RPMD implementation including chemical reaction rate calculations and biased periodic and nonperiodic samplings with a collection of easy to set up potential energy surface (PES) methodologies, thus delivering an all-inclusive approach. Most implemented PESs are based on the QMDFF and EVB-QMDFF methods. Where the quantum mechanically derived force field (QMDFF) can be set up for an arbitrary molecular system in a black-box fashion, the empirical valence bond (EVB)-QMDFF connects two QMDFFs and is able to represent the PES of a chemical reaction. With our previously published flavors of this composite method, PESs for almost arbitrary gas phase thermal ground state reactions can be set up. Given an optimized reaction path, the mechanism of the reaction can be classified and RPMD rate constants can be obtained via umbrella sampling and recrossing calculations on an EVB-QMDFF PES. Further, QMDFFs can be polymerized for the description of liquid systems. In this paper, the internal structure as well as the handling philosophy of Caracal are outlined. Further, examples of the different possible kinds of calculations are given.

Current Status of the X + C2H6 [X ≡ H, F(2P), Cl(2P), O(3P), OH] Hydrogen Abstraction Reactions: A Theoretical Review

This paper is a detailed review of the chemistry of medium-size reactive systems using the following hydrogen abstraction reactions with ethane, X + C₂H₆ → HX + C₂H₅; X ≡ H, F(²P), Cl(²P), O(³P) and OH, and focusing attention mainly on the theoretical developments. These bimolecular reactions range from exothermic to endothermic systems and from barrierless to high classical barriers of activation. Thus, the topography of the reactive systems changes from reaction to reaction with the presence or not of stabilized intermediate complexes in the entrance and exit channels. The review begins with some reflections on the inherent problems in the theory/experiment comparison. When one compares kinetics or dynamics theoretical results with experimental measures, one is testing both the potential energy surface describing the nuclei motion and the kinetics or dynamics method used. Discrepancies in the comparison may be due to inaccuracies of the surface, limitations of the kinetics or dynamics methods, and experimental uncertainties that also cannot be ruled out. The paper continues with a detailed review of some bimolecular reactions with ethane, beginning with the reactions with hydrogen atoms. The reactions with halogens present a challenge owing to the presence of stabilized intermediate complexes in the entrance and exit channels and the influence of the spin-orbit states on reactivity. Reactions with O(³P) atoms lead to three surfaces, which is an additional difficulty in the theoretical study. Finally, the reactions with the hydroxyl radical correspond to a reactive system with ten atoms and twenty-four degrees of freedom. Throughout this review, different strategies in the development of analytical potential energy surfaces describing these bimolecular reactions have been critically analyzed, showing their advantages and limitations. These surfaces are fitted to a large number of ab initio calculations, and we found that a huge number of calculations leads to accurate surfaces, but this information does not guarantee that the kinetics and dynamics results match the experimental measurements.

Effect of (H2O)n (n = 0–3, 13) on the NH3 + OH reaction in the gas and liquid phases

We studied the effect of water clusters on the NH₃ + OH reaction at both the DFT and CCSD(T) levels. The calculated rate constants for the pure reaction are 2.07 × 10⁻¹³ and 1.35 × 10⁻¹³ cm³ molecule⁻¹ s⁻¹ in the gas and liquid phases, respectively, and the gas-phase rate constants are consistent with the corresponding experimental result (1.70 × 10⁻¹³ cm³ molecule⁻¹ s⁻¹), while the liquid-phase rate constants are slightly smaller than the experimental value (5.84 × 10⁻¹² cm³ molecule⁻¹ s⁻¹). In the gas phase, the presence of (H₂O)_n (n = 1–3) decreases the rate constant compared to the pure NH₃ + OH reaction, and these results are in agreement with many reported H₂O-catalyzed reactions. For the liquid phase reaction, compared with the case of n = 0–3, when the size of the water molecule cluster surrounding the OH radical is n = 13, the rate constant of the title reaction increases. Our study also shows that proton transfer is also a factor which accelerates the liquid phase NH₃ + OH reaction.

We used DFT and CCSD(T) methods with PCM solvation model to study the water cluster effect and solvation effect on the NH₃ + OH reaction.

Quantum calculations of the photoelectron spectra of the OH-·NH3 anion: implications for OH + NH3→ H2O + NH2 reaction dynamics

We present the results of quantum dynamics calculations for analyzing the experimentally measured photoelectron spectra of the OH^-·NH₃ anion complex. Detachment of an excess electron of OH^-·NH₃ initially produces a molecular arrangement, which is close to the transition-state structure of the neutral OH + NH₃→ H₂O + NH₂ hydrogen abstraction reaction due to the Franck-Condon principle, and thus finally leads to the OH + NH₃ or H₂O + NH₂ asymptotic channel. We used both the path integral method and the reduced-dimensionality quantum wave packet method to simulate the photoelectron spectra of the OH^-·NH₃ anion. The calculated spectra were found to be in qualitative agreement with the experimental spectra. It was found that the photodetached complex mainly dissociates into the OH + NH₃ channel; however, we found that the hydrogen exchange process also contributes to the photodetachment spectra.

Quasi-Classical Trajectory Study of the CN + NH3 Reaction Based on a Global Potential Energy Surface

Based on a combination of valence-bond and molecular mechanics functions which were fitted to high-level ab initio calculations, we constructed an analytical full-dimensional potential energy surface, named PES-2020, for the hydrogen abstraction title reaction for the first time. This surface is symmetrical with respect to the permutation of the three hydrogens in ammonia, it presents numerical gradients and it improves the description presented by previous theoretical studies. In order to analyze its quality and accuracy, stringent tests were performed, exhaustive kinetics and dynamics studies were carried out using quasi-classical trajectory calculations, and the results were compared with the available experimental evidence. Firstly, the properties (geometry, vibrational frequency and energy) of all stationary points were found to reasonably reproduce the ab initio information used as input; due to the complicated topology with deep wells in the entrance and exit channels and a “submerged” transition state, the description of the intermediate complexes was poorer, although it was adequate to reasonably simulate the kinetics and dynamics of the title reaction. Secondly, in the kinetics study, the rate constants simulated the experimental data in the wide temperature range of 25–700 K, improving the description presented by previous theoretical studies. In addition, while previous studies failed in the description of the kinetic isotope effects, our results reproduced the experimental information. Finally, in the dynamics study, we analyzed the role of the vibrational and rotational excitation of the CN(v,j) reactant and product angular scattering distribution. We found that vibrational excitation by one quantum slightly increased reactivity, thus reproducing the only experimental measurement, while rotational excitation strongly decreased reactivity. The scattering distribution presented a forward-backward shape, associated with the presence of deep wells along the reaction path. These last two findings await experimental confirmation.

Progress and Prospective of Nitrogen-Based Alternative Fuels

Alternative fuels are essential to enable the transition to a sustainable and environmentally friendly energy supply. Synthetic fuels derived from renewable energies can act as energy storage media, thus mitigating the effects of fossil fuels on environment and health. Their economic viability, environmental impact, and compatibility with current infrastructure and technologies are fuel and power source specific. Nitrogen-based fuels pose one possible synthetic fuel pathway. In this review, we discuss the progress and current research on utilization of nitrogen-based fuels in power applications, covering the complete fuel cycle. We cover the production, distribution, and storage of nitrogen-based fuels. We assess much of the existing literature on the reactions involved in the ammonia to nitrogen atom pathway in nitrogen-based fuel combustion. Furthermore, we discuss nitrogen-based fuel applications ranging from combustion engines to gas turbines, as well as their exploitation by suggested end-uses. Thereby, we evaluate the potential opportunities and challenges of expanding the role of nitrogen-based molecules in the energy sector, outlining their use as energy carriers in relevant fields.

Molecular Dynamics Simulations on Relaxed Reduced-Dimensional Potential Energy Surfaces

Molecular dynamics (MD) simulations with full-dimensional potential energy surfaces (PESs) obtained from high-level ab initio calculations are frequently used to model reaction dynamics of small molecules (i.e., molecules with up to 10 atoms). Construction of full-dimensional PESs for larger molecules is, however, not feasible since the number of ab initio calculations required grows rapidly with the increase of dimension. Only a small number of coordinates are often essential for describing the reactivity of even very large systems, and reduced-dimensional PESs with these coordinates can be built for reaction dynamics studies. While analytical methods based on transition-state theory framework are well established for analyzing the reduced-dimensional PESs, MD simulation algorithms capable of generating trajectories on such surfaces are more rare. In this work, we present a new MD implementation that utilizes the relaxed reduced-dimensional PES for standard microcanonical (NVE) and canonical (NVT) MD simulations. The method is applied to the pyramidal inversion of a NH₃ molecule. The results from the MD simulations on a reduced, three-dimensional PES are validated against the ab initio MD simulations, as well as MD simulations on full-dimensional PES and experimental data.

High-level theoretical study of the reaction between hydroxyl and ammonia: Accurate rate constants from 200 to 2500 K

Hydrogen abstraction from NH₃ by OH to produce H₂O and NH₂-an important reaction in combustion of NH₃ fuel-was studied with a theoretical approach that combines high level quantum chemistry and advanced chemical kinetics methods. Thermal rate constants calculated from first principles agree well (within 5%-20%) with available experimental data over a temperature range that extends from 200 to 2500 K. Quantum mechanical tunneling effects were found to be important; they lead to a decided curvature and non-Arrhenius behavior for the rate constant.

F(2P) + C2H6 → HF + C2H5 kinetics study based on a new analytical potential energy surface

An exhaustive kinetics study was performed for the title reaction using two theoretical approaches: variational transition-state theory and quasi-classical trajectory calculations, based on an original new analytical full-dimensional potential energy surface, named PES-2018, which has been fitted to high-level ab initio calculations.

Kinetics study of the CN + CH4 hydrogen abstraction reaction based on a new ab initio analytical full-dimensional potential energy surface

We have developed an analytical full-dimensional potential energy surface, named PES-2017, for the gas-phase hydrogen abstraction reaction between the cyano radical and methane.

Variational transition state theory: theoretical framework and recent developments

This article reviews the fundamentals of variational transition state theory (VTST), its recent theoretical development, and some modern applications.

Dynamics study of the OH + NH3 hydrogen abstraction reaction using QCT calculations based on an analytical potential energy surface

To understand the reactivity and mechanism of the OH + NH3 → H2O + NH2 gas-phase reaction, which evolves through wells in the entrance and exit channels, a detailed dynamics study was carried out using quasi-classical trajectory calculations. The calculations were performed on an analytical potential energy surface (PES) recently developed by our group, PES-2012 [Monge-Palacios et al. J. Chem. Phys. 138, 084305 (2013)]. Most of the available energy appeared as H2O product vibrational energy (54%), reproducing the only experimental evidence, while only the 21% of this energy appeared as NH2 co-product vibrational energy. Both products appeared with cold and broad rotational distributions. The excitation function (constant collision energy in the range 1.0-14.0 kcal mol(-1)) increases smoothly with energy, contrasting with the only theoretical information (reduced-dimensional quantum scattering calculations based on a simplified PES), which presented a peak at low collision energies, related to quantized states. Analysis of the individual reactive trajectories showed that different mechanisms operate depending on the collision energy. Thus, while at high energies (E(coll) ≥ 6 kcal mol(-1)) all trajectories are direct, at low energies about 20%-30% of trajectories are indirect, i.e., with the mediation of a trapping complex, mainly in the product well. Finally, the effect of the zero-point energy constraint on the dynamics properties was analyzed.

Atmospheric formation of the NO3 radical from gas-phase reaction of HNO3 acid with the NH2 radical: proton-coupled electron-transfer versus hydrogen atom transfer mechanisms

The amidogen radical abstracts the hydrogen from nitric acid through a proton coupled electron transfer mechanism rather than by an hydrogen atom transfer process.